MODFLOW-2005 U.S. GEOLOGICAL SURVEY MODULAR FINITE-DIFFERENCE GROUND-WATER FLOW MODEL VERSION 1.8.00 12/18/2009 LIST FILE: imodflow.list UNIT 10 OPENING imodflow.dis FILE TYPE:DIS UNIT 11 STATUS:OLD FORMAT:FORMATTED ACCESS:SEQUENTIAL OPENING imodflow.bas6 FILE TYPE:BAS6 UNIT 12 STATUS:OLD FORMAT:FORMATTED ACCESS:SEQUENTIAL OPENING imodflow.bcf6 FILE TYPE:BCF6 UNIT 13 STATUS:OLD FORMAT:FORMATTED ACCESS:SEQUENTIAL OPENING imodflow.oc FILE TYPE:OC UNIT 14 STATUS:OLD FORMAT:FORMATTED ACCESS:SEQUENTIAL OPENING imodflow.met FILE TYPE:MET UNIT 15 STATUS:OLD FORMAT:FORMATTED ACCESS:SEQUENTIAL OPENING imodflow.riv FILE TYPE:RIV UNIT 16 STATUS:OLD FORMAT:FORMATTED ACCESS:SEQUENTIAL OPENING imodflow.pcg FILE TYPE:PCG UNIT 22 STATUS:OLD FORMAT:FORMATTED ACCESS:SEQUENTIAL # data: BAS -- BASIC PACKAGE, VERSION 7, 5/2/2005 INPUT READ FROM UNIT 12 DISCRETIZATION INPUT DATA READ FROM UNIT 11 # Discretization file 2 LAYERS 100 ROWS 100 COLUMNS 1 STRESS PERIOD(S) IN SIMULATION MODEL TIME UNIT IS DAYS MODEL LENGTH UNIT IS METERS Confining bed flag for each layer: 1 0 DELR = 25.0000 DELC = 25.0000 TOP ELEVATION OF LAYER 1 = 0.00000 MODEL LAYER BOTTOM EL. = 0.00000 FOR LAYER 1 BOT. EL. OF QUASI-3D BED = 0.00000 FOR LAYER 1 MODEL LAYER BOTTOM EL. = 0.00000 FOR LAYER 2 STRESS PERIOD LENGTH TIME STEPS MULTIPLIER FOR DELT SS FLAG ---------------------------------------------------------------------------- 1 0.000000 1 1.000 SS STEADY-STATE SIMULATION # Basic Package file THE FREE FORMAT OPTION HAS BEEN SELECTED OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\BOUNDARY\IMODFLOW_IBOUND1.IDF FILE IDENTIFIED AS IDF BOUNDARY ARRAY = 1 FOR LAYER 2 AQUIFER HEAD WILL BE SET TO -999.99 AT ALL NO-FLOW NODES (IBOUND=0). OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\STARTING_HEADS\IMODFLOW_SHEAD1.IDF FILE IDENTIFIED AS IDF OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\STARTING_HEADS\IMODFLOW_SHEAD2.IDF FILE IDENTIFIED AS IDF OUTPUT CONTROL IS SPECIFIED ONLY AT TIME STEPS FOR WHICH OUTPUT IS DESIRED HEAD PRINT FORMAT CODE IS 0 DRAWDOWN PRINT FORMAT CODE IS 0 HEADS WILL BE SAVED ON UNIT 28 DRAWDOWNS WILL BE SAVED ON UNIT 0 BCF -- BLOCK-CENTERED FLOW PACKAGE, VERSION 7, 5/2/2005 INPUT READ FROM UNIT 13 STEADY-STATE SIMULATION CELL-BY-CELL FLOWS WILL BE SAVED ON UNIT 29 HEAD AT CELLS THAT CONVERT TO DRY= -9999.0 WETTING CAPABILITY IS NOT ACTIVE LAYER LAYER-TYPE CODE INTERBLOCK T -------------------------------------------- 1 0 0 -- HARMONIC 2 0 0 -- HARMONIC COLUMN TO ROW ANISOTROPY = 1.00000 OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\TRANSMISSIVITY\IMODFLOW_KD1.IDF FILE IDENTIFIED AS IDF OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\VERTICAL_RESISTANCE\IMODFLOW_C1.IDF FILE IDENTIFIED AS IDF OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\TRANSMISSIVITY\IMODFLOW_KD2.IDF FILE IDENTIFIED AS IDF RIV -- RIVER PACKAGE, VERSION 7, 5/2/2005 INPUT READ FROM UNIT 16 TNO version v7r1, 28 apr 2009 # River Package file No named parameters MAXIMUM OF 1 ACTIVE RIVER REACHES AT ONE TIME CELL-BY-CELL FLOWS WILL BE SAVED ON UNIT 30 LISTS OF RIVER CELLS WILL NOT BE PRINTED AUXILIARY RIVER VARIABLE: RFCT AUXILIARY RIVER VARIABLE: ISUB RFACT variabel RFCT used as infiltration factor. RSUBSYS variabel ISUB used for sub-system indices. 0 River parameters PCG -- CONJUGATE-GRADIENT SOLUTION PACKAGE, VERSION 7, 5/2/2005 # Preconditioned Conjugate-Gradient Package MAXIMUM OF 150 CALLS OF SOLUTION ROUTINE MAXIMUM OF 20 INTERNAL ITERATIONS PER CALL TO SOLUTION ROUTINE MATRIX PRECONDITIONING TYPE : 1 SOLUTION BY THE CONJUGATE-GRADIENT METHOD ------------------------------------------- MAXIMUM NUMBER OF CALLS TO PCG ROUTINE = 150 MAXIMUM ITERATIONS PER CALL TO PCG = 20 MATRIX PRECONDITIONING TYPE = 1 RELAXATION FACTOR (ONLY USED WITH PRECOND. TYPE 1) = 0.98000E+00 PARAMETER OF POLYNOMIAL PRECOND. = 2 (2) OR IS CALCULATED : 1 HEAD CHANGE CRITERION FOR CLOSURE = 0.10000E-03 RESIDUAL CHANGE CRITERION FOR CLOSURE = 0.10000E+02 PCG HEAD AND RESIDUAL CHANGE PRINTOUT INTERVAL = 1 PRINTING FROM SOLVER IS LIMITED(1) OR SUPPRESSED (>1) = 0 STEADY-STATE DAMPING PARAMETER = 0.10000E+01 TRANSIENT DAMPING PARAMETER = 0.10000E+01 1 STRESS PERIOD NO. 1, LENGTH = 0.000000 ----------------------------------------------- NUMBER OF TIME STEPS = 1 MULTIPLIER FOR DELT = 1.000 INITIAL TIME STEP SIZE = 0.000000 DATA IDENTIFIED AS GCD OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\RIVERS\IMODFLOW_RIV_COND1.IDF FILE IDENTIFIED AS IDF OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\RIVERS\IMODFLOW_RIV_STAGE1.IDF FILE IDENTIFIED AS IDF OPENING FILE ON UNIT 99: P:\1206764-IMOD\IMOD30\TESTBANK\RUNS\TEST_ALL\RIVERS\IMODFLOW_RIV_RBOT1.IDF FILE IDENTIFIED AS IDF REACH NO. LAYER ROW COL STAGE CONDUCTANCE BOTTOM EL. 4 = 1.00000 FOR LAYER 1 100 RIVER REACHES SOLVING FOR HEAD SOLVING FOR HEAD SOLVING FOR HEAD 3 CALLS TO PCG ROUTINE FOR TIME STEP 1 IN STRESS PERIOD 1 41 TOTAL ITERATIONS MAXIMUM HEAD CHANGE FOR EACH ITERATION (1 INDICATES THE FIRST INNER ITERATION): HEAD CHANGE HEAD CHANGE HEAD CHANGE HEAD CHANGE HEAD CHANGE LAYER,ROW,COL LAYER,ROW,COL LAYER,ROW,COL LAYER,ROW,COL LAYER,ROW,COL --------------------------------------------------------------------------- 1 3.728 0 3.289 0 1.283 0 1.056 0 0.7895 ( 1, 5, 99) ( 2, 1, 99) ( 2, 3, 96) ( 2, 2, 87) ( 2, 1, 93) 0 0.7877 0 0.6896 0 0.5541 0 0.8652 0 0.7837 ( 2, 3, 88) ( 2, 3, 81) ( 1, 5, 78) ( 2, 2, 79) ( 2,100, 62) 0 0.6362 0 0.6390 0 0.2760 0 0.1247 0 0.9598E-01 ( 2,100, 56) ( 2, 1, 70) ( 2, 4, 50) ( 2, 1, 55) ( 2, 1, 62) 0 0.4582E-01 0 0.4232E-01 0 0.3952E-01 0 0.3791E-01 0 0.2395E-01 ( 2, 3, 57) ( 2, 1, 61) ( 2,100, 30) ( 1,100, 24) ( 2, 99, 22) 1 -0.1185E-01 0 0.1056E-01 0 0.8742E-02 0 -0.4913E-02 0 0.3680E-02 ( 2, 2, 55) ( 2, 2, 40) ( 2, 1, 46) ( 2, 1, 63) ( 2, 1, 32) 0 -0.3359E-02 0 0.2751E-02 0 -0.1922E-02 0 -0.1095E-02 0 0.7130E-03 ( 2, 4, 66) ( 2, 1, 30) ( 1, 11, 85) ( 2, 2, 91) ( 2, 1, 71) 0 0.4383E-03 0 0.2868E-03 0 0.3123E-03 0 0.2920E-03 0 0.3486E-03 ( 2, 3, 40) ( 1, 5, 43) ( 2, 2, 47) ( 1, 4, 50) ( 2, 2, 54) 0 0.3188E-03 0 0.2623E-03 0 0.2492E-03 0 0.1646E-03 0 0.1135E-03 ( 2, 2, 57) ( 1, 7, 55) ( 2, 1, 45) ( 2, 1, 21) ( 1, 5, 75) 1 0.4091E-04 ( 2, 4, 81) MAXIMUM RESIDUAL FOR EACH ITERATION (1 INDICATES THE FIRST INNER ITERATION): RESIDUAL RESIDUAL RESIDUAL RESIDUAL RESIDUAL LAYER,ROW,COL LAYER,ROW,COL LAYER,ROW,COL LAYER,ROW,COL LAYER,ROW,COL --------------------------------------------------------------------------- 1 -254.2 0 81.22 0 101.9 0 -36.98 0 68.96 ( 2, 1, 99) ( 2, 1, 98) ( 2, 1, 99) ( 2,100,100) ( 2, 2, 98) 0 32.58 0 -44.44 0 -38.68 0 -36.97 0 -48.53 ( 2, 99, 99) ( 2, 2, 98) ( 2, 2, 98) ( 2, 4, 97) ( 2, 3, 97) 0 27.48 0 23.25 0 17.64 0 8.523 0 3.414 ( 2, 99, 99) ( 2, 3, 97) ( 2,100, 99) ( 2,100, 99) ( 2, 4, 96) 0 -3.818 0 -2.487 0 -1.570 0 -2.323 0 -1.234 ( 2,100, 99) ( 2,100, 99) ( 2, 6, 94) ( 2, 99,100) ( 2, 99,100) 1 -0.6545 0 0.3713 0 -0.5793 0 0.3053 0 -0.1802 ( 2, 99,100) ( 2, 8, 93) ( 2,100, 99) ( 2, 99,100) ( 2, 3, 97) 0 0.3048 0 0.1935 0 0.1150 0 0.9613E-01 0 -0.4991E-01 ( 2,100, 99) ( 2, 3, 97) ( 2, 99, 99) ( 2, 99, 99) ( 2,100, 99) 0 -0.3001E-01 0 -0.2932E-01 0 -0.1421E-01 0 0.2132E-01 0 0.1775E-01 ( 2, 99, 99) ( 2, 99, 99) ( 2, 2, 98) ( 2, 99, 99) ( 2, 99, 99) 0 -0.1011E-01 0 -0.1586E-01 0 0.8565E-02 0 -0.9195E-02 0 0.4678E-02 ( 2,100,100) ( 2, 99, 99) ( 2,100, 99) ( 2, 1, 99) ( 2, 99, 99) 1 0.2487E-02 ( 2,100,100) OUTPUT CONTROL FOR STRESS PERIOD 1 TIME STEP 1 SAVE HEAD FOR ALL LAYERS SAVE BUDGET SAVE BUDGET INFO. Sign of cbc-flux FLOW RIGHT FACE is swapped! INFO. Sign of cbc-flux FLOW FRONT FACE is swapped! INFO. Sign of cbc-flux FLOW LOWER FACE is swapped! HEAD WILL BE SAVED ON UNIT 28 AT END OF TIME STEP 1, STRESS PERIOD 1 1 VOLUMETRIC BUDGET FOR ENTIRE MODEL AT END OF TIME STEP 1 IN STRESS PERIOD 1 ------------------------------------------------------------------------------ CUMULATIVE VOLUMES L**3 RATES FOR THIS TIME STEP L**3/T ------------------ ------------------------ IN: IN: --- --- STORAGE = 0.0000 STORAGE = 0.0000 CONSTANT HEAD = 0.0000 CONSTANT HEAD = 1227.5625 RIVER LEAKAGE = 0.0000 RIVER LEAKAGE = 0.0000 TOTAL IN = 0.0000 TOTAL IN = 1227.5625 OUT: OUT: ---- ---- STORAGE = 0.0000 STORAGE = 0.0000 CONSTANT HEAD = 0.0000 CONSTANT HEAD = 0.0000 RIVER LEAKAGE = 0.0000 RIVER LEAKAGE = 1227.1089 TOTAL OUT = 0.0000 TOTAL OUT = 1227.1089 IN - OUT = 0.0000 IN - OUT = 0.4536 PERCENT DISCREPANCY = 0.00 PERCENT DISCREPANCY = 0.04 TIME SUMMARY AT END OF TIME STEP 1 IN STRESS PERIOD 1 SECONDS MINUTES HOURS DAYS YEARS ----------------------------------------------------------- TIME STEP LENGTH 0.0000 0.0000 0.0000 0.0000 0.0000 STRESS PERIOD TIME 0.0000 0.0000 0.0000 0.0000 0.0000 TOTAL TIME 0.0000 0.0000 0.0000 0.0000 0.0000 1 Run end date and time (yyyy/mm/dd hh:mm:ss): 2014/09/25 16:09:07 Elapsed run time: 3.763 Seconds