pictures/imod_logo.png

iMOD User Manual version 4.4 (html)


8.5IPF-FUNCTIONS


8.5.1IPFSTAT-Function

The IPFSTAT function can be used to perform statistical analyses on time series that are defined in IPF files as associated files.

FUNCTION=

IPFSTAT

IPFFILE_IN

Specify the name for the IPF file to be used, e.g. IPFFILE_IN=D:\DATA \OBS.IPF

IPFFILE_OUT

Specify the name for the IPF file to be created and added with additional statistics of the associated time series, e.g. IPFFILE_OUT=D:\DATA \OBS_STAT.IPF

IVARS

Specify the variable(s) that need to be computed, e.g. IVARS=1,3,5 or IVARS=1; this can be a list of the following options:
IVARS=1
This will add the minimal value of the time series to the IPFFILE_OUT file.
IVARS=2
This will add the maximal value of the time series to the IPFFILE_OUT file.
IVARS=3
This will add the average value of the time series to the IPFFILE_OUT file.
IVARS=4
This will add the average-highest (GHG) and average-lowest (GLG) values of the time series to the IPFFILE_OUT file. It is the total average of the three highest and lowest values per hydrological year for the 14\(^{\rm st}\) and 28\(^{\rm st}\) of each month (starting at 1\(^{\rm st}\) of April).
IVARS=5
This will add the specified percentile value (PERC) of the time series to the IPFFILE_OUT file.

SDATE=

Enter the starting date for which the statistics need to be computed, e.g. SDATE=20010101 which is the 1\(^{\rm st}\) of January 2001. This date is included.

EDATE=

Enter the starting date for which the statistics need to be computed, e.g. EDATE=20151115 which is the 15\(^{\rm st}\) of November 2015. This date is included.

PERC= (IVARS=5)

Specify the percentile(s) which are used whenever IVARS=4, e.g. PERC=50.0 or PERC=10.0,50.0,90.0 which will generate the different percentiles and add them as separate columns to the IPF file.

ITXTFILE= (IVARS=5)

Specify whether an associated TXT file need to be created for the different percentiles as specified at PERC. iMOD creates a TXT file with entries per percentile \(z(10.0)\) is the 10-percentile value) as if they are dates, e.g. the PERC=10.0,50.0,90.0 will become:

3
2,1
DATE,-999
OBS,-999
20010110,\(z(10.0)\)
20010219,\(z(50.0)\)
20010310,\(z(100.0)\)

Note
iMOD does not take into account percentiles that have decimals.

MINMEASURE=
(IVARS=4)

Specify the number of minimal measures per year for a complete year will be included in the computation of the GHG and GLG for IVARS=4.

DIFFDAY=
(IVARS=4)

Specify the number of days for which an existing measure might differ from the prescribed 14\(^{\rm st}\) and 28\(^{\rm st}\) of each month, e.g. DIFFDAY=2 means that an observation in between 12-16 and 26-30 is valid. iMOD will select the one closest to the 14\(^{\rm st}\) and 28\(^{\rm st}\) of each month.

IWEIGHT=
(optional)

Specify whether the total measurements \(n\) that satisfied the SDATE and EDATE entries need to be written as a weight value (IWEIGHT=1 thus \(\sqrt (n)\)) or as absolute value IWEIGHT=0 thus \(n\). The default value is IWEIGHT=0.

IACOL=
(optional)

Specify the column in the associated text file that you like to calculate the statistics for. It is not possible to use the date-column (IACOL=1). The default value is IACOL=2.

Example 1

FUNCTION=IPFSTAT
IPFFILE_IN=D:\TESTS\TEST.IPF
IPFFILE_OUT=D:\TESTS\OUT.IPF
IVARS=3

The example above computes the average groundwaterlevels (or equivalent) that are associated with the IPF file TEST.IPF. The result is stored in OUT.IPF.

Example 2

FUNCTION=IPFSTAT
IPFFILE_IN=D:\TESTS\MEASURE.IPF
IPFFILE_OUT=D:\TESTS\MODEL.IPF
VARIABLES=4,5
PERC=50.0
SDATE=20100101
EDATE=20170101
MINMEASURE=18
DIFFDAY=2
IACOL=3

The example above, computes the GHG and GLG for the measurements in the second column associated with the MEASURE.IPF for the 1\(^{\rm st}\) of January 2010 up to the 1\(^{\rm st}\) of January 2017. It takes years into account that have at least 18 measurement for the appropriates dates of 14\(^{\rm st}\) and 28\(^{\rm st}\) per month, which can differ for maximal 2 days (DIFFDAY=2).


8.5.2IPFSPOTIFY-Function

The IPFSPOTIFY function can be used to spotify geological formations in existing model discretisations. It gives per modellayer the fraction that a geological formation exists.

FUNCTION=

IPFSPOTIFY

IPFFILE_IN
(optional)=

Enter the name of an *.IPF file for which the locations need to be spotified in the selected op geological formations, e.g IPFFILE_IN=D:\WELLS.IPF.

IPFFILE_OUT=

Specify an IPF file for which the fraction per geological formation will be saved, e.g. D:\RESULT.IPF.

IXCOL=

Specify the column number in the IPFFILE_IN that defines the X coordinates, by default IXCOL=1.

IYCOL=

Specify the column number in the IPFFILE_IN that defines the Y coordinates, by default IYCOL=2.

IFCOL=

Specify the column number in the IPFFILE_IN that defines the attribute for which the fraction need to be computed, by default IFCOL=3.

IZ1COL=

Specify the column number in the IPFFILE_IN that defines the top elevation to spotify underneath, by default IZ1COL=4.

IZ2COL=

Specify the column number in the IPFFILE_IN that defines the bottom elevation to spotify above, by default IZ2COL=5.

ILCOL=

Specify the column number in the IPFFILE_IN that defines the modellayer, by default ILCOL=6.

OUTPUT
FOLDER
(optional)=

Specify the output folder in which all IDF files will be saved with the individual fractions per geological formation per model layer, e.g OUTPUTFOLDER=D:\AQUIFER

NLAY=

Specify the number of model layers to be spotified, e.g NLAY=\(10\)

TOP_{\(i\)}=

Specify the top elevation for the \(i^{\rm th}\) model layer, e.g TOP_L2=D:\MODEL\TOP_L2.IDF.

BOT_{\(i\)}=

Specify the bottom elevation for the \(i^{\rm th}\) model layer, e.g BOT_L2=D:\MODEL\BOT_L2.IDF.

FORMTOP=

Specify the folder that stores the TOP elevation of geological formations, e.g. FORMTOP=D:\GEOLOGY\*-T-CK.IDF. All files will be used that fit this wildcard definition.

FORMBOT=

Specify the folder that stores the BOT elevation of geological formations, e.g. FORMBOT=D:\GEOLOGY\*-B-CK.IDF. All files will be used that fit this wildcard definition.

Example 1

FUNCTION=IPFSPOTIFY
IPFFILE_IN=D:\WELLS\TEST.IPF
IPFFILE_OUT=D:\SPOTIFIED\OUT.IPF
NLAY=2
TOP_L1=D:\MODEL\TOP_L1.IDF
BOT_L1=D:\MODEL\BOT_L1.IDF
TOP_L2=D:\MODEL\TOP_L2.IDF
BOT_L2=D:\MODEL\BOT_L2.IDF
REGISTOP=D:\REGIS\*-T-CK.IDF
REGISBOT=D:\REGIS\*-B-CK.IDF

The example above computes the fractions for each location in the IPFFILE_IN of all geological formations in REGISTOP and REGISBOT for each model layer.

Example 2

FUNCTION=IPFSPOTIFY
OUTPUTFOLDER=D:\FRACTIONS\AQUIFER
NLAY=2
TOP_L1=D:\MODEL\TOP_L1.IDF
BOT_L1=D:\MODEL\BOT_L1.IDF
TOP_L2=D:\MODEL\TOP_L2.IDF
BOT_L2=D:\MODEL\BOT_L2.IDF
REGISTOP=D:\REGIS\*-T-CK.IDF
REGISBOT=D:\REGIS\*-B-CK.IDF

The example above, computes the fractions for each cell in the model layers for each geological formation defined by the REGISTOP and REGISBOT keywords, the results are stored in the AQUIFER folder. To spotify aquitards in it is neccessary to switch the top and bottom elevations, e.g.

FUNCTION=IPFSPOTIFY
OUTPUTFOLDER=D:\FRACTIONS\AQUITARD
NLAY=1
TOP_L1=D:\MODEL\BOT_L1.IDF
BOT_L1=D:\MODEL\TOP_L2.IDF
REGISTOP=D:\REGIS\*-T-CK.IDF
REGISBOT=D:\REGIS\*-B-CK.IDF


8.5.3IPFSAMPLE-Function

The function IPFSAMPLE samples IDF-files to add values to the points defined in an IPF file.

FUNCTION=

IPFSAMPLE

IPFFILE_IN=

Enter the name of an IPF file with minimal 2 columns that represents x- and y coordinates, e.g. D:\DATA\MEASURE.IPF.

IPFFILE_OUT=

Enter the name of an IPF file that need to be written with the results of the IDF values from the specified IDFFILE, e.g. D:\DATA\CHECK.IPF. Results read from the IDF-files in SOURCEDIR, will be stored as an extra column in IPFFILE_IN, the label will be identical to the name of the IDF-files.

SOURCEDIR=

Enter the name of an IDF-file that needs to be read by the points specified in the IPF file IPFFILE_IN, e.g. D:\DATA\RESULTS\HEAD.IDF.

IXCOL

Enter the column number in the IPF file IPFFILE_IN that represents the X-coordinate, e.g. IXCOL=4. By default IXCOL=1.

IYCOL

Enter the column number in the IPF file IPFFILE_IN that represents the Y-coordinate, e.g. IYCOL=6. By default IYCOL=2.

IACOL

Enter the column number to enter the sampled data from the IDF files, e.g. IACOL=3 which means that the entered starts at columns 3. By default IACOL=0 which means that the sampled data will be added at the end of the IPF file.

Example 1

FUNCTION=IPFSAMPLE
IPFFILE_IN=D:\WELLS.IPF
IPFFILE_OUT=D:\WELLS_KD.IPF
SOURCEDIR=D:\DATA\KD*.IDF

This example, adds values (columns) to all points in the IPF file WELLS.IPF, with the corresponding values from the KD*.IDF-files in the folder D:\DATA.

Example 2

FUNCTION=IPFSAMPLE
IPFFILE_IN=D:\WELLS.IPF
IPFFILE_OUT=D:\WELLS.IPF
SOURCEDIR=D:\DATA\KD*.IDF
IXCOL=4
IYCOL=3

This example, adds values (columns) to all points in the IPF file WELLS.IPF, with the corresponding values from the KD*.IDF-files in the folder D:\DATA. The x- and y coordinates in the IPF file WELLS.IPF, will be read from the fourth and third column, respectively.


8.5.4IPFEDIT-Function

The function IPFEDIT reorganises an existing IPF file.

FUNCTION=

IPFEDIT

IPFFILE_IN=

Enter the name of an IPF file, e.g. D:\DATA\MEASURE.IPF.

IPFFILE_OUT=

Enter the name of an IPF file that needs to be saved with the selected set of columns, e.g. D:\DATA\CHECK.IPF.

NCOLUMNS=

Enter the number of columns from the IPF for IPFFILE_IN to be saved in the new IPF at IPFFILE_OUT, e.g NCOLUMNS=3.

ICOL{i}

Enter the column number from the IPF file IPFFILE_IN to be saved as the \(i^{\rm th}\) column in the IPFFILE_OUT, e.g. ICOL1=3.

Example 1

FUNCTION=IPFEDIT
IPFFILE_IN=D:\WELLS.IPF
IPFFILE_OUT=D:\WELLS_CLEANED.IPF
NCOLUMNS=3
ICOL1=1
ICOL2=2
ICOL3=8

This example, reduces the original D:\WELLS.IPF for 3 columns determined by the columns 1,2 and 8. The results is saved into D:\WELLS_CLEANED.IPF.


8.5.5IPFSUM-Function

The function IPFSUM generates a list of total volumes per specified year for zone (if needed) form IPF files. This can be used to prior check the IPF files before entering the model.

FUNCTION=

IPFSUM

NFILE=

Enter the number if IPF file to be used, e.g. NFILE=3.

SDATE=

Enter the start date for the summation of volumes, e.g. SDATE=20110101 to denote the 1\(^{\rm th}\) of January 2011.

EDATE=

Enter the end date for the summation of volumes, e.g. SDATE=20171231 to denote the 31\(^{\rm st}\) of December 2017.

DDATE=

Enter the steps in between the SDATE and EDATE for the summation of volumes, e.g.

  • 1. Daily;

  • 2. Weekly;

  • 3. Monthly;

    4. Yearly.

IPFFILE{i}=

Enter the name of i\(^{\rm th}\) IPF file, e.g. IPFFILE_1=D:\DATA\WELL_L1.IPF.

ZONEIDF=
optional

Enter the name of an IDF file that needs to be used for the different areas (zones) for which the volumes need to be aggregated, e.g. ZONEIDF=D:\DATA\ZONES.IDF.

OUTPUTFILE=

Enter the name of the output file in which the results are written, e.g. OUTPUTFILE=D:\DATA\OUTPUT.TXT.

Example 1

FUNCTION=IPFSUM
NFILE=1
IPFFILE_1=D:\DATA\WELL_L1.IPF
SDATE=19500101
EDATE=20191231
DDATE=4
OUTPUTFILE=D:\OUTPUT.TXT

This example, summes all volumes from the IPF file on a yearly basis (DDATE=4) in between the 1\(^{\rm th}\) of January 1950 up to the 31\(^{\rm st}\) of December 2019 and saves the volumes in D:\OUTPUT.TXT.