pictures/imod_logo.png

iMOD User Manual version 5.2 (html)


8.11POSTPROCESSING-FUNCTIONS


8.11.1GXG-Function

The GXG function calculates the maximum and minimum groundwaterhead during the hydrological year (from 1\({}^{st}\) April) based on the mean of the three highest, c.q. lowest observed groundwaterheads. The GXG is an indicator used in the Netherlands to indicate the seasonal variation of the groundwaterhead.

FUNCTION=

GXG

ILAYER=

Enter the layer numbers to be used in the GxG computation,
subsequently; e.g. ILAYER=1,3,6.

NDIR=

Enter the number of folders to be processed repeatedly, e.g. NDIR=10.

SOURCEDIR{i}=

Enter the folder and first part of the filename for all files that need to be used, e.g. SOURCEDIR1=C:\DATA\HEAD. This mean that the GXG function will search for IDF-files that meet the name syntax requirement of C:\DATA\HEAD_{yyyymmdd}_L{ILAY}.IDF.

SURFACEIDF=

Enter a name for the surface, e.g. SURFACEIDF=C:\DATA\AHN.IDF. By default a surface elevation of 0.0m+MSL will be considered.

SYEAR=

Enter the start year (yyyy) for which IDF-files are used, e.g. SYEAR=1998.

EYEAR=

Enter the end year (yyyy) for which IDF-files are used, e.g. EYEAR=2011. This keyword will be read whenever SYEAR is included.

IYEAR=

Specify particular years to be used, e.g. IYEAR=2001,2003,2004. This keyword will be read whenever the keyword SYEAR is absent.

STARTMONTH=

Enter the start month from the which the hydrological year starts. Default STARTMONTH=4.

IPERIOD=

Enter two integers (0 or 1) for each month to express the inclusion of the 14\({}^{th}\) and 28\({}^{th}\) of that particular month, e.g. IPERIOD=010101010101010101010101, which mean to use the 14\({}^{th}\) of each month solely. On default IPERIOD=111111111111111111111111.

ISEL=

Enter a code for the area to be processed:
ISEL=1 will compute the entire region
ISEL=2 will compute within given polygons;
ISEL=3 will compute for those cells in the given IDF-file that are not equal to the NoDataValue of that IDF-file.

GENFNAME=

Enter a GEN-filename for polygon(s) for which mean values need to be computed. This keyword is obliged whenever ISEL=2.

IDFNAME=

Enter an IDF-file for which mean values will be computed for those cell in the IDF-file that are not equal to the NoDataValue of that IDF-file. This keyword is compulsory whenever ISEL=3

HGLG3=

Indicates whether or not the HG3 and LG3 needs to be computed output per year. HG3(LG3)=the groundwater stage maps of the 3 days with the highest(lowest) stages per year. Note: Only the first layer indicated by keyword ILAYER, will be taken into account. Its not possible to get the HG3 and LG3 for more layers within one GXG-calculation. Note 2: in case of SYEAR=2011, EYEAR=2013 and startmonth=4, the HG3 and LG3 are computed for year 2011 (=04-2011 till 03-2012) and 2012 (=04-2012 till 03-2013).

Example 1

FUNCTION=GXG
ILAYER=1
NDIR=1
SOURCEDIR1=D:\MODEL\HEAD
SYEAR=1991
EYEAR=2000

This illustrates a simple example of a GxG computation over the years 1991 (actually starts at 14\({}^{th}\) of April 1991) until 2000 (actually 28\({}^{th}\) of March 2000), for all the HEAD* files that are within the folder D:\MODEL. Since the keyword SURFACEIDF is absent, the GxG will be expressed according to 0.0 instead of a true surface level, moreover, ILAY is absent too, but ILAY=1 will be used as default.

Example 2

FUNCTION=GXG
ILAYER=1,2
SURFACEIDF=D:\DATA\AHN.IDF
IYEAR=1994,1995,2000,2001
IPERIOD=000000001111111100000000
ISEL=3
IDFNAME=D:\DATA\ZONE.IDF
NDIR=1
SOURCEDIR1=D:\MODEL\HEAD

This example computes the GxG for the years 1994, 1995, 2000 and 2001 only. This means two hydrological years, namely 14-4-1994-upto 28-3-1995 and 14-4-2000 upto 28-3-2001. In this period the summer months May, June, July, August are included as expressed by the IPERIOD keyword.


8.11.2WBALANCE-Function

The WBALANCE function calculates the water balance based on the model output for the steady state condition or for a specific period and area. The result is a CSV file (Step 1).
Alternatively, this function can create images, IDF files and/or CVS files from aggregation on existing CSV files (Step 2).

FUNCTION=

WBALANCE


Step 1: Use the following keywords to create a waterbalance file

NDIR=

Enter the number of folders to be processed repeatedly, e.g. NDIR=10.

SOURCEDIR{i}=

Enter the main source folder for which underlying files need to be used;
e.g. SOURCEDIR1=C:\DATA\MODEL. It depends on the following keywords: BAL{i}, BAL{i}ISYS, ILAYERS and SDATE/EDATE, what specific files the WBALANCE function will obtain.

ILAYER=

Enter the layer numbers to be included in the waterbalance,
e.g. ILAYER=1,3,6. It is also possible to specify the layers as LAYERS=4:45 to indicate that these layers 4 up to 45 need to used.

SDATE=
(optional)

Enter the starting date (yyyymmdd or yyyymmddhhmmss) for which IDF-files are used, e.g. SDATE=19980201 or SDATE=20141231123015 (the latter expresses the 31\(^{\rm th}\) of December 2014 at 12 hours, 30 minutes and 15 seconds.

EDATE=
(optional)

Enter the ending date (yyyymmdd or yyyymmddhhmmss) for which IDF-files are used, e.g. EDATE=20111231 or EDATE=20180715084500 (the latter expresses the 15\(^{\rm th}\) of July 2018 at 08 hours, 45 minutes and 00 seconds. In case SDATE is specified, EDATE is compulsory as well.

IYEAR=
(optional)

Specify a particular year (within SDATE and EDATE) to be used exclusively, e.g. IYEAR=2001,2003,2005. IYEAR is filled in for all years in-between SYEAR and EYEAR.

NPERIOD=
(optional)

Enter a number of periods to be defined to use IDF-file within these periods solely, e.g. NPERIOD=2. NPERIOD=0 by default.

PERIOD{i}=
(optional)

Enter a period i (ddmm-ddmm), e.g. PERIOD1=1503-3110 to express the period 15\({}^{th}\) of March until the 31\({}^{th}\) of October.

NBAL=

Enter the number of water balance topics, e.g. NBAL=2.

BAL{i}=

Enter for each of NBAL topics one of the folder name, e.g. BDGBND, BDGSFR. iMOD will look for files that are in the folder SOURCEDIR{i}\BAL{i}\BAL{i}_{time}_{layer}.IDF. whenever SDATE/EDATE is absent, {time} is the keyword STEADY_STATE. Repeat BAL{i} for NBAL times. E.g., whenever BAL1=BDGWEL and the simulation is steady state, the following file is appropriate: SOURCEDIR{i}\BDGWEL\BDGWEL_STEADY-STATE_L1.IDF. For transient simulations, the iMOD Batch function will search for SOURCEDIR{i}\BDGWEL\BDGWEL_????????_L1.IDF and SOURCEDIR{i}\BDGWEL\BDGWEL_??????????????_L1.IDF files. After that, the selected set of files will be matched against the given SDATE and EDATE keywords, and if necessary against NPERIOD.

BAL{i}ISYS=
(optional)

Enter the number of systems to be included in the water balance, e.g. BAL1ISYS=1,2,3. This mean to add the systems 1,2 and 3 for the first entered water balance item, specified by BAL1. E.g. iMOD will look for files as BDGDRN_SYS1_STEADY-STATE_L1.IDF. By default, no systems will be distinguished and iMOD will look for files as BDGDRN_STEADY-STATE_L1.IDF.

OUTPUTNAME{i}=

Enter the output filename (*.TXT, *.CSV or *.IPF), e.g. OUTPUTNAME1=C:\DATA\HEAD\WBAL_MIPWA.CSV

ISEL=
(optional)

Enter a code for the area to be processed:
ISEL=1 will compute the entire region
ISEL=2 will compute within given polygons;
ISEL=3 will compute for those cells in the given IDF-file that are not equal to the NoDataValue of that IDF-file. By default ISEL=1.

GENFILE=
(optional)

Enter a GEN-filename for polygon(s) for which mean values need to be computed. This keyword is obliged whenever ISEL=2.

IDFNAME=
(optional)

Enter an IDF-file for which mean values will be computed for those cell in the IDF-file that are not equal to the NoDataValue of that IDF-file. This keyword is compulsory whenever ISEL=3

WBEX=
optional)

Enter WBEX=1 to generate interconnected flux between the zones. This option is only valid whenever the flux terms BDGFRF and BDGFFF are active. By default WBEX=0 and none of the interconnected fluxes is computed.


Step 2: use the following keywords whenever CSVFNAME is available

CSVFNAME=
(optional)

Enter the name of the CSV file that contains a water balance created previously by this function, e.g. CSVFNAME=D:\MODELRESULTS\WBAL.CSV.

DIR=

Enter the name of the output folder that will be used to save the resulting pictures, e.g. DIR=D:\MODELRESULTS\FIGURES.

BDGIACT=
(optional)

Enter BDGIACT=1 to denote that the budget is present in the output file (IDF, CSV and/or time series/graphical representation image). By default BDGIACT=1 for all budget terms.

BDGIGRP=
(optional)

Enter the group numbers for the individual budget terms, e.g. BDGIGRP=1,2,5,5. By default BDGIGRP is a sequence of the numbers 1,2,3 etc.

BDGICLR=
(optional)

Enter the colour number (combination of the individual colour red,green en blue - each ranging between 0-255) of the particular budget term, BDGICLR=2443254.

IAVG=
(optional)

Enter the method for averaging:

  • • 1 = All Time steps: a single value for all entries;

  • • 2 = Year: a single value per year;

  • • 3 = Months: a single value per month, starting at the first month in the series;

  • • 4 = Hydrological Seasons: two seasons are used, 1) April - September and 2) October - March;

  • • 5 = Decade: a single value per 10 days and the remaining days in that month;

  • • 6 = Hydrological Year: four seasons are used 1) December - February 2) March - May 3) June - August 4) September - November;

  • • 7 = Quarters: a single value per 3 months, starting in January;

  • • 8 = None: all entries remain unchanged, this is the default.

NETFLUX=
(optional)

Enter NETFLUX=1 to apply net fluxes in the output, by default NETFLUX=0.

IUNIT=
(optional)

Enter INUIT=1 to apply net fluxes in the output, by default IUNIT=0.

LSUM=
(optional)

Enter LSUM=1 to aggregate all selected model layers into a single water balance, by default LSUM=0.

ZSUM=
(optional)

Enter ZSUM=1 to aggregate all selected zones into a single water balance, by default ZSUM=0.

LAYERS=
(optional)

Enter the number of layer(s) to be used for the water balance, e.g. LAYERS=1,2,5. By default ALL layers are selected.

ZONES=
(optional)

Enter the number of zones(s) to be used for the water balance, e.g. ZONES=10,23. By default ALL zones are selected.

DATES=
(optional)

Enter the number of dates(s) to be used for the water balance, e.g. DATES=20100114,20100128. By default ALL dates are selected.

IOPT=
(optional)

Enter the output option:

  • • 1 = Time Series: select this option to display the selected water balance items in a graph, this is the default;

  • • 2 = Graphical Representation: select this option to present the water balance items in a illustrative image;

  • • 3 = Preview Table: select this option to display each value for the water balance items in a table;

  • • 4 = Export to CSV: select this option to export all water balance items into an CSV file;

  • • 5 = IDF per Layer: select this option to export all water balance items into separate IDF files.

OUTPUTFNAME=
(optional)

Enter the name of the CSV file to be created whenever IOPT=4, e.g. OUTPUTFNAME=D:\OUTPUT\SUMMARY.CSV.

Example 1

FUNCTION=WBALANCE
NBAL=3
BAL1=BDGFRF
BAL2=BDGFFF
BAL3=BDGFLF
ILAYER=3
NDIR=1
ISEL=2
GENFILE=D:\MODEL\zone.gen
SOURCEDIR1=D:\MODEL
OUTPUTNAME1=D:\MODEL\WBAL.TXT

The above mentioned simple example will give a waterbalance for the BDGFRF, BDGFFF and BDGFLF, respectively, for modellayer 3.
The IDF-files will be D:\MODEL\BDGFRF\BDGFRF_STEADY-STATE_L3.IDF;
D:\MODEL\BDGFFF\BDGFFF_STEADY-STATE_L3.IDF; and
D:\MODEL\BDGFLF\BDGFLF_STEADY-STATE_L3.IDF.
The result is written in D:\MODEL\WBAL.TXT.

Example 2

FUNCTION=WBALANCE
NBAL=2
BAL1=BDGRIV
BAL1ISYS=1,2
BAL2=BDGDRN
ILAYER=1,2
SDATE=19900101
EDATE=20001231
IYEAR=1990,1995,1997,2000
NPERIOD=1
PERIOD1=0104-3107
ISEL=2
GENFNAME=D:\DATA\ZONES.GEN
NDIR=2
SOURCEDIR1=D:\MODEL
SOURCEDIR2=D:\SCENARIO
OUTPUTNAME1=D:\OUTPUT\WBAL_MODEL.CSV
OUTPUTNAME2=D:\OUTPUT\WBAL_SCENARIO.CSV

The example above will compute a waterbalance for two modellayers (1,2) for the budgetfiles BDGRIV*SYS1 and BDGRIV*SYS2 and BDGDRN in the period from 1\({}^{th}\) of April until the 31\({}^{th}\) of July for the years 1990,1995,1997,2000. The waterbalance will be summed for the zones that are described by the polygon(s) inside the file ZONES.GEN. Finally, the computation will be executed twice, for those results in D:\MODEL and those in D:\SCENARIO. Results are stored in the folder D:\OUTPUT.


8.11.3MSPNETRCH-Function

The MSPNETRCH function calculates the net Recharch from the MetaSWAP compartment to the watertable of Modflow for a specific period and area. The result is a series of IDF files. The following functions is used to calculate the net recharch.

\begin{equation} NetRCH = \Delta lvgmodf * sc1 - qmodf, \end{equation}

whereby lvgwmodf is the head calculated by Modflow and sc1 a long term calculated mean of the Storage Coefficient. Qmodf is the sum of all MODFLOW stresses on groundwater.
For time step t the \(\Delta \) lvgwmodf is defined as lvgmodf (t) - lvgmodf (t-1).

Note: It is obligatory to provide MSPNETRCH the necessary files for lvgwmodf and qmodf on a daily basis.
Note: To get model results for lvgwmodf and qmodf switch on these keywords in the model input file sel_key_svat_per.inp.

FUNCTION=

MSPNETRCH

SOURCEDIR=

Enter the folder name referring to the model result directory. Underlying MetaSWAP files will to be used;
e.g. SOURCEDIR=C:\IMOD_USER\MODELS\ISLAND.

OUTPUTDIR=
(optional)

Enter the name of the output folder;
e.g. OUTPUTDIR=C:\IMOD_USER\MODELS\ISLAND\METASWAP\MSPNETRCH (by default)

SDATE=
(optional)

Enter the starting date (yyyymmdd) for which IDF-files are used, e.g. SDATE=19980201

EDATE=
(optional)

Enter the ending date (yyyymmdd) for which IDF-files are used, e.g. EDATE=20111231. In case SDATE is specified, EDATE is compulsory as well.

STOAVG=
(optional)

Specify an IDF with a given mean Storage Coefficient, e.g. STOAVG=C:\IMOD_USER\DBASE\STO\STOAVG.IDF. The dimensions of STOAVG (window and cell size) should be the same as the model dimensions. If this is not the case, the calculation is cancelled and a message is given. Automatic rescaling is not included.

Example *.INI file content

FUNCTION=MSPNETRCH
SOURCEDIR=C:\IMOD_USER\MODELS\ISLAND
OUTPUTDIR=C:\IMOD_USER\MODELS\ISLAND\METASWAP\MSPNETRCH

The above given example will calculate the net Recharge IDFs for the whole model period. The IDF-files will be written as C:\IMOD_USER \MODELS \ISLAND \METASWAP \MSPNETRCH \MSPNETRCH_yyyymmdd.IDF.
The MSPNETRCH is calculated using the \(\Delta \) lvgwmodf between 2 time steps (t and t-1). As a result the number of MSPNETRCH file is the number of provided time steps minus 1.


8.11.4MF6TOIDF-Function

Use this iMODFLOW post-processing function to convert standard MODFLOW6 output to IDF files.

FUNCTION=

MF6TOIDF

SDATE=
(optional)

Enter an initial date of the model, e.g. 20121231 to denote the 31\(^{\rm st}\) of December 2012. This is used to include a date string to the converted files, e.g. HEAD_20121231_L1.IDF. MF6 saves the total elapsed time in the output files, these will be added to the SDATE. The first stress period is assumed to be a steady-state and converted files for that stress period are saved as e.g. HEAD_STEADY-STATE_L1.IDF.

IDF=
(optional)

Enter an IDF file for which spatial definition is used to convert the results from the MF6 binary files into IDF-format. Whenever IDF is specified, it is only possible to convert HED-files from MF6 and CBC-file obligates the usage of a GRB-file.

GRB=
(optional)

Specify the GRB of MF6. Whenever the model has been created by iMOD, the GRB file is saved, per (sub)model, in the folder .\GWF_1 \MODELINPUT \{MODELNAME}.DIS6.GRB. This is a binary file that is not editable, it contains all spatial information for the current (sub)model.

HED=
(optional)

Specify the output file (*.HED) for MODFLOW6 that contains the HEAD results, e.g. HED=.\GWF_1 \MODELOUTPUT \HEAD \HEAD.HED. At least HED of BDG need to be specified. Additional two files are created, a PHR.IDF and SUM.IDF. This first is the phreatic-head, that is the value for the first active model layer. Secondly, SUM.IDF represents the total number of active heads found per location.

BDG=
(optional)

Specify the output file (*.CBC) for MODFLOW6 that contains the BUDGET results (this includes all budget from all packages including the cell-by-cell budget fluxes), e.g. BDG=.\GWF_1 \MODELOUTPUT \BUDGET \BUDGET.CBC. At least HED of BDG need to be specified.

SAVESHD=
SAVECHD=
SAVEWEL=
SAVEDRN=
SAVERIV=
SAVEGHB=
SAVERCH=
SAVEEVT=
SAVESTO=
SAVESPY=
(all optional)

Use these keywords to save the hydraulic head per layer or/and results for the CHD, WEL, DRN, RIV, GHB, RCH, EVT, STO and SPY package (latter two for transient models), e.g. SAVESHD=3,4,10 to note that model layers 3, 4 and 10 will be saved only, by default all keyword are 0, meaning no layers will be saved. Specify SAVESHD=-1 to denote that ALL layers will be saved, this is similar for the other packages.

SAVEFLX=
(optional)

Use this keyword to include layers to be saved for the spatial fluxes in x,y and z direction, e.g. SAVEFLX=3,4,10 to note that model layers 3, 4 and 10 will be saved only, by default SAVEFLX=0, meaning no layers will be saved. Specify SAVEFLX=-1 to denote that ALL layers will be saved. Part of this, the BDGFFF, BDGFRF and BDGFLF will be saved.

Example 1

FUNCTION=MF6TOIDF
SDATE=20121231
GRB=.\GWF_1 \MODELINPUT \{MODELNAME}.DIS6.GRB
HED=\GWF_1 \MODELOUTPUT \HEAD \HEAD.HED
BDG=.\GWF_1 \MODELOUTPUT \BUDGET \BUDGET.CBC
SAVESHD=-1
SAVEFLX=-1
SAVERIV=1
SAVEWEL=8,12

The above mentioned examples creates IDF files for all HEADS in the HED-file as well as all budget files for all items in the CBC-file. Multiple folders can be created for BDGRIV, BDGDRN, BDGFLF etc. if these are present in the CBC-file.

Example 2

FUNCTION=MF6TOIDF
IDF=.\GWF_1 \MODELINPUT \DIS \BND_L1.IDF
HED=\GWF_1 \MODELOUTPUT \HEAD \HEAD.HED

This example creates HEAD*.IDF files from the selected HED-file and using the IDF files from that submodel.


8.11.5IDFTIMESERIE-Function

Use this function to generate timeseries out of IDF-files that have the notation {item}_yyyymmdd_l{ilay}.idf. These are IDF-files that yield from a normal iMODFLOW simulation.

FUNCTION=

IDFTIMESERIE

IPF1=

Enter the name of an IPF file that contains the locations of the measurements, e.g. IPF1=D:\DATA\MEASUREMENTS_LOC.IPF.

IPF2=

Enter the name of an IPF file that will be used to store the computed time series, e.g. IPF2=D:\IMOD\HEAD_MODEL_MEASUREMENTS.IPF.

ILAY=

Enter the modellayer, e.g. ILAY=2.

SOURCEDIR=

Enter the directory name of the folder that contains the specific files + the first (similar) part of the name of the files, e.g. SOURCEDIR=D:\RESULT\HEAD\HEAD. This will yield IDF-files that belong to D:\RESULT\HEAD\HEAD_{yyyymmdd}_L{ilay}.IDF.

SDATE=
(optional)

Enter the start date of the time series to be computed, e.g. SDATE=19700803000000 to express the 3\({}^{\rm rd}\) of August 1970. By default SDATE=0 and will not be used, the series starts at the earliest file that confirms the SOURCEDIR.

EDATE=
(optional)

Enter the end date of the time series to be computed, e.g. SDATE=20120601133015 to express the 1\({}^{\rm st}\) of June 2012 at 13 hours, 30 minutes and 15 seconds. By default EDATE=0 and will not be used, the series ends at the latest file that confirms the SOURCEDIR.

LABELCOL=
(optional)

Enter the column to be used for labelling the associated text files. Default LABELCOL=0 and will not be used.

TXTCOL=
(optional)

Enter the column to be used from the associated text files. Default LABELCOL=2 and it takes the second column from associated TXT-files.

NOTE: if IPF1 is an IPF with time series (e.g. in associated TXT files) than be aware: if the new IPF2 is in the same folder, than the original associated TXT files will be overwritten.

Example 1

FUNCTION=IDFTIMESERIE
IPF1=D:\DATA\MEASUREMENTS_LOC.IPF
IPF2=D:\POST_PROCESSING\HEAD_MODEL_MEASUREMENTS.IPF
ILAY=1
SDATE=19500101000000
EDATE=20120101123005
SOURCEDIR=D:\RESULT\HEAD\HEAD

The example above will yield time series from the results in D:\RESULT\HEAD\HEAD_*_L1.IDF for the period between the 1\({}^{\rm st}\) of January 1950 and the 1\({}^{\rm st}\) of January 2012 at 12 hours, 30 minutes and 5 seconds.

Example 2

FUNCTION=IDFTIMESERIE
IPF1=D:\DATA\MEASUREMENTS_LOC.IPF
IPF2=D:\POST_PROCESSING\HEAD_MODEL_MEASUREMENTS.IPF
ILAY=1
SDATE=19500101000000
EDATE=20120101000000
SOURCEDIR=D:\RESULT\HEAD\HEAD
LABELCOL=3

The example above differs for LABELCOL only. The 3\({}^{\rm rd}\) column of IPF1 will be used to generate the name of the text file that stores the time series, most logically this is the ID name of the monitoring location.


8.11.6IPFRESIDUAL-Function

Use this function to calculate residuals based on IPF files containing statistics.

FUNCTION=

IPFRESIDUAL

NIPF=

Enter the number of IPF files to be handled.

IPFFILE{i}=

Enter the name of the i\({}^{th}\) IPF file containing measurement and computed values, or refer to associated TXT files with time series of measurements and computed values.

OUTNAME=

Enter the name of the output filename, the statistics will be written in this file. Also, per IZONE and per LAYER, a separate IPF file will be created called residual_ILAY{\(i\)}_IZONE{\(j\)} .

ILCOL{i}=
(optional)

Enter the column number in the IPF{i} that represents the model layer, by default ILCOL{i}=3.

IMCOL{i}=
(optional)

Enter the column number in the IPF{i} that represents the measurement. In the case that the measurement is given by associated TXT files, it is not necessary to enter this keyword, in the other case the default value is IMCOL{i}=3.

IHCOL{i}=
(optional)

Enter the column number in the IPF{i} that represents the computed head. In the case that the computed head is given by associated TXT files, it is not necessary to enter this keyword, in the other case the default value is IHCOL{i}=3.

W_TYPE{i}=
(optional)

Give whether the load is entered as variance (W_TYPE1=1) or weights (W_TYPE1=2). Whenever variances are entered, a weight is computed internally as: \(w=\frac {1}{\sqrt {v}}\) in which \(w\) is the weight and \(v\) is the variance. By absent of the keyword, so NO weight are expected in that case (in fact the weight will be 1.0 for all locations).

IWCOL{i}=
(optional)

Enter the column number in the IPF{i} that represents the variance or weight, by default IWCOL{i}=3 and is only read whenever W_TYPE{I} \(\ne \) 0.0.

SDATE=
(optional)

Enter the starting date for which statistics need to be gathered, e.g. SDATE=20020726 to express the 26\(^{\rm th}\) of July 2002. By default SDATE=-10\(^{10}\).

EDATE=
(optional)

Enter the starting date for which statistics need to be gathered, e.g. EDATE=20041226 to express the 26\(^{\rm th}\) of December 2004. By default EDATE=10\(^{10}\).

POINTERIDF=
(optional)

Enter the name of an IDF file containing zones, those zones will be used to distinguish in different statistics per zone.

NZONE=

Enter the number of zones in the IDF file to be used, e.g. NZONE=3 will use three zone form the IDF files mentioned at POINTERIDF.

IZONE{i}=

Enter the value in the IDF file for the i\({}^{th}\) zone, e.g. IZONE1=10 denotes that the first zone will be number 10 from the POINTERIDF.

ICOLLECT=
(optional)

Specify ICOLLECT=1 to add the associated TXT files to the created IPF files, by default ICOLLECT=0.

HNODATA=
(optional)

NodataValue of the computed head, by default HNODATA=-999.99. This value will be only used for IPF files with W_TYPE{i}=0 .

Example

FUNCTION=IPFRESIDUAL
NIPF=2
POINTERIDF=D:\MODEL\POINTER.IDF
IPF1=D:\MODEL\HEAD_TSERIES_1.IPF
IPF2=D:\MODEL\HEAD_TSERIES_2.IPF
W_TYPE2=1
NZONE=2
IZONE1=5
IZONE2=6
OUTNAME=D:\RESULTS\STATISTICS.TXT

The example above will give a file containing all the residual values per given IPF-file based on the predefined conditions.


8.11.7PLOTRESIDUAL-Function

Use this function to make a scatter plot or histogram plot of the observations and calculated heads and residuals (calculated-observed) from the iPEST-output data file(s).

FUNCTION=

PLOTRESIDUAL

INPUTFILE=

You can either enter the iPEST output text file name or an IPF file with appropriate information, such as x, y, measurement, computed head and weight values. For IPF files with associated text files, you need to specify ITRANSIENT=1 .

IXCOL=
(optional)

Enter the column number in the IPF file for the X-coordinates, by default this IXCOL=1.

IYCOL=
(optional)

Enter the column number in the IPF file for the Y-coordinates, by default this IYCOL=2.

IMCOL=
(optional)

Enter the column number in the IPF file for the measurement, for TRANSIENT=1, it means the column in the associated text file, by default this IMCOL=2.

IHCOL=
(optional)

Enter the column number in the IPF file for the computed head, for TRANSIENT=1, it means the column in the associated text file, by default this IHCOL=3.

IWCOL=
(optional)

Enter the column number in the IPF file for the weight values, by default this IWCOL=0 and the weight values are assumed to be all equal to 1.0.

ILCOL=
(optional)

Enter the column number in the IPF file for the layer values, by default this ILCOL=0 and the layer is assumed to be all equal to 1.

IPLOT=

Choose the preferred plot type (see for examples further this function description):

  • • 1=Scatter plot;

  • • 2=Histogram plot;

  • • 3=IPF file.

BMPNAME=
(IPLOT=1,2)

Give the name of the output plot file (*.BMP. *.PNG, *.JPG or *.PCX), e.g. d:\SCATTERPLOT_LAYER2.BMP.

IPFNAME=
(IPLOT=3)

Give the name of the output IPF file (*.IPF), e.g. d:\RESIDUAL.IPF.

ITRANSIENT=
(optional)

Choose the type of input file you use:

  • • 0=Steady-state input (this is the default value ITRANSIENT=0);

  • • 1=Transient input.

A transient input file contains date values in the first column whenever the INPUTFILE is an iPEST output file, in the case the INPUTFILE is an IPF-file, the measurement and computed heads are expected to be in associated text files.

IAVERAGE= (optional)
By default IAVERAGE=1 and the average values are compared between the observation and the measurement. If IAVERAGE=0, all timesteps from the associated TXT files will be listed in the final figure. The is keyword is only read for ITRANSIENT=1.

ILAYER=
(optional)

Choose the layers you prefer to plot, e.g. ILAYER=1,4,8 plots data points related to layer 1, 4 and 8. On default all layers are plotted.

IIPF=
(optional)

Choose the IPF-file(s) for which the data needs to be plotted, e.g. IIPF=1,3 plots the data points for the first and 3\(^{\rm rd}\) given IPF-file from the input file. This keyword is not used in case the INPUTFILE is an IPF file.

WC1=
(optional)

Specify the lower limit of a weight value to be included in the selection for the statistics, e.g. WC1=1000.0 means that only points with weight value larger than 1000 will be included in the statistics. By default WC1 is absent, so all will be included.

WC2=
(optional)

Specify the upper limit of a weight value to be included in the selection for the statistics, e.g. WC2=5000.0 means that only points with weight value smaller than 5000 will be included in the statistics. By default WC2 is absent, so all will be included.

HCLASSES=
(optional)

Specify the classes for the histogram, by default this keyword is absent and the following classes are used: -10.0E10,-5.0,-4.5,-4.0,…,4.0,4.5,5.0,10.0E10. However, it is possible to specify a user defined class via e.g. HCLASSES=-10.0,-1.0,-0.5,0.5,1.0,10. In this case the first class is for \(> -10.0\) and \(\le -1.0\), the second for \(> -1.0\) - and \(\le -0.5\) and so on.

IWEIGHT=
(optional)

Choose whether you will include the weight factors in the calculation or not. IWEIGHT=0 the weight factor is not included, IWEIGHT=1 the weight factor is included. In case the weight factor equals 0 the related data point will not be plotted.

IXY=
(optional)

Choose whether you will define the minimum and maximum value of the x,y-axes of the scatter plot. Default IXY=0 and x,y-axis range will be based upon dataset.

XMIN=
(optional, when IXY=1 required)

When IXY=1 give minimum value for x-axis of scatter plot.

XMAX=
(optional, when IXY=1 required)

When IXY=1 give maximum value for x-axis of scatter plot.

YMIN=
(optional, when IXY=1 required)

When IXY=1 give minimum value for y-axis of scatter plot.

YMAX=
(optional, when IXY=1 required)

When IXY=1 give maximum value for y-axis of scatter plot.

Example

FUNCTION=PLOTRESIDUAL
INPUTFILE=D:\LOG_PEST_RESIDUAL_10V2.6.66.TXT
IPLOT=1
BMPNAME=d:\SCATTERPLOT.BMP
ITRANSIENT=1
ILAYER=16
IIPF=1
IWEIGHT=1

This example makes a scatter plot with the name SCATTERPLOT.BMP of the transient data available in the LOG_PEST_RESIDUAL_10V2.6.66.TXT file for a selection of layer 16 and the first IPF-file. The plotted values are multiplied with the given weight factor.

Result of PLOTRESIDUAL, left a scatter plot (IPLOT=1) , right a histogram (IPLOT=2)

pictures/h72-end/residual_2.png pictures/h72-end/residual_1.png


8.11.8PWTCOUNT-Function

Use this iMODFLOW post-processing function to count the number of moments where a PWT situation occurs.

FUNCTION=

PWTCOUNT

SDLIDF=

Enter the name of the IDF file containing the layer number of the first Aquitard.

ILAYIDF=

Enter the name of the IDF file containing the deepest modellayers to be processed.

SDATE=

Enter the starting date (yyyymmdd) for which IDF-files are used.

EDATE=

Enter the ending date (yyyymmdd) for which IDF-files are used.

SOURCEDIR{i}=

Enter the folder and wildcard for all files that need to be used.

OUTPUTIDF=

Enter the name of the IDF file that contains calculated sum of moments where PWT situations occur.

Example 1

FUNCTION=PWTCOUNT
SDLIDF=C:\RESULTS2\PWT\SDL_LAYER.IDF
ILAYIDF=C:\RESULTS2\PWT\PWT_LAYER.IDF
SDATE=19980201
EDATE=20111231
SOURCEDIR=C:\RESULTS\PWT\PWT*.IDF
OUTPUTIDF=D:\RESULTS\PWT_COUNT.IDF

The above mentioned examples creates the IDF file PWT_COUNT.IDF based on a timeseries of PWT result files in directory C:\DATA\