NEWS.md
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.