Software Open Access
Joe Hamman; Bart Nijssen; Ted Bohn; Diana Gergel; Wietse Franssen; Yixin Mao; homefc; Hordur Helgason; Keita Suzuki; Min RK; The Gitter Badger; Tony Craig; Jake Zwart; kiritokun07; Bin Peng
This is a maintanence release that includes numerous bug fixes and a few model enhancements. All code changes were included in the two release candidates (5.1.0.rc1 and 5.1.0.rc2). See the notes below from the two release candidates for more details.Announcements
Fixed datetime handling bug in unit test module (GH#796)
Removed descriptions of forcing disaggregation options from documentation of classic driver global parameter file (GH#831)
Removed descriptions of forcing disaggregation options from documentation of image driver global parameter file (GH#833)
Fixed segmentation fault in lake model caused by incorrect passing of pointer argument to
Relaxed some of the validation of lake parameters, to allow them to be null in cells containing no lakes. (GH#826)
Fixed passing of lake data structures to
Separated the dimensions of arrays related to lake basin shape and of arrays related to the number of lake simulation layers. (GH#826)
Added global parameter option to set the maximum number of (dynamic) lake simulation layers. (GH#826)
Fixed bug causing ET to be 0 in pure bare soil tiles. (GH#823)
Release date: (April 27, 2018)
Source code is available here:
This is a minor update from VIC 5.0.1. The VIC 5.1.0 includes new features, such as a new streamflow routing extension and extended parallelization using OpenMP. The release also includes a number of bug fixes for the CESM driver. See the VIC Github page for more details on the changes included in this release.Model enhancement:
Improved calculation of drainage between soil layers (GH#656)
Drainage from upper layer to adjacent lower layer is calculated according to Brook & Corey curve (where drainage rate is a function of upper-layer soil moisture). In previous versions, a simple numerical solution is applied which uses the timestep-beginning upper-layer soil moisture to calculate drainage rate, and assume this constant rate over the entire timestep. This can cause unreasonably large drainage if the curve has a steep shape and when soil moisture is high. Now, the current version uses exact integral (instead of numerical solution) for layer drainage calculation.
Fixes for the CESM driver
advance_vic_time, respectively, to avoid conflicts with WRF functions with the same names when compiling RFR case.
l2x_Sl_ram1field with aero_resist to agree with the VIC 4 implementation in RASM.
TRUE, and sets defaults for daily and monthly mean output.
GH#880 --Adds option STATENAME_CESM to the option_struct so that statefiles in the CESM driver can be in accordance with RASM naming conventions for model components.
Speed up NetCDF operations in the image/CESM drivers (GH#684)
These changes speed up image driver initialization, forcing reads, and history writes by only opening and closing each input netCDF file once.
Added two new timers to measure time in I/O operations (GH#703)
These two timers count the CPU and WALL time spent in
vic_write. The accumulated time from these timers is printed out at the end of each simulation in the timing table. See also GH#442.
Added gridcell-averaged albedo (STATE_AVG_ALBEDO) as a state file variable (GH#712)
This is for use in the CESM driver for VIC to pass to WRF, but has been implemented in the core structure of VIC (in vic_run) for consistency with the classic and image drivers. Running VIC from a cold start now also includes calculation of gridcell-averaged albedo.
Cleanup of the initialization sections of the
cesm drivers (GH#701)
Codified behavior in the initialization of the
cesm drivers that requires the parameter variables
Cv must sum exactly to 1.0. If using the
SNOW_BAND option, the area weighted
elevation must match the mean grid cell elevation (
elev). VIC will print warnings if any of these criteria are violated.
Added thread parallelization using OPENMP (GH#712)
The VIC image and CESM drivers now may be optionally compiled with OPENMP to enable shared memory thread parallelization. This option should improve the parallel scaling of these drivers by reducing the number of MPI messages and increasing message size.
Added streamflow routing extensions ROUT_STUB and ROUT_RVIC for the VIC image driver (GH#231)
The VIC image driver can be optionally compiled with ROUT_RVIC to enable routing in image mode (ROUT_STUB is the default extension which means no routing). With ROUT_RVIC enabled, the output variable
OUT_DISCHARGE is available, and there will also be an extra state variable
STATE_ROUT_RING stored in the state file.
Moved MAX_ITER_GRND_CANOPY, which controls the maximum number of ground-canopy iterations in CLOSE_ENERGY mode for vegetation types with an overstory, to the parameters struct (GH#771)
Previously this was set in the surface_fluxes.c numerics routine for ground-canopy iterations, which meant that that routine had to be altered to change the maximum number of iterations. It has now been moved to the parameters struct so that it can be overriden in the constants file.
Updated new snow density function by adding a cap to new snow density that is set in the parameters struct by the parameter SNOW_NEW_SNOW_DENS_MAX (GH#776)
Previously the change in cold content of the snowpack term (deltaCC in the snow_data_struct) would get unreasonably large if the Hedstrom and Pomeroy 1998 equation used to calculate snow density, which depends only on air temperature, was calculated with air temperatures above about 2 deg C. We use this term to calculate the ground flux from the snowpack and snow depth, which resulted in extremely small snow depths and unreasonably large ground fluxes from the snowpack (and thus changes in snowpack cold content). Now there is a cap on new snow density with the new parameter SNOW_NEW_SNOW_DENS_MAX as well as a snow depth below which we disregard the ground flux from the snowpack (1.e-8).
Added new option BULK_DENSITY_COMB that enables soil bulk density (mineral and organic) to be read from the parameters file (GH#817)
The option BULK_DENSITY_COMB enables soil bulk density (mineral and organic) to be read in as a parameter when the option is set to true in the global parameter file. Default is false.
ORGANIC_FRACT option that previously existed in VIC 4 and enables the
BULK_DENSITY_COMB option to be used in conjunction with it (GH#837)
ORGANIC_FRACT, if set to True, means that the organic fraction of the soil, soil density of the organic matter, and bulk density of the organic matter (if
BULK_DENSITY_COMB is set to false) will be read in from the parameter file. If
BULK_DENSITY_COMB is set to True, the bulk density of the organic matter will not be read in separately. Default is false.
Added new option
MAX_SNOW_ALBEDO that enables new snow albedo to be read in from the parameters file (GH#835)
MAX_SNOW_ALBEDO, if set to true in the global parameter file, means that new snow albedo will be read in from the parameter file and used in the snow routines for all vegetation types except for bare soil.
|All versions||This version|
|Data volume||1.2 GB||119.4 MB|