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

Initialize transient bottom temperature boundary conditions by performing the following calculations:

1. Calculate the distance from the base of the lithosphere to the bottom of the model.
2. Calculate the final cooler temperature at the bottom of the model temperature_bottom_cooler_final using the temperature at the base of the lithosphere temperature_base_lith, the adiabatic gradient adiabatic_gradient, and the distance from the base of the lithosphere to the bottom of the model.
3. Calculate the initial elevated temperature at the base of the lithosphere temperature_base_lith_warmer_initial using the temperature at the base of the temperature_base_lith' and the elevated temperature anomalydeltatemperaturetransient`.
4. Calculate the initial elevated temperature at the bottom of the model temperature_bottom_warmer_initial using the initial elevated temperature at the base of the lithosphere temperature_base_lith_warmer_initial, the distance from the base of the lithosphere to the bottom of the model, and the adiabatic gradient adiabatic_gradient.

The transient temperature will be the final cooler temperature at the bottom of the model temperature_bottom_cooler_final which will be applied between the start and end times specified by the user. This is used to approximate the effects of a plume that rapidly ascends through the lithosphere, spreads out at the base of the lithosphere, and cools over time.

Arguments

• model::ModelData: Model data object containing model parameters and arrays.

Optional Keyword Arguments

• iuse_bottom_transient::Int64:
  • Flag to use transient temperature at the bottom boundary
• delta_temperature_transient::Float64:
  • Temperature perturbation for initial transient temperature at the bottom boundary
• temperature_base_lith::Float64:
  • Initial temperature at the base of the lithosphere in Kelvin
  • For cases where an transient bottom temperatures this parameter is the final cooler temperature at the base of the lithosphere. If an AnalyticalThreeLayer model is used, the temperature at the base of the lithosphere (temperature_base_lith) is updated to the warmer value when calculating the initial geotherm.
• start_time_bottom_transient::Float64:
  • Start time (Myr) for transient temperature at the bottom boundary
• end_time_bottom_transient::Float64:
  • End time (Myr) for transient temperature at the bottom boundary
• adiabatic_gradient::Float64:
  • Adiabatic gradient in K/km