pictures/imod_logo.png

iMOD User Manual version 4.4 (html)


8.7GEO-FUNCTIONS


8.7.1DINO2IPF-Function

This function will extract from a CSV file exported from DINO (TNO) appropriate data to generate an IPF file with borehole information attached to it. The content of the CSV file is prescribed on the next page.

FUNCTION=

DINO2IPF

CSVFILE=

Enter a CSV file that contains the necessary information from the DINO database, e.g. CSVFILE=D:\DATA\DINO.CSV. The output file (IPF file) will be named after the CSVFILE. Moreover, you can specify a wildcard to transform more CSV files into a single IPF file, e.g. CSVFILE=D:\DATA\*.CSV. In this case you need to specify an IPF filename with the keyword IPFFILE=.

WINDOW=

Specify a window (X1,Y1,X2,Y2) for which the entered RUNFILE will be clipped, WINDOW=125100.0,345000.0,135000.0,355000.0.

GENFILE=

Enter a name for a GEN-file that contains a polygon that determines the area for which the CSV files need to be converted into the IPF files.

IPFFILE=

Enter the name of the IPF file to be created, e.g. IPFFILE=D:\DINO.IPF. This keyword is obliged only whenever the CSVFILE contains a wildcard “*”.

Example 1:

FUNCTION=DINO2IPF
CSVFILE=D:\iMOD-DATA\DINO\*.csv
IPFFILE=D:\iMOD-DATA\DINO\AREA.IPF
WINDOW=130000.0,450000.0,141000.0,461000.0

This example imports all CSVFILES (*.csv) into the IPF file AREA.IPF for a particular window.

Example 2:

FUNCTION=DINO2IPF
CSVFILE=D:\iMOD-DATA\DINO\BOX12.CSV
GENFILE=D:\IMOD-DATA\AREA.GEN
WINDOW=130000.0,450000.0,141000.0,461000.0

The example above imports the boreholes from the BOX12.CSV for the area within the specified polygon(s) in AREA.GEN.


8.7.2GEOTOP-Function

This function will replace the top of a groundwatermodel with a GEOTOP schematization.

FUNCTION=

GEOTOP

RESULTFOLDER=

Enter the name of the folder that will store the merged model

NLAYG=

Enter the number of GEOTOP model layers.

KVG_L{i}=

Enter the IDF for the i\({}^{th}\) modellayer of GEOTOP that represents the KV i, e.g. TOP_L1=D:\GEOTOP\KVG_L1.IDF.

KHG_L{i}=

Enter the IDF for the i\({}^{th}\) modellayer of GEOTOP that represents the KH i, e.g. TOP_L1=D:\GEOTOP\KHG_L1.IDF.

NLAYM=

Enter the number of model layers in the groundwatermodel
Enter the IDF for the j\({}^{th}\) model layer of the model that represents: ...

IBM_L{j}=

... the iBOUND for layer j, e.g. IBM_L1=D:\UTRECHT\iBOUND_L1.IDF

SHM_L{j}=

... the starting head for layer j, e.g. SHM_L1=D:\UTRECHT\SHEAD_L1.ID

TPM_L{j}=

... the top of layer j, e.g. TPM_L1=D:\UTRECHT\TOP_L1.IDF

BTM_L{j}=

... the bot of layer j, e.g. BTM_L1=D:\UTRECHT\BOT_L1.IDF

KHM_L{j}=

... the KH for layer j, e.g. KHM_L1=D:\UTRECHT\KH_L1.IDF

KAM_L{j}=

... the KA for layer j, e.g. KAM_L1=D:\UTRECHT\KA_L1.IDF

KVM_L{j}=

... the KV for layer j, e.g. KVM_L1=D:\UTRECHT\KV_L1.IDF

WINDOW=

Specify a window (X1,Y1,X2,Y2) for which the entered RUNFILE will be clipped, e.g. WINDOW=125100,345000,135000,355000

CELLSIZE=

Enter the cell size (meter) for the IDF-files that will be created, e.g. CELLSIZE=25

Example 1:

FUNCTION=GEOTOP
RESULTFOLDER=D:\MODEL
NLAYG=95
WINDOW=130000,450000,141000,461000
kvg_l1 =D:\MODEL\kv1.idf
khg_l1 =D:\MODEL\kh1.idf
..........
kvg_l95=D:\MODEL\kv95.idf
khg_l95=D:\MODEL\kh95.idf
NLAYM=9
tpm_l1=d:\modelutrecht\GEOHYDROLOGY\VERSION1\TOP_l1.IDF
btm_l1=d:\modelutrecht\GEOHYDROLOGY\VERSION1\BOTTOM_l1.IDF
kdm_l1=d:\modelutrecht\TRANSMISSIVITY\VERSION1\TRANSMISSIVITY_l1.IDF
vcm_l1=d:\modelutrecht\VERTICALRESISTANCE\VERSION1\VERTICALRESISTANCE_l1.IDF
ibm_l1=d:\modelutrecht\IBOUND\VERSION1\IBOUND_l1.IDF
shm_l1=d:\modelutrecht\STARTINGHEADS\VERSION1\SHEAD_l1.IDF
....
tpm_l9=d:\modelutrecht\GEOHYDROLOGY\VERSION1\TOP_l9.IDF
btm_l9=d:\modelutrecht\GEOHYDROLOGY\VERSION1\BOTTOM_l9.IDF
kdm_l9=d:\modelutrecht\TRANSMISSIVITY\VERSION1\TRANSMISSIVITY_l9.IDF
vcm_l9=d:\modelutrecht\VERTICALRESISTANCE\VERSION1\VERTICALRESISTANCE_l9.IDF
ibm_l9=d:\modelutrecht\IBOUND\VERSION1\IBOUND_l9.IDF
shm_l9=d:\modelutrecht\STARTINGHEADS\VERSION1\SHEAD_l9.IDF


8.7.3GEF2IPF-Function

This function will combine information from a series of GEF files to generate an IPF file with borehole or cone penetration test (CPT) information attached to it. The content of the GEF file is prescribed in Section 9.24.1.

FUNCTION=

GEF2IPF

GEFDIR=

Enter a DIR name that contains the GEF files, e.g. GEFDIR=D:\DATA\.

IPFFILE=

Enter the name of the IPF file to be created, e.g. IPFFILE=D:\DINO.IPF.

GEFTYPE=

Enter a number for the type of file you prefer to read in. 1=CPT, 2=Borehole.

WINDOW=

Specify a window (X1,Y1,X2,Y2) for which the GEF files will be selected, WINDOW=125100.0,345000.0,135000.0,355000.0.

GENFILE=

Enter a name for a GEN-file that contains a polygon that determines the area for which GEF files will be selected for conversion.

Example 1:

FUNCTION=GEF2IPF
GEFDIR=D:\iMOD-DATA\DINO\
IPFFILE=D:\iMOD-DATA\DINO\AREA.IPF
GEFTYPE=1

This example imports all GEFFILES (*.GEF) into the IPF file AREA.IPF.

Example 2:

FUNCTION=GEF2IPF
GEFDIR=D:\iMOD-DATA\DINO\*DELFT*.GEF
GENFILE=D:\IMOD-DATA\AREA.GEN
GEFTYPE=1

The example above imports the GEFFILES for the area within the specified polygon(s) in AREA.GEN.

Note: Before the calculation is started, iMOD asks you what type of GEF-file you want to convert to IPF; one that contains CPT (Cone Penetration Test) information or one that contain Borehole information. Below you can find a brief description of both GEF-file types.


8.7.4CUS-Function

Use this function to determine a minimal number of model layers (aquifers) based on layers that describe geological formations (mainly aquitards). The concept idea is that the vertical distribution of aquitards (or other distinguishing layers), determine the minimal number of model layers to represent the aquifers.

The methodology computes the interrelation ship of all the individual parts within a geological formation and those from other geological formations. This interrelation ship is used to compute the minimal number of model layers to capture all of them without any loss of information. In fact is depends mainly on the lateral distribution of the geological formation whether the number of model layers becomes less than the number of geological formations. The more spread in the distribution, probably less model layers are necessary. In figure an overview is given of the methodology as it shows all the interrelation ships that are computed, e.g. \(\Delta F_1-F_5\) represents the distance between a individual element on the the first and fifth geological formation, and moreover, the vertical position of both elements as the first element should always be in a higher model layer than the fifth element. All these interrelation ships are fed into a linear-programming algorithm than find a model layer for each of them such that the total number of model layers is minimized.

The CUS function is fully automatic, which means that a) the vertical order of given geological formations is irrelevant, b) overlapping geological formations may be clipped and c) interrelation ships in vertical and horizontal direction are computed automatically that yield a minimal set of model layers. With some variables it is possible to set the threshold at what distance of interrelation ship geological formation may join together, in this way aquitards can be combined whenever they are separated less than a given distance.

Schematic overview working of CUS

pictures/h72-end/cus_explanation3.png

pictures/h72-end/cus_explanation4.png

FUNCTION=

CUS

ICPOINTERS=

Enter a value whether or not use a given CUS pointer IDF (0 or 1).

ICPOINTERS=1

FDISTANCES=

file containing predefined vertical distances that is created in a former CUS-action.

or

CRIT_THICKNESS=

maximum vertical step size (e.g. CRIT_THICKNESS=25.0 m) to combine elements laterally.

MIN_THICKNESS=

minimal thickness of the element to be included in the final model, e.g. MIN_THICKNESS=0.5 will include elements thicker than 0.5 meter only.

ZCRIT=

Critical vertical distance. Layers will be connected vertically whenever a percentage of their vertical distance is less than ZCRIT, e.g. ZCRIT=0.5 m.

PERCENTAGE=

Give the percentile for which ZCRIT needs to be taken into account, e.g. if a percentage of 90 is given, layers will be connected if 90% of the distance is less than ZCRIT.

ICLIP=

Enter the name of an IDF file (at least at the dimension of the given IDF files are FORMTOP_Li or FORMTOP keyword) that denote the zone for which an entry is not need to be blanked out. E.g. ICLIP=D:\CUS\ZONES.IDF. Whenever a value of 1 is found in the IDF file, all geological formations that refer to this zone 1, will be blanked out for areas not equal to 1.

IEXPZONE=

Enter a number of additional cells around each individual element in each formation (IEXPZONE>1 to have any effect), to be used to include any element laterally to determine the most optimal model layer, e.g. IEXPZONE=2. Adding a value of IEXPZONE will have the effect that elements that are horizontally nearby (less than 2 cells in this case), will be tried to vertically positioned in the same model layer.

ICPOINTERS=0

NLAY=

Enter the fixed number of model layers to be constructed based on the IDf file with pointer values given at PNT_L{i}.

NFORM=

Enter number of geological formations, e.g. NFORM=19.

FORMTOP_L{i}=

Enter an optional zone number after the file name whenever the keyword ICLIP is used, e.g. FORMTOP_L1=D:\INPUT\BEK1_CK.IDF,1.

FORMBOT_L{i}=

Enter an optional zone number after the file name whenever the keyword ICLIP is used, e.g. FORMBOT_L1=D:\INPUT\BEK1_CK.IDF,1.

ICPOINTERS=0

PNT_L{i}=

Enter an IDF for the i\({}^{th}\) formation that gives a pointer values that refers to a modellayer i, e.g. PNT_L1=D:\INPUT\PNT_L1.IDF. This file contains for example the values 1-5. These values serve as a label of a specific aquitard layer (1-5).

or

FORMTOP=

Enter a path and wildcard to specify for a collection of IDF files containing information about the TOP of the geological formations to be included, e.g. FORMTOP_L1=D:\FORMATIONS\*_TOP.IDF.

FORMBOT=

Enter a path and wildcard to specify for a collection of IDF files containing information about the BOT of the geological formations to be included, e.g. FORMBOT_L1=D:\FORMATIONS\*_BOT.IDF.

OUTPUTFOLDER=

Enter the foldername in which the results will be saved, e.g. OUTPUTFOLDER=D:\RESULT.

WINDOW=

Enter the coordinates of the window that needs to be computed. Enter coordinates of the lower-left corner first and then the coordinates of the upper-right corner, e.g. WINDOW=100000, 400000, 200000, 425000. When WINDOW= is absent iMOD will take the WINDOW-extent of the input IDF’s.

CELLSIZE=

Enter the cell size (meter) for the IDF-files that will be created, e.g. CELL_SIZE=25.0.

TOPSYSTEM=

Enter the name of the IDF-file containing the AHN for the specific model area.

BOTSYSTEM=

Enter the name of the IDF-file containing the bottom-depth of the lowest bottom layer of the model.

Example

FUNCTION= CUS
NLAY=2
WINDOW=120000.0,298000.0,240000.0,430000.0
CELLSIZE=100.0
FORMTOP_L1=D:\MODEL\BEK1_T.IDF
FORMTOP_L2=D:\MODEL\BEK2_T.IDF
FORMBOT_L1=D:\MODEL\BEK1_B.IDF
FORMBOT_L2=D:\MODEL\BEK2_B.IDF
OUTPUTFOLDER=D:\OUTPUT
TOPSYSTEM=D:\MODEL\AHN250.IDF
BOTSYSTEM=D:\MODEL\BEDROCK_TOP.IDF

This example corrects the top and bottom IDF-files specified by the FORMTOP_L{i} and FORMBOT_L{i} keywords in a top-bottom consistent manner and scales the IDF-files to the specified WINDOW and CELLSIZE.


8.7.5SOLID-Function

Use this function to generate hypothetical interfaces (i.e. line in between model layers that represent an artificial interface since any resistance layer (e.g. clay). This function computes the transmissivities and vertical resistance between model layers as well. It uses the PCG solver algorithm for the interpolation of the hypothetical interfaces and uses the existence of permeability field and the top- and bottom elevation to compute the nett transmissivity for each a modellayer and vertical resistance between those model layers. By means of masks it is possible to define those areas for which hypothethical interfaces need to be computed.

FUNCTION=

SOLID

NLAY=

Enter the number of modellayers, e.g. NLAY=6.

OUTPUT-
FOLDER=

Enter the foldername in which the results IDF-files will be saved, e.g. OUTPUTFOLDER=D:\RESULT. The following results will be saved:

\MASK

If IMASK=1, for each interface a mask IDF will be created and saved in this folder. Those may be adjusted afterwards, but set IMASK=0 to avoid that those modified mask files will be overwritten

if ICKDC=1

\FFRAC

For each model layer, the fraction of each geological formation will be saved. It represents the fraction (\(0.0-1.0\)) of the geological that is present in the model layer.

\CFRAC

For each in between model layer (aquitard), the fraction of each geological formation will be saved. It represents the fraction (\(0.0-1.0\)) of the geological that is present in the aquitard.

MDL_TOP_{i}

The TOP elevation for each model layer;

MDL_BOT_{i}

The BOT elevation for each model layer;

MDL_KD_{i}

The total transmissivity for each model layer, it becomes zero when the thickness of the aquifer (model layer) is zero;

MDL_VC_{i}

The vertical resistance over aquitards in between each model layer, excluding the resistance due to the vertical resistance in the above- and beneath lying aquifers. Its value becomes zero if the aquitard is absent;

MDL_KHV_{i}

The horizontal permeability for each model layer, it can be zero for layer thicknesses of zero;

MDL_KVA_{i}

The vertical anisotropy for each model layer, it will always have a value larger than 0 and smaller equal to 1. For non existing model layers (aquifers) this parameter will be one;

MDL_KVV_{i}

The vertical permeability for each aquitard in between each model layer, it becomes zero when the aquitard does not exists;

MDL_KDFRAC_{i}

The total fraction of the model layers that has been parameterised by the permeabilities found by the REGISKHV and/or REGISKVV files, whenever the fraction is 1.0 is means that the entire model layer has been filled in correctly, lower values indicate that areas in the model layers have not been filled in properly.

MDL_CFRAC_{i}

The total fraction of the aquitards in between the model layers that has been parameterised by the permeabilities found by the REGISKHV and/or REGISKVV files. So more comment above;

TOP_L{i}=

Enter the IDF for the i\({}^{\rm th}\) modellayer that represents the top of modellayer \(i\), e.g. TOP_L1=D:\INPUT\TOP_L1.IDF.

ICLC_TL{i}=

Enter the option 0 or 1 to define whether this TOP modellayer needs to be interpolated, e.g. ICLC_TL1=1. This is optional, the default is 1.

BOT_L{i}=

Enter the IDF for the i\({}^{\rm th}\) modellayer that represents the bottom of modellayer \(i\), e.g. BOT_L1=D:\INPUT\BOT_L1.IDF.

ICLC_BL{i}=

Enter the option 0 or 1 to define whether this BOTTOM modellayer needs to be interpolated, e.g. ICLC_BL1=1. This is optional, the default is 1.

IMASK=

Specify IMASK=1 to (re)compute masks. Those are IDF files that contain a pointer value that indicates how the interfaces need to be computed. A mask value can have the following values:

0

means that this particular location will be excluded, those locations are initially formed by non-existence of the upper- and lowermost interface;

-1

means that for that particular area no interface will be computed, the original value will be used;

+1

means that the interface will be computed for this locations.

Each mask IDF file will be saved in the MASK folder under the given OUTPUTFOLDER. Whenever IMASK=0, iMOD will look in this particular folder to read the mask IDF files, make sure that those files are in that folder.

if IMASK=1

ZOFFSET=

Specify a vertical offset (meters) for which mask values need to be set af -1. In other words, whenever the difference between the TOP_L{i} and BOT_L{i+1} is larger than ZOFFSET the mask will be put on -1. Small aquitards can be removed in this manner. By default ZOFFSET=0.0 meter.

IHYPO=

Specify IHYPO=1 to compute the hypothetical interfaces, for mask values of 1. As the values for TOP_L|[i} and BOT_L{i+1} need to be identical for mask values of +1, only the interface for TOP_L|[i} will be computed and BOT_L{i+1} will be set equal to that value.

if IHYPO=1

DZ(.)=
(optional)

Specify for each model layer the minimal thickness (meter), e.g. DZ(1)=1.0 means that the minimal thickness will be 1.0 meter. In this way it is possible to have continuous thicknesses for model layers. This minimal thickness requirement can not be met whenever the distance between two aquitards is less than this DZ. In that case a smaller thickness is forced. By default DZ=0.0 for each model layer.

IMIDELEV=
(optional)

Specify IMIDELEV > 0.0 to force te PCG solver to position the hypothetical interface more-or-less such that model layers have uniform thicknesses. By default IMIDELEV=1.0, however, IMIDELEV=0.0 will deactivate this feature. The higher the value for IMIDELEV the more the hypothetical interface is a true average for all appropriate interfaces, the lower the value, the more smooth is the hypothetical interface will be probably, but the constraint of even distributed interfaces is more violated.

IINT_IDF=
(optional)

This keyword can be defined as IINT_IDF=1 (Default) to use upper and lower situated clay layers by the investigation of the hypothetical interfaces.

IBNDCHK=
(optional)

Specify IBNDCHK=1 to check internally for isolated cells that are NOT connected to constant value cells. By default IBNDCHK=0.

HCLOSE=
(optional)

Specify the closure criterion of the PCG solver, e.g. HCLOSE=0.1 m. By default HCLOSE=0.001 meter.

MICNVG=
(optional)

Specify the number of subsequent inner convergences of the PCG solver, e.g. MICNVG=25. Use this whenever the PCG solver does not find a unique solution. By default MICVNG=5.

ICKDC=

Specify ICKDC=1 to compute transmissivities for model layers and vertical resistances for in between model layers.

if ICKDC=1

FNLAY=
(optional)

Number of formation to be specified separately, e.g. FNLAY=10.

if FNLAY specified

FTOP_L{i}=

Specify the IDF file for the i\(^{\rm th}\) TOP elevation of a geological formations, e.g. FTOP_L1=D:\FORMATION\BEK1_TOP.IDF.

FBOT_L{i}=

Specify the IDF file for the i\(^{\rm th}\) BOT elevation of a geological formations, e.g. FBOT_L1=D:\FORMATION\BEK1_BOT.IDF.

FKHV_L{i}=

Specify the IDF file for the i\(^{\rm th}\) horizontal permeability of a geological formations, e.g. FKHV_L1=D:\FORMATION\BEK1_KHV.IDF. If no file is defined, this permeability will be assumed to be 3.0 the given vertical permeability at FKVV_L{i}.

FKVV_L{i}=

Specify the IDF file for the i\(^{\rm th}\) vertical permeability of a geological formations, e.g. FKVV_L1=D:\FORMATION\BEK1_KVV.IDF. If no file is defined, this permeability will be assumed to be 0.3 the horizontal permeability at FKHV_L{i}.

if FNLAY is not specified

FOLDERTOP=

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

FOLDERBOT=

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

FOLDERKHV=

Specify the folder that stores the horizontal permeability of geological formations, e.g. FOLDERKHV=D:\REGIS\*-KH-SK.IDF. All files will be used that fit this wildcard definition. If no file is found for the horizontal permeability for a particular geological formation, this permeability will be assumed to be 3.0 the vertical permeability.

FOLDERKVV=

Specify the folder that stores the vertical permeability of geological formations, e.g. FOLDERKVV=D:\REGIS\*-KV-SK.IDF. All files will be used that fit this wildcard definition. If no file is found for the vertical permeability for a particular geological formation, this permeability will be assumed to be 0.3 the horizontal permeability.

WINDOW=
(optional)

Enter the coordinates of the window that needs to be computed. Enter coordinates of the lower-left corner first and then the coordinates of the upper-right corner, e.g. WINDOW=100000.0, 400000.0, 200000.0, 425000.0. When WINDOW= is absent, the entered IDF-files by TOP_L{i} and BOT_L{i} need to be equally in their dimensions. Otherwise they will be upscaled (mean) or downscaled (interpolation) to the entered CELLSIZE.

CELLSIZE=

Enter the cell size (meter) for the IDF-files that will be created, e.g. CELLSIZE=25.0.

NGEN=
(optional)

Enter a value of the amount of GEN-files containing fault information to be taken into account in the interpolation.

if NGEN \(\ge \) 1

GEN_i=

Related to NGEN with this keyword all the GEN-files are included in the interpolation. Enter first a value for the interface number, followed by the name of the GEN-file. Example: GEN_1=3,faults_l3.GEN.

Example

FUNCTION=SOLID
NLAY=4
TOP_L1=D:\MODEL\TOP_L1.IDF
TOP_L2=D:\MODEL\TOP_L2.IDF
TOP_L3=D:\MODEL\TOP_L3.IDF
TOP_L4=D:\MODEL\TOP_L4.IDF
BOT_L1=D:\MODEL\BOT_L1.IDF
BOT_L2=D:\MODEL\BOT_L2.IDF
BOT_L3=D:\MODEL\BOT_L3.IDF
BOT_L4=D:\MODEL\BOT_L4.IDF
IMASK=1
IHYPO=1
ICKDC=1
FOLDERTOP=D:\REGIS\*-T-CK.IDF.
FOLDERBOT=D:\REGIS\*-B-CK.IDF.
FOLDERKHV=D:\REGIS\*-KH-SK.IDF.
FOLDERKVV=D:\REGIS\*-KV-SK.IDF.
OUTPUTFOLDER=D:\OUTPUT

This example creates a mask files based on the specified top and bottom IDF-files specified by the TOP_L{i} and BOT_L{i} keywords. For those areas that does not contain an aquitard, the hypothetical interfaces will be computed. After that the transmissivities and vertical resistances will be computed and the other output as specified at the keyword OUTPUTFOLDER.

Example of (left) the computed thickness of an aquitard; (right) the corresponding values for the mask IDF (green is +1 and red = -1). The green area will be filled in by hypothetical interfaces.

pictures/h5/example_solidtool_mask.png

Example of computed hypothetical interfaces as TOP and BOTTOM elevation for the model layers.

pictures/h5/cross-section-solidtool.png

Example of computed fractions of a geological formation in three different model layers. The formation has been part of three model layers and the total transmissivity of each model layers depends on the given fraction of the geological formation. Red represents a higher fraction than yellow.

pictures/h5/solidtool_fractions.png


8.7.6FLUMY-Function

Use this function to create Flumy textfiles with borehole information out of IPF-file related textfiles. Flumy textfiles contain information about the fluvial deposits in (former) riverbeds (Flumy is a Geovariances software program).

FUNCTION=

FLUMY

IPFFILE=

Enter the name of the IPF-file that contains the borehole information to be converted to a Flumy readable format.

OFFSET=

Optional variable. Enter a value to elevate the depth of the borehole, e.g. current depth of borehole is -50 m, on OFFSET=50 results in a borehole reference depth of 0 m.

NPARAM=

Enter the number of parameters {i} that needs to be distinguished in the Flumy-textfile(s).

GRAIN{i}=

Enter the name of each grain type as given in IPF related textfile, e.g. GRAIN1=SILT or GRAIN2=’Sandy Clay’

FACIESL{i}=

Enter the verb for the location in the fluvial area where the i\(^{th}\) Grain type will be located, e.g. FACIESL1=OB, in case "SILT" needs to be located in the Overbanks.

FACIESN{i}=

Enter the number related to the facies layer as defined with FACIESL, e.g. FACIESL1=OB, FACIESN1=8.

The created Flumy textfile(s) are(is) stored in the working directory of the FLUMY-batchfile within a new generated folder “FLUMY”.

Example 1

FUNCTION=FLUMY
IPFFILE=’D:\flumy\Boreholes\Borelogs.ipf’
NPARAM=2
GRAIN1=SILT
FACIESL1=OB
FACIESN1=8
GRAIN2=’SANDY Clay’
FACIESL2=PB
FACIESN2=2

Note: Be aware that quotes are obligatory around a “GRAIN”-description containing more than one word, e.g. ’SANDY CLAY’. Otherwise iMOD is not able to read this parameter properly.


8.7.7GEOCONNECT-function

Use this function to:

General Settings
These are the general settings needed in each *.ini file using the GEOCONNECT-Batch function:

FUNCTION=

GEOCONNECT

NLAY=

Enter the amount of model layers, e.g. NLAY=10.

ACTLAYERS=
(optional)

Enter a string of values to include or exclude a specific model layer from the computation; 0=inactive, 1=active, on default all layers are used in de computation (similar to e.g.: ACTLAYERS=1111111111). E.g. in case of the amount of model layers is 10 and it is preferred to only take the first 6 layers into account: ACTLAYERS=1111110000.

REGISFOLDER=

Give the directory and name of the folder where all REGIS-files are stored. Note: subdirectories are not allowed and the filenames need to be of the following format: abbreviation formation name-t/b/ks/kv-ck/sk.idf (’t’ and ’b’ need to be combined with ’ck’, and ’ks’ and ’kv’ with ’sk’), e.g. d:\Model\REGIS\bez1-b-ck.idf.

TOPFOLDER=

Give the directory and name of the folder of the model TOP-files, e.g. d:\Model\TOP\TOP.

BOTFOLDER=

Give the directory and name of the folder of the model BOT-files, e.g. d:\Model\BOT\BOT.

WINDOW=

Enter the coordinates of the window that needs to be computed. Enter coordinates of the lower-left corner first and then the coordinates of the upper-right corner, e.g. WINDOW=100000, 400000, 200000, 425000. When WINDOW= is absent iMOD will take the WINDOW-extent of the input IDF’s.

CELLSIZE=

Enter the cell size (meter) for the IDF-files that will be created, e.g. CELL_SIZE=25.0.

NFORM=

Enter a value for the total amount of formation factors to be read from this *.ini file, e.g. NFORM=127.

FORM{i}=

Give the name of the i’th formation and the corresponding factor, e.g. FORM1=HLC,1.000. The largest i-number needs to correspond with the total amount of factors as defined with NFORM.

OUTPUTFOLDER=

Give the directory and name of the folder to store the results of the preprocessing computation.

IFLAG=

Enter a value for the specific GeoConnect function to be used: 1=Preprocessing, 2=Postprocessing. E.g. IFLAG=1.

Options
In case of IFLAG=1 (Preprocessing option), the *.ini file with specific options needs to contain:

ISAVEK=
(optional)

Use this keyword to save KHV, KVV and KVA values in IDF-format; 1=saved and 0=not saved, on default ISAVEK=1.

ISAVEC=
(optional)

Use this keyword to save KDW and VCW values in IDF-format; 1=saved and 0=not saved, on default ISAVEC=0.

In case of IFLAG=2 (Postprocessing option), the *.ini file with specific options needs to contain:

DBASEFOLDER=

Give the directory and name of the folder containing DBASE model information, e.g. p:\1221301-ibrahym-dld\DBASE_V1

IAGGR=

Give the value related to the aggregation option to be used, IAGGR=1: apply to model output, IAGGR=2: apply to model input, IAGGR=3: apply to ipf-file

if IAGGR=1…

MODELFOLDER=

Give the directory and name of the folder containing the preferred model output information.

MODELTYPE=

Choose the favored variable to apply the aggregation to, e.g. MODELTYPE=2. Options are: 1="HEAD", 2="BDGWEL", 3="BDGRIV", 4="BDGDRN".

if IAGGR=2…

INPUTTYPE=

Choose the preferred variable to apply the aggregation to, e.g. INPUTTYPE=3. Options are: 1="KDW", 2="VCW", 3="KHV", 4="KVV".

if IAGGR=3…

IPFFILE=

Give the directory and name of the IPF-file to apply the aggregation to.

IDUPLICATES=

Give the type of aggregation. When IDUPLICATES=1 the maximum value per grid cell is taken from the files to be aggregated, IDUPLICATES=2: maximum value, IDUPLICATES=3: average value, and IDUPLICATES=4 the sum of the value per grid cell for all files is taken as new value in the aggregated grid cell in the file.

ISAVETB=
(Optional)

Use this option to include TOP- and BOT elevation IDF-files in the aggregation and save these files in the given Output folder. In case ISAVETB=1, TOP and BOT elevation are saved, otherwhise when ISAVETB=0 these files are not saved in the Output folder.

OUTPUTFOLDER=

Give the directory and name of the folder to store the results of the preprocessing computation.

Be sure that the geostratigraphy.txt file is placed in the defined DBASEFOLDER!

Example 1
FUNCTION=GEOCONNECT
NLAY=19
ACTLAYERS=1111111100011100000
REGISFOLDER=d:\Model_Ibrahym2\REGIS21
TOPFOLDER=d:\Model_Ibrahym2\DBASE_V2\TOP\VERSION_1
BOTFOLDER=d:\Model_Ibrahym2\DBASE_V2\BOT\VERSION_1
IFLAG=1
NFORM=8
FORM1=HLC,1.000
FORM2=BXSCK1,1.300
FORM3=BXZ1,1.000
FORM4=BXK1,1.000
FORM5=BXLMK1,2.500
FORM6=BXZ2,1.000
FORM7=BXK2,1.000
FORM8=BXZ3,0.910
OUTPUTFOLDER=d:\Model_Ibrahym2\Results\
ISAVEK=1
ISAVEC=0

With this example the KHV, KVV and KVA grids are (re)calculated (preprocessing) based on a given factor per formation (FORM{i}) for the first 8 formation seen from the top of the Ibrahym model.

Example 2
FUNCTION=GEOCONNECT
NLAY=19
ACTLAYERS=1111111100011100000
REGISFOLDER=d:\Model_Ibrahym2\REGIS21
TOPFOLDER=d:\Model_Ibrahym2\DBASE_V2\TOP\VERSION_1
BOTFOLDER=d:\Model_Ibrahym2\DBASE_V2\BOT\VERSION_1
IFLAG=2
DBASEFOLDER=d:\Model_Ibrahym2\DBASE_V2
NFORM=8
FORM1=HLC,1
FORM2=BXSCK1,4
FORM3=BXZ1,4
FORM4=BXK1,4
FORM5=BXLMK1,4
FORM6=BXZ2,4
FORM7=BXK2,4
FORM8=BXZ3,4
IAGGR=2
INPUTTYPE=3
IDUPLICATES=3
ISAVETB=0
OUTPUTFOLDER=d:\Model_Ibrahym2\DBASE_V4\

In this example the (postprocessing) aggregation is applied to the input model variable "KHV" for which all BX-formations are taken together as one formation. An average value is calculated per grid cell for all BX-formations. Top and Bot elevations are not saved.


8.7.8CREATEIZONE-Function

Use this function to calculate zones and corresponding fractions per model layer based on geologic formations, these IDF files can be used during an iPEST optimization.

FUNCTION=

CREATEIZONE

OFOLDER=

Give an output folder name to store all the fractions per model layer.

PFOLDER=

Give a folder name that contains the fraction per model layers per geological formations. This is the results

NLAY=

Enter the number of modellayers, e.g. NLAY=6

MINF=

Minimum fraction to assign .... to a zone....

IZONEOFFSET=

Enter the number of zones .....

IGROUPOFFSET=

Enter the number of zones .....

NFORMATIONS=

Enter the number of formations.

FORMATION{i}=

Enter the ... for the i\({}^{th}\) formation.

TPARAMETER=

Enter the name of the ...

Example

FUNCTION=CREATEIZONE
OFOLDER=D:\MODEL
PFOLDER=D:\RESULTS
NLAY=2
MINF=...
IZONEOFFSET=...
IGROUPOFFSET=...
NFORMATIONS=2
FORMATION1=...
FORMATION2=...
TPARAMETER=5

The example above will ......