# Presentation: 066109

## What is UCODE_2005?

## Model CalibrationTrial and Error

- Model Design
- Boundary conditions
- Geometry
- Transmissivity
- Recharge

- “Intelligent” mechanism
- for model adjustment

## Model CalibrationUsing UCODE

- Model CalibrationUsing UCODE

- “Intelligent” mechanism
- for model adjustment

- UCODE_2005
- Modified Gauss Newton

- Model Design
- Boundary conditions
- Geometry
- Transmissivity
- Recharge

## How UCODE Works

- UCODE_2005

## How UCODE Works: 1.adjust input

- UCODE_2005

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

- Template file
- jtf !
- 5.1  2.7  50.
- !K                  ! !w                  !
- 9
- 5.
- 10.
- 15.
- 20.
- 25.
- 30.
- 35.
- 40.
- 45.

- model input file
- 5.1  2.7  50.
- .0100417696   2.723437E-4
- 9
- 5.
- 10.
- 15.
- 20.
- 25.
- 30.
- 35.
- 40.
- 45.

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

- Parameters we want to adjust

## How UCODE Works: 2.run model

- UCODE_2005

## How UCODE Works: 3. get output

- UCODE_2005

## Instruction file
jif @
@      x         head@
l1 [h1]17:37
l2 [h2]17:37 
l2 [h3]17:37  
l2 [h4]17:37  
l2 [h5]17:37 
@      q at x=zero@
l1 [qleft]1:25 
@      q at x=width@
l1 [qright]1:25 

- jif @
- @      x         head@
- l1 [h1]17:37
- l2 [h2]17:37
- l2 [h3]17:37
- l2 [h4]17:37
- l2 [h5]17:37
- @      q at x=zero@
- l1 [qleft]1:25
- @      q at x=width@
- l1 [qright]1:25

- Process model output file
-  Hleft          Hright         Width
- 5.10000        2.70000       50.00000
-  K              RechargeRate
-         0.01004        0.00027
-       x         head
-          5.00000        5.49911
-        10.00000        5.75451
-        15.00000        5.88494
-        20.00000        5.89870
-        25.00000        5.79661
-        30.00000        5.57231
-        35.00000        5.21005
-        40.00000        4.67787
-        45.00000        3.90695
-       q at x=zero
-        -0.00493
-       q at x=width
-         0.00869

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

- Communicating with the model:Extract output using Instruction files

## How UCODE_2005 Works

## UCODE_2005 Capabiities

- UCODE_2005 has many capabilities not discussed in class.
- See the documentation.
- Here we go over those needed for class

## Controlling What UCODE does

- Input blocks
- UCODE_2005 input files can have up to 20 input blocks. See Appendix F.
- The input blocks serve 7 purposes
- Control UCODE_2005 operation
- Define parameters
- Define observations and predictions
- Include measurements of parameter values
- Define weight matrices
- Interact with process model
- Parallel execution

## Input block structure

- The input block can be constructed using one of the three possible block formats
- Keywords (a list of keywords)
- Table (a table of input data)
- File (a file from which data is read)

## Input block: Keywords example

- This is one of the input blocks that serve purpose
- 1: Control UCODE_2005 operation
- BEGIN UCODE_CONTROL_DATA KEYWORDS
-   ModelName=Dupuit_Simple.1
-   ModelLengthUnits=meters
-   ModelTimeUnits=days
-   sensitivities=yes
-   optimize=yes
-   DataExchange=yes
- END UCODE_CONTROL_DATA

##  Input block: table example

- This is one of the input blocks that serve purpose
- 3: Define observations and predictions
- BEGIN OBSERVATION_DATA TABLE
- NROW=7  NCOL=5  COLUMNLABELS
- obsname    obsvalue     statistic  statflag    groupname
-      h1            5.505         5.0e-4      sd            head
-      h2            5.897         5.0e-4      sd            head
-      h3            5.812         5.0e-4      sd            head
-      h4            5.223         5.0e-4      sd            head
-      h5            3.915         5.0e-4      sd            head
-      qleft        -0.005         1.0e-8      var           flow
-      qright       0.0088       6.4e-7      var           flow
- END OBSERVATION_DATA TABLE

## Input blocks: example file list

- BEGIN OBSERVATION_DATA FILES
-    ..\ex1a-files\hed.obs
-    ..\ex1a-files\flo.obs
- END OBSERVATION_DATA
- The files would use either Keywords or  Table to define the format of the input block

##  Input

- Defaults available for many things. No input is needed to use the defaults.
- Unknown keywords are IGNORED!!
-    Be careful!
-    In the Options input block use Verbose=5 to check what UCODE_2005 is seeing

## Data-Exchange Files

- There are lots of data-exchange files. Alphabetically listed in appendix B. Listed by topic in tables of Chapters 14-17.
- Use for post-processing
- GW_CHART (Richard Winston)

**UCODE_2005 documentation, **
**Appendix B, p. 231-233**

## UCODE_2005 Modes Use UCODE_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)

##  Weighting

- Weights have to make contributions to the objective function dimensionless
- Observations with different units can be combined
- Weights can be used to emphasize reliable observations compared to less reliable ones
- Can use standard deviations or coefficients of variation (cv) to define the weight.

## SA/PE/UA

- Sensitivity Analysis, Parameter Estimation, and Uncertainty Analysis
- Insight into how data and models interact
- Use quantification provided to communicate data importance and needs to decision makers
- Provides measure of uncertainty