eSTOMP

Observed Data Card Options (e-W)

This card allows the user to input observed values at specific location and time. To have this card take effect, at the first line of the Solution Control Card, the phrase “inverse” is needed before the first comma (e.g., Normal w/inverse). eSTOMP will then generate two of the UCODE (Poeter and Hill, 1998) input files, i.e., the universal file out_uc1.uni and the extraction file out_uc1.ext, and an output file out_uc1.sto. These files follow the convention described in Poeter and Hill (1998). It is noted that, to have the out_uc1.uni and out_uc1.ext files work for UCODE, a UCODE Control Card is also needed. Without the UCODE Control Card, the file out_uc1.uni contains only the data portion and is incomplete.

The out_uc1.sto file is needed by UCODE and may be used for other analyses. It includes all the simulated values corresponding to the observations along with the IDs of the observation data. Each data ID consists of a three-character acronym representing the data type, followed by 6 digits representing the data number. Accepted variable names and their corresponding abbreviations in the output files are summarized in the Observed Data Variable Options section below. The data numbers are assigned consecutively and sequentially by eSTOMP starting from 000001 based on the observation time. For observations at the same time, the data numbers are assigned in the data type order as listed in the Observed Data Card.

Tri-linear interpolation (https://en.wikipedia.org/wiki/Trilinear_interpolation) was used to calculate the simulated value v at the position where the observation was taken based on the values of 8 neighbor cells (see figure below).  The Tri-linear interpolation may be weighted by using different weights (between 0 and 1) to the three directions and time. (The current version of eSTOMP only uses the default weight of 1).

Observation Location Specification Options

By Distance (Field)

Observations can be specified by their coordinate location within the model domain but cannot be at boundary cells that have less than 8 neighbors.

By IJK (Reference)

Observations can be specified by their IJK location within the model domain but cannot be at boundary cells that have less than 8 neighbors.

 

Observed Data Variable Options

Variable Name SI Units Output Variable Abbreviation Qualification
Aqueous Pressure Pa pl_
Aqueous Saturation - sl_
 Aqueous Moisture Content - mcl
 Aqueous Hydraulic Head m hhl
 Matric Potential m mph
 x Aqueous Volumetric Darcy Flux m/s xql
 y Aqueous Volumetric Darcy Flux m/s yql
 z Aqueous Volumetric Darcy Flux m/s zql
 Solute Volumetric Concentration 1/L cv[n]
 Solute Aqueous Concentration 1/L cl[n]
 x Solute Flux 1/m^2 s xc[n]
 y Solute Flux 1/m^2 s yc[n]
 z Solute Flux 1/m^2 s zc[n]

Bold indicates actual keyword(s) that must be specified in the Observed Data Card.

[n] indicates solute number

 

Reference: 

Poeter EP, and MC Hill. 1998. Documentation of UCODE, a computer code for universal inverse modeling. U.S. Geological survey, water-resources investigations report 98-4080. Available at https://pubs.er.usgs.gov/publication/wri984080 (Last accessed on July 10, 2019).

eSTOMP User Guide Home

Modules