iMOD User Manual version 5.2 (html)

A.3Grid2MetaSWAP manual


The Grid2MetaSWAP program is incorporated in the iMODFLOW and iMOD framework to be able to couple iMODFLOW with MetaSWAP on a grid level. In essence, the program combines table information, grid data and predefined conditions into information readable for MetaSWAP. The following paragraphs describe the structure of Grid2MetaSWAP, scaling methods used within, and translation of all available data to MetaSWAP input files.

A.3.2Structure of the Grid2MetaSWAP program

As already mentioned briefly, the Grid2MetaSWAP program couples the MetaSWAP information to an iMODFLOW grid cell. Such a grid cell can be covered by, a selection of, three coverage types; surface water, paved area (urban and roads) and unpaved areas. Usually surface water is not explicitly taken into account by MetaSWAP but is an important parameter for iMODFLOW. In case of the urban (paved) areas some specific characteristics are linked to the grid cell, like soil intrusion, evapotranspiration and runoff to the sewage system. Grid2MetaSWAP reads the listed files and constant values from the CAP block in the iMODFLOW RUNFILE. As shown in Figure A.2, in total 18 parameters can be distinguished of which 4 are split up in a parameter for “urban” and one for “rural” (see also Table 1); ponding depth, runoff resistance, runon resistance, and infiltration capacity. An extensive description of these parameters is given in Chapter 10 of the iMOD manual (“Data set 10: Number of files”, paragraph 10.11). Besides these compulsory parameters that are used to create the needed MetaSWAP input files, additional input files can be added to the RUNFILE. After Grid2MetaSWAP created the coupling files and necessary input files, it makes a copy of the additional input files and places them in the MetaSWAP model folder . An exception on this is made for the files "para_sim.inp" and "mete_grid.inp". These files are copied and rewritten based on certain conditions; this is explained in paragraph "Parameter processing and adjustments".


Figure A.2: Overview of the Grid2MetaSWAP general structure and workflows

Table A.1 gives an overview of the compulsory iMODFLOW-MetaSWAP parameters, the used abbreviations in the iMOD manual and examples of file names or default values to be used in the iMODFLOW RUNFILE. The purpose of a value that is given for a certain parameter may differ. A value of "999.0" in case of the artificial recharge capacity means that this variable is turned off during the computation. Similar for ponding depth, a value of 999.0 [m] means that the overland flow will only take place when groundwater level exceeds this; which in The Netherlands will never happen and the ponding depth is technically turned off. A default value for other parameters is set to -9999.0. The soil moisture factor and conductivity factor are set to 1.0 on default. For more information about the specific parameters and what type of values or files to be used, the iMOD manual can be consulted.

Table A.1: Overview of iMODFLOW-MetaSWAP parameters

iMODFLOW parameter Beregen_laag.idf or beregen_locatie.ipf, VXMU_SOPP.idf, VXMU_ROPP.idf, SOF_HYDRO.IDF or
abbreviation iMODFLOW parameter Example file / default value
BND Boundary setting ibound_l1.idf
LGN Landuse code LGN250.idf
RTZ Rootzone thinkness Rootzone.idf
SFU Soil physical unit Bofek2012.idf
MET Meteo station number Metestat.idf
SEV Surface elevation AHN.idf
ART Artificial recharge type Beregen.idf
ARL Artificial recharge layer default=1.0
ARC Artificial recharge capacity Default = 999.0
WTA Wetted area NOPP.idf
UBA Urban area SOPP.idf
PDU + PDR Ponding depth default = 999.0
OFU + OFR Runoff resistance Default = 1.0
ONU + ONR Runon resistance Default = 1.0 or -9999 .0
QUI + QIR Infiltration capacity QUI =0.0 and QIR=0.5
PWT Perched water table -9999.0
SMF Soil moisture factor Default = 1.0
CFC Conductivity factor Default = 1.0

As briefly mentioned before, seen from the Grid2MetaSWAp code a distinction is made between urban and rural areas. Grid2MetaSWAP processes the parameters based on this division. However, only one of each MetaSWAP table input file (with extension *.inp) is created, containing both data for the urban and for the rural area as shown by Figure A.3. Within the input file set 3 input files are written with specific coupling information; 2 related to the meteorological forcing (precipitation and evapotranspiration grids) and 1 to couple the iMODFLOW cells to the MetaSWAP units.


Figure A.3: MetaSWAP input files created by Grid2MetaSWAP, split up in 3 categories: 1. urban, 2. rural and 3. coupling

A.3.3Creation of MetaSWAP input files and content

As already shown in Figure A.3, 11 MetaSWAP input files are created by the Grid2MetaSWAP program. From Table A.2 it can be read which MetaSWAP parameters are included in which input file. In case of the meteorological coupling files ("svat2precgrid.inp" and "svat2etrefgrid.inp") information from the wetted area and urban area parameters are used to fill these coupling files. In the meteorological coupling files, each svat number is coupled to the row and column numbers of the grids in case there is urban or rural area present in the cell. In the iMODFLOW coupling file ("mod2svat.inp") the iMODFLOW cell number is coupled to the svat number and combined with the artificial recharge layer of the specific cell.

Table A.2: Overview of MetaSWAP parameters included or used for the creation of the MetaSWAP input files


In case of the artificial recharge (only an option for cells with rural area) that is stored in the ’scap_svat.inp’ a couple of choices are made in to code before storage. The flow diagram in Figure A.4 shows the possibility of using a grid file (IDF-file) as input or a point file (IPF-file) containing the necessary information. The IPF-file only can contain information about artificial recharge from a groundwater resource.


Figure A.4: Flow diagram of the creation of "scap_svat.inp" containing artificial recharge information


Each grid file given in the iMODFLOW RUNFILE is scaled towards the right extent using a standardized routine. Different scaling methods are used depending on the type of parameter as can be seen in Table A.3. There are 4 upscaling methods used and 2 downscaling methods.

Upscaling methods (blue coloured in Table A.3):

Downscaling methods (orange coloured in Table A.3):

Table A.3: Overview of upscale (blue) and downscale (orange) methods used for rescaling input grid files


A.3.5Parameter processing and adjustments

Within the Grid2MetaSWAP program code a number of adjustments are made to be able to create MetaSWAP input files that are readable within the iMODFLOW-MetaSWAP framework. In all cases the following conditions are processed:

When looking at cells that are completely our partly covered by urban area, the following specific conditions are taken into account:

A similar set of conditions is in place for the cells containing rural area:

As given in Figure A.2, the "para_sim.inp" and the "mete_grid.inp" are rewritten. Changes in the "para_sim.inp" include:

Nothing is changed in the "mete_grid.inp". Related to the settings in de "para_sim.inp" and "mete_grid.inp", you need to keep in mind that the following settings are not touched:

Finally, some tips for running you model with the iMODFLOW-MetaSWAP framework: