Syntax
>>> ELEMENT: eleme (eleme_i ...) (+ iplus ([islotb-islote]))
or
>>> ELEMENT COORDINATE: X Y Z
or
>>> ELEMENT COORDINATE BOX (ROTATE) (OUTSIDE): Xmin Ymin Zmin Xmax Ymax Zmax
or
>>> ELEMENT COORDINATE ELLIPSOID (ROTATE) (OUTSIDE): Xcenter Ycenter Zcenter Xsemiaxis Ysemiaxis Zsemiaxis
or
>>> ELEMENT COORDINATE CYLINDER (ROTATE) (OUTSIDE): Xstart Ystart Zstart Xend Yend Zend Radius
or
>>> ELEMENT COORDINATE CUBE (ROTATE) (OUTSIDE): Xcenter Ycenter Zcenter Xsemilength Ysemilength Zsemilength
or
>>> ELEMENT COORDINATE PRISM (OUTSIDE) AXIS: ijkAXIS POLYGON: NP : Axismin Axismax C1,1 C1,2 ... CNP,1 CNP,2
or
>>> ELEMENT PROFILE/MAP/CROSS-SECTION
or
>>> MATERIAL: material (material2 ...)
Parent Command
all second-level commands in block > OBSERVATION requiring element names.
Subcommand
all fourth-level commands in block > OBSERVATION
Description
This command reads one or more element names. Most observation types refer to a variable that is associated with a grid block (as opposed to a connection or sink/source name). Element names are designated by a three-character/two-integer (FORTRAN format: AAAII) code name. Blanks in the element names as printed in the TOUGH2 output file must be replaced by underscores (e.g. an element name specified in the TOUGH2 input file as ‘B 007’ is printed as ‘B 0 7’ in the TOUGH2 output file. Therefore, it must be addressed in the iTOUGH2 input file as ‘B_0_7’). Multiple elements can be specified, and iTOUGH2 calculates the sum or mean of the corresponding output variable (see subcommands >>>> SUM and >>>> MEAN, respectively). A sequence of iplus elements can be generated by increasing the number found in substring [islotb-islote] (default: [4-last]) of the last element. The following two command lines are identical:
>>> ELEMENT: AA__1 BB_15 +3 [4-5]
>>> ELEMENT: AA__1 BB_15 BB_16 BB_17 BB_18
If keyword COORDINATES is present, coordinates can be specified, and the element closest to these coordinates will be identified.
If keywords COORDINATE BOX, COORDINATE ELLIPSOID, COORDINATE CYLINDER, COORDINATE CUBE, or COORDINATE PRISM POLYGON: NP AXIS: ijkAXIS (where NP is the number of points defining the polygon, and ijkAXIS = 1, 2, or 3 defines the alignement of the prism axis with the X, Y, or Z coordinate axis, respectively) are present, all elements within the specified region will be included in the element list. The region is defined according to the following table. If keyword COMPLEMENT or OUTSIDE is present, the complementary region is selected (i.e., all elements outside the defined geometry). Element coordinates must be present in columns 51-80 of the TOUGH2 ELEME block.
| KEYWORD | i=1 | i=2 | i=3 | i=4 | i=5 | i=6 | i=7 | i=8 | i=9 |
|---|---|---|---|---|---|---|---|---|---|
| (none) | X | Y | Z | – | – | – | – | – | – |
| BOX | X-min | Y-min | Z-min | X-max | Y-max | Z-max | azimuth | dip | plunge |
| ELLIPSOID | X-center | Y-center | Z-center | X-semi-axis | Y-semi-axis | Z-semi-axis | azimuth | dip | plunge |
| CYLINDER | X-start | Y-start | Z-start | X-end | Y-end | Z-end | radius | – | – |
| CUBE | X-center | Y-center | Z-center | X-half-length | Y-half-length | Z-half-length | azimuth | dip | plunge |
| PRISM | Axismin | Axismax | C1,1 | C1,2 | Ci,1 | Ci,2 | CNP,1 | CNP,2 |
Add keyword ANGLE or ROTATE if the region is not aligned with the coordinate axes, and provide azimuth, dip, and plunge correction angles.
Keywords CONSTANT (default), LINEAR, SPHERICAL, EXPONENTIAL, or RAMP may be added to weigh the contribution of each observation according to its distance from the center of the region.
If keyword PROFILE, CROSS-SECTION, or MAP is present, the data (which must be provided in an external file, see command >>>> DATA FILE: filename) are organized spatially by lists of coordinates followed by columns of measured values for each survey time. The data will be assigned to the nearest element (element coordinates must be present in columns 51-80 of the TOUGH2 ELEME block). Keyword PROFILE supports one-dimensional, vertical data profiles (i.e., Z-coordinates), whereas CROSS-SECTION and MAP both require X-, Y-, and Z-coordinates. Profile data are given as follows:
NZ NT
XCoord YCoord
Time(1) Time(2) ... Time(NT)
Z(1) Data(1,1) Data(1,2) ... Data(1,NT)
Z(2) Data(2,1) Data(2,2) ... Data(2,NT)
... ... ... ... ...
Z(NZ) Data(NZ,1) Data(NZ,2) ... Data(NZ,NT)
Two- or three-dimensional data are given as follows:
NP NT
Time(1) Time(2) ... Time(NT)
X(1) Y(1) Z(1) Data(1,1) Data(1,2) ... Data(1,NT)
X(2) Y(2) Z(2) Data(2,1) Data(2,2) ... Data(2,NT)
... ... ... ... ... ... ...
X(NP) Y(NP) Z(NP) Data(NP,1) Data(NP,2) ... Data(NP,NT)
If command >>> MATERIAL is used, all elements assigned to the specified material(s) will be selected.
Example
> OBSERVATION
>> GAS SATURATION
>>> ELEMENTS: ELM_0 + 99
>>>> ANNOTATION : Mean saturation
>>>> take the MEAN of the saturation in all 100 elements
>>>> DATA on FILE : Sg.dat
>>>> DEVIATION : 0.05
<<<<
<<<
>> CAPILLARY pressure
>>> ELEMENT COORDINATE BOX: 10. 0.4 -5.0 40.0 0.5 -1.0
>>>> ANNOTATION : Mean Pcap
>>>> take the VOLUMEtric MEAN of all elements in box
>>>> DATA on FILE : Pc.dat
>>>> DEVIATION : 1.0E4
<<<<
<<<
>> MASS FRACTION
>>> MATERIAL: CLAY1 CLAY6
>>>> COMPONENT: 4
>>>> PHASE: 2
>>>> AVERAGE
<<<<
<<<
<<
See Also
>>> CONNECTION | >>> SINK | >>>> MEAN | >>>> SUM