DATA SET 4: PRINTER OUTPUT, DISK STORAGE CONTROL, AND TIMES FOR RESETTING STEP SIZE

The number of lines here depends on the number of time increments, NTI, and the number of times to reset the time step, NDTCHG. The number of lines is [(NTI/80+1)*2] + NLINE: (NTI/80+1) lines for printer output control, (NTI/80+1) lines for storage control, and NLINE lines for time-step-size resetting control, which depends on NDTCHG. NLINE lines are required to enter NDTCHG real numbers with free format.

Line 1 to Line [(NTI/80+1)]: FORMAT(80I1)

KPR0, KPR(I), I=1,NTI, where

KPR0 = Printer control for initial and steady state conditions.

KPR(I) = Printer control for the I-th time step:

0 = Print nothing.

1 = Print chemical component data from within biogeochemical subroutine.

2 = Print above plus chemical component species data.

3 = Print above plus product species data.

4 = Print above plus stoichiometric coefficients and reaction data.

5 = Print above plus mass flow rates (FLOW, FRATE, and TFLOW).

6 = Print above plus velocities.

Note: summary concentration values at each node are printed when a convergent solution is reached for each time step. Options 1 through 4 provide more detailed output on the chemical speciation in the system from within the biogeochemical portion of the program every NHGCI iterations. NHGCI is specified in Data Set 1.

Line [(NTI/80+1)+1] to Line [(NTI/80+1)*2]: FORMAT(80I1)

KDSK0, KDSK(I), I=1,NTI, where

KDSK0 = Auxiliary storage control for steady-state and initial conditions.

KDSK(I) = Auxiliary storage control for the I-th time-step, the same allowable values as for KDSK0;

0 = No storage.

1 = Store in Logical Unit 12, in FORTRAN binary.
(This file can be processed by the postprocessor PostHBGC.)

Line greater than line (NTI/80+1)*2: Free-field format input contains NDTCHG real numbers:

TDTCH(I), I=1, NDTCHG, where TDTCH(I) is the time (T) at which the time-step size is to be reset to the initial value for the I-th time.

v1.0 - 1/14/1999

Jin-Ping Gwo, email: jgwo@umbc.edu