para.RdThe input parameter must have the necessary items shown in the example.
para(parameter)
| parameter | list of the input parameters |
|---|
list of the calculated parameters
parameter <- para(list( D_vector = c(1e-5, 1e-2), # m2 s-1, diffusion coefficient in Fick's law k2_vector = c(0, 10 ^ seq(log10(1.0 / (3600.0 * 24.0 * 30.0)), log10(1 / 671.1), len = 2)), # reaction coefficient depth = 10, # m, measurement height difference layer = 20, # how many layers in depth flux = 10, # flux from the bottom layer, mol m-2 s-1 c0 = 0, # initial ambient concentration t_max = 3.0 * 24.0 * 3600, # simulation time, seconds t_bin = 20 # time bins ))