# Load seismic traces from .mat files *call trcmat nmodes.mat # matrix variable placed in .mat file # Each trace produces a row in this matrix: matrix Skin # first columns in 'matrix' taken from trace headers: column norder column lorder column freq # add columns by using projections onto basis functions defined in 'radial_model': basis Skin radial_model