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:
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) |
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. |
Field | Type | Units | Description |
---|---|---|---|
criticalrays | struct | - | This struct contains all new information in sub-fields as described below. It is generated on the fly if needed by calling MKIMPROVEMODEL, but you can also provide models that already contain it. |
.z | array | [km] | Subfield of the criticalryas field. This is a list of critical depths in the velocity model, e.g. depths of discontinuities, sudden changes in velocity gradient, Low velocity tops and bottoms and others. A critical depth is understood as a depth where the behaviour of rays changes. |
.p | array | s/° | Subfield of the criticalryas field. These are the ray parameters of P waves that have their turning point depths at the depths given in .z. With these ray parameters, rays can be sent exactly to the critical depths - these rays e.g. delimit of the various branches of the PKP phase. |
.s | array | s/° | Subfield of the criticalryas field. These are the ray parameters of S waves that have their turning point depths at the depths given in .z. With these ray parameters, rays can be sent exactly to the critical depths - these rays e.g. delimit of the various branches of the SKS phase. |
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. |