param_b90$water_table_depth
). In this way, capillary rise from a water table can be simulated. Default is -9999, meaning no groundwater influence (i.e. unit gradient flow at the bottom of the soil profile). Any depth can be specified, also a depth below the lowest soil layer. If the water table is within the soil profile, soil layers below the water table are saturated with water.layer_output
, and 'snowlq' (liquid water content of snow on the ground) and 'cc' (cold content of snow pack) to output
. The latter two are important initial variables for continuing a simulation.param_b90
). These can now be passed as initial variables to snowlq and cc, for simulation continuationBudgMon
and BudgAnn
using miscellaneous function process_output_LWFB90()
rgroper
) was set to 0, to prevent root growth. The parameter defines the period during which the final root length density distribution is reached. The user now needs to actively switch on root growth by choosing a value > 0. It as advised to carefully investigate layer-wise root water uptake (tran
) when using the root growth module for young stands.Basic output is now at the precipitation interval level:
options_b90$prec_interval
> 1, the output of run_LWFB90()
now contains one row for each day
and precipitation interval of the simulation.options_b90$prec_interval = 1
).daily_output
of the return was renamed to output
.output
-argument and exported the function process_outputs_LWFB90()
that does the job.output
: sthr
.slvp
from layer_output
(soil evaporation).psiti
from layer_output
(total potential of soil layer): can be calculated by summing the layers' matrix potententials (psimi
) and the gravity potentials at the soil layers' midpoints.run_multisite_LWFB90()
(Thanks to Henrik Bengtsson)slrad
), net solar (solnet
), net longwave (lngnet
) and net radiation above (aa
) and below canopy (asubs
).make_rootden()
with method = 'betamodel'
led to negative
root density of the top layer when specified maxrootdepth
was lower than soil depth.hydpar_wessolek_tab()
caused NAs in return for texture = 'fSms'
soil_to_param()
: some unique soil materials were omitted in the return, if they had multiple occurences in the input soil
-data.frame.run_LWFB90()
:
options_b90$root_method = 'table'
gave an error in make_rootden()
because of a wrong argument name.verbose = T
.run_LWFB90()
(before runLWFB90()
):
output
, two tables
including all available output variables are returned: general daily outputs
and layer outputs. A proper description of model output variables was added to the help pages.output
) are now in lower case.climate
-argument, instead of a data.frame
.run_multisite_LWFB90()
(before msiterunLWFB90()
):
param_b90
(before param.b90
) input parameter objects for
individual climate/soil combinations (i.e. individual locations). The option to
provide a list of options_b90
(before options.b90
) input objects was disabled.climate
-data.frames for a multisite-simulation, a function
can be provided for on-the-fly creation of climate
-data.frames
. Arguments
for the function have to specified via the new climate_args
-argument of run_multisite_LWFB90()
.climate
, soil
and param_b90
objects are automatically
passed from run_multisite_LWFB90()
to run_LWFB90()
and thus become available to output_fun
.?ptfs
.run_multisite_LWFB90()
and run_multi_LWFB90()
from superseded packages 'snow' and 'doSNOW' to 'future', 'doFuture' and 'progressr' for parallel computation and progress reporting thereof. Pacifies a check note and is more future-proof (thanks @rnuske).MakeRelRootDens()
: The root depth distribution provided in a table is now redistributed to the soil nodes under preservation of the total root mass.MakeRelRootDens()
. See ?MakeRelRootDens
.relawat
values in MISCDAY.ASC
-output item.MakeRelRootDens()
returned the increment of the cumulative root proportion at the soil nodes, instead of the relative root density, which actually is the former value divided by the layer thickness.msiterunLWFB90()
runLWFB90()
: model input (param.b90
, options.b90
, standprop_daily
) is appended
to the return value before evaluating the output_fun
-argument. In this way,
on-the-fly post-processing of model results including model-input is possible now.output_fun
argument replaces gof_fun
-argument in runLWFB90()
for more flexible output.