Open Source Software


Normal-mode analysis of gravity due to earthquakes

eqpotential Computes the first-order Eulerian potential perturbation due to an earthquake
findrad Nearest-neighbor interpolation of a sorted sequence with repeats
getmodel Reads an Earth model from the output format of minosy
getspheroidal Reads spheroidal-mode eigenfunctions computed by minosy from a binary table
modestrain Computes spheroidal-mode strain eigenfunctions
modesum Performs normal-mode summation efficiently
prem800 Link an ASCII file output by minosy, the Earth model
fct_prem800S.direct Link to a second binary file output by minosy, the spheroidal modes file for this Earth model
smoment Calculates a scalar seismic moment

Spherical harmonics and localization: the definitive suite

Supported by the U. S. National Science Foundation under grant EAR-0105387.
abelpoisson Calculates the Abel-Poisson kernel (Amirbekyan, Michel & Simons 13-14)
bcoupling Calculates the boxcar coupling matrix (Dahlen & Simons 56)
bpboxcap Calculates the power spectrum of zonal boxcar caps or belts (Dahlen & Simons 48-50)
galpha Spatial eigensolutions of the localization problem (Simons, Dahlen & Wieczorek 3.15/5.7)
galphapto Spatial eigensolutions of the localization problem for Euler-angle rotated functions
gammap Calculates a matrix important to calculate the variance of multitaper estimates (Dahlen & Simons 167)
gammab Another matrix important to calculate the variance of multitaper estimates (Dahlen & Simons 164)
geoboxcap Calculates the power spectrum of geographically defined boxcars
glmalpha Spectral eigensolutions of the localization problem (Simons, Dahlen & Wieczorek 4.6/5.2)
glmalphapto Spectral eigensolutions of the localization problem for Euler-angle rotated functions
mcoupling Calculates the individual multitaper coupling matrix (Dahlen & Simons 131)
mcouplings Calculates the eigenvalue-weighted multitaper coupling matrix (Dahlen & Simons 145)
multivar Calculates the multitaper estimation variance, the fast way (Dahlen & Simons 163)
mvarratios Calculates the eigenvalue-weighted multitaper estimation variance (Dahlen & Simons 178)
periodovar Calculates the periodogram estimation variance (Dahlen & Simons 174-175)
plm Associated Legendre functions (Dahlen & Tromp B.48)
plm2map Makes an array of spherical harmonic expansion coefficients nicely mappable
plm2slep Performs a Slepian expansion of a function of known harmonics
ptoslep Calculates single-order single cap Slepian functions rotated about three Euler angles
slep2map Makes an array of Slepian expansions coefficients nicely mappable
slep2plm Performs an inverse Slepian expansion of a set of known Slepian coefficients
ulmmp Single-degree transformation matrix that maps real to complex harmonics
ummp Block-transformation matrix that maps real to complex harmonics
plm Associated Legendre functions (Dahlen & Tromp B.48)
rcenter Computes the centroid or center of gravity of a closed domain on the surface of the unit sphere
xlm Normalized associated Legendre polynomials (Dahlen & Tromp B.58)
xyz2slep Local Slepian-basis expansion of scattered data points
xyz2spl Spline interpolation using the Abel-Poisson kernel
ylm Normalized real spherical harmonics (Dahlen & Tromp B.72)
pxyerh Error handling for PLM, XLM and YLM
pxynrm Normalization verification (Dahlen & Tromp B.49/B.59/B.73)
randpatch Calculates points randomly distributed on a spherical patch
randsphere Calculates points randomly distributed on a spherical surface
whitevariance Calculates the multitaper estimation variance, the slow way (Dahlen & Simons 162-163)

Spatiospectral localization on the sphere

africa Coordinates and splining of Africa
amazon Coordinates and splining of the Amazon basin
antarctica Coordinates and splining of a rotated Antarctica
australia Coordinates and splining of Australia
bezier Smooths a geograhic curve by B-spline fitting
boxcap Calculates zonal polar boxcar tapers
boxcaptarget Targets a zonal polar boxcar tapers to a point on the sphere
caploc Calculates the location of a spherical cap
cdlkernel Calculates a Christoffel-Darboux kernel with Legendre functions
christoffeldarboux Verifies various Christoffel-Darboux theorems (Simons, Dahlen & Wieczorek 3.10/5.15)
curvecheck Checks that a closed geographic curve makes sense
dlmlmp Calculates the axisymmetric spatiospectral concentration kernel
dphpatch Calculates the integration domain of a spherical patch
dphregion Calculates the integration domain of a closed spherical surface
england Coordinates and splining of England
eurasia Coordinates and splining of Eurasia
greenland Coordinates and splining of Greenland
grunbaum Calculates circular single-cap tapers of a single angular order
grunbaum2 Calculates circular double-cap tapers of a single angular order
kernelc Finds the kernel whose eigenfunctions are optimally concentrated
klmlmp2rot Rotates a kernel so that its eigenfunctions are rotated versions of the unrotated kernel
localization Finds bandlimited functions localized to a closed spherical domain
namerica Coordinates and splining of North America
nsubm Calculates partial Shannon number at a single angular order
orinoco Coordinates and splining of the Orinoco basin
orthocheck Checks orthogonality, determines integration points, calculates eigenvalues
phicurve Finds longitude crossings at query latitudes for closed spherical curves
samerica Coordinates and splining of South America
sdupcont Shows the effect of upward/downward continuation of the Slepian functions
sdwallcons Calculates localization solutions to all of Earth's continents
sdwcap Calculates circular bandlimited single-cap tapers of a single angular order
sdwcap2 Calculates circular bandlimited double-cap tapers of a single angular order
sdwcapt Calculates circular spacelimited single-cap tapers of a single angular order
sdwcapt2 Calculates circular spacelimited double-cap tapers of a single angular order
sdwelm Computes global ranking of eigenvalues for the single-cap concentration problem
sdwtarget Targets a single-order circular taper to a point on the sphere
sinsin Performs longitudinal integral with (co)sine functions
sincos Performs longitudinal integral with (co)sine functions
coscos Performs longitudinal integral with (co)sine functions
coscosexp Performs longitudinal integral with (co)sine functions and complex exponentials
sincosexp Performs longitudinal integral with (co)sine functions and complex exponentials
spharea Computes areas of closed domains on the unit sphere

Spatiospectral localization on a sphere: Figures

sdwdiagram Figure 3.1 from Simons et al., SIAM Review (2006)
sdwspace Figure 5.1 from Simons et al., SIAM Review (2006)
sdwspectral Figure 5.2 from Simons et al., SIAM Review (2006)
sdwvals Figure 5.3 from Simons et al., SIAM Review (2006)
sdwfried Figure 5.4 from Simons et al., SIAM Review (2006)
sdwsumall Figure 5.5 from Simons et al., SIAM Review (2006)
sdwcmb Figure 5.6 from Simons et al., SIAM Review (2006)
sdwregions Figure 6.1, 6.2 and 6.3 from Simons et al., SIAM Review (2006)
sdwconsum Figure 6.4 from Simons et al., SIAM Review (2006)

The spherical harmonics transform and its inverse

addmabout Find the running index in a block array with SH degrees and orders
addmin Repeats the degrees in a vector l+1 times
addmoff Calculates total number of orders belonging to an alternative SH expansion
addmon Lists degrees and orders belonging to a SH expansion
addmout Lists degrees and orders belonging to an alternative SH expansion
addmup Calculates total number of orders belonging to a SH expansion
blanco Calculates spherical harmonics rotation (Wigner-D) matrix, method I
cpx2rsh Converts complex to real spherical harmonics coefficients
datafit Least-squares solutions to a linear system with (un)known data covariance
dlmb Calculates spherical harmonics rotation (Wigner-D) matrix, method II
gaunt Calculates Gaunt coefficients using three different algorithms
gausslegendre Performs Gauss-Legendre quadrature
igrf10 Illustrates some of these routines to plot Earth's geomagnetic field
gausslegendrecof Calculates weights and abscissas for Gauss-Legendre integration
geopotential Plots geopotential fields for Earth, Moon and Mars
guseinov Calculates Clebsch-Gordan and Gaunt coefficients using binomial coefficients
legendrecof Calculates roots and coefficients of Legendre polynomials
legendrediff Computes derivatives of Legendre functions
legendreprodint Integrates products of Legendre functions
legendreroots Calculates roots of Legendre polynomials
libbrecht Libbrecht's algorithm to calculate high-order Legendre functions
plmfilt Filters a set of spherical harmonic expansion coeffficients
plm2pot Reads in and scales geopotential coeffficients
plm2rot Rotation of spherical harmonics coefficients using Euler angles
plm2spec Calculates the power spectrum of real spherical harmonic coefficients
plm2th Inverse single-order spherical harmonic transform
plm2xyz Inverse spherical harmonic transform
plotonsphere Plots data on a sphere
plotonearth Plots data on an Earth sphere with continents
plotplm Plots spherical harmonics in various ways
plotslep Plots Slepian functions in various ways
pl2th Inverse Legendre transform
rotenberg Calculates Wigner symbol values from Rotenberg's peculiar notation
rsh2cpx Converts real to complex spherical harmonics coefficients
rotcof Rotation of spherical harmonics coefficients around z-axis
shcos Retrieves spherical harmonics cosine coefficients from matrix listing
shsin Retrieves spherical harmonics sine coefficients from matrix listing
simpson Performs Simpson quadrature
sinsinexp Performs longitudinal integral with (co)sine functions and complex exponentials
sixj Calculates Wigner 6j coefficients, from lookup table
th2pl Forward Legendre transform
th2plm Forward single-order spherical harmonic transform
threej Calculates Wigner 3j coefficients, from lookup table
topography Plots topography fields for Earth, Moon and Mars
trapeze Performs trapezoidal quadrature
triangle Checks the triangle inequality
wignercycle Calculates Wigner 3j coefficients, for lookup table
wigner0j Calculates zero-bottom-row Wigner 3j coefficients, by recursion
wigner3j Calculates Wigner 3j coefficients, by recursion... so-so algorithm
wigner3jm Calculates Wigner 3j coefficients, by recursion... better algorithm
wigner6j Calculates Wigner 6j coefficients, by recursion
xyz2plm Forward spherical harmonic transform
zeroj Calculates zero-bottom-row Wigner 3j coefficients, from lookup table

Spatiospectral localization in the Cartesian plane

Supported by the U. S. National Science Foundation under grant EAR-0710860.
devilliers Builds a matrix whose eigenvectors are Jacobi expansion coefficients for radial Slepian functions
jacobi Calculates Jacobi polynomials in an explicit way
kernelc2D Finds the kernel for bandlimited concentration to a closed Cartesian domain
localization2D Finds bandlimited functions localized to a closed Cartesian domain
swconsum2d Plots the cumulative energy of geographical Slepian functions in the spatial and spectral domains
swdisk Constructs radial Slepian functions for concentration to a disk in the plane
swfried2d Plots Slepian functions for the disk in the spatial domain
swregions2d Plots Slepian functions for geographical areas in the spatial domain
swspace2d Plots radial Slepian functions for the disk in the spatial domain
swspectral2d Plots radial Slepian functions for the disk in the spectral domain
swvals2d Eigenvalue plot of the Slepian functions for the disk

Earthquake early warning

wtALL Analyzes a whole directory of data in SAC form and prints results to a file
wtALL_diagnostics The actual results of the above analysis used and interpreted by Simons et al., EPSL (2006)
allen1 One of a sequence of programs related to Simons et al., EPSL (2006)
allen2 One of a sequence of programs related to Simons et al., EPSL (2006)
allen3 One of a sequence of programs related to Simons et al., EPSL (2006)
allen4 One of a sequence of programs related to Simons et al., EPSL (2006)
allen5 One of a sequence of programs related to Simons et al., EPSL (2006)
allen6 One of a sequence of programs related to Simons et al., EPSL (2006)
allen7 One of a sequence of programs related to Simons et al., EPSL (2006)
allen8 One of a sequence of programs related to Simons et al., EPSL (2006)
allen9 One of a sequence of programs related to Simons et al., EPSL (2006)
locident An auxiliary program to the above codes used to Simons et al., EPSL (2006)
statident An auxiliary program to the above codes used to Simons et al., EPSL (2006)
installen Computes the instantaneous frequency of a signal

The wavelet transform and its inverse in various incarnations

abank Analysis filter bank
compres Figures out compression level of the wavelet approximation
decay Figures out the decay of the wavelet coefficients
denoise Denoises a signal by wavelet thresholding
dnums Figures out the number of scaling and wavelet coefficients
dyadplot Plots wavelet and scaling coefficients on dyadic grid
even Returns logical downsampling array
fbank Filter bank illustration of perfect reconstruction
graphs Computes graphs of wavelets and scaling functions
intersub Interpolation subdivision for interpolating boundary filters
iwt Inverse wavelet transform (this is the main routine)
landd Returns lengths and downsampling arrays
liftco Wavelet coefficients and lifting steps
makewc Makes and saves data base with wavelet coefficients
plotnorm Plots data on unit interval
polyphase Returns the polyphase forms for a given filter pair
prodco Checks the product filter
sbank Synthesis filter bank
succapp Successive approximation of a signal by wavelet analysis
succaplot Plots successive wavelet approximation of a signal
tenlectures63 Reproduces Figure 6.3 of Daubechies' "Ten Lectures" book
threshold Thresholds wavelet or scaling coefficients
threshold2 Thresholds wavelet or scaling coefficients slightly differently
wavetour710 Reproduces Figure 7.10 of Mallat's "Wavelet Tour" book
wavetour714 Reproduces Figure 7.14 of Mallat's "Wavelet Tour" book
wc Lists wavelet filter coefficients
wigner Computes the alias-free Wigner-Ville distribution
wt Forward wavelet transform (this is the main routine)
wtimax Makes a wavelet time axis for polyphase and lifting implementations
wtlift Fast discrete wavelet transform by (integer) lifting

Some very particular wavelet transforms

Supported by the U. S. National Science Foundation under grant DMS-0530865.
angularD4WT Forward and inverse 4-tap wavelet transform for the angular parts of a three-dimensional box
angularpart513to128 Dimensionality reduction of the angular functions from 513 to 128 specific pixels
cubejac Produces a single face of the cubed sphere of Ronchi et al.
cubemats Produces rotation matrices for the cubed sphere of Ronchi et al.
cube2sphere A special realization of the cubed sphere of Ronchi et al.
d4box Forward 4-tap wavelet transform in a finite three-dimensional box
d4boxcof Collects the coefficients for the 4-tap forward and inverse wavelet transform on the interval
d4boxi Inverse 4-tap wavelet transform in a finite three-dimensional box
d4boxstep One iteration in the forward 4-tap wavelet transform on the interval
d4boxstepi One iteration in the inverse 4-tap wavelet transform on the interval
imageletter A code to help make plots of (flattened) cubed spheres
interpS40RTS Makes maps at certain depths of Jeroen Ritsema's model s40_prelim.sph
plm2cube Transforms spherical harmonics to the faces of the cubed sphere
plotonchunk Plots data on the faces of the cubed sphere in the plane
plotconcube A code to make plots of (flattened) cubed spheres
preconD4 Preconditioning for the 4-tap wavelet transform for a three-dimensional box
radialHaarWT Forward and inverse Haar transform for the radial part of a three-dimensional box
radialpart129to37 Dimensionality reduction of the radial functions from 129 to 37 specific layers
raydyntrace An interface to the FORTRAN program raydyntrace.f
raykernels An interface to the MATLAB program raydyntrace.m
readGNmodel Read Raffa Montelli's model PRI-05 at one particular depth
readS40RTS Read a preliminary version of Jeroen Ritsema's model s40_prelim.sph
rtssplines Get spline expansion coefficients for Jeroen Ritsema's model s40_prelim.sph
skernelplot Makes a two-dimensional plot of a kernel at depth
sphere2cube Converts spherical coordinates to the equivalent standard cubed sphere coordinates

Robust power estimation and 1-D spectral analysis

bracewell Illustrates the real/imaginary properties of the Fourier transform
fhanning Calculates a Hann window for spectral analysis
jackknife Returns delete-one sampling matrix for a set of values
pchave Chave's method of robust spectral estimation
pchavedemo Illustration of Chave's method of spectral estimation
spectrogram Calculation of the spectrogram
specdensplot Calculates and then nicely plots the spectral density using PCHAVE
timdomplot Reads and then nicely plots seismic data using READSAC
timspecplot Calculates and then nicely plots the spectrogram using SPECTROGRAM

Lithospheric flexure, loading and 2-D spectral analysis

annulus Calculates a 2-dimensional annulus
bayram Calculates a Hermite-multiple window spectral estimate
choice Picks a set of Hermite multitapers
cohen Calculates the smoothing kernel for a Cohen's class spectral estimator
forsyth Predicts coherence and admittance using Forsyth's (1985) method
geninv Calculates the generalized inverse
halfmazi Calculates maximum and half maximum along an azimuth in a 2-D data set
halfmiso Calculates maximum and half maximum of an azimuthally averaged 2-D data set
hermite Computes Hermite-polynomial multitapers
hermitian Makes a matrix Hermitian so its inverse Fourier transform is real
knum2 Makes a Hermitian two-dimensional wavenumber scale
laguerre Computes Laguerre polynomials
loading Computes topography, gravity and their relation for an arbitrary loading model
loading1 A demo version of loading
maxazi Calculates mean and maximum of mean azimuthal line integral of a data set.
mckenzie Predicts coherence and admittance using McKenzie's (2003) method
mtm Two-dimensional multi-taper analysis
mtm3 Two-dimensional multi-taper analysis with Slepian tapers
mtm4 Two-dimensional multi-taper analysis with Hermite tapers
parker Computes terrain/Bouguer corrections using Parker's (1972) method
planef Calculates and/or draws a plane
planefit Calculates a best-fitting plane through a data set
randphase Returns random phase angles
syncoh Creates a data set with a known coherence to another one
transl Predicts transitional coherence wavelength for a loading model

Image analysis bits and pieces

blockcorr Calculates correlation coefficients between tiles of two images
blockcorr2 Uses BLOCKCORR to find best-fitting shifts to align two images
blockcorr3 Uses BLOCKCORR2 to successively correlate a series of images
blockisolate Retrieves a submatrix from within a larger matrix
blockisolate2 Retrieves a subregion from within a larger area, in physical coordinates
blockmean Calculates means of submatrices within a matrix
blocktile Calculates (the number of) overlapping squares tiling a matrix
innermerge Writes newly aligned images from parameters calculated by BLOCKCORR3

Seismic processing bits and pieces

readsac Reads SAC-formatted binary seismograms
writesac Writes SAC-formatted binary seismograms
lowpass Low pass filtering
bandpass Bandpass filtering
hipass High pass filtering
stalta STA/LTA triggering algorithm
wtevents Makes a diagnostic wavelet analysis plot of a seismogram
cmtday PERL to extract events by date from a `cat *dek` CMT master
cmtloc PERL to extract event attributes from a `cat *dek` CMT master
cmtsol PERL to extract an event by code from a `cat *dek` CMT master and write SPECFEM output

Surface-wave ray tracing and azimuthal anisotropy

Supported by the U. S. National Science Foundation under grant EAR-9614341.
acor2ind Matches a physical coordinate to a sequential equal-area grid cell index
aind2sub Returns row and 'column' index given a sequential equal-area grid cell number
apex Calculates highest points on a great circle
asub2ind Returns a sequential equal-area grid cell number given a row and 'column' index
authalic Calculates an equal-area grid for a regional tomography experiment
clutser Simplified clustering algorithm to equalize coverage in discrete cells
distgrc Calculates coordinates of points at a certain distance on a great circle
eqplot Plots an equal-area grid as calculated by AUTHALIC
fillauth Fills the cells of an equal-area grid from AUTHALIC with prescribed colors
fridplot Plots a regular grid
grcazim Calculates points lying on a great circle of specified azimuth
grcircle Calculates great circle points through a pair of points
grcdist Calculates uncorrected great-circle distance
gridcheck Check if grid dimensions are up to containing all event/station pairs
gridlod2 Loads the three-dimensionsal grid defined by GRIDINFO and INAMAK3D
gridsec Performs surface-wave ray tracing and calculates sensitivity matrix elements - with demo
lcsg Loads all matrices with information written by and for SURFTRACE
lcsd Loads some matrices with information written by SURFTRACE
neqplot Plots a clustered grid as calculated by CLUTSER
surfaz Calculates azimuthal coverage from the matrices written by SURFTRACE
surfdens Calculates ray path density from the matrices written by SURFTRACE
surftrace Uses GRIDSEC to calculate and write the matrices for inversion (this is the main routine)
truecourse Calculates the true course at the first point between great-circle points

Auxiliary functions, for graphics, indexing, etc.

addcb Adds a colorbar for use with multiple color maps
adrc Adds a row and a column to a matrix suitable for plotting
bindens Density plot for data where histograms are made for the y-data only
binsearch Performs a binary search through a presorted array
boxmid Finds the middle of a box
boxtex Makes rectangular boxes filled with text
cbarticks Forces agreement of colorbar tick marks and labels with figure
cellnan Initializes a cell array
circ Plots circles and circle segments
cor2ind Switches between physical coordinates and matrix indices
dat2col Converts data values to indices in one of multiple color maps
decibel Converts a range of values to decibel
defval Supplies default value to a named variable
degamini Undoes gamini which replicates elements of a vector
deggies Put degree signs on a graph
demmap Applies colormap suitable to plot Earth's topography
difer Evaluates the difference between two quantities
dindeks Returns indices to the ith plane perpendicular to the dth dimension of a 3D array
equistat Expresses a geographical grid in terms of its projected E-W and N-S distances
errorfill Plots error bars as shaded regions
errorxy Plots symmetric error bars on data in the the x and y direction
errorxy2 Plots asymmetric error bars on data in the the x and y direction
ext2lrtb Transforms 'extent' to 'left-right-top-bottom' coordinates for use with FILLBOX
fftaxis Makes a frequency axis for Fourier transformed data
fftaxis1D Makes a frequency axis for Fourier transformed time series
figdisp Suggests print command for figures
fig2print Makes the figure window wysiwyg
fillbox Makes rectangular boxes filled with color
flipflop Flips and duplicates a matrix suitable for spectral analysis
fralmanac An almanac with Earth and planetary data
fsize Determines a file size
gamini Replicates elements of a vector
gamini2 Replicates elements of a matrix
gauss Makes a 1-D gaussian
gauss2 Makes a 2-D gaussian
guess Returns a random integer
getcol Converts a color symbol into RGB values
getkids Gets axis children
getpos Gets axis handle positions
grey Defines a scaled grey color
halverange Sets color axis limits to reasonable values
imagef Plots a map using pixel-centered coordinates
imagefdir Plots a map using pixel-centered coordinates using direct indexing
imagefnan Plots a map with NaN's using pixel-centered coordinates
indeks Returns the ith entry of a matrix
ind2cor Switches between matrix indices and physical coordinates
insert Insert elements at certain positions into a vector
isav Returns azimuthally isotropic average of a two-dimensional array
kelicol Uses the color scheme designed by Hrafnkell Kárason
krijetem Creates figure panels and handles
kindeks Returns the ith column of a matrix
latitude Yields coordinates for a line of latitude
label Puts smart labels on figure panels
laxis Positions an axis on top of another one for easy annotation
lcs2spr Converts the format written by SURFTRACE to a conventionally sparse matrix
ldetrend Removes best-fitting linear trend from data
letter Returns letters from the alphabet
linazim Plots directional data on a map
longitude Yields coordinates for a meridian
loadb Loads binary data
longitude Yields coordinates for a meridian
lrtb2ext Transforms 'left-right-top-bottom' to 'extent' coordinates for use with FILLBOX
ls2cell Collects directory contents in a cell array
longticks Produces outward facing tickmarks
matranges A supremely useful function to make index vectors
minmax Finds the minimum and maximum value of a matrix
mollweide Projects data into a Mollweide projection
movev Moves a figure object vertically
moveh Moves a figure object horizontally
moving Computes moving average of a signal
nolabels Dispenses with axis labels
noticks Dispenses with axis ticks
nounder Gets rid of a filename underscore
number Adds numbers to the entries in a vector
openup Makes the axis slightly more than "tight"
osdep Changes flag according to the operating system
parse Performs some simple string parsing
pauli Returns a Pauli-type matrix useful in the generation of Toeplitz matrices
plotcont Plots continental outlines from a database
plotplates Plots plate boundaries from a database
plotprop Plots geographically distributed properties as color-scaled dots
plotprops Plots x, y, z scatter plots as color-scaled dots
polargrid Creates a geographical map-projection grid centered on the North pole, equally spaced in E-W and N-S
pref Gets a filename prefix
project Projects a simple coordinate transformation on a map
recindex Returns the mapping of one array into another
revolve Transforms a radial profile into a cylindrically symmetric data set
rindeks Returns the ith row of a matrix
rotdom Rotates a two-dimensional data domain over a specified angle
rots Rotates the Slepian polar-cap basis around the z-axis
rottp Rotates three-dimensional points about three Euler angles
rottick Rotates tick mark labels on a graph
rotx Rotates vector in three dimensions around the x-axis
roty Rotates vector in three dimensions around the y-axis
rotz Rotates vector in three dimensions around the z-axis
row2stats Find statistics of data presented belonging to different groups
scale Scales data values to fall within some range
seemax Equalizes axis ranges on plots
serre Brings subplot panels closer together
setnans Prevents color plotting values below a certain threshold
shrink Reduces width and height of a figure object
shuffle Shuffles the rows or entries in a matrix or vector
skip Reduces a vector by skipping certain positins
spr2lcs Converts a conventionally sparse matrix to the format written by SURFTRACE
star69 Returns calling function name it it exists
subnum Returns a figure panel subnumber
suf Gets a filename suffix
supertit Finds space for a title spanning several panels
tindeks Returns the ith third dimension of a 3D array
top Moves a graphics handle to the top of the stack
topofilt Makes filtered maps of EGM2008 topography and its gradients
tridiag Constructs a tridiagonal matrix
undeggies Removes degree signs from a graph
unproject Unprojects a simple coordinate transformation on a map
writeb Writes binary data
xpand Expands a graphics axis by a percentage of the plotted data range
xtraxis Adds a second axis to a plot
xtraxis1d Adds a logarithmic wavelength axis to a wavenumber plot
xtraxis1dlin Adds a linear wavelength axis to a wavenumber plot

Data sets, coefficients, colormaps, etc

conm.mat Continental data for use in PLOTCONT
cont.mtl Continental data for use in PLOTCONT
cost.mtl Continental data for use in PLOTCONT
Daubechies.mat Daubechies' wavelet coefficients for use in WC
CDF.mat Cohen, Daubechies, Feauveau wavelet coefficients for use in WC
CDFI.mat Cohen, Daubechies, Feauveau integer wavelet coefficients for use in WC
Earth.mat Earth data for use in FRALMANAC
gji2002lonlat Simons et al. (GJI, 2002) data set for use in the GRIDSEC demo
generic.m An example program with information pertaining to the depth grid for GRIDLOD2
gridinfo An example file with information pertaining to a surface grid for GRIDLOD2
inamak3d An example file with information pertaining to a depth grid for GRIDLOD2
kelim.mat A color scheme designed by Hrafnkell Kárason for use in KELICOL
Mars.mat Mars data for use in FRALMANAC
Moon.mat Moon data for use in FRALMANAC
plates.mtl Plate boundary data for use in PLOTPLATES
platm.mat Plate boundary data for use in PLOTPLATES
PRI5w.gz The coefficients of Raffa Montelli's PRI-05 model
SHM.mat Spherical harmonic data for use in FRALMANAC
staltademo.zip Sample seismograms for use in the STALTA demo
topo.mat A color scheme good to plot Earth's topography for use in DEMMAP
SHM.mat Spherical harmonic data for use in FRALMANAC
tapestry.mat A tectonic regionalization for North America
xpall.mat Spline-basis interpolation coefficients for 'all' depths used by RTSSPLINES
xpcube.mat Spline-basis interpolation coefficients for the 'cubed-sphere' depths used by RTSSPLINES
XYZ.mat Spatial data for use in FRALMANAC


Frederik Simons home
Last modified: Tue Nov 17 14:28:34 EST 2009
Valid HTML 4.0 Transitional