pictures/imod_logo.png

iMOD User Manual version 5.2 (html)


8.3ISG-FUNCTIONS


8.3.1ISGEXPORT-Function

The ISGEXPORT function reads an ISG file and export the entire content into readable files.

FUNCTION=

ISGEXPORT

ISGFILE=

Give an ISG file to be exported.

IEXPORT=

Enter the type of export.

IEXPORT=0
This type exports a complete dump of all individual files of an ISG file structure, see section 9.9. These files can be modified and imported again via the iMOD Batch function ISGIMPORT, see section 8.3.1.
*_ISG.TXT
1,0,“Date”,“Water level”,“Bottom level”,“Resistance”,“Inf.factor”
“SEGMENT 1”,1,4,1,2,1,1,1,1,1,1
*_ISP.TXT
RECORD X Y
1 200250.000 402250.000
2 201250.000 401875.000
3 201250.000 400625.000
4 202250.000 400250.000
*_ISD1.TXT
RECORD N IREF DIST CNAME
1 2 1 0.000 CALCPOINTID2
2 2 3 3386.001 CALCPOINTID2
*_ISD2.TXT
RECORD IDATE WLVL BTML RESIS INFF
1 20040101 7.000 6.000 1.000 0.330
2 20040701 7.000 6.000 1.000 0.330
3 20040101 9.000 8.000 1.000 0.330
4 20040701 8.500 8.000 1.000 0.330
*_IST1.TXT
RECORD N IREF DIST CNAME
1 2 1 500.000 WEIR ON SEGM
*_IST2.TXT
RECORD IDATE WLVL_UP WLVL_DOWN
1 0 1.200 1.000
2 0 1.000 0.800
*_ISC1.TXT
RECORD N IREF DIST CNAME
1 3 1 150.000 CROSSSECTION
*_ISC2.TXT
RECORD DISTANCE BOTTOM MRC ZP
1 -5.000 5.000 25.000 0.002
2 0.000 0.000 25.000 0.002
3 5.000 5.000 25.000 0.002
*_ISQ1.TXT
RECORD N IREF DIST CNAME
1 2 1 750.000 QH RELATION
*_ISC2.TXT
RECORD Q W D F
1 10.000 1.000 7.500 1.000
2 20.000 1.500 15.250 1.200

IEXPORT=1
This type exports all cross-sections; an example is given below:
ISEGMENT,SEGNAME ,ICROSS,DISTANCE, X, Y, CROSS.NAME
1,SEGMENT 1, 1, 150.00,200390,402197,CROSSSECTION
IEXPORT=2
This type exports all calculation points; an example is given below:
ISEGMENT,SEGNAME ,ICALC,DISTANCE, X, Y, CALC.NAME
1,SEGMENT 1, 1, 0.00, 200250,402250,CALCPOINTID1
1,SEGMENT 1, 2, 3386.00, 202250,400250,CALCPOINTID2

EXPORT-
FNAME=

Specify the output name for the file created whenever IEXPORT=1 or IEXPORT=2

EXPORT-
MAP=

Specify the output folder for all the files whenever IEXPORT=0

Example 1

FUNCTION=ISGEXPORT
ISGFILE=D:\MODEL\RIVER.ISG
IEXPORT=0


8.3.2ISGIMPORT-Function

The ISGIMPORT function generates an ISG file from files created via the IEXPORT=0 option on the iMODBatch Function ISGEXPORT, see section 8.3.1.

FUNCTION=

ISGIMPORT

ISGEXPORT=

Give one of the exported ISG equivalent files, e.g. ISGEXPORT=D:\RIVER_ISG.TXT. All the other corresponding files will be assigned automatically, i.e. RIVER_ISP.TXT, RIVER_ISD1.TXT, RIVER_ISD2.TXT, RIVER_ISC1.TXT, RIVER_ISC2.TXT, RIVER_IST1.TXT, RIVER_IST2.TXT, RIVER_ISQ1.TXT and RIVER_ISQ2.TXT.

ISGFILE=

Specify the ISG file that need to be created, e.g. ISGFILE=D:\RIVER.ISG.

ISGDOUBLE=
(optional)

Specify whether the ISG needs to be a single precision (ISGDOUBLE=0) or double precision (ISGDOUBLE=1) files. BY default ISGDOUBLE=0.

Example 1

FUNCTION=ISGIMPORT
ISGEXPORT=D:\MODEL\RIVER_ISG.TXT
ISGFILE=D:\RIVER.ISG
ISGDOUBLE=1


8.3.3GEN2ISG-Function

The GEN2ISG function reads a GEN file and creates a ISG file. There are two ways to use this function; IUSEDAT=1 and IUSDAT=0. These will be explained below.

Use IUSEDAT=1 to force using the *.DAT file (with GEN file information).

FUNCTION=

GEN2ISG

GENFNAME=

Give a GEN file containing x and y coordinates of GEN segments.

IUSEDAT=1

Enter the column number in the DAT file that represents:

STAGE1_COLUMN=

Stage on Startpoint of the Segment, e.g. STAGE1_COLUMN=1

STAGE2_COLUMN=

Stage on Endpoint of the Segment, e.g. STAGE2_COLUMN=2

BOTL1_COLUMN=

Bottom on Startpoint of the Segment, e.g. BOTL1_COLUMN=3

BOTL2_COLUMN=

Bottom on Endpoint of then Segment, e.g. BOTL2_COLUMN=4

shape of the Cross-section of the Segment:

SLOPE_L_COLUMN=

Left slope, e.g. SLOPE_L_COLUMN=5

SLOPE_R_COLUMN=

Right slope, e.g. SLOPE_R_COLUMN=6

BWIDTH_COLUMN=

Bottom width, e.g. BWIDTH_COLUMN=7

Use IUSEDAT=0 if you do not have a *.DAT file (with GEN file information).

FUNCTION=

GEN2ISG

GENFAME=

Give a GEN file containing x and y coordinates of GEN segments.

IUSEDAT=0

the following keywords must be given:

IDFSTAGE=

Enter the column number in the DAT file that represents the y coordinate, e.g. STAGE1_COLUMN=1.

IDFSUMMER=

Enter the filename of the IDF-file containing summer heads.

IDFSUMMER_BACKUP=

Enter the filename of the IDF-file containing summer heads. This file is created as backup-file in case no values are found in IDFSUMMER or by using the SAMPLE_SEARCH function.

IDFWINTER=

Enter the filename of the IDF-file containing winter heads.

IDFWINTER_BACKUP=

Enter the filename of the IDF-file containing winter heads. This file is created as backup-file in case no values are found in IDFWINTER or by using the SAMPLE_SEARCH function.

SUMMERPERIOD=
(optional)

Enter the day and month at which the summer period starts (e.g. default value is SUMMERPERIOD=0104 in case summer starts at the 1\(^{\rm st}\) of April each year).

WINTERPERIOD=
(optional)

Enter the day and month at which the winter period starts (e.g. default value is WINTERPERIOD=0110 in case summer starts at the 1\(^{\rm st}\) of October each year)

START_YEAR=

Enter the starting year of the calculation period, e.g START_YEAR=2014.

END_YEAR=

Enter the end year of the calculation period, e.g. END_YEAR=2020.

IDFBOTTOM=

Enter the name of the IDF-file containing all bottom values of the area.

BOTTOMVALUE=

If there’s no IDFBOTTOM-file available give a constant value for the bottom height, e.g. BOTTOMVALUE=12.32.

SAMPLE_DISTANCE
(optional)=

Enter a value (in meters) that accounts for the distance iMOD needs to sample the IDF-file on, e.g. 100 m. On default SAMPLE_DISTANCE=250 m.

CCFFNAME=

Enter the name of a CCF-file (See also section 9.10) that describes the cross-section that will be used to insert cross-section for all segments.

SEARCH_DISTANCE=
(optional)

Enter a value (in meters) that accounts for the distance iMOD needs to resample the IDF-file on, e.g. 50 m. On default SEARCH_DISTANCE=250 m. This value will only be used if there were no values found by making use of the SAMPLE_DISTANCE.

IPUZZLE=

Enter a 0.0 or 1.0. If a value of 1.0 is given, iMOD tries to find all possible places where segments can be joined and join them if needed.

IDFRESISTANCE=

Enter the name of IDF-file that contains resistance values of the area.

RESISTANCE=

Enter a resistance value that will be applied to the whole area. Can only be used if IDFRESISTANCE is not available.

IDFINFILTRATIONFACTOR=

Enter the name of IDF-file that contains infiltration factors for the area.

INFILTRATIONFACTOR=

Enter an infiltration factor that will be applied to the whole area. Can only be used if IDFINFILTRATIONFACTOR is not available.

ICORDIR=
(optional)

Enter ICORDIR=1 to correct the direction of flow (order of coordinates) based upon the levels read on the first- and last calculation point for a segment. By default ICORDIR=0 and no correction is applied.

OUTFILE=

Enter the name of the ISG output-file, e.g. NEWFILE.ISG.

Example 1

FUNCTION=GEN2ISG
FUNCTION=GEN2ISG
GENFNAME=d:\Model\Data\Shape_data\river_lines.gen
IUSEDAT=0
IDFSUMMER=d:\Model\Basic_data\SUMMER_LEVEL_RIVER.IDF
IDFWINTER=d:\Model\Data\Basic_data\WINTER_LEVEL_RIVER.IDF
IDFBOTTOM=d:\Model\Data\Basic_data\BODEMHOOGTE_RIVER.IDF
SAMPLE_DISTANCE=25.0
CCFFNAME=d:\Model\Data\Basic_data\CROSS-SECTION.CCF
IDFRESISTANCE=d:\Model\Data\Basic_data\RIV_RESISTANCE.idf
IDFINFILTRATIONFACTOR=d:\Model\Data\Basic_data\INFFACTOR_RIVER.IDF
OUTFILE=d:\Model\Data\ISG_data\River.isg

Above an example is given how to use the GEN2ISG. This example will generate and ISG file based on river levels, bottom elevation, resistance, infiltration factor and the location of the river segments. The ISG will contain the following information: as many as cross-sections per river segment as there are unique ID-numbers in CCFFNAME related file, calculation nodes on each segment intersection and segment nodes as much as there are coordinate points defined in the GEN file.


8.3.4ISGGRID-Function

Use this function to rasterize the selected ISG-files into IDF-files that can be used by iMODFLOW in a runfile.

FUNCTION=

ISGGRID

ISGFILE_IN=

Enter an ISG-file that need to be simplified, e.g. ISGFILE_IN=D:\PO.ISG.

CELL_SIZE=

Enter the cell size (meter) for the IDF-files that will be created from the ISG-file mentioned by ISGFILE_IN, e.g. CELL_SIZE=25.0.

MINDEPTH=
(optional)

Enter the minimum water depth (meter) used for the calculation of the conductance of the stream bed, e.g. MINDEPTH=1.0. The default value is MINDEPTH=0.1.

MAXWIDTH=
(optional)

Enter the maximal width of a stream (meter) used for the calculation of the conductance of the stream bed. Introducing this parameter limits any stream width larger than MAXWIDTH, e.g. MAXWIDTH=1000. The default value is MAXWIDTH=250.

OUTPUT
FOLDER=

Enter a foldername to save all rasters into, e.g.
OUTPUTFOLDER=D:\OUTPUT. The following rasters will be saved:

1

COND

The computed stream bed conductance (m\({}^{2}\)/day)

2

STAGE

The interpolated riverlevel (m+MSL)

3

BOTTOM

The interpolated riverbed height (m+MSL)

4

INFFCT

The interpolated river infiltration factor (-)

5

TOTAL_LENGTH

Total length if existing river segments in a single rastercell (meter)

6

MEAN_
WPERIMETER

Mean wetted perimeter within a river segment

7

MEAN_WIDTH

Mean stream bed width

8

RESISTANCE

Interpolated river resistance (days). IMPORTANT to note is that a minimal resistance is applied of 0.001 days to avoid extraordinary conductance (COND) values.

9

EROSION

Erosion matrix to be used to extent the riverbed existence over more rastercells

The following will be created only whenever ICDIST=1

10

EFFECT

The computed water level that are influenced by the weirs.

11

CUR_ID

Identification of structures for current segment.

12

NEX_ID

Identification of following structure for current segment.

POSTFIX=
(optional)

Enter a postfix to be used to add to the end of the IDF-file names mentioned above, e.g. POSTFIX=_SUMMER yields STAGE_SUMMER.IDF instead of STAGE.IDF.

NODATA=

Enter a NodataValue for which water levels will be skipped in determining the waterlevels along profiles, e.g. NODATA=-999.

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 entire ISG will be gridded for its maximum extent.

ISAVE=

Enter the numbers to be saved solely, e.g. ISAVE=1,1,1,0,0,0,0,0,0,0,0,0 to identify that the IDF-files COND, STAGE and BOTTOM need to be saved.

IPERIOD=
(optional)

Specify whether the water levels need to be calculated for a specific period, by default IPERIOD=1 which means that water levels will be computed as the mean over the entire existing periods within the ISG-file (which can be different among the segments). Specify IPERIOD=2 to enter a date over which the water levels will be averaged.

Enter the following keyword only for IPERIOD=2.

SDATE=

Enter a starting date to compute averaged water levels for, e.g. SDATE=19910101 to represent the 1\({}^{st}\) of January 1991.

EDATE=

Enter a starting date to compute averaged water levels for, e.g. EDATE=19911231 to represent the 31\({}^{st}\) of December 1991.

DDATE=
(optional)

Enter a date-difference to be used to compute more rasters for different periods, e.g. DDATE=14 means that a sequence between SDATE and EDATE will be computed with length of 14 days. By default DDATE=0 which will ignore any time steps in-between the SDATE and EDATE variables. The names of the IDF-file will be extended to include a date notification, e.g. STAGE{POSTFIX}_19910101.IDF

ICDIST=
(optional)

Set this value to \(1\) to compute the effects of weir as stored in the ISG file. By default ICDIST=0. See keyword OUTPUTFOLDER to get the names of the extra IDF files that will be created.

ISIMGRO=
(optional)

Set this value to 1 to export the gridded values for the ISG into a MetaSWAP file svat_swnr_drng.inp.
SVAT2SWNR_DRNG=

Enter the name for the svatswnr_drng.inp file.

SEGMENTCSVFNAME=

Enter the CSV that contains the list of ...

THIESSENFNAME=

Enter an IDF file that represents the SVAT-id for MetaSWAP.

AHNFNAME=

Enter an IDF file with the surface level.

SYSID=

Enter a single value for the system identification.

WDEPTH=

Enter a water depth that will be used to define the appropriate trapezia for MetaSWAP.

IAVERAGE=
(optional)

Enter IAVERAGE=1 to apply an arithmetic mean for stages, bottomlevels, resistances and infiltrationfactor over time. Enter IAVERAGE=2 to apply the median value for those parameters. The default value is IAVERAGE=1.

IEXPORT=
(optional)

Set this value to 1 to export the gridded ISG into a MODFLOW river file, important to notice is that it yield a single value for each gridded cell. The export river file will be called OUTPUTFOLDER \modflow.riv. By default IEXPORT=0 and IDF files will be created.

NLAY=

Enter the number of model layers for which the gridded ISG file need to be assigned vertically, e.g. NLAY=3. This option is only valid whenever IEXPORT=1.

TOP_Li=

Enter an IDF file that represents the TOP elevation of the i\(^th\) layer, e.g. TOP_L1=D:\TOP_L1.IDF.

BOT_Li=

Enter an IDF file that represents the BOT elevation of the i\(^th\) layer, e.g. BOT_L2=D:\BOT_L2.IDF.

KHV_Li=

Enter an IDF file that represents the horizontal permeability of the i\(^th\) layer, e.g. KHV_L2=D:\KHV_L2.IDF.

BND_Li=

Enter an IDF file that represents the boundary condition of the i\(^th\) layer, e.g. BND_L4=D:\BND_L4.IDF.

DMM
FILE=
(optional)

Enter DMMFILE is 1 to export the coupling files DFM1DTOMFRIV_WL.DMM and MFRIVTODFM1D_Q.DMM used by the coupling with DFLOW-FM.

Example:

FUNCTION=ISGGRID
ISGFILE_IN=D:\PO.ISG
CELL_SIZE=100.0
NODATA=-999.99
ISAVE=1,1,1,1,0,0,0,0,0,0,0,0
IPERIOD=2
SDATE=19980101
EDATE=19980131
OUTPUTFOLDER=D:\PO_GRIDS

The example above will rasterize the entire ISG for the period of the 1\(^{\rm th}\) of January up to the 31\(^{\rm th}\) of January 1998 on a 100x100 meter grid.


8.3.5ISGADDCROSSSECTION-Function

Use this function to add cross-sections to an existing ISG-file (see section 9.9.3 for more information about the content of an ISG-file and storage of cross-sections). The methodology is twofold:

FUNCTION=

ISGADDCROSSSECTION

ISGFILE_IN

Enter an ISG-file for which cross-sections need to be added, e.g. ISGFILE_IN=D:\DATA\MAAS.ISG.

ISGFILE_OUT=

Enter an ISG-file to save the renewed ISG for, e.g.
ISGFILE_OUT=D:\DATA\MAAS_NEWCROSSSECTIONS.IDF.


2-D Cross-sections

CROSS_PNTR=

Enter the name of an IDF-file describes the spatial distribution of two-dimensional cross-sections, e.g. CROSS_PNTR=D:\DATA\PNTR.IDF. This pointer file is used to denote areas with equal values as the pointer value at the location of the cross-section on the segmenet. The bathymetry for those areas will be read from CROSS_BATH and applied as a 2-D cross-section on the segment.

CROSS_BATH=

Enter the name of an IDF-file that describes the bathymetry for the riverbed at the locations where the values for CROSS_PNTR are \(\ne \) to the pointer value at the corresponding cross-section at the segment, e.g. CROSS_BATH=D:\DATA\RIVERBED.IDF.

CROSS_ZCHK
(optional)=

Enter the name of an IDF-file that describes the Reference Height to be used to distinguish between areas with positive and negative values for CROSS_PNTR, e.g. CROSS_ZCHK=D:\DATA\REFHEIGHT.IDF.

CROSS_CVAL
(optional)=

Enter the name of an IDF-file that describes the Resistance values to be used to distinguish different resistance in inundation areas, e.g. CROSS_CVAL=D:\DATA\RESISTANCE.IDF. The values in this IDF act as a multiplication factor to the given resistance (attribute RESIS in the ISD2-file, see 9.9.2) at the nearest calculation point. Bear in mind that the multiplication factor will be saved in the ISG as an integer with a maximal value of 256.


1-D Cross-sections

CROSSSEC-
TION_IN=

Enter the filename that stores the renewed cross-sections, e.g. CROSSSECTION_IN=D:\DATA\CROSS.TXT. The syntax of the CROSSSECTION_IN file is a free-formatted, comma-separated-values file with for which each row is defined as follows (be aware that you do not include a header in the text-file):
XC,YC,LABEL,N,X\({}_{1}\),X\({}_{2}\),..,X\({}_{N}\),Z\({}_{1}\),Z\({}_{2}\),..,Z\({}_{N}\)

XC

X-coordinate (meter) for the cross-section;

YC

Y-coordinate (meter) for the cross-section

LABEL

Label for the cross-section, maximum length is 32 characters.

N

Number of cross-sections points.

X\({}_{i}\)

Specify as many distances as needed to define the bathymetry of the riverbed. The amount of definitions N should be \(>\)3.

Z\({}_{i}\)

Specify as many elevations as needed to define the bathymetry of the riverbed. The amount of definitions N would be equal the number of definitions used for X\({}_{i}\).

ICLEAN=
(optional)

Enter ICLEAN=1 to clean ALL cross-sections before adding new ones, apply ICLEAN=2 to remove existing cross-section only for those segments where an update of the cross-section will be applied. By default ICLEAN=1.

WIDTH_IDF=
(optional)

Specify an IDF that represents the width of default cross-sections to be placed on all segments and a default water depth of 5 meter where no cross-section will be placed based on the entered data by the following keywords, e.g. WIDTH_IDF=D:\DIST.IDF. This keyword is necessary only whenever ICLEAN=1.

MAXDIST=
(optional)

Specify a distance (meter) over which the cross-section will be snapped to the segment, e.g. MAXDIST=5.0. By default MAXDIST=0.0 meter.

Example 1 (one-dimensional cross-sections):

FUNCTION=ISGADDCROSSSECTION
ISGFILE_IN=D:\iMOD-DATA\MAAS.ISG
CROSSSECTION_IN=D:\DATA\CROSS.TXT
WIDTH_IDF=D:\DATA\WIDTH.IDF
MAXDIST=2.5
ISGFILE_OUT=D:\iMOD-DATA\MAAS_RENEWEDCROSSSECTIONS.ISG

The example above will add cross-sections based on the entered CROSS.TXT file that specifies a cross-section for “New Cross” as follows:

12000.0,45300.0,”New Cross”,-10.0,-5.0,-2.5,2.5,7.5,12.0,5.0,3.0,2.0,1.0,2.5,5.0

the results will be saved in MAAS_RENEWEDCROSSSECTION.ISG.

Example 2 (two-dimensional cross-sections):

FUNCTION=ISGADDCROSSSECTION
CROSS_PNTR=D:\DATA\PNTR.IDF
CROSS_BATH=D:\DATA\BATHEMETRY.IDF
ISGFILE_OUT=D:\iMOD-DATA\MAAS_RENEWEDCROSSSECTIONS.ISG

The example above will transform the existing cross sections with two-dimensional definitions based on the pointerfile read in CROSS_PNTR and the corresponding bathymetry read in BATHEMETRY.IDF.


8.3.6ISGSIMPLIFY-Function

Use this function to reduce the amount of calculation points in a ISD file (part of the ISG-files, see section 9.9.2). iMOD will eliminate calculation points that do not add significant information to the declination of waterlevels, in other words, whenever the gradient of the waterlevel can be described by less calculation points, iMOD will locate those calculation points that are able to represent the original waterlevel most optimally. iMOD will use the mean waterlevels for all calculation nodes to determine a mean descent of waterlevels along a segment. Simplification will be carried out for segments as a whole. Whenever segments will be very short, this function will have a minor effect.

FUNCTION=

ISGSIMPLIFY

ISGFILE_IN=

Enter an ISG-file that need to be simplified, e.g. ISGFILE_IN=D:\DATA\MAAS.ISG.

ZTOLERANCE=

Specify a distance (meter) for which the simplified waterlevel along a profile may differ from the original one, e.g. ZTOLERANCE=0.10.

NODATA=

Enter a NoDataValue for which waterlevels will be skipped in determining the waterlevels along profiles, e.g. NODATA=-999.

ISGFILE_OUT=

Enter an ISG-file to save the simplified ISG for, e.g. ISGFILE_OUT=D:\DATA\MAAS_SIMPLIFIED.IDF.

Example:

FUNCTION=ISGSIMPLIFY
ISGFILE_IN=D:\iMOD-DATA\MAAS.ISG
ZTOLERANCE=0.10
NODATA=-999.99
ISGFILE_OUT=D:\iMOD-DATA\MAAS_SIMPLIFIED.ISG

The example above will reduce the amount of calculation points such that the simplified waterlevel will be differ more than 0.10 from the original one, the results will be saved in MAAS_SIMPLIFIED.ISG.


8.3.7ISGADJUST-Function

Use this function to perform changes to an existing ISG.

FUNCTION=

ISGADJUST

SESFILE=

Enter the name of the SES file,
e.g. D:\iMOD-DATA\ISGEDIT\ISG-change-stage.SES.

LOGFILE=

Enter a name for the logfile showing all changes by listing both old and new parameter values , e.g. D:\iMOD-DATA\ISGEDIT\ISG-LOG.TXT
default = .\log_ses.txt

OUTNAME=

Foldername for new location to save all ISG related files (*.isg, *.isp, *.isd etc),
e.g. OUTNAME=D:\RIV\ISG_new

Example:

FUNCTION=ISGADJUST
SESFILE=D:\iMOD-DATA\ISGEDIT\ISG-change-stage.SES
LOGFILE=D:\iMOD-DATA\ISGEDIT\ISG-LOG.TXT
OUTNAME=D:\iMOD-DATA\RIV\ISG_new

The example above will produce new ISG files based on new data in the SES file.


8.3.8ISGADDSTRUCTURES-Function

Use this function will add weirs to an ISG.

FUNCTION=

ISGADDSTRUCTURES

ISGFILE_IN=

Enter the name of the ISG file.

IPFFILE_IN=

Enter the name of the IPF file containing weir data

IXCOL=
(optional)

Column number in the specified IPF for the X coordinates, by default IXCOL=1.

IYCOL=
(optional)

Column number in the specified IPF for the Y coordinates, by default IYCOL=2.

IDCOL=
(optional)

Column number in the specified IPF for the ID of the structure, by default IDCOL=3.

IOCOL=
(optional)

Column number in the specified IPF for the orientation of the structure (east is 0.0 degrees, north is 90.0 degrees etc.), by default IOCOL=4. The orientation is used to compute the downstream level of the structure as this not measured. iMOD computes this from the nearest calculation point downstream of the structure. iMOD computes the flow direction from the nearest up- and downstream levels and compares this with the specified angle.

ISCOL=
(optional)

Column number in the specified IPF for the summer upstream level of the structure, by default ISCOL=5.

IWCOL=
(optional)

Column number in the specified IPF for the winter upstream level of the structure, by default IWCOL=6.

MAXDIST=
(optional)

Specify the maximum distance to snap a structure to a nearby segment, by default MAXDIST=1000.

START_YEAR=
(optional)

Specify the start year, by default START_YEAR=1980.

END_YEAR=
(optional)

Specify the end year, by default END_YEAR=2012.

START_PERIOD_
SUMMER=
(optional)

Specify the start date of summer, by default START_PERIOD_SUMMER=’01-04’.

END_PERIOD_
SUMMER=
(optional)

Specify the end date of summer, by default END_PERIOD_SUMMER=’30-09’.

START_PERIOD_
WINTER=
(optional)

Specify the start date of winter, by default START_PERIOD_WINTER=’01-10’.

END_PERIOD_
WINTER=
(optional)

Specify the end date of winter, by default END_PERIO_WINTER=’31-03’.

DATE_WLEVEL=
(optional)

Date of measure to be used to compute undisturbed water level. Give 0-0-0 to compute the mean of all values or 28-02-1994 for a fixed date, by default DATE_WLEVEL=0-0-0.

IPFLOGFILE=
(optional)

Enter a name for the log file, by default IPFLOGFILE=LOG.IPF. In the logfile, the added structure will be logged with an column IOKAY that denotes the status of the addition, IOKAY=1 means the structures is added as is, IOKAY=-1 the angle of structures differs more than 90 degrees with the angle of the segment (might indicate a wrong placement), IOKAY=2 means that the structure existed already and is updated, IOKAY=-2 means that the structures it not positioned as it could not find a nearby segment within MAXDIST.

ICORDIR=
(optional)

Enter ICORDIR=1 to correct the direction of flow (order of coordinates) based upon the angle in the IPF file. If the difference between the specified and computed angle (from the GEN-file) is larger than 90 degrees, the coordinates are flipped. By default ICORDIR=0 and no correction is applied.

ISGFILE_OUT=

Enter the name of the ISG output file.

Example:

FUNCTION=ISGADDSTRUCTURES
ISGFILE_IN=D:\RIV-DATA\LEGGER.ISG
IPFFILE_IN=D:\RIV-DATA\WEIR.IPF
ISGFILE_OUT=D:\RIV-DATA\LEGGER_V2.ISG

The example above will produce new ISG files based on new data in WEIR.IPF.


8.3.9ISGADDSTAGES-Function

Use this function will add water levels to an ISG. Existing water levels wil be overwritten if the entered IPF with water levels has water levels for identical dates as already mentioned in the ISG file.

FUNCTION=

ISGADDSTAGES

ISGFILE_IN=

Enter the name of the ISG file.

IPFFILE=

Enter the name of the IPF file containing water level data via associated TXT-files. The construction of the IPFFILE_IN file must be:

Column 1

X coordinate.

Column 2

Y coordinate.

Column \(i\)

refers to the IEXT variable in the IPF-file to indicate the column with the associated TXT-files, see section 9.7 for more detailed description of the IPF-files.

ISGFILE_OUT=

Enter the name of the ISG output file.

STAGETYPE=
(optional)

Specify STAGETYPE=1 to enter river stages in the IPFFILE in depth values. By default STAGETYPE=0 and stages are entered in absolute values (m+MSL).

ICLEAN=
(optional)

Specify ICLEAN=1 to remove all existing entries for stages in the ISG file prior to adding new entries. By default ICLEAN=0 and existing entries are retained.

Example:

FUNCTION=ISGADDSTAGES
ISGFILE_IN=D:\RIV-DATA\LEGGER.ISG
IPFFILE=D:\RIV-DATA\WATERLEVEL.IPF
ISGFILE_OUT=D:\RIV-DATA\LEGGER_V2.ISG

The example above will produce new ISG files based on new data in WATERLEVEL.IPF.


8.3.10SFRTOISG-Function

This function can be used to convert the output of the SFR package onto an ISG file. In this way the existing functionalities of ISG Edit can be used to inspect and analyse the output of the SFR package.

FUNCTION=

SFRTOISG

ISGFILE_IN=

Enter the name of the ISG file, e.g. SFRFILE_IN=D\MODEL\SFR.ISG.

ISGFILE_OUT=

Enter the name of the ISG file that need to be created, e.g. SFRFILE_OUT=D\MODEL\SFR_RESULT.ISG.

SFRFILE_IN=

Enter the name of the output file of the SFR package, e.g. SFRFILE_IN=D\MODEL\TUT_FSFR.TXT. The file need to contain exactly the amount of calculation points as the ISG files given at ISGFILE_IN. You should always use the ISG files that iMOD created during the export of your model to MF2005 files. This file contains per reach and per time step the following columns:

Column 1:

Layer number.

Column 2:

Row number.

Column 3:

Column number.

Column 4:

Stream number.

Column 5:

Segment number.

Column 6:

Flow into Stream Reach, a positive number means flow in the stream, a negative number would not be feasible.

Column 7:

Flow to Aquifer, negative means that groundwater migrates to surface water (stream acts as a drain), a positive value denotes that surface water infiltrates.

Column 8:

Flow out of Stream Reach, a positive number means flow out of the stream, a negative number would not be feasible.

Column 9:

Overland Runoff.

Column 10:

Precipitation.

Column 11:

Evaporation.

Column 12:

Stream Head, the stage in the river.

Column 13:

Stream Depth, the depth of the river.

Column 14:

Stream Width, the width of the river.

Column 15:

Conductance, the area of contact between the surface water and the groundwater, divided by the resistance of the riverbed.

Column 16:

Gradient, the hydraulic gradient across the stream bed, negative gradient means that the stream acts as a drain. Large hydraulic gradients might indicate some conceptual errors. There is a difference in how the SFR package deals with unsaturated flow underneath the stream bed. As the conventional RIV package applies a hydraulic gradient as the difference between the surface water level and the river bed level, the SFR package takes the thickness of the river bed material instead.

Not all of the above mentioned columns are converted to the ISG file, the average stream discharge (half of the sum of the Flow into Stream Reach plus the Flow out of Stream Reach) is converted from m\(^3\)/d to m\(^3\)/s. Also the stream head, depth and width are transferred to the ISG file. Those four attributes are considered to be most valuable.

Example:

FUNCTION=SFRTOISG
SFRFILE_IN=D\MODEL\SFR.ISG.
SFRFILE_OUT=D\MODEL\SFR_RESULT.ISG.
SFRFILE_IN=D\MODEL\TUT_FSFR.TXT.

The example above will produce a new ISG files based on the results of the SFR package.


8.3.11IPFTOISG-Function

This function can be used to convert an IPF file with appropriate columns into an ISG suitable for the SFR package. Two following pairs of coordinates will form a single segment in the ISG file. Also, only one calculation point and cross-section will appear on each segment. The cross-section will be rectangular (4 points) based on the given stream width (see below).

FUNCTION=

IPFTOISG

IPFFILE=

Enter the name of the IPF file, e.g. IPFFILE=D\DATA\STREAM.IPF.

IXCOL=
(optional)

Enter the column number in the IPF file that represents the X-coordinate (m), e.g. IXCOL=2, by default IXCOL=1.

IYCOL=
(optional)

Enter the column number in the IPF file that represents the Y-coordinate (m), e.g. IYCOL=3, by default IYCOL=2.

ILABELCOL=
(optional)

Enter the column number in the IPF file that represents the stream label, e.g. ILABELCOL=5, by default ILABELCOL=3. This will be used to label the stream in the ISG file. If the label appeared to be empty (not available), the stream segment identification will be used as given by the keyword ISEGMCOL. In the end the segment name in the ISG file becomes “S_{stream label}_R_{\(i\)}” where \(i\) is the sequential order number.

ISEGMCOL=
(optional)

Enter the column number in the IPF file that represents the stream segment identification, e.g. ISEGMCOL=5, by default ISEGMCOL=4. If two following pairs of coordinates in the IPF have similar stream segment identification, they become a segment in the ISG file.

IWIDTHCOL=
(optional)

Enter the column number in the IPF file that represents the stream total width (m), e.g. IWIDTHCOL=6, by default IWIDTHCOL=5.

IBOTTOMCOL=
(optional)

Enter the column number in the IPF file that represents the stream bottom height (m+MSL), e.g. IBOTTOMCOL=8, by default IBOTTOMCOL=6.

ISTAGECOL=
(optional)

Enter the column number in the IPF file that represents the stream stage (m+MSL), e.g. ISTAGECOL=12, by default ISTAGECOL=7.

IPERMCOL=
(optional)

Enter the column number in the IPF file that represents the stream bed permeability (m/d), e.g. IPERMCOL=2, by default IPERMCOL=8.

SDATE=
(optional)

Enter the date for which the data will be recorde in the ISG file, e.g. SDATE=20101231123030 to represent the 31\(^{\rm st}\) of December 2010 at 12:30:30, by default SDATE=20000101000000.

ISGFILE=

Enter the name of the ISG file that need to be created, e.g. SFRFILE=D\MODEL\SFR_RESULT.ISG.

Example:

FUNCTION=IPFTOISG
IPFFILE=D\DATA\STREAM.IPF.
ISGFILE=D\MODEL\SFR.ISG.
ILABELCOL=4
ISEGM=3

The example above will produce a new ISG files compatible with the SFR package, for more clarification an example is given for the IPF file as well. This file will create 2 ISG segment called S_2391526_R_1 and S_2391526_R_2. An example of the IPF file is given below:

8
X
Y
SECTION
NAME
WIDTH
WATERLEVEL
BOTTOMLEVEL
PERM
0,TXT
459420.449184,5495001.214598,252-2391526,"2391526",4.0,90.72400,89.724000,100.0
459327.814399,5495038.775708,252-2391526,"2391526",4.0,90.27100,89.271000,100.0
459235.179516,5495076.336717,252-2391526,"2391526",2.0,89.14100,88.141000,100.0