Page contents:
Program and Source
- Source with makefiles (approx. 1.25 MB)
Known Issues
There are no known issues with this version of MESH.
General Updates
MESH 1.4.1149 adds new functionality to the new inline option for routing (RTE). The option uses the same physics as in the standalone WATROUTE (RPN_watroute) code, but is not completely implemented and does not include streamflow insertion (from the streamflow input file), resume functionality, temporally-driven parameterization, or diversions.
WF_R2
, WF_A1
, WF_A2
, WF_A3
, and WF_A4
parameters for WF_ROUTE do not transfer.Migrating to MESH 1.4
See the migration notes for the MESH 1.4.1022. All the same formatting changes are required for this version, in addition to the following.
XSLP
and DD
(or DDEN
) parameters from the MESH_parameters.r2c input file directly (activated with the INPUTPARAMSFORMFLAG control flag).Components of the reservoir output location must now be separated by one or more spaces in MESH_input_reservoir.txt.
Previously:
2940-5185
Now required:
2940 -5185
Revision History
September 22, 2017
General Updates
- Added printing diagnostic information about reservoir outlet locations, including
RANK
, if DIAGNOSEMODE is enabled - Added "point mode" with option 2 on SHDFILEFLAG to run single column configurations of the land surface scheme (more details here)
- Added checks to flag where the derived X/Y location is missing for active grids in the domain (where a
RANK
has no association to an active grid in the stride of 1 toNA
, which can occur from manual creation of modification of the drainage database file) - Disabled index of "water tiles" where the mosaic ID
MID
is less than zero (unused and incompatible with the revised indexing for MPI parallel processing) - Disabled parameter checks using the minmax_parameters.txt input file as the code no longer works as originally intended since the modular reorganization of the code, addition of new process modules, removal of hard-coded maximum limits for domain size and GRUs, and restructuring of internal variables in the driver
- Deactivated LOCATIONFLAG (location components are read with decimal precision from MESH_input_streamflow.txt and MESH_input_reservoir.txt regardless, since removing the fixed-formatting requirement)
- Modified the allocation of driver variables to only allocate variables associated with active process modules (e.g., only CLASS variables when CLASS is active, and only SVS variables when SVS is active)
- Removed the association of
FLZ
andPWR
with RTE to ensure they are allocated if RTE is disabled (WF_LZS is a standalone process that is activated independently to RTE) - Updated diagnostic and error messages throughout
- Added support to skip records in the reservoir input file (for controlled reservoirs/streamflow insertion) if using the tb0 format file (MESH_input_reservoir.tb0); it is not possible to skip records in the text format file (MESH_input_reservoir.txt), as the file contains no fields to specify the date corresponding to the first record
- Updated how grid and basin area are accumulated for basin totals, which can cause minor rounding differences in accumulated variables (particularly evident in large basins)
- Updated the version to "1149"
Climate Forcing Files
- Added the BASINFORCINGFLAG control flag with support for the "met" option to read the CLASS MET format forcing mode (e.g., for "point mode")
Inputs and Outputs (general)
- Added STREAMFLOWFILEFLAG with the "tb0" and "txt" options to specify the format of the streamflow input file (the default remains "txt")
- Added support for the MESH_input_streamflow.tb0 input file (can be used in place of MESH_input_streamflow.txt)
- Added printing diagnostic information about streamflow gauge locations, including derived drainage area
DA
andRANK
, if DIAGNOSEMODE is enabled - Added RESERVOIRFILEFLAG with the "tb0" and "txt" options to specify the format of the reservoir input file (the default remains "txt")
- Added support for the MESH_input_reservoir.tb0 input file (can be used in place of MESH_input_reservoir.txt)
- Removed the fixed-formatting requirement for reading values in MESH_input_reservoir.txt
- Activated the "csv" option for INPUTPARAMSFORMFLAG to support the MESH_parameters.csv input file
- Added support for
PWR
andFLZ
("BASEFLOWFLAG wf_lzs") andR2N
,R1N
,MNDR
,AA2
,AA3
,AA4
, andWIDEP
(RTE) fields in MESH_parameters.csv - Added
IWF
as a GRU-dependent parameter (to be assigned by GRU) in MESH_parameters_hydrology.ini (version 2.0 and higher) - Added
FLZ
andPWR
as GRU-dependent parameters (to be assigned by GRU) in MESH_parameters_hydrology.ini (version 2.0 and higher) for WF_LZS - Added support for
FCAN
andLNZ0
(both requiring the "nl", "bl", "cr", "gr", or "ur" specifiers),SDEP
,XSLP
,DD
(orDDEN
),SAND
(optionally with a layer specifier),CLAY
(optionally with a layer specifier), andIWF
in MESH_parameters.r2c - Added columns for
LZS
(used if using BASEFLOWFLAG) andDZS
(not used), and included these components inSTG
, in the Basin_average_water_balance.csv output file - Changed the
NO_DATA
value printed to screen and file (when no more records exist in the streamflow input file) from -999.0 to -1.0
Prairie Blowing Snow Model (PBSM)
- Re-enabled the Prairie Blowing Snow Model (PBSM) after updating the code to work with the tile-based indexing
- Added support for
fetch
,ht
,n_s
,a_s
, anddistrib
in MESH_parameters.r2c (the option to parameterize these as GRU-dependent parameters is still possible with MESH_parameters_hydrology.ini)
CLASS
- Reactivated the
DEGLAT
andDEGLON
fields from the CLASS input file if "point mode" is active (still ignored if using a drainage database instead) - Removed
IWF
as a constant and replaced it as a tile-based vector so different tiles can use different IWF routines (updated the WATROF/WATDRN, PDMROF, LATFLOW, ICEBAL, CLASSW, CLASSA and related code throughout) - Removed the experimental WATDRN3 code
SVS
- Updated the SVS code
- Added support to map properties for more than 3 soil layers in SVS
- Fixed passing
EVAP
andROFS
to the driver from SVS
RTE
- Updated the RPN_watroute (RTE) code
- Activating reservoir routing in RTE
- Implemented the same "avr_qo2" averaging as found in RPN_watroute in RTE
- Added the "cap_shd" option to the "rte" option flag in MESH_parameters_hydrology.ini to use the bankfull capacity read from the drainage database; without this option, channel geometry is determined using the inline functions
WF_ROUTE
- Revised the reservoir routing to set values to zero, instead of to arbitrary values greater than zero
- Added the "cap_shd" option to the "wf_route" option flag in MESH_parameters_hydrology.ini to use the bankfull capacity read from the drainage database (internally set as the product of this value and channel length, as previously in the code); without this option, channel geometry is determined using
WF_A2
,WF_A3
, andWF_A4
as before - Added the "rl_shd" option to the "wf_route" option flag in MESH_parameters_hydrology.ini to use the channel length read from the drainage database; without this option, channel geometry is determined using
WF_A1
as before
BASEFLOWFLAG
- Added word options for BASEFLOWFLAG
- Modified WF_LZS to run at the grid-scale (like in WATROUTE); the "luo_2012" option still runs at the tile-scale
- Added
LZS
to the MPI exchange and updated the BASEFLOWFLAG routines to be compatible with MPI parallel processing
Standalone WATROUTE (RPN_watroute)
- Updated the RPN_watroute (RTE) code