Initialize Earth layering geometry.

Arguments

• model::ModelData: The model data container used to store model parameters and arrays.

Keyword arguments:

• thick_lith::Float64:
  • Thickness in meters of lithosphere
• thick_crust::Float64:
  • Thickness in meters of crust
• thick_upper_crust::Float64:
  • Thickness in meters of upper crust
• thick_lower_crust::Float64:
  • Thickness in meters of lower crust
• thick_upper_lith::Float64:
  • Thickness in meters of upper mantle lithosphere
• thick_middle_lith::Float64:
  • Thickness in meters of middle lithosphere
• thick_lower_lith::Float64:
  • Thickness in meters of lower mantle lithosphere

Returns

• Nothing

Initialization Scenarios

Scenario 1: Layer thicknesses provided as input are used as-is for the upper continental crust (thick_upper_crust), lower continental crust (thick_lower_crust), upper mantle lithosphere (thick_upper_lith), middle mantle lithosphere (thick_middle_lith) and lower mantle lithosphere (thick_lower_lith).

Scenario 2: If the total thickness of the lithosphere (thick_lith), the total thickness of the continental crust (thick_crust), and the thickness of the upper continental crust (thick_upper_crust) are provided as inputs, remaining thicknesses are derived as follows:

• The thickness of the lower continental crust (thick_lower_crust) will be calculated as the difference between the total thickness of the continental crust (thick_crust) and the thickness of the upper continental crust (thick_upper_crust).
• The thickness of the mantle lithosphere (thick_mantle_lith) will be calculated as the difference between the thickness of the total lithosphere (thick_lith) and the thickness of the continental crust (thick_crust).
• The thickness of the upper mantle lithosphere (thick_upper_lith) will be set to one-third of the thickness of the mantle lithosphere (thick_mantle_lith).
• The thickness of the middle mantle lithosphere (thick_middle_lith) will be set to one-third of the thickness of the mantle lithosphere (thick_mantle_lith).
• The thickness of the lower mantle lithosphere (thick_lower_lith) will be set to one-third of the thickness of the mantle lithosphere (thick_mantle_lith).

Example 1: Direct input of layer thicknesses

eb = EarthBoxState(
    xnum = 200, ynum = 100, xsize = 500000.0, ysize = 120000.0,
    dx_marker = 100.0, dy_marker = 100.0
    )
model = eb.model_manager.model # ModelData container
EarthBox.MaterialGeometry.EarthLayering.initialize!(
    model, 
    thick_upper_crust = 22_000.00, 
    thick_lower_crust = 18_000.0, 
    thick_upper_lith  = 28_333.3, 
    thick_middle_lith = 28_333.3, 
    thick_lower_lith  = 22_000.0
)

Example 2: Calculating layer thicknesses from total lithosphere and crust thicknesses

eb = EarthBoxState(
    xnum = 200, ynum = 100, xsize = 500000.0, ysize = 120000.0,
    dx_marker = 100.0, dy_marker = 100.0
)
model = eb.model_manager.model # ModelData container
EarthBox.MaterialGeometry.EarthLayering.initialize!(
    model, 
    thick_lith        = 100_000.0, 
    thick_crust       = 40_000.0, 
    thick_upper_crust = 22_000.0
)