toymodel: Toy model

View source: R/toymodel.R

toymodelR Documentation

Toy model

Description

The function solves the toy model equation for the emission heigth making use of the energy balance equation.

\frac{z_e}{dt} \approx \frac{-1}{\gamma c_s} \left[ \frac{S_0}{4} (1-A) - \sigma (T_E - \gamma z_E)^4 + \int_0^H e^{-\tau z} \rho \sigma (T_E - \gamma (z_E - z))^4 dz - \frac{\rho 10^{(11.40-2353/[T_E - \gamma z_E])}}{2} \eta_z(t) \right]

Usage

toymodel(t0 = 1, t1 = 10, N = 1000, S0 = 1361, A = 0.3, gamma = -5, eta = 1+0.5*cos(seq(0,30,length=1000)), sigma = 5.67 * 10^-8, cs = 1000, tau = 10, rho = 0.1, plot = TRUE)

Arguments

t0

start time

t1

end time

N

number of steps

S0

The 'solar constant' (W/m**2)

A

The albedo

gamma

The lapse rate (C/km)

eta

The overturning rate

sigma
cs

specific heat of the ground

tau

optic depth

rho

density of the air at surface

plot

TRUE - plot the results

Examples

toymodel()

brasmus/GHE documentation built on Jan. 3, 2024, 6:03 p.m.