# INPUT FOR MAIN OPTIMIZATION FILE # Additional comments possible by starting line with # sign # All entry beyond a ! are comments that are not read #------------------------------------------------------- # DATA GROUP 1 - FILE SPECIFICATIONS # INPUT_CHTEST.INP.txt 9 !! Composite MODFLOW-MODOPTIM file OUTPUT_CHTEST.MAIN-OPT.txt 8 !! Main MODOPTIM output file OUTPUT_CHTEST.Resi&Sens.txt 85 !! MODOPTIM residuals and sensitivities file associated with lowest SS error. OUTPUT_CHTEST.BAS.txt 5 !! MODFLOW-BAS file is generated by MODOPTIM OUTPUT_CHTEST.MODFLOW-List.txt 6 !! Standard MODFLOW output OUTPUT_CHTEST.Parameters.txt 86 !! MODOPTIM file with parameter estimates and statistics associated with lowest SS error. # #------------------------------------------------------- # DATA GROUP 2 - CLOSURE CRITERIA FOR OPTIMIZATION # maxit=15 smin= 0.0001 dermin= 0.0100 CCmin > 0.950 # # maxit is maximum number of parameter estimation iterations # smin is the minimum overall parameter change # dermin is the minimum reduction in Sum-of-Squares objective fuction # CCmin is the minimum correlation coefficient shown for parameter pairs #------------------------------------------------------- # DATA GROUP 3 - OPTIMIZATION CONTROL VARIABLES # smax= 0.500 stfrc= 0.500 frcsec= 0.010 rkmin= 0.005 cfi= 0.010 LIMIT:Keep # # smax is the maximum step length for all parameters, typical value is 0.5 # stfrc is the maximum fraction any given parameter can # change between its lower and upper limits, typical value is 0.25 # frcsec determines LM value for damping Hessian after BFGS # update fails, typical value is 0.1 for low damping and 1.0 for high damping # rkmin is the threshold below which a diagonal value is # treated as zero, typical value is 0.0001 # cfi is the pertubation factor for calculating gradients, typical value is 0.01 # #------------------------------------------------------- # DATA GROUP 4 - OPTIONAL - DEFINITION OF WEIGHT MATRICIES # # If a weight matrix is not read, no weight matrices will be defined # and estimated parameters cannot be modified by a weight matrix # #------------------------------------------------------- # DATA GROUP 5 - DEFINITION OF PARAMETERS # # Begin defining Parameters to be optimized # # BEST ESTIMATES # Parameter Unique Tag Type Parameter Estimate # 1 kICU PERM 5.000000E-02 # 2 kUFA PERM 3.800000E+01 # 3 kMSU PERM 0.300000E+00 # 4 kWEL PERM 0.300000E+10 # 5 aUFA ANIS 0.100000E+00 # 6 sICU STOR 1.000000E-06 # 7 sUFA STOR 5.000000E-06 # 8 sMSU STOR 3.000000E-06 # 9 sWEL STOR 0.100000E+01 # # ______________ Hydraulic Conductivity ________________ # Parameter:kICU Tag:PERM Status:Modify Type:Cnst,Vol,Exp,NON Lower:5.000000E-03 Initial:5.000000E-02 upper:5.000000E-01 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 2 )[All] !! Intermediate Confining Unit # Parameter:kUFA Tag:PERM Status:Modify Type:Cnst,Vol,Exp,NON Lower:3.800000E+00 Initial:3.800000E+01 upper:3.800000E+02 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 3 )[All] !! Upper Floridan Aquifer # Parameter:kMSU Tag:PERM Status:Modify Type:Cnst,Vol,Exp,NON Lower:0.300000E-01 Initial:0.300000E+00 upper:0.300000E+01 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 4 )[All] !! Middle Semi-Confining Unit # Parameter:kWEL Tag:PERM Status:Modify Type:Cnst,Vol,Exp,NON Lower:0.100000E-01 Initial:0.100000E+01 upper:0.100000E+03 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 0.300000E+10*( 6 )[All] !! WELLBORE # # _________ Vertical-to-Horizontal Anisotropy __________ # Parameter:aUFA Tag:ANIS Status:Modify Type:Cnst,Vol,Exp,NON Lower:0.100000E-02 Initial:0.100000E+00 upper:0.100000E+01 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 3 )[All] !! Anisotropy in Penetrated Aquifer # # _________________ Specific Storage ___________________ # Parameter:sICU Tag:STOR Status:Modify Type:Cnst,Vol,Exp,NON Lower:1.000000E-07 Initial:1.000000E-06 upper:1.000000E-05 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 2 )[All] !! Intermediate Confining Unit # Parameter:sUFA Tag:STOR Status:Modify Type:Cnst,Vol,Exp,NON Lower:5.000000E-07 Initial:5.000000E-06 upper:5.000000E-05 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 3 )[All] !! Upper Floridan Aquifer # Parameter:sMSU Tag:STOR Status:Modify Type:Cnst,Vol,Exp,NON Lower:3.000000E-07 Initial:3.000000E-06 upper:3.000000E-05 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 4 )[All] !! Middle Semi-Confining Unit # Parameter:sWEL Tag:STOR Status:Modify Type:Parallel,None,Exp,NON Lower:0.100000E+00 Initial:0.100000E+01 upper:0.100000E+02 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 0.387597E-02*( 6 )[All] !! WELLBORE Multiplier is inverse of well bore length = 1/(SCBOT-SCTOP) # ! End Parameter Input # #------------------------------------------------------- # DATA GROUP 6 - GENERAL OBSERVATION MODIFIERS # Time interpolation # # Flag that allows simulated water levels to be interpolated in time to the # observation time # Insert Blanks between observations # # Flag that indicates blank lines should be printed between sites in the ouput. # default is no blank lines # START:0.0 SCALE:1.000 # # Indicates time 0 of simulation in time units of data input. Default is 0. # Indiates scale unit 1 of data input used within MODFLOW simulation. Default is 1. # Z-Interpolate:No # # Simulated water levels can either be or not be interpolated between layers # to observed location #------------------------------------------------------- # DATA GROUP 7 - OBSERVATIONS - BEGIN READING OBSERVATIONS # # HEAD is the 4 character discriptor that identifies values to be compared # IJK: is the discrete location of the observations. Choices are: # XYZ: Cart-X, Cart-Y,Cart-Z # XYK: Cart-X, Cart-Y, Layer # XJZ: Cart-X, Row, Cart-Z # XJK: Cart-X, Row, Layer # IYZ: Column, Cart-Y, Cart-Z # IYK: Column, Cart_Y. Layer # IJZ: Column, Row, Cart-Z # IJK: Column, Row, Layer # SITE: is a descriptor used to identify an observation site # Add BLANKS - Insert a blank line between time-drawdown series from each observation well. # UFA1 XYK: 430.000 -592.500 1 SITE:UFA1_OBS_wel # Time, Days s, ft Weight 0.00278 0.11 1.000 W-Default !! 0 HOURS, 4 MINs., 0 SECs AFTER PUMPING Redirect to file:INPUT_UFA1.OBS.txt # ICU1 XYK: 100.000 -252.500 1 SITE:ICU1_OBS_wel # Time, Days s, ft Weight 0.01250 0.14 1.000 W-Default !! 0 HOURS, 18 MINs., 0 SECs AFTER PUMPING Redirect to file:INPUT_ICU1.OBS.txt # MSC1 XYK: 200.000 -787.500 1 SITE:MSC1_OBS_wel # Time, Days s, ft Weight 0.00694 0.11 1.000 W-Default !! 0 HOURS, 10 MINs., 0 SECs AFTER PUMPING Redirect to file:INPUT_MSC1.OBS.txt #