Release notes - iMOD
Release Notes
Release notes - current
The latest PDF version of the iMOD release notes can be downloaded here:
Release notes - history
The release notes of iMOD 3.0, 3.1, 3.2 and 3.4 can be found as the last 2 chapters in the iMOD User Manual.
- Technical release notes iMOD-GUI, releases 3.0 - 3.4
- Technical release notes iMODFLOW, releases 3.0 - 3.4
For iMOD 3.4 both technical release notes (see above) and a descriptive summary of all new, changed, extended and fixes functionalities (see below) is available.
- Release notes iMOD 3.4:
Starting from iMOD 3.6 we provide a descriptive summary of the release notes only, see below.
---------------------------------------------------------------------------------------------------------------------------------------------------
The main reason of this release is the new GeoConnect Tool functionality. iMOD 3.4 also contains improved features and bugfixes.
Below please find:
- A summary of new and improved features, and a summary of most important bugfixes.
- A summary of current developments.
Ad 1) Summary of new and improved features and bugfixes in iMOD 3.4.
In the Graphical User Interface of iMOD (iMOD-GUI):
- The new GeoConnect Tool:
- In many cases groundwater model layer parameters are based on assembled (e.g. summed) detailed geological aquifer- and aquitard-properties (elevation, thickness, geometry (extent), hydraulic conductivity, etc.). The GeoConnect Tool allows the modeller to define and utilize permanent links between 1) the (unassembled) geologic layers (incl. its properties) and 2) the aggregated model layers. With the GeoConnect Tool you can re-calculate the hydraulic conductivities of a model layer after adapting the individual weights of each contributing geological layer. For more info, see section 7.8 of the iMOD User Manual.
- Redesign of the 3D Tool:
- The 3D Tool window and the 3D-Tool menu-options have been merged to one window that optionally synchronizes with the zooming in the 2D-window.
- Output of Pathline calculations:
- The number of the deepest model layer the pathline passes through is now recorded (as in the past).
- Contouring and legends:
- Improved contour algorithm for delineation of flat areas.
- A legend-label now uses a proper number of decimals.
- Legends can now be generated using a user-specified interval.
- Solid Tool:
- During interpolation, separate IPF-files can now be included, with associated txt-files that specify the elevation of each interface.
- IMF-files:
- In IMF-files, relative directory names can now be used.
- Thanks to the iMOD-forum discussions, several bugfixes were identified and solved:
- Assignment of the appropriate output filename when starting iMODPATH.
- PST-module now mentioned correctly in the runfile.
- The Water Balance Tool now recognizes filenames containing *sys*.
In iMODFLOW:
- As was the case in iMOD 3.3, iMOD 3.4 includes MetaSWAP SVN-version 1047. Compared with the previous version of MetaSWAP only some in- and output statements for post-processing purposes have been improved.
- Automatic assignment to layers of faults in HFB (Horizontal Flow Barrier) package:
- When using GEN-files as a representation of the location of faults, it was already possible to have horizontal flow barriers assigned automatically to the nearest cell faces resulting in interconnected series of HFB-cell faces. We have extended this functionality with automatic assignment to layers; it requires a GEN-file including elevation data (a so-called 3D-GEN-file) and TOP and BOT elevations of the model layers. For input instructions, see section 10.14 "Data Set 11: Input ?le assignment" ; for theory and examples see section 12.15 in the iMOD User Manual.
- New value for HNOFLOW cells:
- To avoid possible confusion when modelling deep systems the HNOFLOW value is now set to 1.43E+37 (the old value was -999).
- On-the-fly correction of input data:
- Starting from iMOD 3.0 the iMODFLOW computational core is based on MODFLOW2005. Although we prefer to keep the pre-processing and processing separate, and thus not deviate from the standard MODFLOW2005-code, we have implemented 2 on-the-fly corrections on input data:
- If the elevation of the riverbed is higher than the waterlevel, the riverbed elevation is set to the river waterlevel. We intend to migrate this on-the-fly correction to a to be developed iMOD batch function that allows for the mentioned correction and similar checks on input data as an optional pre-processing step.
- MODFLOW2005 cannot handle active cells with zero thickness. In iMODFLOW these cells can optionally be assigned a minimum transmissivity value; for input instructions see section 10.6 "Data Set 4: Simulation mode".
- Starting from iMOD 3.0 the iMODFLOW computational core is based on MODFLOW2005. Although we prefer to keep the pre-processing and processing separate, and thus not deviate from the standard MODFLOW2005-code, we have implemented 2 on-the-fly corrections on input data:
- Saving OLF output separately:
- OverLand Flow is now saved in separate budget files. For more info, see section 7.15 "Compute Waterbalance" in the iMOD User Manual.
- Bug fix in waterbalance post-processing regarding diagonal ANI-fluxes
- Diagonal fluxes between 'ANI-cells' and 'diagonal-adjacent' fixed head cells are now computed correctly and processed in the (cumulative) water balance.
- Bug fix in ANI, incorrect nodata values:
- When the ANI package was activated, this could cause errors near closed model boundaries (IBOUND=0) for anisotropy factors < 1. The reason for this is that for stability reasons small minimal horizontal conductivities (1E-07) were initially imposed by this package for all MODFLOW cells, causing a) arti?cial diagonal ?uxes being computed near closed boundaries from the small horizontal conductivities and large default HNOFLO (1.43E+37), and b) incorrect coef?cient matrix cleaning in the PCG solver due to a bug in this package (see below).
- Bug fix in PCG, cleaning matrix coef?cients:
- The cleaning of the matrix coef?cients was corrected. The cleaning routine of CC/CR/CV was a sequence-dependent method for closed model boundaries
(IBOUND=0). Also cleaning was done in the same loop as the initialization of the coef?cients matrix. In case of closed East and South boundaries, a certain initial residual was included in the calculation. The issue only occurred with models that included the ANI package, because with this package it is possible to initialize CC and CR with very small values.
- The cleaning of the matrix coef?cients was corrected. The cleaning routine of CC/CR/CV was a sequence-dependent method for closed model boundaries
- Bug fix regarding inactive package handling:
- In the runfile it is possible to temporarily inactive packages. Sometimes this option was erroneously interpreted as "re-use"; this is now corrected.
- Bug fix in SFT/ISG:
- When utilizing the StreamFlow Thickness package a correction regarding the calculated total length of the river reaches in a cell is implemented.
Current developments for upcoming releases:
- Additional runfile-options:
- We are implementing a number of additional runfile-options; for a specification of these runfile-options scheduled for a next iMOD-release, please have a look in chapter 10 of the iMOD User Manual.
- Implementation of the SFR2- and LAK-packages:
- SFR2- and LAK-package are being implemented in iMOD as part of ongoing projects and will become available in one of the upcoming iMOD-releases.
- A new Parallel Krylov Solver package (PKS):
- We are currently developing a new parallel solver package for MODFLOW called PKS (Parallel Krylov Solver). It is based on overlapping domain decomposition combining both the techniques of MPI and OpenMP. This development is done together with the USGS and Wageningen Environmental Research (Alterra) in cooperation with Utrecht University and Delft University of Technology. This year we are piloting this new PKS package in iMOD to test running large iMOD models on powerful multi-core computers allowing a drastic reduction in runtimes. We hope we will be able to make this new solver available in one of the upcoming iMOD-releases.
Other ongoing developments (status on the release date of iMOD 3.4):
- A new SUB-CR package (including 'creep'):
- The existing SUBsidence-package has been extended with 'creep' for clay and peat; the implemented concept is consistent with D-Settlement (geotechnical software). We hope we will be able to include this extension in an upcoming iMOD-release.
- Coupling of iMOD and DgFlow:
- Proof of Concept: iMODFLOW has been off-line coupled to DgFlow, a finite-element code developed by Deltares for simulating the consolidation process of dikes/embankments. An explicit head-flux boundary coupling has been realised, demonstrating the added value of incorporating regional, transient, groundwater flow boundary conditions for dike stability analysis. For more info, please contact Deltares.
- We also plan to make iMOD-MT3DMS/SEAWAT open source and include it in the iMOD-framework. This modified version of MT3DMS/SEAWAT uses an easy-to-use runfile for specifying model data for many layers and supports the common iMOD file-formats. We hope to finalize this open source process at the end of 2017, depending on the availability of resources.
- We are working on extending the list of iMOD-supported MODFLOW2005-packages, among others the MNW2-, UZF- and SWR1-packages; a fixed schedule of implementation/prioritization is not available yet.
---------------------------------------------------------------------------------------------------------------------------------------------------
This release contains the iMOD-implementation of five existing USGS MODFLOW-packages: the SFR-package (Surface Flow Routing), the LAK-package (Lake), the UZF-package (Unsaturated Zone Flow), the FHB-package (Flow and Head Boundary) and the MNW-package (Multi Node Well). Implementation of these packages was realized within the framework of foreign projects.
This iMOD 3.6 release also contains a new version of MetaSWAP which facilitates the use of an extensive database containing 370 soil profiles; the new version of MetaSWAP was made possible by the The Netherlands Hydrological Instrument - project.
iMOD 3.6 also contains some significant improvements of existing functionalities.
Below please find:
- A summary of new and improved features, and a summary of most important bug fixes.
- A summary of current developments.
New features:
- SFR-package:
The existing Surface Flow Routing package of MODFLOW2005 has been implemented in iMOD. With this package it is possible to route surface water through a network of segments. These segments are stored in an additional type of ISG file that contains more records than the conventional ISG. These additional records are specific for the SFR package and describe the relationship between individual segments and how the water depth in the segments needs to be computed. The water depth is a function of the computed hydraulic head, the SFR is therefore a significant improvement in the interaction between surface- and groundwater. - LAK-package:
The existing Lake package of MODFLOW2005 has been implemented in iMOD. With this package it is possible to compute lake stages in relation to the groundwater head. Given a bathymetry of the lake, iMOD will assign the lake to the appropriate model layers and assign the correct conductances to the lake. The SFR can be connected to locations where a lake is defined to route water in- and from a lake. - UZF-package:
The existing Unsaturated Zone Flow package of MODFLOW2005 has been implemented in iMOD. Percolation of precipitation through unsaturated zones is important for recharge of ground water. Rain and snow melt at land surface are partitioned into different pathways including runoff, infiltration, evapotranspiration, unsaturated-zone storage, and recharge. The package was developed to simulate water flow and storage in the unsaturated zone and to partition flow into evapotranspiration and recharge. The package also accounts for land surface runoff to streams and lakes. A kinematic wave approximation to Richards' equation is solved by the method of characteristics to simulate vertical unsaturated flow. The approach assumes that unsaturated flow occurs in response to gravity potential gradients only and ignores negative potential gradients (upward capillary flow); the approach further assumes uniform hydraulic properties in the unsaturated zone for each vertical column of model cells. - FHB-package:
The existing Flow and Head Boundary package of MODFLOW2005 was implemented in iMOD. With this package it is possible to combine constant head- and constant flux boundaries for all model layers independently. For each model cell it is possible to define whether it becomes a constant head or constant flux boundary condition. - MNW-package:
The existing Multi Node Well package of MODFLOW2005 was implemented in iMOD, among others to be able to couple with RIBASIM (Deltares-software for planning and analysis of water resources). The Multi-Node Well package is used to simulate "long" wells that are connected to more than one model layer; the abstraction rate is vertically distributed proportional to the transmissivity adjacent to the well screen; e.g. when a hydraulic head gradually drops below the top of a well screen the yield of this shallow part of the well will also gradually drop. This package also connects to overlying aquifers via the well screen; an example of this behaviour is given in the new MNW-tutorial in the iMOD User Manual. - In the iMOD-GUI:
- Configuration of the SFR-, LAK-, MNW- and UZF-packages has been implemented in the iMOD Project Manager. These packages are not used through by the so-called 'runfile' (*.RUN) but by an iMOD-project file (*.PRJ) only. Based on the iMOD-project file-definitions iMOD generates/exports the input files for iMODFLOW.
- An iMODBatch function had been added that allows converting the model output to ISG-format thus facilitating visualizing results vector-wise.
- We now also facilitate the use of unconfined layers.
- In iMODFLOW:
- An option to sum flux-output of sub-systems (e.g. summing bdgisg_sys1_*.idf + bdgisg_sys2_*.idf to bdgisg_sys_*.idf) has been added for RIV, DRN and ISG.
- The volumetric water balance term 'DXC FLUXES' (summed data exchange fluxes for coupled non-USGS-packages) can now be stored as separate idf-output (e.g. MetaSWAP-fluxes).
- The debug option IDEBUG has been reactivated.
- In MetaSWAP:
The MetaSWAP update is compatible with the new distributed format of the unsaturated zone database. The distributed format has a separate database for each of the soil types. This format has the advantage that new soil types can more easily be added, and the files are kept a manageable size. Extracting the zip archive does yield a large number of files in a directory tree; these files form part of the documentation of how the database was generated. In the para_sim.inp file the path to the root of the tree should be specified.
Compatibility with the old database format has technically been maintained. But one should be aware of the fact that exhaustive testing of all combinations with all old databases is not feasible. In the case of problems we will first advise to update to the latest database. Apart from containing the latest available soil information the database also generates storage coefficients that are not smaller than 0.01, which speeds up the running of the model.
The MetaSWAP code has been made more robust for situations with extremely sensitive relationships in the MODFLOW model. Since the model no longer crashes in such a situation, this makes it easier to diagnose where the sensitivity is located and thus address the issue.
The most recent database with 72 soil types can be downloaded here: ftp://ftp.wur.nl/simgro/tests/PreMetaSWAP_V7_3_13_BOFEK2012_Fits2_a/out/LHM2016_v01.zip, info on the database with 370 soil type can be found in the MetaSWAP-release notes.
Improvements, bug fixes and extension of existing functionalities in iMOD 3.6:
- In the iMOD-GUI:
- Mainly improvement of legend-features, e.g. to use legends for multiple maps simultaneously more easily.
- ISG Edit has been extended significantly. For figures containing water levels and cross-sections now axes have been added. For use in the SFR-context (SFR uses the extended ISG-file type) the set of SFR-definitions has been extended, e.g. interactively or automatically generating the definition of the connection of segments in a downstream direction.
- In iMODFLOW:
- The accuracy of the gridding of GEN-files for the definition of Horizontal Flow Barriers has been improved and is now similar to common GIS-gridding methods. Applying this improved gridding can cause significant differences in calculated heads in the vicinity of the Horizontal Flow Boundary.
- The ISG-gridding has been made more robust, avoiding a divide-by-zero in rare cases where the water level coincides with a buckling point elevation of a specified cross-section.
- A number of bugs in iPEST have been fixed, e.g. during conductance optimization of rivers and drains water levels were optimized instead, this is fixed now. Also the statistical information on the optimized parameter has been extended.
- Names of output files can now also include the starting date of a time step as opposed to the end date of a time step; in a PRJ-file this is now default, the RUN-file still used the 'old' method.
- In the MODFLOW-MetaSWAP-coupling:
- In the N:1 coupling of iMODFLOW and MetaSWAP the second MetaSWAP-column coupled to a MODFLOW-cell was not included correctly. These second MetaSWAP-column usually represents the paved/urban fraction of a cell area which results in a relatively small share of the water balance. This is corrected in the current version.
- The standard parameterization for urban fraction of a cell has been made more realistic including a surface elevation comparable with the rest of the cell and including a small threshold (2 mm) to allow runoff calculated by MetaSWAP. Deltares and Wageningen Environmental Research (Alterra) will communicate more recommendations regarding runoff- and evaporation parameter values for urban areas in the near future.
- In MetaSWAP:
- A 'MSWP-error: update rootzone - floating invalid'-issue seems to be fixed, by doing some dynamic array allocations more strictly than required by the Intel compiler.
- Known issues: sprinkling using groundwater from multiple model layers of the same cell simultaneously does not work properly. Also a water balance problem may occur when using groundwater for sprinkling purposes from a different location than where it is sprinkled. In case you encounter the latter situation, a temporary workaround is available via the iMOD help desk. For both issues we are working on a solution to be made available in a next release.
- For a the full list of release notes, please have a look at the MetaSWAP-release notes.
Ad 2) Developments (status on the release date of iMOD 3.6):
Current developments for upcoming releases:
- A new Parallel Krylov Solver package (PKS):
- We are currently developing a new parallel solver package for MODFLOW called PKS (Parallel Krylov Solver). It is based on overlapping domain decomposition combining both the techniques of MPI and OpenMP. This development is done together with the USGS and Wageningen Environmental Research (Alterra) in cooperation with Utrecht University and Delft University of Technology. This year we are piloting this new PKS package in iMOD to test running large iMOD models on powerful multi-core computers allowing a drastic reduction in runtimes.
- We made considerable progress; preliminary results have been presented during the 2016 AGU Fall Meeting (San Francisco, USA) and the MODFLOW and more 2017 conference (Golden, USA). See also Deltares R&D Highlights 2016, page 64. During the Dutch Deltares Software Dagen we organize a beta-workshop 'iMOD Workshop Parallel Krylov Solver package' (in Dutch); participants can experiment with a beta-version of the PKS-package on their own laptop. The use of PKS-package will most likely also be programmed during the iMOD International User Day (to be scheduled in the last week of October / first week of November 2017).
- We hope to make the PKS-package available in a next iMOD-release.
- Extension of the Waterbalance Tool:
- The current water balance tool is extended significantly to allow constructing water balances for sub-regions, and selecting / de-selecting MODFLOW- and MetaSWAP-water balance components interactively. The new water balance tool also allows specifying only a subset of the output time-series and contains several time-related aggregation possibilities (e.g. averages per year, month, season). Visualization includes interactive stack-bars time-series plots and schematic vertical cross-sectional overviews. This development is done together with Tauw and financed by the Dutch iMOD-CGO group.
- Movie tool:
- Also part of the above mentioned 'Water balance Tool extension'-project a movie tool is developed. It facilitates visualizing a time-series of maps (e.g. calculated groundwater levels) and record it as a playable file.
- 3D fault assignment tool:
- As part of a consortium-project led by the Province of Limburg, the automated 3D-assignment of fault-lines to the appropriate model layers is being improved. We are currently looking for opportunities to make releasing this beta-functionality possible.
- Fence-diagrams
- The current 3D Tool is extended with an option to interactively create so-called fence-diagrams by using clipping planes along the major eucledian axes. This development is funded by the Alberta Geological Survey.
- Deviated wells:
- The current 3D Tool is extended with the option to visualize deviated wells. This development is funded by the Alberta Geological Survey.
- Additional runfile-options:
- We are implementing a number of additional runfile-options; for a specification of these runfile-options scheduled for a next iMOD-release, please have a look in chapter 10 of the iMOD User Manual.
- iMODBatch:
- Further improvement of the iMODBatch function XYZTOIDF including the computation of voxel models from bore hole data.
Other ongoing developments (status on the release date of iMOD 3.6):
- A new SUB-CR package (including 'creep'):
- The existing SUBsidence-package has been extended with 'creep' for clay and peat; the implemented concept is consistent with D-Settlement (geotechnical software). We hope we will be able to include this extension in an upcoming iMOD-release.
- Coupling of iMOD and DgFlow:
- Proof of Concept: iMODFLOW has been off-line coupled to DgFlow, a finite-element code developed by Deltares for simulating the consolidation process of dikes/embankments. An explicit head-flux boundary coupling has been realized, demonstrating the added value of incorporating regional, transient, groundwater flow boundary conditions for dike stability analysis. For more info, please contact Deltares.
- We also plan to make iMOD-MT3DMS/SEAWAT open source and include it in the iMOD-framework. This modified version of MT3DMS/SEAWAT uses an easy-to-use runfile for specifying model data for many layers and supports the common iMOD file-formats. We are making progress but slower than we hoped for; we will update you on the progress later.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.0, July 14th 2017:
iMOD 4.0 features the new Parallel Krylov Solver (PKS) package, allowing running your iMOD-models much faster or with much higher resolutions on multi-core computers. The Parallel Krylov Solver package was developed by Deltares, USGS and WEnR (Alterra) in close collaboration with Utrecht University and Technical University Delft (TUDelft).The new Parallel Krylov Solver can be used by doing the following 3 steps:
- Install free 'third party' MPI-software MPICH.
- Specify the PKS-solver settings.
- Start your parallel run, either from within the iMOD-GUI or from the command prompt.
- The partitioning process of your model domain into subdomains is fully automated; it results in a number of subdomains equal to the assigned number of CPU's (or threads).
- Your model input does not need any restructuring (thanks to iMOD's subdomain modelling techniques).
- The PKS solver settings can be specified in the iMOD-GUI interactively (step 2).
- Existing iMOD-runfiles can be fully reused, only the solver settings need a few small modifications. Solver settings can simply be edited manually too, for an example see the Runfile-section 10.20.2 'Updating from iMOD 3.6 to iMOD 4.0' of the iMOD User Manual.
- Tutorial 6, 'Model Simulation' now contains an exercise on how to approach load balancing your model on a multi-core machine.
- We have included a test program in the iMOD 4.0 allowing you to check your MPI-installation prior to using the PKS-package.
Please find below:
Ad 1) Summary of new and improved features and bug fixes.
New features:
- A new Parallel Krylov Solver package (PKS):
Deltares, USGS and WEnR (Alterra) together with Utrecht University and Technical University Delft have developed the new parallel solver package for MODFLOW called PKS (Parallel Krylov Solver). It is based on overlapping domain decomposition combining both the techniques of MPI and OpenMP. Preliminary results were presented during the 2016 AGU Fall Meeting (San Francisco, USA) and the MODFLOW and More 2017 conference (Golden, USA). Deltares R&D Highlights 2016 (page 64) summarizes first results for the Dutch National Hydrological Instrument and the global groundwater model PCR-GLOBWB. During the annual Dutch Deltares Software Days we organized a beta-release workshop 'iMOD Workshop Parallel Krylov Solver package' (June 14th, 2017, in Dutch); feedback from participants experimenting with a beta-version of the PKS-package on their own laptop was integrated in the current release.
Scope, limitations and expectations:- The PKS package was implemented for the 64-bit version of iMODFLOW only. So the PKS package can be used on Windows-based multi-core 64-bit computer systems, starting for example from Windows-based laptops and desktops (containing e.g. Intel Core i5 or i7 processors) to dedicated Windows-based multi-core 64-bit supercomputers.
- The PKS package uses MPI-techniques, therefore prior to using the PKS package so-called 'third party' MPI software (available free of charge) should be installed, see iMOD 4.0 installation instructions for details. We have included a test program in the iMOD 4.0 distribution allowing you to check your MPI-installation.
- The expected overall speedup strongly depends on the available hardware, load balancing of the partitions and the data input/output intensity:
- On a typical home computer / laptop with an Intel Core i5 processor inside (maximum 4 threads) the speedup will usually be limited. When assigning 3 (or 4) threads one may expect an acceleration factor of around 1.5 to 2; since MPI creates extra overhead communicating between the subdomains typically large and computational intensive models will potentially benefit from a parallel approach; also limiting the model output (frequency) contributes to the potential speedup. Relatively small models most likely will not benefit from a parallel approach due to the necessary extra communication between the subdomains.
- Test runs for the Dutch National Hydrological Model were carried out on a Windows server consisting of 2 Haswell 16-core CPUs with 128GB RAM resulting in a speedup of approximately a factor 5; in this case overall speedup is limited because only the MODFLOW-MetaSWAP part could be run in parallel mode; the coupled surface water (MOZART/DM) effectively still runs in serial mode (see Deltares R&D Highlights 2016, page 64).
- A test run for an artificial groundwater model covering the area of California (23897 column, 27974 rows, 335 million active cells, cell size 50 meters) using 128 cores on the SURFsara Cartesius Dutch supercomputer showed a reduction of an estimated single core runtime of 12 hours to approximately 17 minutes (a speedup factor of around 42) prior to load balancing efforts.
- Another test on the Cartesius supercomputer for the global PCR-GLOBWB model showed speedup factors ranging from ~17 to ~45 using 48 cores, depending on the partioning method (see Deltares R&D Highlights 2016, page 64).
- To play it safe, for large models, expect speedup factors starting from 1/3 of the number of the assigned number of cores. To get maximum speedup for large high resolution models, tuning is necessary; for more practical considerations, see section 12.32.3 of the iMOD User Manual.
- The PKS package can be used in combination with most iMOD-packages except the LAK-, MNW-, PST-, SFR- and UZF-packages; when an iMOD-model contains one or more of these latter package the single-core PCG-solver has to be used. The PKS package has been implemented for models including the MetaSWAP concept.
- The iMOD-GUI (i.e. the 'Model Simulation' window) has been extended for the use of the PKS package. In addition, a parallel run can also be started from the command prompt.
- The PKS package is not yet supported in the Project Manager.
- Tutorial 6, 'Model Simulation' of the iMOD User Manual has been extended with an exercise on:
- how to configure the PKS package.
- how to switch between the PCG- and PKS-solver.
- how to approach load balancing your model on a multi-core machine.
Improvements, bug fixes and extension of existing functionalities:
- In the iMOD-GUI:
- The 'Solver Settings' tab of the 'Model Simulation' window has been extended for the PKS-package. Here you can specify in how many subdomains your model domain should be partitioned by specifying how many CPU's should be used. It is recommended not to assign all available cores since your computer should be left some computational power for overall operating purposes.
- In the cross-section tool (ideal to compare single core and subdomain model results) the synchronization of overlying idf's has been tuned.
- In iMODFLOW:
- At downscaling:
- Sometimes for interpolation there was a wrong supporting point for the edge of a subdomain, this has been fixed.
- Sometimes cells with nodata-values were ignored and interpolated with surrounding data values, this has been fixed.
- Sometimes values of the first column were wrongly copied to the last column if the window of the fine grid not exactly coincided with the grid-lines of the coarser grid, this has been fixed.
- Sometimes for cells without precipitation an interpolation with surrounding cells values was performed resulting in precipitation awarded to cells originally having no precipitation, this has been fixed.
- At upscaling: KHV en KVV were upscaled arithmetically instead of geometrically, this has been fixed.
- At downscaling:
Ad 2) Developments (status on the release date of iMOD 4.0):
Current developments for upcoming releases:
- Extension of the Waterbalance Tool:
- The current water balance tool is extended significantly to allow constructing water balances for sub-regions, and selecting / deselecting MODFLOW- and MetaSWAP-water balance components interactively. The new water balance tool also allows specifying only a subset of the output time series and contains several time-related aggregation possibilities (e.g. averages per year, month and season). Visualization includes interactive stack-bars time series plots and schematic vertical cross-sectional overviews. This development is done together with Tauw and financed by the Dutch iMOD-CGO group.
- Movie tool:
- Also part of the above mentioned 'Water balance Tool extension'-project a movie tool is developed. It facilitates visualizing a time series of maps (e.g. calculated groundwater levels) and records it as a playable file.
- 3D fault assignment tool:
- As part of a consortium-project led by the Province of Limburg, the automated 3D-assignment of fault-lines to the appropriate model layers is being improved. We are currently looking for opportunities to make releasing this beta-functionality possible.
- Fence diagrams
- The current 3D Tool is extended with an option to interactively create so-called fence-diagrams by using clipping planes along the major Euclidian axes. This development is funded by the Alberta Geological Survey.
- Deviated wells:
- The current 3D Tool is extended with the option to visualize deviated wells. This development is funded by the Alberta Geological Survey.
- Additional runfile-options:
- We are implementing a number of additional runfile-options; for a specification of these runfile-options scheduled for a next iMOD-release, please have a look in chapter 10 of the iMOD User Manual.
- iMODBatch:
- Further improvement of the iMODBatch-function XYZTOIDF including the computation of voxel models from bore hole data.
Other ongoing developments (status on the release date of iMOD 4.0):
- A new SUB-CR package (including 'creep'):
- The existing SUBsidence-package has been extended with 'creep' for clay and peat; the implemented concept is consistent with D-Settlement (geotechnical software). We hope we will be able to include this extension in an upcoming iMOD-release.
- Coupling of iMOD and DgFlow:
- Proof of Concept: iMODFLOW has been offline coupled to DgFlow, a finite-element code developed by Deltares for simulating the consolidation process of dikes/embankments. An explicit head-flux boundary coupling has been realized, demonstrating the added value of incorporating regional, transient, groundwater flow boundary conditions for dike stability analysis. For more information, please contact Deltares.
- We also plan to make iMOD-MT3DMS/SEAWAT open source and include it in the iMOD-framework. This modified version of MT3DMS/SEAWAT uses an easy-to-use runfile for specifying model data for many layers and supports the common iMOD file-formats. We are making progress but slower than we hoped for; we will update you on the progress later.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.0.1, September 8th 2017:
iMOD 4.0.1 is a bugfix-release: in iMOD 4.0 output from MetaSWAP to idf-files was accidentally blocked, this has been repaired in this release.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.1, October 27th 2017:
iMOD 4.1 was released because the automated 3D-assignment of fault-lines to the appropriate model layers in the Horizontal Flow Barrier package has been improved. This release was financed by the Dutch IBRAHYM-consortium through a consortium-project led by the Province of Limburg. iMOD 4.1 also faciltates the use of non-uniform grids (not in combination with the PKS-package and in combination with ISG containing cross sections wider than the cell size).
Please find below:
- A summary of new and improved features.
- A summary of current developments.
- Other ongoing developments.
Ad 1) Summary of new and improved features and bug fixes.
- Extension and improvement of the 3D-parameterization of faults in the Horizontal Flow Barrier package:
- The automatic 3D-fault-parameterization has been extended: it is now possible to use 3D GEN-files to assign fault to the appropriate model layers. Also the concept of how to calculate the fault resistances has been refined (see section 12.15 in the iMOD User Manual).
- When using 3D GEN-files each intersection of a cell face will result in additional resistance harmonically. When all 3D line segments have been processed the resistance is corrected for situations the fault is partly unconfined; a default resistance value equal to 1 is assigned to that unconfined part of the fault and included in the harmonic calculation of the overall fault resistance.
- Contrary to using 3D GEN-files: when using 2D GEN-files the resistance is now solely determined by the first line segment intersecting the cell face. In the earlier release resistance was added when the GEN-file contained an additional line segment intersecting the same cell face, this is not the case anymore; as a result final resistances may be significant lower.
- To allow inspection the resulting parameterization is exported in GEN- and DAT-files.
- Non-uniform grids
- Earlier 2.x versions of iMOD featured the use of non-uniform grids. During the transition to open source in 2013 this feature was not included. In iMOD 4.1 this functionality has been revived. Around the area of interest a buffer can be defined where cell sizes gradually increase in the direction of the outer boundary of the model area. This is also facilitated in the iMOD-GUI. Optionally pre-defined irregular finite difference grids can be used. Limitation: non-uniform grids cannot be combined with 1) the PKS-package and 2) the ISG-package containing 1D- or 2D-cross sections wider than the cell sizes.
- Refinement of automated layer-assignments:
- In earlier versions a minimum of 5% of the total penetration depth of a river reach or well screen was used to assign a river reach or well to a particular model layer. In iMOD 4.1 this minimum criterion is not used anymore because it results in more realistic schematizations when rivers or wells penetrate very thin layers.
- For areas containing many thin shallow lying aquifers containing clayey sediments this may result in significant different river results. For wells this effect will be smaller because the vertical discharge distribution was and still is proportional to the penetrated thickness times the conductivity. The conductivities values are now based on the user input transmissivities, not on the average of the CC- and CR-terms anymore; this may cause a vertical redistribution of the assigned well rates.
- Correction of the ‘brush'-option of the ISG-gridding:
- The ‘brush'-option of the ISG-gridding has been corrected for cases 1D cross-sections are wider than a model cell by maximizing the conductance according to the cell width. This may result in lower overall conductances causing less interaction between ground- and surface water.
- When using the ‘brush'-option now only cells covered with a minimum fraction of 50% will be assigned a river reach.
- iPEST:
- iPEST output parameters are now generated with additional significant digits.
- The contents of the logfile have been extended.
- Before iPEST stopped as soon the steepest gradient and the gauss-newton direction conflicted, now iPEST continues ignoring the generated parameter values of that cyclus.
- KHV, KVV en KVA are exported for re-use purposes.
- iPEST can now be invoked through the additional argument ‘–ipest *.pst‘
- Now the Brook-Coley parameter of the UZF-package can be included in iPEST.
- SFT-package:
- Constant input values are now processed correctly.
- The SFT-package can now be ignored when IFVDL is active; in that case river conductances are calculated using the ‘De Lange formula' using thicknesses and transmissivities of model layer 1 instead of the CC- and CR-terms.
- Starting conditions:
- Now a non-steady state run can start with one steady-state ‘stress period'; This is particularly useful if you want to start a non-steady state simulation where the groundwater system is in full equilibrium with its initial boundary conditions.
- A run is set as transient based on the definition of the stress periods in the DIS-file (SS or TR), not on the length of the first stress period anymore.
- Generation of NAM-files:
- The value of MAXWBALERROR is equal to 0 by default. If a value > 0 is specified this value is also copied to the STOPERROR keyword in the BAS-file.
- The OCD-file now contains correct output layer numbers.
- Extension of output filenames:
- The keyword 'ISAVEEND' has been added to the MODFLOW MET file. This facilitates the option to use the end time of a time step in the output filenames.
- Improvement of the optional drain level correction:
- Via the ICONCHK keyword DRN-levels can be corrected each time new DRN-input is being read. However, this correction was not performed if only new RIV- or ISG-input was processed (in case no new DRN-input was read for that particular time step); this has been corrected.
- If ICONCHK is active this is logged in the LST output file.
- ICONCHK has been reprogrammed for memory efficiency.
- Well data:
- Associated TXT files containing well data now optionally can contain additional 2 columns without causing iMOD to generate an error. iMOD continues using the first 2 column only.
- Time discretisation:
- The iMOD-Batch function RUNFILE now supports ‘hours', ‘minutes' and ‘seconds'.
- The Time Series Tool in the iMOD-GUI has been extended to support processing folders containing filenames with and without ‘hours', ‘minutes' and ‘seconds'.
- Downscaling (going from coarse data to a finer model network) of conductances contained an error, this has been corrected.
- The combination of a constant + factor + a to be added impuls was processed correctly only when an IDF file was used; this has been corrected.
- Water balance terms of inactive cells are now assigned the HNOFLO-value (‘huge number').
- The maximum value of a resistance has now been increased form 1.0E6 to 1.0E8 days.
- A constant head cell becomes an inactive cell if fully surrounded by inactive cells or constant head cells.
- The iMOD Manager has been extended with the option to (de-)select files using wildcards.
Ad 2) Developments (status on the release date of iMOD 4.1):
Current developments for upcoming releases:
- Extension of the Waterbalance Tool:
- The water balance tool is extended significantly to allow constructing water balances for sub-regions, and selecting / deselecting MODFLOW- and MetaSWAP-water balance components interactively. The new water balance tool also allows specifying only a subset of the output time series and contains several time-related aggregation possibilities (e.g. averages per year, month and season). Visualization includes interactive stack-bars time series plots and schematic vertical cross-sectional overviews. This development is done together with Tauw and financed by the Dutch iMOD-CGO group.
- Movie tool:
- Also part of the above mentioned 'Water balance Tool extension'-project a movie tool is developed. It facilitates visualizing a time series of maps (e.g. calculated groundwater levels) and records it as a playable file.
- Fence diagrams
- The current 3D Tool is extended with an option to interactively create so-called fence-diagrams by using clipping planes along the major Euclidian axes. This development is funded by the Alberta Geological Survey.
- Deviated wells:
- The current 3D Tool is extended with the option to visualize deviated wells. This development is funded by the Alberta Geological Survey.
- Additional runfile-options:
- We are implementing a number of additional runfile-options; for a specification of these runfile-options scheduled for a next iMOD-release, please have a look in chapter 10 of the iMOD User Manual.
- iMODBatch:
- Further improvement of the iMODBatch-function XYZTOIDF including the computation of voxel models from bore hole data.
Ad 3) Other ongoing developments (status on the release date of iMOD 4.1)::
- A new SUB-CR package (including 'creep'):
- The existing SUBsidence-package has been extended with 'creep' for clay and peat; the implemented concept is consistent with D-Settlement (geotechnical software). We hope we will be able to include this extension in an upcoming iMOD-release.
- Coupling of iMOD and DgFlow:
- Proof of Concept: iMODFLOW has been offline coupled to DgFlow, a finite-element code developed by Deltares for simulating the consolidation process of dikes/embankments. An explicit head-flux boundary coupling has been realized, demonstrating the added value of incorporating regional, transient, groundwater flow boundary conditions for dike stability analysis. For more information, please contact Deltares.
- We also plan to make iMOD-MT3DMS/SEAWAT open source and include it in the iMOD-framework. This modified version of MT3DMS/SEAWAT uses an easy-to-use runfile for specifying model data for many layers and supports the common iMOD file-formats. We are making progress but slower than we hoped for; we will update you on the progress later.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.1.1, November 24th 2017:
iMOD 4.1.1 is a bugfix-release: an error in the GENSNAPTOGRID-functionality of iMOD 4.1 used to parameterize the Horizontal Flow Barrier has been fixed. The iMOD 4.1 and iMOD 4.1.1 releases were financed by the Dutch IBRAHYM-consortium through a consortium-project led by the Province of Limburg. iMOD 4.1.1 also contains a number of other improved features and minor bugfixes.
Ad 1) Summary of new and improved features and bug fixes.
Improvements, bug fixes and extension of existing functionalities:
In iMOD-GUI:
- GENSNAPTOGRID triggered a crash when used without top/bot IDF files, this has been fixed.
- Visualizing length-profiles for ISG-files in combination with location of structures / cross-sections etc. reactivated sorting calculations points along a segment for correct visualisation on a profile, this has been fixed.
- Sometimes usage of NaN's inside an IDF file gave a crash; this has been resolved by skipping these values.
- The package definition of the project manager gave wrong results in case files were used with time information but needed to be treated as time contstant in the file selection procedure. Also there was a crash when no appropriate files were found and the user continued; this has been fixed.
- Invoking the iMODbatch PLOT function for only one IPF, GEN of IFF-file without specifying a window resulted in an improper plot, this has been corrected.
- iMOD crashed when trying to create a legend with a single class, this has been fixed.
- Displaying IPF-points with 7 significant numbers triggered an infinite loop when trying to generate a scale bar; this has been fixed.
- Plotting of IPF files went wrong in case the number of IPF files sequentially drawn was less and in a different order than of a previous plotting sequence, this has been fixed.
- In the ‘Start iMOD'-window the user can now sort the available IMF-files alphabetically, by size or by date.
- The user cannot start iMOD accidentally anymore via the option ‘Open an existing iMOD project' without selecting an existing IMF-file.
- Plotting of time series:
- In the status bar the space for displaying the mouse position coordinates has been increased;
- two digits for displaying dates are added whenever the mouse is moved;
- the to be plotted time series need to be sorted; iMOD stops plotting the remainder of the time series when encountering dates earlier than already plotted.
- nodata lines in TXT files will be skipped.
- Whenever the Model Export Configuration is 'MODFLOW2005' (NAM file for simulation) and the simulation time extends the input of the txt files for wells, the well rate remains equal to the last entry.
- When plotting ISG profiles the plot-sequence has been revised: first the left parameter (a) is plotted followed by the right parameter (b). Also the assignment of the parameters colours has been corrected.
- The function to import CSV's in ISGEDIT with segment/calculation number has been extended to include remaining attributes. It allows to insert time series of data sequentially to different number of selected segments/calculation nodes. Futhermore ISGEDIT can now also export to CSV including these additional attributes.
- Further improvement of the iMODBatch-function XYZTOIDF including the computation of voxel models from bore hole data.
- IPF Analyse now doesn't crash when the user visualizes an IPF accidentally containing associated TXT files definitions that do not exist.
- The Interactive Pathline Simulator (IPS) now allows the use of the $DBASE$ keyword in iMODPATH.RUN files.
- The PKS-settings in iMOD-GUI 4.1 were inconsistent, this has been fixed.
- Shapefiles containing POLYLINEZ shapes are now converted correctly.
- The 3D Tool has been extended with the optioin 'Deactivate Associated Files' to make visualizing IPF's containing secondary Z-coordinates in 3D possible.
- In the Project Manager in the 'Define Characteristics' window for the RCH and EVT packages the value 'Assign parameter to modellayer' = 0 (assign to modellayers automatically) is not allowed anymore.
In iMODFLOW:
- Whenever the simulation time extends the input of the txt files for wells, the well rate remains equal to the last entry.
- The use of NSCL=0 (network specified by an IDF file) generated an error, this has been fixed.
- When ISAVE_DATE = 1 (save results for enddate) and hours, minutes and seconds are not specified (zero's) the names of the IDF's will not contain zero's; the filenames will contain e.g. '20171117' instead of '21071117000000'.
- The multiplication of RCH in iPEST was also done in case the recharge was not reread from file, causing an erroneous repetition of multiplications, this has been fixed.
- Now by default the switch "simgro_opt = -1" is included in the MetaSWAP inputfile para_sim.inp; "simgro_opt= -1" is the switch to have MetaSWAP generate IDF-output.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.2, December 21st, 2017:
iMOD 4.2 contains the following new functionalities:
- Extension of the water balance tool, funded by the Dutch iMOD-CGO group and implemented in collaboration with Tauw.
- A new movie tool, also funded by the Dutch iMOD-CGO group.
- Extension of the 3D Tool with Fence diagrams & deviated wells, funded by Alberta Energy Regulator | Alberta Geological Survey.
Please find below:
- A summary of new and improved features.
- A summary of current developments.
- Other ongoing developments.
New and extended functionalities in iMOD-GUI:
- Extension of the Waterbalance Tool
The current water balance tool was extended significantly to allow constructing water balances for sub-regions, and selecting / deselecting MODFLOW- and MetaSWAP-water balance components interactively. The new water balance tool also allows specifying only a subset of the output time series and contains several time-related aggregation possibilities (e.g. averages per year, month and season). Visualization includes interactive stack-bars time series plots and schematic vertical cross-sectional overviews. This extension was implemented in collaboration with Tauw and financed by the Dutch iMOD-CGO group. - A new movie tool
Also part of the above mentioned 'Water balance Tool extension'-project a movie tool was developed. It facilitates visualizing a time series of maps (e.g. calculated groundwater levels) and records it as a playable file. - Extension of the 3D Tool, funded by the Alberta Energy Regulator | Alberta Geological Survey:
- Fence diagrams
It is now possible to create 3D fence-diagrams. The fence-diagrams can both be defined interactively using clipping planes along the major Euclidian axes and can be generated automatically. Clipping can also be activated for IPF- and IFF data. Section 7.3 '3D Tool' of the iMOD User Manual has been upgraded accordingly. - Deviated wells
The current 3D Tool was extended with the option to visualize deviated wells. An additional format in associated txt files was added to support the use of deviated wells; the values can vary in depth by x and y. - An 'Explode IDF'-option has been added. It allows to increase the vertical distance between selected IDF's in a 3D view interactively. For each individual IDF ?le this vertical distance can be specified by setting a slider interactively. See section 7.3.3 of the iMOD User Manual.
- It is now possible to Import PETREL pillar-files and convert them to 3D-GEN files.
- iMOD now displays boreholes in the 3D Tool in a fancy way by default when less than 1000 points are read, previously this was 500.
- Fence diagrams
- The keyword IBLNTYPE was added to the IMODBATCH function XYZTOIDF to allow using BLNFILES not only as a fault but also as a polygon.
- When the ‘CheckRun' is performed in the ‘Start Model Simulation' window now also a result-message is written to the logfile when all files listed in the runfile exist (so not only when files do not exist).
- Now it is also possible to search for time dependent files without timesteps using wildcards; those files will be assigned a steady-state timestamp.
- Non-uniform grids can now be used with the ISG-package containing 1D- or 2D-cross sections wider than the cell sizes.
- Prior to closing the ISGEDIT window iMOD now askes whether to save changes.
- When exporting a model to a namfile, iMOD assigns a NODATA value to the VCW and KVV parameters in the LPF package in case the cell underneath is inactive.
- In SFR now a warning is given if the upstream bottom level is lower/equal to the downstream bottom level since this is not allowed by the SFR package.
- Now you cannot assign a layer number equal to 0 anymore in the RCH- and EVT-package since this is not a logical thing to do.
- The keyword IDATE in the PLOTRESIDUAL IMODBATCH function has been removed since a transient model yields an average not involving a particular date anymore.
- The waterbalance analyse tool crashed when the Modflow/Metaswap buttons were used and no files were available, this has been fixed.
- An error in 'search for segment number' in ISGEDIT was corrected.
- Reading an MDF file with non-existing IDF file(s) generated a crash in the profile-tool, this has been fixed.
- The IMODBATCH function IPFSPOTIFY generated results from non-initialised variables when a geological layer was not present in a model layer, this has been resolved by initialising those fractions as 0.0.
- Plotting of IUPSEG and IUDWN for SFR was mixed up, this has been fixed.
- Labels can be displayed in 3D again.
- An error in creating Metaswap input files for non-uniform grids (available since the iMOD 4.1 release) has been fixed.
- A more efficient allocation of zones to parameters has been implemented in iPEST.
- When using iPEST the RCH parameter was modified for each stress period regardless it was reread from disc, this has been fixed.
Current developments for upcoming releases:
- Larger coordinates
As a spin-off in a modelling project funded by Landesamt für Landwirtschaft, Umwelt und ländliche Räume, Geologischer Dienst Schleswig-Holstein, Germany, iMOD will be extended for the use of larger coordinates. This extension will become very handy with respect to the planned European transition to the ETRS89 coordinate system including an extra two-digit zonal prefix. The extension for larger coordinates will be part of a release scheduled for the spring of 2018.
- Additional runfile-options
We are implementing a number of additional runfile-options; for a specification of these runfile-options scheduled for a next iMOD-release, please have a look in chapter 10 of the iMOD User Manual. - A new SUB-CR package (including 'creep')
The existing SUBsidence-package has been extended with 'creep' for clay and peat; the implemented concept is consistent with D-Settlement (geotechnical software). We hope we will be able to include this extension in an upcoming iMOD-release. - Coupling of iMOD and DgFlow
Proof of Concept: iMODFLOW has been offline coupled to DgFlow, a finite-element code developed by Deltares for simulating the consolidation process of dikes/embankments. An explicit head-flux boundary coupling has been realized, demonstrating the added value of incorporating regional, transient, groundwater flow boundary conditions for dike stability analysis. For more information, please contact Deltares. - We also are working on making iMOD-MT3DMS/SEAWAT open source and include it in the iMOD-framework. This modified version of MT3DMS/SEAWAT uses an easy-to-use runfile for specifying model data for many layers and supports the common iMOD file-formats. We are making progress but slower than hoped for; we will update you on the progress later.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.2.1, January 31st, 2018:
iMOD 4.2.1 is a bugfix-release, it fixes a rare OpenGL-error present in the GUI of iMOD 4.2 only.
Current iMOD 4.2 users not experiencing OpenGL-errors don't have to upgrade to iMOD 4.2.1 because the source codes of iMOD 4.2.1 and iMOD 4.2 are exactly the same.
Compared to earlier releases the compiler-settings used for the GUI of iMOD 4.2 were a bit more strict. Recently in a couple of individual cases these stricter settings cause an OpenGL-error. For the iMOD 4.2.1-release the pre-iMOD 4.2-compiler settings have been used.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.3, December 19th, 2018:
iMOD 4.3 contains the following major new functionalities:
- Large coordinates, as a spin-off in a modelling project funded by Landesamt für Landwirtschaft, Umwelt und ländliche Räume, Geologischer Dienst Schleswig-Holstein, Germany.
- Extended queries and artificial boreholes in the 3D-Tool, funded by Alberta Energy Regulator | Alberta Geological Survey.
- Pumping Tool, initially funded by HDSR, now reactivated by Deltares.
Please find below: Ad 1) Summary of new and improved features and bug fixes.
Major new functionalities:
- Large coordinates
As a spin-off in a modelling project funded by Landesamt für Landwirtschaft, Umwelt und ländliche Räume, Geologischer Dienst Schleswig-Holstein, Germany, iMOD has been extended for the use of large coordinates. This extension is very handy with respect to the planned European transition to the ETRS89 coordinate system including an extra two-digit zonal prefix. iMOD now supports large coordinates with 12 significant numbers for displaying IDF, IPF, ISG and GEN files. Also IDF's can be stored in a double precision format. - Extended query-possibilities in the 3D-tool
This option allows querying on specific items of boreholes to define a subset. Tutorial 3 has been extended accordingly. The development of this extension was funded by Alberta Energy Regulator | Alberta Geological Survey. - Artificial boreholes in the 3D-tool
This option allows specifying artificial boreholes through a selected set of IDF files. Tutorial 1 has been extended accordingly. The development of this option was funded by Alberta Energy Regulator | Alberta Geological Survey. - Pumping Tool
The Pumping Tool (previously known as Scenario Tool) has been reactivated and updated. It is now also possible to save the configuration of the PUMPING TOOL in IPF files for the observations and/or extractions. The iMOD User Manual has been extended with a separate tutorial on the Pumping Tool (tutorial 12). An initial version of the Pumping Tool was funded by HDSR, the update by Deltares.
- Binary GEN files
GEN files can now also have a binary format, resulting in much faster processing performance. The new binary GEN file format supports the use of circles, rectangles, polygon, points and lines. The binary GEN files have been fully incorporated in the iMOD-GUI and iMODFLOW. A binary GEN file can be used in the HFB package. Limitations: It is not possible to use ILAY=0 with a binary GEN file. Currently the binary GEN files is not replacing the 3D ASCII GEN file; these 3D BINARY GEN files will be part of a next release of iMOD. - Convert ESRI shapefile to binary GEN file
There is now a separate option to convert ESRI Shapefiles to the new iMOD binary GEN file format, and to convert a BINARY GEN file to the ASCII GEN format. - WINDOWS system variables in PRF file
This option allows using Windows system variables, such as %USERNAME% in PRF files (for an example see section 9.1 of the iMOD User Manual). - Selecting and sorting in iMOD Manager
It is now possible to select files in the iMOD Manager via a query. Also it is now possible to sort the selected files using keywords; this allows to sort e.g. a large list of IDF files according to desired TOP- and BOT-sequences. - Moveable axes and scale bar
In the 2D canvas of iMOD it is now possible to change the position of the scale bar (by dragging the inner upper left corner) and axis (dragging the x- or y-axis) interactively (see section 2.6 in the iMOD User Manual), and turn them on and off as well (via View, Layout). - Moveable 2D survey-plot in cross-section tool
In the cross-section tool it is now possible to resize and re-position the plot feature 'Plot Survey Cross-section' interactively (see iMOD User Manual page, section 7.1.1.5, 'Plot Survey Cross-section'). - Pilot points
The concept of pilot points is now supported in iPEST and steered via the iMOD PRJ-file when generating RUN- and NAM-files. Also iMODFLOW has been extended to support the use of pilot points. It is advised to use iPEST in combination with NAM-files as this is more time-efficient. A dedicated tutorial on model calibration in iMOD using iPEST is scheduled for a next iMOD release. - Extended Define Startpoint options
- It was already possible to use a circle, rectangle, line, polygon and point to define start points, this now has been extended with the option to use a grid.
- Start points can now be snapped to the centres of the cells of the underlying model. Remark: since the start point definition is independent of the underlying model grid, the Preview-button does not include the effect of 'snap-to-modelgrid'.
- 3D IPF files usable in Pathline Simulation
In addition to the ISD-files (iMOD Start point Definition files) for Pathline Simulation it is now possible to use existing 3D IPF files (containing x-, y- and z-coordinates) to specify start points in the 'Input'-tab of the Start Pathline Simulation window. - ISG option to export to CSV-format
It is now possible to save an entire selected segment including dates as a CSV file when exiting the ISG Edit window. This makes it easier to modify and adjust time series for multiply calculation points along (several) segments. See iMODBATCH functions below. - ARR-files added to list of 'Known Files'
If a model is converted to a default MF2005 file format, the yielding ARR-files can now be loaded into iMOD as they are converted to IDF files. In that manner it is easy to check the actual (modified) simulation data. The ARR files can be loaded for all packages and modules, excluding SFR. - The following iMODBATCH functions have been added and/or extended:
- ISGEXPORT / ISGIMPORT: it is now possible to convert (export) an entire ISG file to TXT-files and convert (import) them back into a binary ISG. This allows a) to correct a corrupt ISG file and/or b) to modify it outside iMOD.
- XYZTOIDF: the iMOD BATCH function XYZTOIDF is extended with functionalities to generate VOXELS model out of IPF files. Also the use of GEN files is supported to delineate surfaces in the interpolation.
- RUNFILE:
- it is now possible to overrule an export for the WEL/ISG package for each time step (DWEL/DISG).
- it is possible to save the modelling results in double-precision (IDOUBLE).
- IMODPATH has been extended with the post-processing option to thin out particles that pass interfaces. It also contains the keyword WINDOW to allow for a simulation within a part of the model.
- IDFSCALE: has been modified to support the upscaling of hydraulic gradients in x- and y-direction.
- New exercise on Movie Tool
Tutorial 3 has been extended with an exercise on the Movie Tool. - IMF files
- The position of axes and the scale bar is saved in the IMF file.
- Whenever opening IMF files, iMOD informatively lists missing files.
- The position of axes and the scale bar is saved in the IMF file.
- Transparency
Turning on the transparency-option for more than 10 IDF files now generates a warning, allowing the user to avoid undesired results. - IDF Calculator
- the LOG10() function has been added to IDF CALC.
- if an entered formulae is not correct, now a warning message appears.
- IDF Edit
- When exiting the IDF EDIT Calculation window without having used the 'Calculate'-button, iMOD now issues a warning message.
- For the selected cells now statistics can be generated.
- Pathline Simulation in strong sinks
The computation of the duration of a streamline inside a strong-sink is computed with the velocity the particle had at the entrance of the cell. Previously, the velocity was set to zero as standard MODPATH does not compute particles inside a cell with a strong sink. - 'Map value' coordinates fixed
Now also the 'Current Location'-coordinates in the 'Map Value' window are fixed as soon as the left-mouse button has been pressed once; when hovering over the canvas the x- and y-coordinates are still updated in the lower left corner of the canvas. - Improved algorithm for computing tic-intervals in graphs and plots
- Locked iMOD Manager when 3D Tool is active
An active 3D tool does not synchronise with the iMOD Manager when new or other files are selected in the iMOD Manager. To avoid unjust expected synchronisation it is not allowed anymore to add files to the iMOD Manager whenever the 3D tool is active. - Extra SFR checks
- A check has been added whether the stage declines along the stream-direction
- SFR simplifies cross-sections to maximum 8 points; a comparison check between the length of the original cross-section and the 8-point cross-section has been added issuing a warning when the two deviate too much.
- Maximum value of vertical resistance removed
There is not a maximum vertical resistance in iMODFLOW anymore, this was 1.0E8 days. - iPEST reviewed and updated
- The source code of the build-in parameter estimation program has been reviewed and updated completely to be more robust.
- The output files are commented and restructured as well as the algorithm has been improved to be more robust.
- Importing SOBEK models
The import of SOBEK models is now supported in the 64-bit version of iMOD-GUI instead of the 32-bit version.
- Order-dependency in runfile for RIV / ISG and DRN / OLF
It appeared that rivers and drains were omitted when in the runfile the positions of RIV and ISG were switched, or when the positions of DRN and OLF were switched. This has been corrected. - A combination of FACTOR and IMPULSE did not work for BND, this has been corrected.
- Smoothing of downscaled starting heads
In case the modelling grid was not aligned or smaller than the head grid, fixed boundary heads were not properly smoothed when downscaled, this has been corrected. - Up- and downscaling of KVA-factor corrected
For cases the KVA-distribution is heterogeneous the up- and downscaling has been corrected (in most existing regional iMOD-CGO-models this bugfix has no consequences because KVA is 1/3 x KH, or KVA is not used). KVA is now a mandatory parameter for 3-D model that use TOP, BOT and KHV parameters. - Upscaling of ROOTDEPTH and LANDUSE now consistent
Upscaling the ROOTDEPTH-variable is now consistent with the upscaling algorithm of the LANDUSE-variable (both MetaSWAP-variables), namely by selecting the most frequent LANDUSE and ROOTDEPTH-class. Previously the most frequent LANDUSE-class was selected and an average value for the ROOTDEPTH was used, causing inconsistencies between LANDUSE and ROOTDEPTH. - Assignment of rivers and drains to model layers corrected
In the previous version of iMOD the rivers the RIV's and ISG's were assigned to model layers similar to the method used for the WEL package. As a consequence, the RIV's and ISG's were not assigned to model layers with low permeability, which is not correct. The assignment for RIV's and ISG's, and therefore the conductance distribution over model layers, is now based on the thickness of intersected model layers. - ISG Edit
Gridding of an ISG with a single transient entry yielded a crash, this has been resolved. - IBDG / ISUMSAVE
The global option to sum all fluxes per package (IBDG) was misprogrammed, it has been fixed now. Per stress period this can be overruled via ISUMSAVE. - ICONCHK
The ICONCHK-check was only performed when transient DRN- and RIV-inputs were read simultaneously; wrongly no check was performed when only new transient RIV-input was processed, this has been corrected. - IDEBUG
This option wrote the IDF file to the input folders instead of the output folder, this has been corrected.
Current developments for upcoming releases:
- The 64-bit version of iMODFLOW contains a coupling between MODFLOW and the MetaSWAP module; in a project funded by the iMOD-CGO group iMOD is extended to facilitate inspection and visualisation of these couplings, and results in an iMOD release scheduled in May 2019.
- We are working on making iMOD-MT3DMS/SEAWAT open source and include it in the iMOD-framework. This modified version of MT3DMS/SEAWAT uses an easy-to-use runfile for specifying model data for many layers and supports the common iMOD file-formats. A release is now scheduled mid 2019.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 4.4, July 17th, 2019:
iMOD 4.4 contains the following major new functionalities:
- Use of the 3D tool for 4D visualization.
- iPEST (parallel)
- iPEST analyser
- METASWAP analyser
Please find below: Ad 1) Summary of new and improved features and bug fixes.
New functionalities:
- 3D Tool:
- 4D-FILES; the 3D TOOL can display time-variant VOXELS;
- Colour legends can be modified from the 3D TOOL on IDF files, it is not needed anymore to leave the 3D TOOL to modify the legends.
- iPEST
- Added Conductance of General Head Boundary as parameter (GC);
- Added Specific Yield as parameter (SY);
- Added Corey Brooks as parameter (EP);
- Parameters for iPESTP are present in the input window for iPEST in iMOD;
- Prior estimates are supported in combination with a Tikhonov Regularisation.
-
iPEST Analyser
Functionality to analyse the results of iPEST quickly. -
METASWAP Analyser
Functionality to analyse the converted input files for MetaSWAP. -
BACKGROUND IMAGES
It is possible to display background images with different percentages of transparency and/or grey tones. -
ISG EDIT
-
Flipping of coordinates can be done on all selected segments;
- Calculation Points, Structures, Cross-sections are sorted internally based on the distances.
-
-
UZF Analyser
-
Functionality to analyse the moisture-profiles of the UZF-package;
- A new TUTORIAL on model calibration using iPEST and iPESTP.
-
-
Changed HFB data handling in an iMODFLOW RUN file!
For reasons, the resistance of HFB input in a RUN file is now defined as the value FCT * IMP (see Data Set 10 and 11). In precious iMOD versions the resistance was defined by FCT only. As a result:
- In case IMP=1.0, the value of the resistance does not change compared to previous versions.
- In case IMP=0.0, the value of the resistance will become 0.0, representing an imperial barrier!
-
IPF PLOT
Option to display labels with a format, especially for numbers this is an improvement; labels will be plotted in black always. -
iMODFLOW
-
Measurements inside the buffer are not excluded anymore;
- Upscaling of RCH and EVT is taking NODATA as 0.0 value instead of discarding them;
- Improved modification to the gridding of faults, small caps might occur in between nodes of segments, this has been resolved
-
- iMODPATH
- Option to enter a folder name to save results for the particle tracking;
- Inclusion of the fraction parameter (as present in the GUI) in the iMODPATH runfile.
-
$DBASE$
Usage of $DBASE$ generates an error message whenever DBASE is not specified in the PRF file. -
IMODBATCH
- IDFMEAN can be used without the keyword ILAYER, SDATE and EDATE such that all IDF’s can be SUMMED, regardless of ILAYER, SDATE and EDATE;
- RUNFILE supports:
- a time discretisation of MINUTES at the keyword IDT. Be aware that IDT=1 was hourly timesteps previously, which is now MINUTES. IDT=2 is hourly now as IDT=1 is minutes. The rest has to be shifted accordingly. Improved modification to the gridding of faults, small caps might occur in between nodes of segments, this is resolved;
- the option to ignore polynomial interpolation of input data if the modelled network is finer than the input data. This can be invoked via the options INTKHV, INTTOP, INTVCW etc.
- IPFSTAT added the option to cluster points.
-
IPEST
Optimization of RCH in combination with unsaturated conditions did not allow recharge to be changed for layers other than 1. -
iMODFLOW
- Usage of IMP for RCH and EVT added m instead of mm to the package;
- Export of IDF in DEBUG mode did not include any modifications from FCT and/or IMP values;
- Optimization of the parameter RE was not correct in case of unsaturated conditions and the recharge may be directed to a model layer other than 1;
- Writing of the SCAT_SVAT.INP did not take into account the location of an extraction well for artificial recharge correctly, this only affected configurations using IPF files for artificial recharge.
-
iMODPATH
Flowlines are saved with 3 digits (F15.3). -
WATERBALANCE TOOL
The computation of the BDGFRF and BDGFFF were incorrect. -
3D HFB PACKAGE
-
The thickness at the location of a fault was wrongly defined, now it is the minimum (bottom) and maximum value of both cells on both sides of the fault. The z-value presented in the 3D GEN files that come out of iMODFLOW and/or iMOD express this value now as well;
- Output GEN files are written in appropriate number of digits.
-
-
IDFCALC
The function MIN in the SPECIAL tab did not work properly and has been corrected. -
PROJECT MANAGER
Starting a model with iPEST did not initiate iPEST correctly. Via the iMOD Batch Function RUNFILE this worked correctly via the keyword IPEST. -
PROFILE TOOL
-
Returning from the Profile Tool forced to redraw the graphical canvas which is not necessary;
- Labels and points for IPF points were not plotted.
-
-
IDFDRAW
Vector drawing of direction in IDF-files via the ADJUST LEGEND option was not converted to double precision. -
CREATEIPF
Saving an IPF yielded a BINARY GEN file, this has been resolved. -
ISGEDIT
- Structures were not read in correctly;
- Flipping of the coordinates did not modify all the attributes correctly.
-
IMODBATCH
- IPF2ISG/GENTOISG crashed due to undefined setting for single/double precision;
- GENTOISG crashed importing a BINARY GEN file;
- SOLID yielded incorrect results for the hypothetical interfaces;
- XYZTOIDF
- added the option to apply individual weights to points in Kriging;
- improved the SOURCEDIR option in combination with MEAN;
- ISGGRID keyword ICDIST did not work;
- GENTOISG flips coordinates based on water levels from specified IDF files;
-
3D TOOL
- Labels were not plotted for IPF points without associated files;
- If all files are outside visible extent, iMOD does not crash anymore.
-
IMPORTSOBEK
Usage of the option to convert a SOBEK model into an ISG yielded a crash due to the renewed double precision format. By default, the converted ISG will be a double-precision ISG. -
GENERAL
Sorting of files from the iMOD Manager in alphabetical order.
Ad 2) Developments:
End 2019 Deltares will release iMOD 5.0, which includes the new iMOD Water Quality (iMOD-WQ) component. iMOD-WQ integrates the functionalities of SEAWAT (density-dependent groundwater flow) and MT3DMS (multi-species reactive transport calculations) and brings these codes into the iMOD environment by using iMOD formats for input and output. The new executable of iMOD-WQ will become a separate executable with its own runfile format, and will be described in a separate chapter in the iMOD manual. A Beta version of iMOD 5.0 is available for the iMOD community and can be requested via the helpdesk imod.support@deltares.nl. The iMOD community is encouraged to try and test the new functionalities and communicate the findings with the iMOD team, to enable further improvement of the code towards de definitive iMOD 5.0 release.
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
Release notes iMOD 5.0, December 24th, 2019:
iMOD 5.0 contains the following major new functionalities:
- New Simulation Manager
- Support the usage of MT3D/SEAWAT
- Use of the 3D tool for 4D visualization
- support the usage of MODFLOW6.
Please find below:
Ad 1) Summary of new and improved features and bug fixes.
New functionalities:
- iMOD WQ
iMOD Water Quality (WQ) is a new main feature within iMOD. iMOD-WQ integrates the functionalities of SEAWAT (density-dependent groundwater flow) and MT3DMS (multi-species reactive transport calculations) and brings these codes into the iMOD environment by using iMOD formats for input and output.
- IPEST
- Option to wait before starting another process using the keyword NSWAIT in combination with iPESTP;
- Option to save the Jacobian Matrix and Eigenvalues/-vectors in a Matlab *.m file for debugging. Use the keyword PDEBUG=1 to activate that;
- Usage of the DRES in iPEST (acceptable error) is modified, if the absolute difference was more than the entered DRES, the difference became measured minus observation but need to be corrected for DRES;
- Usage of layer number of 0.0 to get the first active model layers for measurements; needed for unconfined simulations. Also, if a layer falls dry and the layer number is not 0, iPEST looks for the next one;
- The truncation level for the SVD-Eigenvalue truncation can be specified via SVD_EIGV. The default value is 99.0 (equal to the original value, it might be encouraged to increase this number).
-
IMODBATCH
-
CUS-FUNCTION: added option to specify the keyword MINEXTENT to discard elements with spatial extent less than entered at the keyword;
-
IDFSCALE, improved the option #12 to remove all cells above- and below the entered threshold for numerical robustness.
-
Improved functionalities:
-
PROJECT MANAGER
-
Complete reconstructed dialogs to support all functionalities from the MODELTOOL and the keywords present in the IMODBATCH function RUNFILE;
-
Order of modules and packages are sorted differently to be more logical;
-
Support the usage of MT3D/SEAWAT, packages can be extended with information on concentrations;
-
For the computation of KVA combined with values for FCT and IMP not equal 1.0 and 0.0, respectively, needs to be inversed in the internal computation;
-
Preparation of future usage of MODFLOW6.
-
-
MF2005
-
Always save dates as-is as written in the RUNFILE, those names are written the DIS file, MF2005 is adjusted accordingly to use those;
-
MNW package automatically correct for filters above- and beyond the vertical limits of the model (minimal 0.10 m);
-
MNW package does not support the automatic correction for partial penetration due to numerical issues in combination with unconfined, dry layers, PPFLAG=0 instead of PPFLAG=1.
-
-
IMODBATCH
-
WBALANCE TOOL-FUNCTION, usage does not give a pop-up window after completion, also output folder is now created;
-
IDFTIMESERIES-FUNCTION, modified such that timeseries are added as second column in the associated TXT files instead of first column. This is done as iPEST needs measurements to be at the first column in those files always.
-
-
MetaSWAP
A new version of MetaSWAP is used. Additional communication will be sent to Dutch users how this new version can be used. (download here the Dutch Release Notes MetaSWAP_iMOD5).
Bugfixes:
-
PARTICLE TRACKING
Usage of the options to repeatedly apply fluxes in a transient model is improved to cope with simulation windows that are shorter than the entered budget files.
-
ROTOOL
Made this functionality compatible with the double precision IDF file format.
-
MF2005
-
MNW package of MF2005 had a bug in determining the start position of a well in combination with inactive cells above;
-
The timestamps from the RUNFILE and/or during the conversion from the PRJ files to the NAMFILE are now echoed in the DIS-files for MF2005. These are used in the construction of the output names. More steering for this is now possible.
-
-
GENERAL
Conversion of an ASCII GEN to a BINARY GEN crashed if the number of columns in the associated DAT file was not equal for all rows.
-
TIME SERIES TOOL
Filenames with {name}{date}{ext} are recognized by the iMOD Time Series Tool, before only {name}_{date}_{text}.{ext} were recognized, e.g. RCH_20110101.IDF is now applicable as it was not before.
-
ISGEDIT
There was a crash if two ISG files were present for which one ISG was active in ISG Edit and the other was used to zoom to full extent.
-
IMPORTSOBEK
Aggregation to daily values has been improved.
-
IMODBATCH
-
CUS-FUNCTION determining the WINDOW from given IDF files instead of specifying the WINDOW priory was incorrect;
-
IDFMEAN function was not working correct for those configurations where ILAYER was defined;
-
RUNFILE function to write MF2005 files the combination KVA with FCT and/or IMP values not equal 1.0 and 0.0, respectively, needs to be inverse.
-
-
CROSS-SECTION TOOL
Display of boreholes and/or seismic-measurements went wrong on last point of a cross-section while snapping the points to the boreholes/seismic locations.
Ad 2) Developments:
In June 2020 Deltares will release iMOD 5.1, which enables more advanced usage of the computational kernel of MODFLOW 6 within iMOD.
In December 2020 Deltares will release iMOD 5.2, with improved functionalities of the use of MODFLOW 6 within iMOD.