1 @@@@@ @@ @ @ @@@ @ @ @@ @@@ @ @ @ @ @ @ @@ @@@@@ @ @@ @ @ @@@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @@ @@ @ @ @ @ @ @ @ @ @ @@ @ @ @ @ @ @@ @ @ @ @ @ @ @ @@ @@@@ @ @@ @ @ @ @ @ @ @ @@@@ @ @ @ @ @ @ @ @@@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @@ @ @ @ @ @ @@ @ @@ @@ @@@ @ @ @@@@ @@@ @ @ @ @@ @@@@ @ @ @ @ @@ @ @ @ @ @@ @ @ TOUGH2 IS A PROGRAM FOR MULTIPHASE MULTICOMPONENT FLOW IN PERMEABLE MEDIA, INCLUDING HEAT FLOW. IT IS A MEMBER OF THE MULKOM FAMILY OF CODES, DEVELOPED BY KARSTEN PRUESS AT LAWRENCE BERKELEY NATIONAL LABORATORY. COPYRIGHT (C) 1999, THE REGENTS OF THE UNIVERSITY OF CALIFORNIA. ******************************************************************************** ********* THIS IS TOUGH2 RUN NO. 0 CALLED BY iTOUGH2 RUN NO. 1 ********* ************* iTOUGH2 V6.5 (OCTOBER, 2010) FOR PC, S. FINSTERLE ************** ******************** EQUATION OF STATE MODULE:777 ******************* ******************** TODAY'S DATE: Provide date in ******************* ******************************************************************************** PROBLEM TITLE : *SAM7C3* ... one-d radial injection into saturated region. =================================================================================================================================== PARAMETERS FOR FLEXIBLE DIMENSIONING OF MAJOR ARRAYS (SEE FILE maxsize.inc) ARE SUMMARIZED AT THE END. =================================================================================================================================== *********************************************************************************************************************************** * * * M A T R I X S O L V E R * * * *********************************************************************************************************************************** THE SOLVER IS DETERMINED FROM MOP(21) THE SOLUTION METHOD INDICATOR MATSLV = 0 - RESET INTERNALLY TO THE DEFAULT, MATSLV = 3 THE SOLUTION METHOD INDICATOR MATSLV = 0 MATSLV = 1: MATB - (PROPRIETARY - MAY NOT BE AVAILABLE) MATSLV = 2: DSLUBC - BI-CONJUGATE GRADIENT SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 3: DSLUCS - LANCZOS-TYPE CONJUGATE GRADIENT SQUARED SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 4: DSLUGM - GENERALIZED MINIMUM RESIDUAL CONJUGATE GRADIENT SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 5: DLUSTB - STABILIZED BI-CONJUGATE GRADIENT SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 6: LUBAND - DIRECT SOLVER USING LU DECOMPOSITION RITMAX: MAXIMUM # OF CG ITERATIONS AS FRACTION OF THE TOTAL NUMBER OF EQUATIONS = 1.00000E-01 (0.0 < RITMAX <= 1.0, DEFAULT = 0.1) CLOSUR: CONVERGENCE CRITERION FOR THE CG ITERATIONS = 1.00000E-12 (1.0E-12 <= CLOSUR <= 1.0E-6, DEFAULT = 1.0E-6) NMAXIT: MAXIMUM # OF CG ITERATIONS - NOT TO EXCEED THE TOTAL NUMBER OF EQUATIONS NELA*NEQ = 20 (20 < NMAXIT <= NREDM) THE MATRIX Z-PREPROCESSING SYSTEM IS ZPROCS = Z1 ZPROCS = Z0: NO Z-PREPROCESSING; DEFAULT FOR NEQ = 1 AND FOR MATSLV = 6 ZPROCS = Z1: REPLACEMENT OF ZEROS ON THE MAIN-DIAGONAL BY A SMALL NUMBER; DEFAULT FOR NEQ > 1 AND FOR 2 < MATSLV < 6 ZPROCS = Z2: LINEAR COMBINATION Of EQUATIONS IN EACH ELEMENT TO PRODUCE NON-ZERO MAIN DIAGONAL ENTRIES ZPROCS = Z3: NORMALIZATION OF EQUATIONS, FOLLOWED BY Z2 ZPROCS = Z4: SAME AS IN OPROCS = O4 THE MATRIX O-PREPROCESSING SYSTEM IS OPROCS = O0 OPROCS = O0: NO O-PREPROCESSING; DEFAULT (AND ONLY OPTION FOR NEQ=1) OPROCS = O1: ELIMINATION OF LOWER HALF OF THE MAIN-DIAGONAL SUBMATRIX WITH CENTER PIVOTING OPROCS = O2: O1+ELIMINATION OF UPPER HALF OF THE MAIN-DIAGONAL SUBMATRIX WITH CENTER PIVOTING OPROCS = O3: O2+NORMALIZATION - RESULTS IN UNIT MAIN-DIAGONAL SUBMATRICES OPROCS = O4: PRE-PROCESSING WHICH RESULTS IN UNIT MAIN-DIAGONAL SUBMATRICES WITHOUT CENTER PIVOTING VALUE OF INCREMENT FACTOR FOR NUMERICAL DERIVATIVES: 0.00000E+00 TIME STEPPING PARAMETERS ------------------------ MINIMUM TIME STEP VALUE....................(DTMIN) = 0.1000E-09 MAX. ALLOWABLE SAT. CHANGE/TIME-STEP......(DSTNOM) = 0.2000E+00 MAXIMUM TIME STEP MULTIPLIER............. (TMULFC) = 0.5000E+01 END OF TOUGH2 INPUT JOB --- ELAPSED TIME = 0.000 SECONDS *********************************************************************************************************************************** * ARRAY *MOP* ALLOWS TO GENERATE MORE PRINTOUT IN VARIOUS SUBROUTINES, AND TO MAKE SOME CALCULATIONAL CHOICES. * *********************************************************************************************************************************** MOP(1) = 0 *** ALLOWS TO GENERATE A SHORT PRINTOUT FOR EACH NEWTON-RAPHSON ITERATION = 0, 1, OR 2: GENERATE 0, 1, OR 2 LINES OF PRINTOUT MORE PRINTOUT IS GENERATED FOR MOP(I) > 0 IN THE FOLLOWING SUBROUTINES (THE LARGER MOP IS, THE MORE WILL BE PRINTED). MOP(2) = 0 *** CYCIT MOP(3) = 0 *** MULTI MOP(4) = 0 *** QU MOP(5) = 0 *** EOS MOP(6) = 0 *** LINEQ MOP(7) = 0 *** IF UNEQUAL ZERO, WILL GENERATE A PRINTOUT OF INPUT DATA CALCULATIONAL CHOICES OFFERED BY MOP ARE AS FOLLOWS: MOP(9) = 0 *** CHOOSES FLUID COMPOSITION ON WITHDRAWAL (PRODUCTION). = 0: ACCORDING TO RELATIVE MOBILITIES. = 1: ACCORDING TO COMPOSITION IN PRODUCING ELEMENT. MOP(10) = 0 *** CHOOSES INTERPOLATION FORMULA FOR DEPENDENCE OF THERMAL CONDUCTIVITY ON LIQUID SATURATION (SL). = 0: K = KDRY + SQRT(SL)*(KWET-KDRY) = 1: K = KDRY + SL*(KWET-KDRY) = 2: K = C0 + C1*T + C2*SL + C3*PHI (RAUTMAN) MOP(11) = 0 *** CHOOSES EVALUATION OF MOBILITY AND ABSOLUTE PERMEABILITY AT INTERFACES. = 0: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS UPSTREAM WEIGHTED. = 1: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS UPSTREAM WEIGHTED. = 2: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS HARMONIC WEIGHTED. = 3: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS HARMONIC WEIGHTED. = 4: MOBILITY * PERMEABILITY PRODUCT IS HARMONIC WEIGHTED. MOP(12) = 0 *** CHOOSES PROCEDURE FOR INTERPOLATING GENERATION RATES FROM A TIME TABLE. = 0: TRIPLE LINEAR INTERPOLATION. = 1: "STEP FUNCTION" OPTION. = 2: RIGOROUS STEP RATE OPTION. MOP2(3) > 0: TIME STEPS ADJUSTED TO MATCH GENERATION TIMES. GENER * : READS TIMES AND RATES [AND ENTHALPIES] IN FREE FORMAT. GENER T : READS TIMES AND RATES [AND ENTHALPIES] IN FREE TABULAR FORM (T_i,Q_i,[E_i], i=1,LTAB). MOP(13) = 0 *** DEFINES CONTENT OF INCON AND SAVE FILE. = 0: STANDARD CONTENT. = 1: WRITES USER-SPECIFIED INITIAL CONDITIONS TO FILE SAVE. = 2: READS PARAMETERS OF HYSTERESIS MODEL FROM FILE INCON. MOP(15) = 0 *** ALLOWS TO SELECT A SEMI-ANALYTICAL HEAT EXCHANGE CALCULATION WITH CONFINING BEDS. = 0: NO SEMI-ANALYTICAL HEAT EXCHANGE > 0: SEMI-ANALYTICAL HEAT EXCHANGE ENGAGED (WHEN A SPECIAL SUBROUTINE MODULE *QLOSS* IS PRESENT) MOP(16) = 0 *** PERMITS TO CHOOSE TIME STEP SELECTION OPTION = 0: AUTOMATIC TIME STEPPING BASED ON MAXIMUM CHANGE IN SATURATION (RECOMMENDED). = 1: AUTOMATIC TIME STEPPING BASED ON NUMBER OF ITERATIONS NEEDED FOR CONVERGENCE > 1: INCREASE TIME STEP BY AT LEAST A FACTOR OF TWO, IF CONVERGENCE OCCURS IN .LE. MOP(16) ITERATIONS. MOP(17) = 0 *** HANDLES SCALING OPTIONS. = 0: NO SCALING. = 7: SCALING. MOP(18) = 0 *** ALLOWS TO SELECT HANDLING OF INTERFACE DENSITY. = 0: PERFORM UPSTREAM WEIGHTING FOR INTERFACE DENSITY. > 0: COMPUTE INTERFACE DENSITY AS AVERAGE OF THE TWO GRID BLOCK DENSITIES. HOWEVER, WHEN ONE OF THE TWO PHASE SATURATIONS IS ZERO, DO UPSTREAM WEIGHTING. MOP(21) = 3 *** PERMITS TO SELECT LINEAR EQUATION SOLVER FROM PACKAGE t2cg22 = 0: DEFAULTS TO MOP(21) = 3 = 1: MATB (MAY NOT BE AVAIALBALE) = 2: DSLUBC: BI-CONJUGATE GRADIENT SOLVER; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 3: DSLUCS: BI-CONJUGATE GRADIENT SOLVER - LANCZOS TYPE; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 4: DSLUGM: GENERALIZED MINIMUM RESIDUAL CONJUGATE GRADIENTS; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 5: DLUSTB: STABILIZED BI-CONJUGATE GRADIENT SOLVER; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 6: LUBAND: DIRECT SOLVER USING LU DECOMPOSITION = 7: TOGGLES BETWEEN 3, AND 5 = 8: TOGGLES BETWEEN 2, 3, 4, AND 5 = 9: TOGGLES BETWEEN 5, AND 1 MOP(22) = 0 *** SUBROUTINE USERBC = 0: DON'T CALL USERBC = 1: CALL USERBC = 2: CALL USERBC AND EOS AFTER CONVERGENCE MOP(24) = 0 *** MAKES CHOICE OF 5- OR 9-POINT MESH = 0: 5-POINT = 1: 9-POINT; USE ONLY WITH ROTATION ANGLE DEG = 0 OR 90 = 2: 9-POINT; SYMMETRY AT I=1 *********************************************************************************************************************************** ADDITIONAL OPTIONS CAN BE PROVIDED ON A LINE IN FORMAT 80I1, FOLLOWING BLOCK IDENTIFIER "MOMOP" MOP2(1) = 0 *** MINIMUM NUMBER OF NEWTON-RAPHSON ITERATIONS 0 or 1: ALLOW CONVERGENCE IN ONE ITERATION (DEFAULT) 2: PERFORM AT LEAST 2 ITERATIONS; PRIMARY VARIABLES ARE ALWAYS UPDATED 3: ALLOW CONVERGENCE IN ONE ITERATION FOR NEGATIVE SIMULATION TIMES BUT REQUEST TWO ITERATIONS FOR POSITIVE TIMES 4: ALLOW CONVERGENCE IN ONE ITERATION FOR POSITIVE SIMULATION TIMES BUT REQUEST TWO ITERATIONS FOR NEGATIVE TIMES MOP2(2) = 0 *** LENGTH OF ELEMENT NAMES = 0: 5-CHARACTER ELEMENT NAMES (DEFAULT) >= 5: 5 <= ELEMENT NAME LENGTH <= 9 FORMAT IN TOUGH2 INPUT BLOCK LENGTH ELEME CONNE INCON GENER 5: (A3,I2,I5,I5,A2,A3,6E10.4) (2(A3,I2),I5,2I5,I5,4E10.4) (A3,I2,I5,I5,E15.8,4E12.4) (A3,I2,A3,I2,I5,2I5,I5,5X,A4,A1,3E10.4) 6: (A3,I3,I5,I4,A2,A3,6E10.4) (2(A3,I3),I5,2I4,I5,4E10.4) (A3,I3,I5,I4,E15.8,4E12.4) (A3,I3,A3,I2,I6,2I4,I5,5X,A4,A1,3E10.4) 7: (A3,I4,I4,I4,A2,A3,6E10.4) (2(A3,I4),I5,2I3,I5,4E10.4) (A3,I4,I4,I4,E15.8,4E12.4) (A3,I4,A3,I2,I5,2I4,I5,5X,A4,A1,3E10.4) 8: (A3,I5,I4,I3,A2,A3,6E10.4) (2(A3,I5),I3,2I3,I5,4E10.4) (A3,I5,I4,I3,E15.8,4E12.4) (A3,I5,A3,I2,I4,2I4,I5,5X,A4,A1,3E10.4) 9: (A3,I6,I3,I3,A2,A3,6E10.4) (2(A3,I6),I3,2I2,I5,4E10.4) (A3,I6,I3,I3,E15.8,4E12.4) (A3,I6,A3,I2,I5,2I3,I5,5X,A4,A1,3E10.4) MOP2(3) = 0 *** HONORING GENERATION TIMES = 0: GENERATION TIMES IGNORED > 0: TIME STEPS ADJUSTED TO MATCH GENERATION TIMES MOP2(4) = 0 *** REDUCE VAPOR PRESSURE = 0: NO VAPOR PRESSURE REDUCTION AT LOW SATURATION > 0: REDUCES VAPOR PRESSURE FOR SL < 0.02 TO PREVENT LIQUID DISAPPEARANCE (ONLY CERTAIN EOS MODULES) MOP2(5) = 0 *** ACTIVE FRACTURE MODEL = 0: ACTIVE FRACTURE MODEL APPLIED TO LIQUID PHASE ONLY > 0: ACTIVE FRACTURE MODEL APPLIED TO ALL PHASES MOP2(6) = 0 *** LEVERETT SCALING > 0: RESCALE CAPILLARY PRESSURE: PC=PC_REF*SQRT(K_REF/K) IF ELEMENT-BY-ELEMENT PERMEABILITIES ARE DEFINED MOP2(7) = 0 *** ZERO NODAL DISTANCE = 0: TAKE ABSOLUTE PERMEABILITY FROM OTHER GRID BLOCK > 0: TAKE ABSOLUTE AND RELATIVE PERMEABILITY FROM OTHER GRID BLOCK. MOP2(8) = 0 *** VERSION OF SINK/SOURCE SUBROUTINE = 0: USE ORIGINAL VERSION OF QU > 0: USE QU VERSION 2 MOP2(9) = 0 *** TIME STEP AFTER FORCED TIME-STEP REDUCTION TO HONOR PRINTOUT TIME = 0: CONTINUE WITH TIME STEP USED BEFORE FORCED TIME-STEP REDUCTION > 0: CONTINUE WITH TIME STEP IMPOSED BY FORCED TIME-STEP REDUCTION MOP2(10)= 0 *** PRINT SAVE FILE AFTER EACH OUTPUT *********************************************************************************************************************************** DOMAIN NO. 1 MATERIAL NAME -- cly01 DOMAIN NO. 2 MATERIAL NAME -- snd01 DOMAIN NO. 3 MATERIAL NAME -- BC01 HAVE READ UNKNOWN BLOCK LABEL "----*" --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA WRITE FILE *MESH* FROM INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA WRITE FILE *GENER* FROM INPUT DATA WRITE FILE *INCON* FROM INPUT DATA HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA *********************************************************************************************************************************** * * * M A T R I X S O L V E R * * * *********************************************************************************************************************************** THE SOLVER IS DETERMINED FROM THE SOLVR DATA BLOCK - MOP(21) IS OVERRIDEN! THE SOLUTION METHOD INDICATOR MATSLV = 5 MATSLV = 1: MATB - (PROPRIETARY - MAY NOT BE AVAILABLE) MATSLV = 2: DSLUBC - BI-CONJUGATE GRADIENT SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 3: DSLUCS - LANCZOS-TYPE CONJUGATE GRADIENT SQUARED SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 4: DSLUGM - GENERALIZED MINIMUM RESIDUAL CONJUGATE GRADIENT SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 5: DLUSTB - STABILIZED BI-CONJUGATE GRADIENT SOLVER INCOMPLETE LU FACTORIZATION PRECONDITIONING MATSLV = 6: LUBAND - DIRECT SOLVER USING LU DECOMPOSITION RITMAX: MAXIMUM # OF CG ITERATIONS AS FRACTION OF THE TOTAL NUMBER OF EQUATIONS = 8.00000E-01 (0.0 < RITMAX <= 1.0, DEFAULT = 0.1) CLOSUR: CONVERGENCE CRITERION FOR THE CG ITERATIONS = 1.00000E-07 (1.0E-12 <= CLOSUR <= 1.0E-6, DEFAULT = 1.0E-6) NMAXIT: MAXIMUM # OF CG ITERATIONS - NOT TO EXCEED THE TOTAL NUMBER OF EQUATIONS NELA*NEQ = 20 (20 < NMAXIT <= NREDM) THE MATRIX Z-PREPROCESSING SYSTEM IS ZPROCS = Z1 ZPROCS = Z0: NO Z-PREPROCESSING; DEFAULT FOR NEQ = 1 AND FOR MATSLV = 6 ZPROCS = Z1: REPLACEMENT OF ZEROS ON THE MAIN-DIAGONAL BY A SMALL NUMBER; DEFAULT FOR NEQ > 1 AND FOR 2 < MATSLV < 6 ZPROCS = Z2: LINEAR COMBINATION Of EQUATIONS IN EACH ELEMENT TO PRODUCE NON-ZERO MAIN DIAGONAL ENTRIES ZPROCS = Z3: NORMALIZATION OF EQUATIONS, FOLLOWED BY Z2 ZPROCS = Z4: SAME AS IN OPROCS = O4 THE MATRIX O-PREPROCESSING SYSTEM IS OPROCS = O0 OPROCS = O0: NO O-PREPROCESSING; DEFAULT (AND ONLY OPTION FOR NEQ=1) OPROCS = O1: ELIMINATION OF LOWER HALF OF THE MAIN-DIAGONAL SUBMATRIX WITH CENTER PIVOTING OPROCS = O2: O1+ELIMINATION OF UPPER HALF OF THE MAIN-DIAGONAL SUBMATRIX WITH CENTER PIVOTING OPROCS = O3: O2+NORMALIZATION - RESULTS IN UNIT MAIN-DIAGONAL SUBMATRICES OPROCS = O4: PRE-PROCESSING WHICH RESULTS IN UNIT MAIN-DIAGONAL SUBMATRICES WITHOUT CENTER PIVOTING VALUE OF INCREMENT FACTOR FOR NUMERICAL DERIVATIVES: 0.00000E+00 END OF TOUGH2 INPUT JOB --- ELAPSED TIME = 0.000 SECONDS *********************************************************************************************************************************** * ARRAY *MOP* ALLOWS TO GENERATE MORE PRINTOUT IN VARIOUS SUBROUTINES, AND TO MAKE SOME CALCULATIONAL CHOICES. * *********************************************************************************************************************************** MOP(1) = 1 *** ALLOWS TO GENERATE A SHORT PRINTOUT FOR EACH NEWTON-RAPHSON ITERATION = 0, 1, OR 2: GENERATE 0, 1, OR 2 LINES OF PRINTOUT MORE PRINTOUT IS GENERATED FOR MOP(I) > 0 IN THE FOLLOWING SUBROUTINES (THE LARGER MOP IS, THE MORE WILL BE PRINTED). MOP(2) = 0 *** CYCIT MOP(3) = 0 *** MULTI MOP(4) = 0 *** QU MOP(5) = 3 *** EOS MOP(6) = 0 *** LINEQ MOP(7) = 0 *** IF UNEQUAL ZERO, WILL GENERATE A PRINTOUT OF INPUT DATA CALCULATIONAL CHOICES OFFERED BY MOP ARE AS FOLLOWS: MOP(9) = 0 *** CHOOSES FLUID COMPOSITION ON WITHDRAWAL (PRODUCTION). = 0: ACCORDING TO RELATIVE MOBILITIES. = 1: ACCORDING TO COMPOSITION IN PRODUCING ELEMENT. MOP(10) = 0 *** CHOOSES INTERPOLATION FORMULA FOR DEPENDENCE OF THERMAL CONDUCTIVITY ON LIQUID SATURATION (SL). = 0: K = KDRY + SQRT(SL)*(KWET-KDRY) = 1: K = KDRY + SL*(KWET-KDRY) = 2: K = C0 + C1*T + C2*SL + C3*PHI (RAUTMAN) MOP(11) = 2 *** CHOOSES EVALUATION OF MOBILITY AND ABSOLUTE PERMEABILITY AT INTERFACES. = 0: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS UPSTREAM WEIGHTED. = 1: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS UPSTREAM WEIGHTED. = 2: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS HARMONIC WEIGHTED. = 3: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS HARMONIC WEIGHTED. = 4: MOBILITY * PERMEABILITY PRODUCT IS HARMONIC WEIGHTED. MOP(12) = 0 *** CHOOSES PROCEDURE FOR INTERPOLATING GENERATION RATES FROM A TIME TABLE. = 0: TRIPLE LINEAR INTERPOLATION. = 1: "STEP FUNCTION" OPTION. = 2: RIGOROUS STEP RATE OPTION. MOP2(3) > 0: TIME STEPS ADJUSTED TO MATCH GENERATION TIMES. GENER * : READS TIMES AND RATES [AND ENTHALPIES] IN FREE FORMAT. GENER T : READS TIMES AND RATES [AND ENTHALPIES] IN FREE TABULAR FORM (T_i,Q_i,[E_i], i=1,LTAB). MOP(13) = 0 *** DEFINES CONTENT OF INCON AND SAVE FILE. = 0: STANDARD CONTENT. = 1: WRITES USER-SPECIFIED INITIAL CONDITIONS TO FILE SAVE. = 2: READS PARAMETERS OF HYSTERESIS MODEL FROM FILE INCON. MOP(15) = 0 *** ALLOWS TO SELECT A SEMI-ANALYTICAL HEAT EXCHANGE CALCULATION WITH CONFINING BEDS. = 0: NO SEMI-ANALYTICAL HEAT EXCHANGE > 0: SEMI-ANALYTICAL HEAT EXCHANGE ENGAGED (WHEN A SPECIAL SUBROUTINE MODULE *QLOSS* IS PRESENT) MOP(16) = 4 *** PERMITS TO CHOOSE TIME STEP SELECTION OPTION = 0: AUTOMATIC TIME STEPPING BASED ON MAXIMUM CHANGE IN SATURATION (RECOMMENDED). = 1: AUTOMATIC TIME STEPPING BASED ON NUMBER OF ITERATIONS NEEDED FOR CONVERGENCE > 1: INCREASE TIME STEP BY AT LEAST A FACTOR OF TWO, IF CONVERGENCE OCCURS IN .LE. MOP(16) ITERATIONS. MOP(17) = 0 *** HANDLES SCALING OPTIONS. = 0: NO SCALING. = 7: SCALING. MOP(18) = 0 *** ALLOWS TO SELECT HANDLING OF INTERFACE DENSITY. = 0: PERFORM UPSTREAM WEIGHTING FOR INTERFACE DENSITY. > 0: COMPUTE INTERFACE DENSITY AS AVERAGE OF THE TWO GRID BLOCK DENSITIES. HOWEVER, WHEN ONE OF THE TWO PHASE SATURATIONS IS ZERO, DO UPSTREAM WEIGHTING. MOP(21) = 3 *** PERMITS TO SELECT LINEAR EQUATION SOLVER FROM PACKAGE t2cg22 = 0: DEFAULTS TO MOP(21) = 3 = 1: MATB (MAY NOT BE AVAIALBALE) = 2: DSLUBC: BI-CONJUGATE GRADIENT SOLVER; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 3: DSLUCS: BI-CONJUGATE GRADIENT SOLVER - LANCZOS TYPE; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 4: DSLUGM: GENERALIZED MINIMUM RESIDUAL CONJUGATE GRADIENTS; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 5: DLUSTB: STABILIZED BI-CONJUGATE GRADIENT SOLVER; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 6: LUBAND: DIRECT SOLVER USING LU DECOMPOSITION = 7: TOGGLES BETWEEN 3, AND 5 = 8: TOGGLES BETWEEN 2, 3, 4, AND 5 = 9: TOGGLES BETWEEN 5, AND 1 MOP(22) = 0 *** SUBROUTINE USERBC = 0: DON'T CALL USERBC = 1: CALL USERBC = 2: CALL USERBC AND EOS AFTER CONVERGENCE MOP(24) = 0 *** MAKES CHOICE OF 5- OR 9-POINT MESH = 0: 5-POINT = 1: 9-POINT; USE ONLY WITH ROTATION ANGLE DEG = 0 OR 90 = 2: 9-POINT; SYMMETRY AT I=1 *********************************************************************************************************************************** ADDITIONAL OPTIONS CAN BE PROVIDED ON A LINE IN FORMAT 80I1, FOLLOWING BLOCK IDENTIFIER "MOMOP" MOP2(1) = 0 *** MINIMUM NUMBER OF NEWTON-RAPHSON ITERATIONS 0 or 1: ALLOW CONVERGENCE IN ONE ITERATION (DEFAULT) 2: PERFORM AT LEAST 2 ITERATIONS; PRIMARY VARIABLES ARE ALWAYS UPDATED 3: ALLOW CONVERGENCE IN ONE ITERATION FOR NEGATIVE SIMULATION TIMES BUT REQUEST TWO ITERATIONS FOR POSITIVE TIMES 4: ALLOW CONVERGENCE IN ONE ITERATION FOR POSITIVE SIMULATION TIMES BUT REQUEST TWO ITERATIONS FOR NEGATIVE TIMES MOP2(2) = 0 *** LENGTH OF ELEMENT NAMES = 0: 5-CHARACTER ELEMENT NAMES (DEFAULT) >= 5: 5 <= ELEMENT NAME LENGTH <= 9 FORMAT IN TOUGH2 INPUT BLOCK LENGTH ELEME CONNE INCON GENER 5: (A3,I2,I5,I5,A2,A3,6E10.4) (2(A3,I2),I5,2I5,I5,4E10.4) (A3,I2,I5,I5,E15.8,4E12.4) (A3,I2,A3,I2,I5,2I5,I5,5X,A4,A1,3E10.4) 6: (A3,I3,I5,I4,A2,A3,6E10.4) (2(A3,I3),I5,2I4,I5,4E10.4) (A3,I3,I5,I4,E15.8,4E12.4) (A3,I3,A3,I2,I6,2I4,I5,5X,A4,A1,3E10.4) 7: (A3,I4,I4,I4,A2,A3,6E10.4) (2(A3,I4),I5,2I3,I5,4E10.4) (A3,I4,I4,I4,E15.8,4E12.4) (A3,I4,A3,I2,I5,2I4,I5,5X,A4,A1,3E10.4) 8: (A3,I5,I4,I3,A2,A3,6E10.4) (2(A3,I5),I3,2I3,I5,4E10.4) (A3,I5,I4,I3,E15.8,4E12.4) (A3,I5,A3,I2,I4,2I4,I5,5X,A4,A1,3E10.4) 9: (A3,I6,I3,I3,A2,A3,6E10.4) (2(A3,I6),I3,2I2,I5,4E10.4) (A3,I6,I3,I3,E15.8,4E12.4) (A3,I6,A3,I2,I5,2I3,I5,5X,A4,A1,3E10.4) MOP2(3) = 0 *** HONORING GENERATION TIMES = 0: GENERATION TIMES IGNORED > 0: TIME STEPS ADJUSTED TO MATCH GENERATION TIMES MOP2(4) = 0 *** REDUCE VAPOR PRESSURE = 0: NO VAPOR PRESSURE REDUCTION AT LOW SATURATION > 0: REDUCES VAPOR PRESSURE FOR SL < 0.02 TO PREVENT LIQUID DISAPPEARANCE (ONLY CERTAIN EOS MODULES) MOP2(5) = 0 *** ACTIVE FRACTURE MODEL = 0: ACTIVE FRACTURE MODEL APPLIED TO LIQUID PHASE ONLY > 0: ACTIVE FRACTURE MODEL APPLIED TO ALL PHASES MOP2(6) = 0 *** LEVERETT SCALING > 0: RESCALE CAPILLARY PRESSURE: PC=PC_REF*SQRT(K_REF/K) IF ELEMENT-BY-ELEMENT PERMEABILITIES ARE DEFINED MOP2(7) = 0 *** ZERO NODAL DISTANCE = 0: TAKE ABSOLUTE PERMEABILITY FROM OTHER GRID BLOCK > 0: TAKE ABSOLUTE AND RELATIVE PERMEABILITY FROM OTHER GRID BLOCK. MOP2(8) = 0 *** VERSION OF SINK/SOURCE SUBROUTINE = 0: USE ORIGINAL VERSION OF QU > 0: USE QU VERSION 2 MOP2(9) = 0 *** TIME STEP AFTER FORCED TIME-STEP REDUCTION TO HONOR PRINTOUT TIME = 0: CONTINUE WITH TIME STEP USED BEFORE FORCED TIME-STEP REDUCTION > 0: CONTINUE WITH TIME STEP IMPOSED BY FORCED TIME-STEP REDUCTION MOP2(10)= 0 *** PRINT SAVE FILE AFTER EACH OUTPUT *********************************************************************************************************************************** ================================================================================================================================== ARRAY DIMENSIONS (SEE FILE maxsize.inc) ---------------------------------------------------------------------------------------------------------------------------------- MAXEL = 10000 Maximum number of elements MAXCON = 50000 Maximum number of connections MAXK = 2 Maximum number of components MAXEQ = 3 Maximum number of equations MAXPH = 2 Maximum number of phases MAXB = 8 Maximum number of phase-dependent secondary variables MAXSS = 300 Maximum number of sinks/sources MAVTAB = 100 Maximum average number of table entries per sink/source MAXROC = 200 Maximum number of rock types MAXTSP = 5 Maximum number of specified time steps, divided by eight MAXLAY = 10 Maximum number of reservoir layers for wells on deliverability MXRPCP = 14 Maximum number of parameters for relative permeability and capillary pressure functions MXPCTB = 5 Maximum number of points in table for ECM capillary pressure MXTBC = 5 Maximum number of elements with time vs. boundary condition MXTBCT = 100 Maximum number of time vs. pressure data MAXTIM = 1000 Maximum number of calibration times MAXN = 50 Maximum number of parameters to be estimated MAXO = 800 Maximum number of datasets MAXM = 1400 Maximum number of calibration points MAXPD = 1000 Maximum number of paired data MAXR = 500 Maximum number of elements or indices of each parameter or observation MAXBRK = 20 Maximum number of points in time at which SAVE file is written for restart MAXEBRK = 50 Maximum number of elements with new initial conditions after restart MAXCOEFF = 5 Maximum number of coefficients for data modeling functions MAXMCS = 500 Maximum number of Monte Carlo simulations MAXCURVE = 500 Maximum number of curves to be plotted MAXXGR = 3 Dimension of third index of array XGUESSR MTYPE = 25 Number of observation types MPFMT = 6 Number of plot file formats ---------------------------------------------------------------------------------------------------------------------------------- ---------------------------------------------------------------------------------------------------------------------------------- PROGRAM VERSION DATE COMMENT ---------------------------------------------------------------------------------------------------------------------------------- iTOUGH2 Current version iTOUGH2 V6.5 (OCTOBER, 2010) ---------------------------------------------------------------------------------------------------------------------------------- iTOUGH 1.0 1 AUGUST 1992 ITOUGH User's Guide, Version 1.0, Report NIB 92-99 iTOUGH2 2.2 1 FEBRUARY 1994 iTOUGH2 User's Guide, Version 2.2, Report LBL-34581 iTOUGH2 3.0 12 JULY 1996 YMP Software qualification, Report LBNL-39489 iTOUGH2 3.1 1 APRIL 1997 iTOUGH2 Command Reference, Version 3.1, Report LBNL-40041 iTOUGH2 3.2 30 JUNE 1998 YMP Software Qualification, Report LBNL-42002 iTOUGH2 3.3 1 OCTOBER 1998 Parallelization using PVM, Report LBNL-42261 iTOUGH2 4.0 19 JANUARY 1999 Released by Energy Science and Technology Software Center iTOUGH2 5.0 31 JULY 2002 Qualified for use within Yucca Mountain Project ---------------------------------------------------------------------------------------------------------------------------------- WHATCOM 1.0 5 MARCH 1996 Q: WHAT COMPUTER IS USED? A: PC CALLSIG 1.0 5 MARCH 1996 #112: SIGNAL HANDLER (NOT INSTALLED) CPUSEC 1.0 5 MARCH 1996 RETURNS CPU-TIME (VERSION PC) OPENFILE 6.5 15 APRIL 2010 OPENS MOST OF THE FILES LENOS 1.0 1 MARCH 1992 RETURNS LENGTH OF LINE PREC 1.0 1 AUGUST 1992 CALCULATE MACHINE DEPENDENT CONSTANTS ITHEADER 3.2 27 MAY 1998 PRINTS iTOUGH2 HEADER DAYTIM 1.0 5 MARCH 1996 RETURNS DATE AND TIME (VERSION PC) THEADER 6.0 29 JUNE 2007 PRINTS TOUGH2 HEADER SOLVTYPE 1.0 1 OCTOBER 1999 INITIALIZE PARAMETERS FOR THE SOLVER PACKAGE INPUT 6.4 5 MARCH 2009 READ ALL DATA PROVIDED THROUGH FILE *INPUT* CHECKMAX 6.0 4 FEBRUARY 2008 CHECK KEY DIMENSIONS ERRORMSG 5.1 19 SEPTEMBER 2005 PRINTS ERROR MESSAGES TERMINAT 6.5 27 JANUARY 2010 PERFORM ERROR ANALYSIS AND TERMINATE iTOUGH2 ================================================================================================================================== --- 1170th iTOUGH2 simulation job stopped fatally on Provide date in --- CPU time used = 0.00 sec. ********************************************************************************************************************************** * iTOUGH2 stopped due to fatal error. Check all output files for error messages! * **********************************************************************************************************************************