# 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= .0010 dermin= .0100 CCmin > .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= .500 stfrc= .500 frcsec= .010 rkmin= .005 cfi= .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 .486686E-01 # 2 kUFA PERM .364524E+02 # 3 kMSU PERM .336445E+00 # 4 kWEL PERM .300000E+10 # 5 aUFA ANIS .100000E+00 # 6 sICU STOR .977969E-06 # 7 sUFA STOR .499982E-05 # 8 sMSU STOR .326127E-05 # 9 sWEL STOR .100000E+01 # # ______________ Hydraulic Conductivity ________________ # Parameter:kICU Tag:PERM Status:Active Type:Cnst,Vol,Exp,NON Lower: .486686E-03 Initial: .486686E-01 upper: .486686E+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:Active Type:Cnst,Vol,Exp,NON Lower: .304524E+02 Initial: .384524E+02 upper: .504524E+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:Active Type:Cnst,Vol,Exp,NON Lower: .336445E-02 Initial: .336445E+00 upper: .336445E+02 # 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: .100000E-01 Initial: .100000E+01 upper: .100000E+03 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis .300000E+10*( 6 )[All] !! WELLBORE # # _________ Vertical-to-Horizontal Anisotropy __________ # Parameter:aUFA Tag:ANIS Status:Active Type:Cnst,Vol,Exp,NON Lower: .100000E-02 Initial: .100000E+00 upper: .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:Active Type:Cnst,Vol,Exp,NON Lower: .977969E-07 Initial: .977969E-06 upper: .977969E-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:Active Type:Cnst,Vol,Exp,NON Lower: .499982E-06 Initial: .499982E-05 upper: .499982E-04 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 3 )[All] !! Upper Floridan Aquifer # Parameter:sMSU Tag:STOR Status:Active Type:Cnst,Vol,Exp,NON Lower: .326127E-06 Initial: .326127E-05 upper: .326127E-04 # 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: .100000E+00 Initial: .100000E+01 upper: .100000E+02 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis .392157E-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 .00278 .12 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 .06944 .11 1.000 W-Default !! 1 HOURS, 39 MINs., 59 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 .00833 .19 1.000 W-Default !! 0 HOURS, 11 MINs., 59 SECs AFTER PUMPING Redirect to file:INPUT_MSC1.OBS.txt #