MODEL structure

The MODEL structure type defines the seismological structure of a planet in terms of discrete layers and named discontinuities. It allows to keep all necessary data within a single input parameter and thus a quick change between models and planets.

The structure fields can be divided into three groups:

  1. Layer Definition Stuff
  2. Discontinuity Definition Stuff
  3. Miscellaneous Stuff
Layer Definition Stuff
Field Type Units Description
z numeric vector km specifies the depth of the upper boundaries of all layers in model. z(1) is the top of the uppermost layer.
vp numeric vector km/s P wave velocities at depths specified in z. vp(1) corresponds to z(1)
vs numeric vector km/s S wave velocities at depths specified in z. vs(1) corresponds to z(1)
rho numeric vector g/cm3 Density at depth specified in z. rho(1) corresponds to z(1)
qp numeric vector dimensionless P wave quality factor at depth specified in z. qp(1) corresponds to z(1)
qs numeric vector dimensionless S wave quality factor at depth specified in z. qp(1) corresponds to z(1)
Discontinuity Definition Stuff
IMPORTANT: TTBOX needs the core mantle boundary and the inner core boundary being defined. Without defining these discontinuities, MKX4P and MKTIM4P will crash sooner or later, since CMB and ICB depths are used to restrict ray parameter ranges.
Field Type Units Description
conr numeric scalar km Depth to Conrad Discontinuity. Must be NaN if not specified.
moho numeric scalar km Depth to Mohorovicic Discontinuity. Must be NaN if not specified.
d410 numeric scalar km Depth to Olivine-α - β-spinel phase transition. Must be NaN if not specified. This is named after the mean depth on earth in order to get a short name.
d520 numeric scalar km Depth to β-spinel-&gamma-spinel phase transition. Must be NaN if not specified. This is named after the mean depth on earth in order to get a short name.
d660 numeric scalar km Depth to &gamma-spinel-Perovskite Boundary. Must be NaN if not specified. This is named after the mean depth on earth in order to get a short name.
cmb numeric scalar km Depth to the Core-Mantle-Boundary. Must be NaN if not specified.
icb numeric scalar km Depth to the Core-Mantle-Boundary. Must be NaN if not specified.
dz numeric vector km Depths of non-standard discontinuites. Must be empty if none specified.
dname string matrix - Names of non-standard discontinuities. dname(i,:) is the name string for the i-th non-standard discontinuity, the depth of which is given in dz(i). Must be empty if none specified.
Miscellaneuos Stuff
Field Type Units Description
rp numeric scalar km Radius of Planet. If not specified in .nd file, the largest depth value given in z is used.
name string - Model name, used in plot annotations.
year numeric year Year of model publication.

eof.