Syntax
>>> EMPIRICAL ORTHOGONAL FUNCTIONS (MATRIX: ndim) (iTOUGH2) (CORRELATION)
Parent Command
>> ERROR
Subcommand

Description
This command must be used in combination with command >>> MONTE CARLO to invoke stochastic simulations of correlated parameters using Empirical Orthogonal Functions (EOF). EOF is a variant of Monte Carlo simulations to quantify the uncertainty of model predictions as a result of parameter uncertainty. Many parameter sets are generated, following the predefined covariance matrix C_{pp}. The ith parameter set Y_{i} is obtained by linear combination of the eigenvectors u_{k} of C_{pp} and stochastic coefficients phi(x_{i}):
where xsi is a standard normal distributed random variable, and a_{k} is the kth eigenvalue of C_{pp}. For more details, see KitterÃ¸d and Gottschalk, [1997].
The elements of matrix C_{pp} can be supplied using one of the following options:
(1) Provide indices and elements of C_{pp}. Example:
>>> EMIRICAL ORTHOGONAL FUNCTIONS
1 1 0.80643E04
2 2 0.71921E04
2 1 0.64412E04
Use keyword CORRELATION if offdiagonal term is correlation coefficient instead of covariance. Example:
>>> EMPIRICAL ORTHOGONAL FUNCTIONS, CORRELATION
1 1 0.80643E04
2 2 0.71921E04
2 1 0.864
(2) Provide keyword MATRIX, followed by a colon and the dimension ndim of the square matrix C_{pp}. The lower triangle of the covariance matrix is then provided on exactly ndim additional lines. If keyword CORRELATION is present, the off diagonal terms represent correlation coefficients rather than covariances. Example:
>>> EMPIRICAL ORTHOGONAL FUNCTIONS, dim. of CORRELATION MATRIX: 3
.80643E04  
.864  .71921E04  
.253  .500  .53843E05 
(3) If calculated during a previous iTOUGH2 inversion, the covariance matrix can be taken from the iTOUGH2 output file and directly copied after the command line. This option is invoked by keyword iTOUGH2. The matrix will be read by formatted input, so it is crucial that the correct format is maintained. If ndim is greater than 6, the matrix is split in multiple submatrices. All submatrices must be copied exactly as they were printed to the iTOUGH2 output file. Example:
>>> EOF error analysis, read MATRIX of dim.: 3 in iTOUGH2 format
log(abs. perm.)  POROSITY SAND  Gas entrapped  
log(abs. perm.)  .80643E04  .846  .253 
POROSITY SAND  .64412E04  .71921E04  .500 
Gas entrapped  .52623E05  .98296E05  .53843E05 
Example
> COMPUTATION
>> STOP
>>> number of Monte Carlo SIMULATIONS: 250
<<<
>> ERROR propagation analysis
>>> MONTE CARLO simulations, SEED number: 777
>>> EMPIRICAL ORTHOGONAL FUNCTIONS, read MATRIX of dim.: 3 in iTOUGH2 format
log(abs. perm.)  POROSITY SAND  Gas entrapped  
log(abs. perm.)  .80643E04  .846  .253 
POROSITY SAND  .64412E04  .71921E04  .500 
Gas entrapped  .52623E05  .98296E05  .53843E05 
<<<
<<
<
See Also
>>> FOSM  >>> MONTE CARLO