****************************************************
                                    *                                                  *
                                    *    PHAST: A Three Dimensional Flow and Solute    *
                                    *         Reactive Transport Simulator             *
                                    *         Version:                       *
                                    *                                                  *
                                    ****************************************************






                          . Central Oklahoma aquifer,                                                     
                          . demonstration of PHAST                                                        

                         *** Fundamental Information ***
                    Isothermal simulation
                    Solute transport simulation
                    Cartesian coordinates
                    Data is in metric units
                    Time unit selected is years   
                    Solute concentration is expressed as mass fraction

          Iterative generalized conjugate gradient solver with d4 zig-zag reduction is selected

                         *** Problem Dimension Information ***
          Number of nodes in x-direction .................................. NX ...     16
          Number of nodes in y-direction .................................. NY ...      9
          Number of nodes in z-direction .................................. NZ ...      5
          Total number of nodes ........................................... NXYZ .    720

          Number of porous media zones .................................... NPMZ .    424
          Number of specified pressure or mass fraction b.c. .............. NSBC .      2
          Number of specified flux b.c. cells ............................. NFBC .     76
          Number of leakage cells ......................................... NLBC .     25
          Number of river leakage cells ................................... NRBC .     34
          Number of wells ................................................. NWEL .      1

          A4 array dimension requirement (D4ZGCG iterative solver)......... NRAL .        0 elements
 ------------------------------------------------------------------------------------------------------------------------


                              ***  Static Data ***


                              X-Direction Node Coordinates    (m)

           1           2           3           4           5           6           7           8           9          10
           0.00     6000.00    12000.00    18000.00    24000.00    30000.00    36000.00    42000.00    48000.00    54000.00



          11          12          13          14          15          16
       60000.00    66000.00    72000.00    78000.00    84000.00    90000.00




                              Y-Direction Node Coordinates   (m)

           1           2           3           4           5           6           7           8           9
           0.00     6000.00    12000.00    18000.00    24000.00    30000.00    36000.00    42000.00    48000.00




                              Z-Direction Node Coordinates   (m)

           1           2           3           4           5
           0.00      100.00      200.00      300.00      400.00




                              Z-Axis is Positive Vertically Upward


                                        *** Fluid Properties ***
                                             Physical
                         Fluid compressibility ...................................... BP ...   0.00E+00  (1/Pa )

                         Reference pressure for density ............................. P0 ..         0.0  (Pa )
                         Reference temperature for density .......................... T0 ..        15.0  (Deg.C)
                         Fluid density .............................................. DENF0  1000.0      (kg/m^3)
                         Fluid density at solute  mass fraction of  0.0500........... DENF1  1000.0      (kg/m^3)

                         Viscosity factor ........................................... VISFAC -1.000E-03


                                        *** Well Data ***
     -------------------------------------------------------------------------------------------------------------------
     Well       Location            Screen Interval          Screen Interval              Calculation     Well Diameter
     No.      X          Y           Z1          Z2       Top Depth  Bottom Depth             Type             (m)
            (m)       (m)        (m)        (m)        (m)         (m)
     -------------------------------------------------------------------------------------------------------------------
       1  1.200E+04  3.600E+04   90.0       110.                                              11               5.080E-02

  Well No.   1
    Specified Flow Rate                              
    Allocation by Mobility                           
    Explicit Layer Rates                    

                                   Layer     Effective Ambient      Well Flow Factor
                                   No.         Permeability
                                                  (m^2)                (m^3)
                  ---------------------------------------------------------------------------

                    Well No.   1
                                     2           1.400E-12             1.557E-11
                                                 1.400E-12

                                   A free-surface water table is specified for this simulation

                                        *** Calculation Information ***
          Backwards-in-time (implicit) differencing for temporal derivative
          Backwards-in-space (upstream) differencing for advective terms
          The cross-derivative solute flux terms will NOT BE calculated

          Direction index for d4 zig-zag renumbering ...................... IDIR..    1
          Incomplete LU [f] or modified ILU [t] factorization ............. MILU..     T
          Number of search directions before restart ...................... NSDR..    5
          Tolerance on iterative solution ................................. EPSSLV 1.0E-10

 ------------------------------------------------------------------------------------------------------------------------

                    Automatic Time Step Control Parameters for Steady Flow
          Maximum pressure change allowed per time step.................... DPTAS    1.18E+06  (Pa )
          Minimum time step required ...................................... DTIMMN    2.0      (yr )
          Maximum time step allowed ....................................... DTIMMX   2.00E+06  (yr )


                    Initial fluid in region ...............................  2.692531E+14  (kg) ;  2.692531E+11  (m^3)

                    Initial solute in region:H         ....................  1.475961E+09  (kg)

                    Initial solute in region:O         ....................  1.910320E+11  (kg)

                    Initial solute in region:As        ....................  2.917440E+05  (kg)

                    Initial solute in region:C         ....................  1.852146E+10  (kg)

                    Initial solute in region:Ca        ....................  7.263407E+12  (kg)

                    Initial solute in region:Cl        ....................  9.168961E+13  (kg)

                    Initial solute in region:Mg        ....................  1.524286E+12  (kg)

                    Initial solute in region:Na        ....................  4.835982E+13  (kg)

                    Initial solute in region:S         ....................  5.899490E+10  (kg)

                    Initial solute in region:Charge    .................... -4.281038E+00  (kg)

 ------------------------------------------------------------------------------------------------------------------------
 ------------------------------------------------------------------------------------------------------------------------


                                        *** Transient Data ***

     Time ..............................   0.000      (yr )


                                        *** Calculation Information ***

          Fixed time step length .......................................... DELTIM  2.000E+03  (yr )
          Time at which next set of transient
               parameters will be read .................................... TIMCHG  1.000E+05  (yr)

 ------------------------------------------------------------------------------------------------------------------------


                                       ***** Simulation Completed *****
               Last time value calculated ...................................... 1.0000E+05 (yr)
               Last time step index ............................................      50