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.