API

# LadeGCM.RCP3PD — Constant.

Alias for the RCP3PDPathway constructor.

source

# LadeGCM.RCP45 — Constant.

Alias for the RCP45Pathway constructor.

source

# LadeGCM.RCP6 — Constant.

Alias for the RCP6Pathway constructor.

source

# LadeGCM.RCP85 — Constant.

Alias for the RCP85Pathway constructor.

source

# LadeGCM.calculate — Method.

calculate(rcp; [c=constants(), tspan=(1765., 2100.), solve_args...])

Entry point to the module. Will construct all required parameters and solve the DAE.

The tspan here is set to mimic the Lade et al. paper, but the emission parameters in each of the RCP files project out to 2500; meaning the 1765–2500 range is valid here.

Any flags you wish to pass to the solver can be added to the end of the call.

Examples

results_6 = calculate(RCP6);

results_45_grid = calculate(RCP45, reltol=1e-10, abstol=1e-10, saveat=1765:0.1:2100);

source

# LadeGCM.constants — Method.

constants(<keyword arguments>)

A constructor to generate all required constants for the model, any of which can be overridden via a keyword argument.

Arguments

cₐ₀::Float64=589.0: PgC, Pre-industrial atmospheric carbon
cₜ₀::Float64=1875.0: PgC, Pre-industrial soil and vegetation carbon
cₘ₀::Float64=900.0: PgC, Pre-industrial ocean mixed layer carbon
λ::Float64=1.8: K, Climate sensitivity (TCR)
τ::Float64=4.0: yr, Climate lag
D::Float64=1.0: yr⁻¹, Atmosphere–ocean mixed layer CO₂ equilibration rate
r::Float64=12.5: Revelle (buffer) factor
DT::Float64=0.0423: 4.23 % K⁻¹, Solubility temperature effect
B₀::Float64=13.0: PgC yr⁻¹, Pre-industrial biological pump
BT::Float64=0.032: 3.2 % K⁻¹, Temperature dependence of biological pump
w₀::Float64=0.1: yr⁻¹, Solubility pump rate
wT::Float64=0.1: 10 % K⁻¹, Weakening of overturning circulation with climate change
QR::Float64=1.72: Terrestrial respiration temperature dependence
NPP₀::Float64=55.0: PgC yr⁻¹, Pre-industrial NPP
KC::Float64=0.3: Fertilisation effect

Examples

Obtaining the default set of constants

julia> constants()
LadeGCM.Constants(589.0, 1875.0, 900.0, 1.8, 4.0, 1.0, 12.5, 0.0423, 13.0, 0.032, 0.1, 0.1, 1.72, 55.0, 0.3)

Override $c_{a0}$ and $KC$ with custom values.

julia> constants(cₐ₀=600.0, KC=0.5)
LadeGCM.Constants(600.0, 1875.0, 900.0, 1.8, 4.0, 1.0, 12.5, 0.0423, 13.0, 0.032, 0.1, 0.1, 1.72, 55.0, 0.5)

source

# LadeGCM.Constants — Type.

Houses all constants for the model. See: constants.

source

# LadeGCM.Pathway — Type.

Representation of a Representative Concentration Pathway.

source

# LadeGCM.RCP3PDPathway — Type.

RCP3PD Representation.

source

# LadeGCM.RCP45Pathway — Type.

RCP45 Representation.

source

# LadeGCM.RCP6Pathway — Type.

RCP6 Representation.

source

# LadeGCM.RCP85Pathway — Type.

RCP85 Representation.

source

# LadeGCM.Results — Type.

Storage for all model output. See results.

source

# LadeGCM.generate_emission_parameters — Method.

(E, LUC) = generate_emission_parameters(rcp)

For a given concentration pathway, generate continuous functions for fossil fuel emissions E(t) and land use emissions LUC(t).

Data comes from files on disk and is linearly interpolated to provide the continuous output.

source

# LadeGCM.load_pathway_data — Method.

load_pathway_data(rcp)

Loads csv file for a given RCP scenario into a dataframe for processing. The information in these files comes from Meinshausen et al. (2011), which were generated using MAGICC6.

source

# LadeGCM.results — Method.

results(sol, c)

Collects all required outputs from the DAE solution.

source