Program and Source
The program and source code may be updated if critical bugs have been found and are fixed. The version written to screen by the program may still show '1022'. If possible, sign in to the Wiki and 'watch' this page to receive an email when this content is updated.
These issues are corrected in the most recent program and source file links.
The LQWS, FRWS, and ALWS fields for individual soil layers are printed to the water balance files in CSV format in an order that does not match the order of the fields in the file header. The order of these fields in the header is (('LQWS', 'FRWS', 'ALWS'), j = 1, NSL), whereas the order of the values printed is ((LQWS, j = 1, NSL), (FRWS, j = 1, NSL), (ALWS, j = 1, NSL)). A bug-fix in r1024 corrects this, which prints these fields in the same order as in the file header. The values of these fields are unaffected by this bug.
The fields for total LQWS, FRWS, and ALWS (of the entire soil profile, and which contain no integer label in the header), STG, and DSTG are correctly labelled and unaffected by this bug.
Bug-fixes to the in-line routing and time-stepping
- In previous versions of MESH, a counter was incremented before a call to the in-line routing. As a consequence, the output of precipitation, evaporation, and runoff at each day to the screen was offset by a day, relative to the observed and simulated streamflow values.
- In all prior versions of MESH, a variable type-set bug in the in-line routing caused channel capacity to be determined using incorrect parameters. Moreover, in MESH 1.3, an incorrect coefficient was determined for the overbank roughness term, R1. The implications of these bug-fixes are minimal in well calibrated setups, and these bug-fixes may in fact improve the fitness of the model. Poorly calibrated setups may see variations in simulated streamflow.
- Continued code clean-up and optimization have improved the general run-time of MESH by up to 20% and reduced memory usage by over 50%. MESH 1.4 also supports the MPI (message passing interface) protocol for parallelized runs, which can further improve run-time by over 80%, depending on the number of grid cells and GRUs in the setup and the number of processes utilized.
- Optional control flags have been introduced to minimize write operations to disk, to further reduce the general run-time of MESH. Further, the default output options have been reverted to write only daily basin total and streamflow output files. Additional outputs, including CLASS output and general outputs at various time intervals, can still be saved by enabling control flags.
- MESH 1.4 introduces LATFLOW with IWF option 3, which combines the PDMROF (overland ponding model) and WATROF (sloped soil column interflow) algorithms. Additionally, PDMROF, WATROF, and "flat-CLASS" proper can still be run with IWF options 2, 1, and 0, respectively.
The formatting of several CSV format output files, as well as many input files, has changed.
- All input files, with the exception of MESH_input_reservoir.txt are now read in free format. All fields in a single line should be separated by one or more spaces. The length and precision of the fields is no longer fixed.
MESH_input_streamflow.txt: Components of the gauge coordinate should now be separated by one or more spaces.
MESH_parameters_CLASS.ini: Individual fields should now be separated by one or more spaces.
Input forcing file time-stepping
Input forcing files are no longer required to use the same time-stepping. The HOURLYFLAG option should be set before all BASIN*FLAG's; however, the 'hf=' option on any BASIN*FLAG can be used to override the HOURLYFLAG value. If HOURLYFLAG is not defined, the default stepping is half-hourly (e.g., HOURLYFLAG 30).
Updated Output control flags
A number of new control flags can be added to MESH_input_run_options.ini to enable or disable saving certain output files at various time intervals.