plot multiple figures based on the output of the km_ode2() function. not done yet.

plot multiple figures based on the output of the km_ode2() function. not done yet.

plot_km(grdnt, parameter)

Arguments

grdnt	the output of the km_ode2() function
parameter	a list of the parameters

Value

figures

Examples

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 = 3)),  # reaction coefficient
 depth = 10, # m, measurement height difference
 layer = 200,
 flux = 10, # flux from the bottom layer, mol m-2 s-1
 t_max = 3.0 * 24.0 * 3600, # simulation time, seconds
 t_bin = 1000
))

# calculation -------------------------------------------------------------
grdnt <- km_ode2(parameter, savedata = c('RDS', 'csv'), path = 'data')
#> 2019-01-16 13:44:16 1e-05 0 ......
#> Error in checkInput(y, times, func, rtol, atol, jacfunc, tcrit, hmin,     hmax, hini, dllname): `y' must be numeric

# plot -------------------------------------------------
plot_km(grdnt, parameter)
#> Error: object 'grdnt' not found