# 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= .000100 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 .153115E-01 # 2 kUFA PERM .485236E+02 # 3 kMSU PERM .105588E+00 # 4 kLFA PERM .413734E+02 # 5 kWEL PERM .300000E+10 # 6 aLFA ANIS .100000E+00 # 7 sICU STOR .360099E-06 # 8 sUFA STOR .278746E-05 # 9 sMSU STOR .151061E-05 # 10 sLFA STOR .668044E-05 # 11 sWEL STOR .100000E+01 # # ______________ Hydraulic Conductivity ________________ # Parameter:kICU Tag:PERM Status:Active Type:Cnst,Vol,Exp,NON Lower: .153115E-03 Initial: .153115E-01 upper: .153115E+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: .485236E+00 Initial: .485236E+02 upper: .485236E+04 # 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: .105588E-02 Initial: .105588E+00 upper: .105588E+02 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 4 )[All] !! Middle Semi-Confining Unit # Parameter:kLFA Tag:PERM Status:Active Type:Cnst,Vol,Exp,NON Lower: .413734E+00 Initial: .413734E+02 upper: .413734E+04 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 5 )[All] !! Lower Floridan Aquifer # 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:aLFA 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*( 5 )[All] !! Anisotropy in Penetrated Aquifer # # _________________ Specific Storage ___________________ # Parameter:sICU Tag:STOR Status:Active Type:Cnst,Vol,Exp,NON Lower: .360099E-07 Initial: .360099E-06 upper: .360099E-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: .278746E-06 Initial: .278746E-05 upper: .278746E-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: .151061E-06 Initial: .151061E-05 upper: .151061E-04 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 4 )[All] !! Middle Semi-Confining Unit # Parameter:sLFA Tag:STOR Status:Active Type:Cnst,Vol,Exp,NON Lower: .668044E-06 Initial: .668044E-05 upper: .668044E-04 # Multiplier*(Zones)[SP] - it is required to indicate zone input between parenthesis 1.00*( 5 )[All] !! Lower Floridan Aquifer # 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 .129032E-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: .750 -560.000 1 SITE:UFA1_OBS_wel # Time, Days s, ft Weight .00069 13.32 1.000 W-Default !! 0 HOURS, 0 MINs., 59 SECs AFTER PUMPING Redirect to file:INPUT_UFA1.OBS.txt # UFA2 XYK: 2430.000 -560.000 1 SITE:UFA2_OBS_wel # Time, Days s, ft Weight .13889 .11 1.000 W-Default !! 3 HOURS, 20 MINs., 0 SECs AFTER PUMPING Redirect to file:INPUT_UFA2.OBS.txt # ICU1 XYK: 100.000 -235.000 1 SITE:ICU1_OBS_wel # Time, Days s, ft Weight .08333 .13 1.000 W-Default !! 1 HOURS, 59 MINs., 59 SECs AFTER PUMPING Redirect to file:INPUT_ICU1.OBS.txt # MSC1 XYK: 200.000 -772.000 1 SITE:MSC1_OBS_wel # Time, Days s, ft Weight .01111 .14 1.000 W-Default !! 0 HOURS, 15 MINs., 59 SECs AFTER PUMPING Redirect to file:INPUT_MSC1.OBS.txt # LFA1 XYK: 230.000 -1130.000 1 SITE:LFA1_OBS_wel # Time, Days s, ft Weight .00139 .11 1.000 W-Default !! 0 HOURS, 2 MINs., 0 SECs AFTER PUMPING Redirect to file:INPUT_LFA1.OBS.txt #