Lisflood OS - LISFLOOD
This release adds the computation and reporting of total water storage maps.
use of numba instead of cython for parallel computations in the routing module: numba was already used in other parallelized functions in soil modules. Now this version extends the use of numba to the kinematic wave parallel routing module. The change allows to completely exclude the use of Cython in the whole Lisflood OS software with no performance loss and skipping further compilation steps and potential run issues when porting kinematic wave compiled library to other systems/hardware. Numba use a Just in time (JIT) compilation that compiles python routines at runtime, generating optimized code and allowing parallelization by skipping the python Global Interpreter Lock (GIL) of specific parallelizable computation functions.
xml files and documentation files are now updated removing the variables used by kinematic wave parallel routing module and skipping the steps for the compilation of Cython routing functions. Numba parameters in xml files should be used instead of the old kinematic wave ones to allow or disable parallelization on multiple CPUs hardware systems.
This release includes the following fixes:
Small fix to use the netcdf template to set the coordinates in the netcdf outputs rather than using a pcraster clone. This avoids inconsistencies between the netcdf inputs and outputs.
This release features the following changes:
This release features the following changes:
The volume of lakes and reservoirs is correctly updated when water is abstracted for anthropogenic use.
The initialization run (or prerun; activated by switching on the "InitLisflood" option) implements a consistent solution throughout the code for the water abstraction module. As in v.4.0.1, the warm start of the initialization run is fully functional.
A hot fix to prevent the warm start (i.e. restarting from a previous state) to occasionally fail due to small approximation errors (order of magnitude: 10exp-12) when writing and reading the state files.
The changes (1) and (2) of this release have impacts on the results of LISFLOOD-OS when the following three optional modules are simultaneously activated: reservoirs, lakes, water use. The magnitude of the impact of the above changes mainly depends on the magnitude of the water demand for energy use, and on the extent of the rice fields within the modelled area.
The changes of this release allow the warm start of the initialization run (or prerun; <setoption choice="1" name="InitLisflood"/>
).
More specifically, the changes of this release allow to print the end step of all the state files required to correctly perform the warm start of the initilization run.
Before this release, the warm start could be correctly performed only for the run (<setoption choice="0" name="InitLisflood"/>
).
Furthermore, a unit test for waterbalance has been added.
LISFLOOD-OS v4.0.0 includes the following changes:
Moreover, to ensure compatibility with LISFLOOD-EPIC, modules source code is compatible with XArray variables, and XArray variables will be allocated when EPIC module is ON.
LISFLOOD-EPIC has been developed by JRC-D2. When compared to LISFLOOD v3.2.0, LISFLOOD-EPIC enables a more detailed representation of crop growth and irrigation. References: ASR - Assessing groundwater irrigation sustainability in the Euro-Mediterranean region with an integrated agro-hydrologic model (copernicus.org) , Gelati, E. et al.: Integrated model of irrigation, crop growth and water resources (LISFLOOD-EPIC): description and application in the Euro-Mediterranean area, in preparation, 2020.
Important note: EPIC modules are not yet available from the OS-LISFLOOD repository.
computation of water infiltration into the soil within the soilloop.py module.
LISFLOOD v3.2.0: all the pixels of a land use fraction are allocated to a vector and the number of loops is defined by the most critical pixel. More specifically, all the pixels within the computational domain are allocated the number of iterations that are necessary to allow the numerical stability of the most critical pixel. The number of iterations consequently depends on the specific computational domain. This can lead to different results when modelling sub-catchmennts and entire basins with highly heterogenous soils.
LISFLOOD v4.0.0: pixel by pixel computation and parallelization of the computations by using the python package numba.
Benefits of v4.0.0: each pixel is allocated the number of iterations that allow its numerical stability, this guarantees the consistency of results when modelling sub-catchments and the entire basin with highly heterogeneous soils. Decreased computational time due to parallel computations.
changes in the allocation of water abstraction.
LISFLOOD v3.2.0: direct abstraction of the consumptive water use.
LISFLOOD v4.0.0: abstraction of the demanded water volume for each use (industrial, domestic, energy, livestock); the consumptive water use leaves the system; the unused water volume is returned to the channels.
Minor note: small differences in the computation of the leakages.
patches in the modules riceirrigation.py and frost.py to avoid non-realistic results in specific conditions.
riceirrigation.py: LISFLOOD v3.2.0: in presence of a very thick third soil layer (10^2 m) the drainage of the rice fields causes non realistic large flow discharge values in the channels.
LISFLOOD v4.0.0: only soil layers 1a and 1b are used by the rice module (for both the saturation and the drainage phases). By definition, the second soil layer includes the roots depth.
frost.py: LISFLOOD v3.2.0: there is no upper boundary to the frost index value. Rain events on deeply frozen soils are totally converted into runoff and can cause false alarms.
LISFLOOD v4.0.0: the maximum frost index value is set to 57.
improvements in the computation of the mass balance. LISFLOOD v3.2.0: the mass balance is computed by considering the input and output volumes up to the current computational step. Small numerical errors accumulate over time and give the erroneous impression of large mass balance errors. LISFLOOD v4.0.0: computation of the mass balance by considering the single computational step.
bug fixes in the use of the transientlandusechange option, in the use of the option simulatePF, in the writing of the monthly outputs.
Correct use of the rice fraction by landusechange.py DYNAMIC; use of the correct soil water content variable for the computation of pF1 (ws1a replaced ws1); correct reporting of the monthly time steps by the optional module indicatorcalc.py.
Lisflood 4.0.0 is the version used for EFAS 5 and GloFAS 4 calibration.
This release includes the following changes:
All of the changes above improve the code capability to handle input and outputs. The modelling of the hydrological processes has not been changed.
IMPORTANT NOTE: the results of the unit tests of this release are different from the results of the unit tests of the release 3.1.1. The differences are due exclusively to a different use of of the optional moduels (e.g. groundwatersmooth) within the .xml setttings file of the unit tests. When using the same settings of .xml file, and the same dataset, v3.1.1 and v3.2.0 provide the same results.