pictures/imod_logo.png

iMOD User Manual version 5.2 (html)


13.3Output

This chapter describes the output files that are, by default or optionally, produced by iMOD-WQ. The output files can be divided into 6 categories:

These 6 categories are described in the subsequent section.


13.3.1Standard output

The set of standard output files that iMOD-WQ produces consists of the following files:

Other standard output files produced by Seawat, MT3DMS and RT3D are the unformatted heads output file (.HDS), the Observation File (.OBS) and the Model Grid Configuration File (MT3D.CNF). In iMOD-WQ, the unformatted heads output file is not produced anymore as it is replaced by grid output (see section 3.2). The Observation File is replaced by IPF-files (see Section 3.3). The writing of the MT3D.CNF file is suppressed in iMOD-WQ as it is considered obsolete in the iMOD environment.


13.3.2Grid output

Fluxes, heads and concentrations can be saved in iMOD grid format (.IDF) per model layer. iMOD-WQ produces subfolders for the different output types within the output folder specified in the runfile (keyword RESULT_DIR). Table 13.2 presents an overview of the different subfolders.

Folder Data Unit Nodata Produced when Remark
bdgfrf Fluxes over the \(L^3/T\) -9999. ILPFCB > 0 or Negative when
right face of all IBCFCB > 0 flux is directed
model cells in out of the cell
the layer
bdgflf Fluxes over the \(L^3/T\) -9999. ILPFCB > 0 or Negative when
left face of all IBCFCB > 0 flux is directed
model cells in out of the cell
the layer
bdgsto Storage flux in \(L^3/T\) -9999. ILPFCB > 0 or Negative when
all model cells IBCFCB > 0 flux goes into
in the layer the storage
bdgwel Well package \(L^3/T\) -9999. IWELCB > 0 Negative when
flux in all model directed out of
cells in the layer the cell
bdgdrn Drain package \(L^3/T\) -9999. IDRNCB > 0 Negative when
flux in all model directed out of
cell in the layer the cell
bdgriv River package \(L^3/T\) -9999. IRIVCB > 0 Negative when
flux in all model directed out of
cells in the layer the cell
bdgevt Evapotranspiration \(L^3/T\) -9999. IEVTCB > 0 Negative when
package flux in all directed out of
model cells in the the cell
layer
bdgghb GHB package flux \(L^3/T\) -9999. IGHBCB > 0 Negative when
in all model cells in directed out of
the layer the cell
bdgrch Recharge package \(L^3/T\) -9999. IRCHCB > 0 Negative when
flux in all model directed out of
cells in the layer the cell
bdgbnd Constant head flux \(L^3/T\) -9999. Always Negative when
in all model cells in directed out of
the layer the cell
head Groundwater head \(L\) -9999. Always in SEAWAT -
in all model cells mode never in
in the layer MT3DMS mode
conc Concentration in \(L\) -9999. Always -
all model cells in
the layer
dcdt Change in \((M/L^3)/T\) -9999. Always -
concentration
per unit of time

Table 13.2: IDF folders.

Timing/frequency of the grid output is controlled by keywords in the OC package. The table below explains the relevant keywords. The keywords offer the possibility to select stress periods and layers for which output is created. Be aware that output files for budget terms are only written if their corresponding flags are switched on in the respective packages (see also the column named "Produced when" in Table 13.2).

KEY TYPE DEFAULT DESCRIPTION
SAVEHEAD_P?_L? logical F Save head at the end
of the stress period
SAVEBUDGET_P?_L? logical F Save budgets at the end
of the stress period
SAVECONCLAYER_P?_L? logical F Save concentrations at
the end of the stress period

13.3.3Point output

The Observation File type (.OBS) as present in MT3MDS and RT3D is replaced in iMOD-WQ by the iMOD Point File type (.IPF). The output .IPF file is placed in the conc results folder. It comes in the form of an IPF file (max. one per layer) with associated .txt files (one per observation location in the layer) with the time-variant information on concentration at that observation point. The associated .txt files are also places in the conc results folder. The resulting time-series can be visualized in the iMOD GUI (see Figure 13.26). Observation types other than concentrations are not (yet) available.

pictures/Ch-imod_water_quality/pointOutput.png

Figure 13.26: Example of visualizing point output in the iMOD GUI using an output IPF with an associated.txt file that contains the time-dependent output concentrations for four different species.


13.3.4Tecplot output

Tecplot (www.tecplot.com) is a popular 3D visualisation and analysis tool. iMOD-WQ supports writing model output in the format required by Tecplot. Writing heads, concentrations and velocities in the three principal directions is facilitated. Currently, Tecplot output is only available if RUNMODE = SEAWAT. The Tecplot output is managed in the OC package. The table below explains the relevant keywords.

KEY TYPE DEFAULT DESCRIPTION
SAVEHEADTEC_P?_L? logical F Save head at the end of
the stress period (Tecplot)
SAVECONCTEC_P?_L? logical F Save concentration at the
end of the stress period
(Tecplot)
SAVEVXTEC_P?_L? logical F Save velocity in x-direction
at the end of the stress
period (Tecplot)
SAVEVXTEC_P?_L? logical F Save velocity in y-direction
at the end of the stress
period (Tecplot)
SAVEVXTEC_P?_L? logical F Save velocity in x-direction
at the end of the stress
period (Tecplot)
TECFILE character concvelo.tec Name of the output file
for Tecplot
TECITMUN character Y Desired time unit of the output
file. S =seconds, M = minutes,
H = hours, D = days, Y = years

Example:

[OC]

SAVEVXTEC_P?_L? = TRUE

SAVEVYTEC_P?_L? = TRUE

SAVEVZTEC_P?_L? = TRUE

SAVEHEADOUT_P?_L? = TRUE

SAVECONCOUT_P?_L? = TRUE

TECFILE = "SCENARIO_1.TEC"

TECITMUN = D

Figure 13.27: Example settings in the OC package


13.3.5Summary output

iMOD-WQ produces .FLX output files, one per species, which gives for each transport time both the total (water) flux (\(L^3/T\)) and the total load (\(M/T\)) to each type of boundary condition (currently implemented are CHD, WEL, DRN, RIV and GHB). This is particularly useful when RUNTYPE = MT3DMS, as in that case water fluxes are not readily available in other MT3DMS output files. The information is written separately for incoming and outgoing fluxes. The naming of these output files is as follows: MT3D.FLX, and cannot be chosen differently. An example of the contents of such a .FLX file is given below.

TIME Q_CHD_IN M_CHD_IN Q_CHD_OUT M_CHD_OUT Q_WEL_IN M_WEL_IN Q_WEL_OUT ...

0.100000E+01 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.301978E+01 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.503956E+01 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.705934E+01 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.907912E+01 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.110989E+02 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.131187E+02 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.151385E+02 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.171582E+02 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.191780E+02 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

0.211978E+02 0.0E+00 0.0E+00 0.0E+00 0.0E+00 0.313500E+04 0.0E+00 -0.481780E+04 ...

Figure 13.28: Example of the contents of a .FLX file.


13.3.6Standard Seawat input

The runfile translator of iMOD-WQ creates input files for Seawat and MT3DMS that contain iMOD-specific formatting for efficient input handling. This makes that these files cannot be read by standard executables of Seawat and MT3DMS. For checking purposes, however, iMOD-WQ can on demand also produce standard Seawat input files in the formats accepted by Seawat v4. Currently, this functionality has only been implemented for the drn, riv, and .ghb file. The files are filled per stress period during runtime and are located in the subdirectory seawat_ori in the imod-wq_tmp directory.