toymodel | R Documentation |
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]
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)
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 |
toymodel()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.