initialize!(model::ModelData; kwargs...)::Nothing

Initialize reference lithosphere parameters for relative base level model. The reference lithosphere is used to define global sea level (reference base level) as the top of the lithosphere column in isostatic equilibrium with the average pressure at the base of the model column (Kneller et al., 2023).

Arguments

• model::ModelData: The model data container containing the model parameters and arrays.

Keyword Arguments

• thickness_upper_continental_crust_ref::Float64
  • Thickness of upper continental crust of the reference lithosphere in meters
• thickness_lower_continental_crust_ref::Float64
  • Thickness of lower continental crust of the reference lithosphere in meters
• thickness_lithosphere_ref::Float64
  • Thickness of lithosphere for the reference lithosphere in meters
• gridy_spacing_ref::Float64
  • Spacing of grid cells in y-direction for the reference lithosphere in meters
• temperature_top_ref::Float64
  • Temperature at top of lithosphere for the reference lithosphere in Celsius
• temperature_moho_ref::Float64
  • Temperature at Moho of the reference lithosphere in Celsius
• temperature_base_lith_ref::Float64
  • Temperature at base of lithosphere for the reference lithosphere in Celsius
• adiabatic_gradient_ref::Float64
  • Adiabatic gradient below the reference lithosphere in K/km
• iuse_linear_segments::Int
  • Integer flag that controls the type of temperature profile used in the reference lithosphere model. If iuse_linear_segments = 1, a temperature profile with four linear segments is used that is controlled by the following user defined parameters: temperature_top_ref, temperature_moho_ref, and temperature_base_lith_ref. If iuse_linear_segments = 0, the temperature profile is calculated using an analytical 3-layer model defined during marker temperature. See the AnalyticalThreeLayer model description at MarkerTemperature). Default is 1.