What is UCODE_2005?What is UCODE_2005?

Model CalibrationModel CalibrationTrial and ErrorTrial and Error

Data

Model OutputModel DesignBoundary conditions GeometryTransmissivityRecharge

“Intelligent” mechanism for model adjustment

Compare

Model CalibrationModel CalibrationUsing UCODEUsing UCODE

Data

Model Output

“Intelligent” mechanism for model adjustment

UCODE_2005 Modified Gauss Newton

Model DesignBoundary conditions GeometryTransmissivityRecharge

Compare

How UCODE WorksHow UCODE Works

UCODE_2005

Model Input Files

TemplateFile

creates model input files

Model executes executes the model

Model Output Files

InstructionFile

reads model output files

How UCODE Works: 1.adjust inputHow UCODE Works: 1.adjust input

UCODE_2005

Model Input Files

TemplateFile

creates model input filesTemplate file

Model executes executes the model

Model Output Files

InstructionFile

reads model output files

Communicating with the model:Communicating with the model:Creating input using Template filesCreating input using Template files

Template fileTemplate filejtf !jtf ! 5.1 2.7 505.1 2.7 50. .

!K ! !w !!K ! !w !9 9 5. 5. 10. 10. 15. 15. 20. 20. 25. 25. 30. 30. 35. 35. 40. 40. 45.45.

model input filemodel input file

5.1 2.7 50. 5.1 2.7 50.

.0100417696 .0100417696 2.723437E-42.723437E-49 9 5. 5. 10. 10. 15. 15. 20. 20. 25. 25. 30. 30. 35. 35. 40. 40. 45.45.

Same as in PEST, via John Doherty’s involvement in the JUPITER API

Parameters we want to adjust

How UCODE Works: 2.run modelHow UCODE Works: 2.run model

UCODE_2005

Model Input Files

TemplateFile

creates model input filesTemplate file

Model executes executes the model

Model Output Files

InstructionFile

reads model output files

How UCODE Works: 3. get outputHow UCODE Works: 3. get output

UCODE_2005

Model Input Files

TemplateFile

creates model input files

Model executes executes the model

Model Output Files

InstructionFile

reads model output filesTemplate file

Instruction fileInstruction filejif @jif @

@ x head@@ x head@

l1 [h1]17:37l1 [h1]17:37

l2 [h2]17:37 l2 [h2]17:37

l2 [h3]17:37 l2 [h3]17:37

l2 [h4]17:37 l2 [h4]17:37

l2 [h5]17:37 l2 [h5]17:37

@ q at x=zero@@ q at x=zero@

l1 [qleft]1:25 l1 [qleft]1:25

@ q at x=width@@ q at x=width@

l1 [qright]1:25 l1 [qright]1:25

Process model output fileProcess model output file Hleft Hright WidthHleft Hright Width5.10000 2.70000 50.000005.10000 2.70000 50.00000 K RechargeRateK RechargeRate 0.01004 0.000270.01004 0.00027 x head x head 5.00000 5.00000 5.499115.49911 10.00000 5.7545110.00000 5.75451 15.00000 15.00000 5.884945.88494 20.00000 5.8987020.00000 5.89870 25.00000 25.00000 5.796615.79661 30.00000 5.5723130.00000 5.57231 35.00000 35.00000 5.210055.21005 40.00000 4.6778740.00000 4.67787 45.00000 45.00000 3.906953.90695 q at x=zeroq at x=zero -0.00493-0.00493 q at x=widthq at x=width 0.008690.00869

Same as in PEST, plus Standard File for values in a continuous column

Communicating with the model:Communicating with the model:Extract output using Instruction filesExtract output using Instruction files

START

Evaluate,Cleanup

2. Generate parameter values

6. Use extracted values

1. Define job

Con

trol

Loo

p

UCODE_2005 output, including data-exchange files

Main input files

Process model

Process-model

input files 3. Adapt parameter values

5. Extract model-calculated values

4. Execute process model

Initialize

Process-model

output files

Instruction files +

Template files

END

Use

Input files.

Communication links with the process model and its input and output files.

UCODE_2005 tasks. Numbered tasks form the control loop.

Bold italic letters identify items with programming unique to an UCODE_2005. Others are accomplished by modules of the JUPITER API or are part of the process model.

How UCODE_2005 How UCODE_2005 WorksWorks

UCODE_2005 CapabiitiesUCODE_2005 Capabiities

UCODE_2005 has many capabilities not UCODE_2005 has many capabilities not discussed in class.discussed in class.

See the documentation.See the documentation.

Here we go over those needed for classHere we go over those needed for class

Controlling What UCODE doesControlling What UCODE does

Input blocksInput blocks– UCODE_2005 input files can have up to 20 UCODE_2005 input files can have up to 20

input blocks. See Appendix F.input blocks. See Appendix F.– The input blocks serve 7 purposesThe input blocks serve 7 purposes

1.1. Control UCODE_2005 operationControl UCODE_2005 operation

2.2. Define parametersDefine parameters

3.3. Define observations and predictionsDefine observations and predictions

4.4. Include measurements of parameter valuesInclude measurements of parameter values

5.5. Define weight matricesDefine weight matrices

6.6. Interact with process modelInteract with process model

7.7. Parallel executionParallel execution

Input block structureInput block structure

The input block can be constructed using The input block can be constructed using one of the three possible block formatsone of the three possible block formats

Keywords (a list of keywords)Keywords (a list of keywords)

Table (a table of input data)Table (a table of input data)

File (a file from which data is read)File (a file from which data is read)

Input block: Keywords exampleInput block: Keywords example

This is one of the input blocks that serve purpose This is one of the input blocks that serve purpose 1: Control UCODE_2005 operation1: Control UCODE_2005 operation

BEGINBEGIN UCODE_CONTROL_DATAUCODE_CONTROL_DATA KEYWORDSKEYWORDS ModelName=Dupuit_Simple.1 ModelName=Dupuit_Simple.1 ModelLengthUnits=metersModelLengthUnits=meters ModelTimeUnits=daysModelTimeUnits=days sensitivities=yessensitivities=yes optimize=yesoptimize=yes DataExchange=yesDataExchange=yesENDEND UCODE_CONTROL_DATAUCODE_CONTROL_DATA

Input block: table exampleInput block: table exampleThis is one of the input blocks that serve purpose This is one of the input blocks that serve purpose 3: Define observations and predictions3: Define observations and predictions

BEGINBEGIN OBSERVATION_DATAOBSERVATION_DATA TABLETABLENROW=7 NCOL=5 COLUMNLABELS NROW=7 NCOL=5 COLUMNLABELS obsname obsvalue statistic statflag groupnameobsname obsvalue statistic statflag groupname h1 5.505 5.0e-4 sd head h1 5.505 5.0e-4 sd head h2 5.897 5.0e-4 sd head h2 5.897 5.0e-4 sd head h3 5.812 5.0e-4 sd head h3 5.812 5.0e-4 sd head h4 5.223 5.0e-4 sd head h4 5.223 5.0e-4 sd head h5 3.915 5.0e-4 sd head h5 3.915 5.0e-4 sd head qleft -0.005 1.0e-8 var flow qleft -0.005 1.0e-8 var flow qright 0.0088 6.4e-7 var flow qright 0.0088 6.4e-7 var flow ENDEND OBSERVATION_DATAOBSERVATION_DATA TABLETABLE

Input blocks: example file listInput blocks: example file list

BEGINBEGIN OBSERVATION_DATA OBSERVATION_DATA FILESFILES

..\ex1a-files\hed.obs..\ex1a-files\hed.obs

..\ex1a-files\flo.obs..\ex1a-files\flo.obs

ENDEND OBSERVATION_DATAOBSERVATION_DATA

The files would use either Keywords or The files would use either Keywords or Table to define the format of the input Table to define the format of the input blockblock

InputInput

Defaults available for many things. No Defaults available for many things. No input is needed to use the defaults.input is needed to use the defaults.

Unknown keywords are IGNORED!! Unknown keywords are IGNORED!!

Be careful! Be careful!

In the In the Options input blockOptions input block use use Verbose=5Verbose=5 to check what UCODE_2005 is seeingto check what UCODE_2005 is seeing

Data-Exchange FilesData-Exchange Files

There are lots of data-exchange files. There are lots of data-exchange files. Alphabetically listed in appendix B. Listed Alphabetically listed in appendix B. Listed by topic in tables of Chapters 14-17.by topic in tables of Chapters 14-17.

Use for post-processingUse for post-processing– GW_CHART (Richard Winston)GW_CHART (Richard Winston)

Table B- 1. Files produced by UCODE_2005, RESIDUAL_ANALYSIS, LINEAR_UNCERTAINTY, and MODEL_LINEARITY named using the fn prefix specified on the command line, in alphabetic order by letter in the file extension.

[File extensions that begin with an underscore identify data-exchange files and are listed first. File extensions that begin with # are main output files and are listed at the end of the table. Shading identifies files not produced by UCODE_2005. ‘iteration’ refers to parameter-estimation iteration. * is replaced by ‘conf’ for results related to confidence intervals or ‘pred’ for results related to prediction intervals.]

Extension Brief description (see also tables 16-22, 24-28, 34, 35) Content1 Use1

_b1 Parameter sets for _b2 19 -- _b1adv* Parameter sets for _b2adv* 35 34 _b2 Values simulated using parameter sets from _b1 19 --

_b2adv* Values simulated using parameter sets from _b1adv* 35 34 _b3* Parameter sets for _b4* 35 34 _b4* Values simulated using parameter sets from _b3* 35 34 _cf* Correction factors 35 34 _cfsu Sensitivities used by CORFAC_PLUS 35 34 _dm Information about model structure, fit and parsimony 21 -- _dmp Information about predictions 21 -- _gm Observation groups 16 -- _gmp Prediction groups 20 -- _init As for _paopt evaluated at other parameter values 15 -- _init._mv As for the _mv file evaluated at other parameter values 15 -- _init._su As for the _su file evaluated at other parameter values 15 -- _init._supri As for the _supri file evaluated at other parameter values 15 -- _int* Nonlinear intervals 35 34 _int*par Parameter values at nonlinear interval limits 35 34 _int*wr Weighted residuals at nonlinear interval limits _linp Predictions and their linear confidence intervals 26 -- #linunc LINEAR_UNCERTAINTY main output file 26 -- _mc Parameter correlation coefficient matrix 19 30 #modlin MODEL_LINEARITY main output file 27 -- #modlinadv MODEL_LINEARITY_ADV main output file _mv Parameter variance-covariance matrix 19 30 _nm Weighted residuals, probability plotting positions 16 31 _os Unweighted simulated and observed or prior values 16 31 _p Predictions 20 --

UCODE_2005 documentation, Appendix B, p. 231-233

UCODE_2005 ModesUCODE_2005 Modes Use Use UCODE_Control_Data input blockUCODE_Control_Data input block

Modes run Independently. Typically proceed in this order

Forward

Sensitivity-analysis

Parameter-estimation

Modes that need to follow a successful parameter estimation run

Test-model-linearity (modified Beale’s measure)

Prediction

Modes that need to be coordinated with runs of other modes

Advanced-test-model-linearity

Nonlinear-uncertainty

Mode run independently

Investigate-objective-function (ModelViewer. Winston, Hsieh)

WeightingWeightingWeights have to make contributions to the objective Weights have to make contributions to the objective function dimensionlessfunction dimensionless– Observations with different units can be combinedObservations with different units can be combined

Weights can be used to emphasize reliable Weights can be used to emphasize reliable observations compared to less reliable onesobservations compared to less reliable ones– Can use standard deviations or coefficients of variation Can use standard deviations or coefficients of variation

(cv) to define the weight.(cv) to define the weight.

SA/PE/UASA/PE/UA

Sensitivity Analysis, Parameter Estimation, Sensitivity Analysis, Parameter Estimation, and Uncertainty Analysisand Uncertainty Analysis– Insight into how data and models interactInsight into how data and models interact– Use quantification provided to communicate Use quantification provided to communicate

data importance and needs to decision data importance and needs to decision makersmakers

– Provides measure of uncertaintyProvides measure of uncertainty