Began at      --> Wed Oct  6 16:40:30 EDT 1999
  Proc            0 reporting for duty
  
  
 ***************************************************
                                                    
         aaaa    dddd    ppppp    aaaa    ccccc     
        a   a   d   d   p   p    a   a   c          
       aaaaa   d   d   ppppp    aaaaa   c           
      a   a   d   d   p        a   a   c            
     a   a   dddd    p        a   a   ccccc         
                                                    
            version ---->  ADPAC V1.0               
                                                    
 adpac 3-d multiple-block euler/navier-stokes solver
  (a)dvanced (d)ucted (p)ropfan (a)nalysis (c)odes  
      Serial/Parallel Production Code ADPAC V1.0    
              release date Feb., 1998               
                                                    
 ***************************************************
  
  
  version V1.0 notes:
  *******************
  
  Available turbulence models:
  Mixing-Length (algebraic)   
  Baldwin-Lomax (algebraic)   
  Spalart-Allmaras (1-Equation)
  Goldberg k-R     (2-Equation)
  (wall functions are available)
  wall functions now ON in separated flow
  
  thermodynamic properties are for air
  with variable or constant gamma
  2-d or 3-d mesh blocks supported
  serial/parallel from same source
  implicit algorithm is available
  FIXED B.C. is OK in parallel
  PINT B.C. is OK with multigrid/parallel
  2-D and 3-D yield identical results
  (corrects previous diffs due to turbulence models)
  Other bug corections up to 2/02/98 EJH
  
  Searching Input File for ADPAC Case Name...
  *******************************************
  data item identified --->CASENAME             = box                 
  data item identified --->FCART                =    1.000000    
  data item identified --->FMULTI               =    1.000000    
  
  case name defined as   :                                                                              box
  multigrid level set as :            1
  grid type set as      : Cartesian
  
  Solution type is explicit: IIMPLIC =            0
  1-Equation turbulence model is disabled: I1EQ =            0
  2-Equation turbulence model is disabled: I2EQ =            0
  
 **************************************************
                begin grid input                   
 **************************************************
  
  
  mesh file name: 
 box.mesh                                                                        
  ***************
  
  number of grids: (nblks) =            1
  ****************
  
  grid block dimensions :
  ***********************
  
      block       imx       jmx       kmx
      -----       ---       ---       ---
          1        25        13        13
  
  block/processor assignment information
  **************************************

  Block/Processor Assignment file 
 box.blkproc                                                                     
  does not exist
  Default block/processor assignment used

  IPR(   Block #  )     Processor
  -----------------     ---------
  ipr(           1) =           0
  
  initialize storage arrays and pointer locations 
  *********************************************** 
  
     block     3-d ijk      2-d ij     1-d iii      length  total length
         1           1           1           1        5096        5096      26
  
  array sizes required for this mesh (No multi-grid)
  **********************************
  
  nra3d (3-d array size) >=         5096
  nra1d (1-d array size) >=           26
  nbl2d (2-d block size) >=          676
  
  
  program array dimensions: 
  *************************
  
  nra3d (3-d array size) =       600001
  nra1d (1-d array size) =         3201
  nbl2d (2-d block size) =       100001
  
  
  begin reading grid block coordinates
  ************************************
  
  grid block            1 size           26          14          14
  grid block            1 coordinates successfully read
  
 **************************************************
               begin flow input                    
 **************************************************
  


    input file read on unit   5
    ***********************
  data item identified --->CASENAME             = box                 
  data item identified --->DIAM                 =    1.000000    
  data item identified --->RMACH                =   0.5000000    
  data item identified --->CFL                  =   -5.000000    
  data item identified --->FINVVI               =   0.0000000E+00
  data item identified --->FCART                =    1.000000    
  data item identified --->FREST                =   0.0000000E+00
  data item identified --->FFULMG               =   0.0000000E+00
  data item identified --->FMULTI               =    1.000000    
  data item identified --->FNCMAX               =    10.00000    
  
  inviscid flow calculation selected
  **********************************
  
 4-stage runge_kutta algorithm selected
  
  
  
  The following debug triggers have been enabled
  **********************************************
  


    input file data    
    ***************
     rm     ( reference mach number          ) =      0.5000
     gamma  ( specific heat ratio            ) =      1.4000
     cfl    ( cfl number (+ time accurate    ) =      5.0000
     cflmax ( max cfl number time-step scheme) =      2.5000
     zetarat( exponent for diss/res.sm. ratio) =      0.6000
     fvtsfac( viscous time step factor       ) =      2.5000
     vt     ( =1, time acc., =0, local time  ) =      0.0000
     vis2   ( 2nd order dissipation factor   ) =      0.5000
     vis4   ( 4th order dissipation factor   ) =      0.0156
     fncmax ( maximum iterations             ) =     10.0000
     frest  ( 1 = restart solution           ) =      0.0000
     save   ( 1 = save solution              ) =      1.0000
     p3dprt ( 0 = no plot3d output           ) =      1.0000
     fitchk ( checkpoint iteration interval  ) =    100.0000
     viscg2 ( coarse grid 2nd order diss coef) =      0.1250
     fistep ( recalculate time step          ) =      1.0000
     fnprnt ( 1 = print flowfield            ) =      1.0000
     fiprnt ( 1 = print partial flowfield    ) =      1.0000
     finvvi ( 0 = inviscid flow solution     ) =      0.0000
            ( 1 =  viscous flow solution     ) =
     epsx   ( i implicit smoothing factor    ) =      1.0000
     epsy   ( j implicit smoothing factor    ) =      1.0000
     epsz   ( k implicit smoothing factor    ) =      1.0000
     ftotsm ( 1 = post multigrid smoothing   ) =      0.0000
     epstot ( post multigrid smoothing coef. ) =      0.2000
     iunint ( # it. for unsteady plot3d file ) =      0.0000
     itimei ( #it between time step update   ) =      1.0000
     iturbi ( #it between turb. model update ) =      1.0000
     iturbb ( #it before starting turb model ) =     10.0000
     iturbf ( #it before freezing turb model ) =      0.0000
     ibconf ( #it before freezing bound. cond) =      0.0000
     iwallf ( wall function trigger (on=1)   ) =      1.0000
     imassav( mxng pln avrgng (area=1,mass=2)) =      1.0000
     itimerm( 1= quit if too little time left) =      0.0000
     fdeltat( implicit time step (sec.)      ) =      0.0000
     fsolve ( solution trigger               ) =      1.0000
     iimplicit (implicit algorithm trigger   ) =           0
     fturbcht(   1) (C/O Turb. Model Trigger ) =      0.0000
     igrafix( trigger for graphics display   ) =      0.0000
     igrafint(#it between graphics update    ) =      0.0000
     iimgsav( trigger for image capture      ) =      0.0000
     iimgint( #it between image captures     ) =      0.0000
     fcart  ( 0 = Cylindrical, 1 = Cartesian ) =      1.0000
     diam   ( blade diameter (ft)            ) =      1.0000
     tref   ( total temperature (deg. r)     ) =    518.6700
     pref   ( total pressure    (lbf/ft**2)  ) =   2116.2200
     rgas   ( gas constant ((ft/s)**2/deg r) ) =   1716.2600
     pr     ( prandtl number                 ) =      0.7000
     prt    ( turbulent prandtl number       ) =      0.9000

non-dimensional initial values calculated as:
********************************************
     rho0   ( initial density              ) =      0.8852
     u0     ( initial axial velocity       ) =      0.5774
     v0     ( initial radial velocity      ) =      0.0000
     w0     ( initial theta velocity       ) =      0.0000
     ei0    ( initial internal energy      ) =      2.5476
     h0     ( initial enthalpy             ) =      3.5000
     p0     ( initial pressure             ) =      0.8430
     t0     ( initial temperature          ) =      0.9524
     dmu0   ( initial viscosity            ) =      0.0000

dimensional reference values calculated as:
*******************************************
     rho0   ( density         (slugs/ft**3)  ) = 0.210433E-02
     u0     ( axial velocity    (ft/s)       ) =  544.724    
     v0     ( radial velocity   (ft/s)       ) = 0.000000E+00
     w0     ( theta velocity    (ft/s)       ) = 0.000000E+00
     ei0    ( static int.energy (ft-lbf/lbm) ) = 0.226782E+07
     h0     ( static enthalpy   (ft-lbf/lbm) ) = 0.311560E+07
     p0     ( static pressure   (lbf/ft**2)  ) =  1784.01    
     t0     ( static temperature(deg. r)     ) =  493.971    
     dmu0   ( viscosity         (lbf-s/ft**2)) = 0.359781E-06

values used for nondimensionalization:
**************************************
     rho    ( density         (slugs/ft**3)  ) = 0.237731E-02
     u      ( axial velocity    (ft/s)       ) =  943.490    
     ei     ( total int. energy (ft-lbf/lbm) ) =  890173.    
     p      ( total pressure    (lbf/ft**2)  ) =  2116.22    
     t      ( total temperature (deg. r)     ) =  518.670    
     dmu    ( viscosity         (lbf-s/ft**2 ) =  2.24297    
     omega  ( rotational speed  (rad/s)      ) =  943.490    
     length ( length            (ft)         ) =  1.00000    
     time   ( time              (seconds)    ) = 0.105990E-02

multigrid solution parameters specified as:
*******************************************
     nmulti ( number of multi-grid grid level) =    1
     nsubit ( multigrid coarse grid subiterat) =    1
     nfulmg ( full multigrid solution trigger) =    0
     ncoag1 ( full multigrid initial level   ) =    1
     ncoag2 ( full multigrid final   level   ) =    2
     nitfmg ( full multigrid coarse mesh iter) =    0
     nrdmul ( multigrid bc read trigger (1-y)) =    0

 Block Solution Information Table
 ********************************

Block           RPM   Dimensional    Non-Dimen. Cyl/Cart.    C-Grid Processor
  Num   Revolutions         Omega         Omega   Trigger Turb Trig    Number
         Per Minute       (Rad/s)                  1-Cart 1-On0-Off      MYID
----- ------------- ------------- ------------- --------- --------- ---------
    1  0.000000E+00  0.000000E+00  0.000000E+00         1         0         0
  
  input data successfully entered 
  ******************************* 
  
  
 **************************************************
               begin boundary data input           
 **************************************************
  
Opening casename.boundata, file is box.boundata                  
INLETG     1   1  I  I  P  P  J  K   1   1   1  13   1  13   1  13   1  13     
  b.c. data item identified --->INLETG              
  
    PTOT   TTOT                                                                 
    1.0    1.0                                                                  
    1.000000       1.000000      0.0000000E+00  1.0000000E-03  0.0000000E+00
   0.0000000E+00  0.0000000E+00  0.0000000E+00
  Stored emass =   0.0000000E+00
  Stored relax =   0.0000000E+00
                                                                               
# j=1 is y=1.0                                                                 
# outward normal is +y                                                         
SSIN       1   1  J  J  P  P  I  K   1   1   1  25   1  13   1  25   1  13     
  b.c. data item identified --->SSIN                
# if this were a 'free' bc, you'd need this                                    
    PTOT   TTOT  EMINF   ALPHA                                                 
    1.0    1.0    0.5   0.0                                                    
                                                                               
# j=13 is y=0.0                                                                
# outward normal is -y                                                         
FREE       1   1  J  J  M  M  I  K  13  13   1  25   1  13   1  25   1  13     
  b.c. data item identified --->FREE                
  
    PTOT   TTOT  EMINF   ALPHA                                                  
    1.000000       1.000000      0.5000000      0.0000000E+00  0.0000000E+00
   1.0000000E-03
  
                                                                               
# k=1 is z=0.0                                                                 
# outward normal is -z                                                         
SSIN       1   1  K  K  P  P  I  J   1   1   1  25   1  13   1  25   1  13     
  b.c. data item identified --->SSIN                
# if this were a 'free' bc, you'd need this                                    
    PTOT   TTOT  EMINF   ALPHA                                                 
    1.0    1.0    0.5   0.0                                                    
                                                                               
# k=13 is z=1.0                                                                
# outward normal is +z                                                         
FREE       1   1  K  K  M  M  I  J  13  13   1  25   1  13   1  25   1  13     
  b.c. data item identified --->FREE                
  
    PTOT   TTOT  EMINF   ALPHA                                                  
    1.000000       1.000000      0.5000000      0.0000000E+00  0.0000000E+00
   1.0000000E-03
  
                                                                               
EXITG      1   1  I  I  M  M  J  K  25  25   1  13   1  13   1  13   1  13     
  b.c. data item identified --->EXITG               
  
    PEXIT    EMDOT   PRELAX                                                     
  0.843019 40.000000  0.000000
  
  b.c. data item identified --->ENDDATA             
  
  boundary data successfully entered 
  ********************************** 
  
 
N    B       L  L  L  L  L  L  L  L   L   L   M   M   N   N   M   M   N   N 
U    C       B  B  F  F  D  D  S  S   1   2   1   1   1   1   2   2   2   2 
M    T       L  L  A  A  I  I  P  P   L   L   L   L   L   L   L   L   L   L 
B    Y       O  O  C  C  R  R  E  E   I   I   I   I   I   I   I   I   I   I 
E    P       C  C  E  E  1  2  C  C   M   M   M   M   M   M   M   M   M   M 
R    E       K  K  1  2        1  2           1   2   1   2   1   2   1   2 
             1  2 
     ------ -- -- -- -- -- -- -- -- --- --- --- --- --- --- --- --- --- --- 
   1 INLETG  1  1  I  I  P  P  J  K   1   1   1  13   1  13   1  13   1  13
   2 SSIN    1  1  J  J  P  P  I  K   1   1   1  25   1  13   1  25   1  13
   3 FREE    1  1  J  J  M  M  I  K  13  13   1  25   1  13   1  25   1  13
   4 SSIN    1  1  K  K  P  P  I  J   1   1   1  25   1  13   1  25   1  13
   5 FREE    1  1  K  K  M  M  I  J  13  13   1  25   1  13   1  25   1  13
   6 EXITG   1  1  I  I  M  M  J  K  25  25   1  13   1  13   1  13   1  13
  
  
  *****************************************
     Examining all Block Outer Boundary    
       Conditions for Completeness         
  *****************************************
  
  
  
  set up block data:
  ******************
  
  cell area terms calculated for grid block   #           1
  calculating cell volumes for grid block   #           1
  

Mesh Cell Expansion  Ratio Table
********************************

Block I Index   I    J    K   J Index   I    J    K   K Index   I    J    K    |
Numbr  Expans.  Max  Max  Max  Expans.  Max  Max  Max  Expans.  Max  Max  Max  |
    N   Ratio                   Ratio                   Ratio                  V
----- -------- ---- ---- ---- -------- ---- ---- ---- -------- ---- ---- ---- --
    1  1.00000   15    2    5  1.00000    9    2    2  1.00000    4    2    8   
  
 * - WARNING!! 
     ********* 

 Solution Accuracy and Convergence May Be Adversely 
 Affected By Large Mesh Cell Expansion Ratios ---
 Mesh Blocks indicated by a * in the final column 
 Possess Coordinate Mesh Cell Expansion Ratios Which
 Exceed the Recommended Maximum Value of  1.300000    

  starting imploop            1
  calling impupd


 convergence history ---> mesh level    1
 ****************************************

  cycle      log10  <--- max ---->     log10  #ss      mass          mass      
 number       max       error           rms   pts.     flow          flow
             error bl#   i   j   k     error            in           out      
 ------  --------- --- --- --- --- --------- ----- --------------  ------------

      1    -6.9550   1   2  13  13   -9.2393     0   36.880405     36.880394
      2    -6.8251   1  24  12  12   -9.1769     0   36.880405     36.880383
      3    -6.8189   1   2  10  13   -9.1783     0   36.880394     36.880390
      4    -6.6871   1  25  12  11   -9.1450     0   36.880390     36.880383
      5    -6.8182   1   7  12  10   -9.1347     0   36.880394     36.880383
      6    -6.6993   1   2   5   8   -9.1566     0   36.880398     36.880379
      7    -6.8091   1  25   9  10   -9.1531     0   36.880398     36.880379
      8    -6.8251   1  15   5  12   -9.1446     0   36.880398     36.880375
      9    -6.6988   1   2  13   5   -9.1315     0   36.880398     36.880383
     10    -6.8172   1  23   6   6   -9.1218     0   36.880402     36.880379
  
  new restart file name 
 box.restart.new                                                                 
  
  number of blocks: (nblks) =            1
  


*** plot3d absolute output on unit   9 ***
*** file box.p3dabs                    


*** plot3d relative output on unit   8 ***
*** file box.p3drel                    
  
 Approximate CPU Time Breakdown (Parallel Runs Only)
 ---------------------------------------------------
 TIMING INIT   tinit  0.0000000000000000E+00 Initialization Time 
 TIMING GRID   tgrid  0.0000000000000000E+00 Grid Input Time
 TIMING COMM   tcomm  0.0000000000000000E+00 Boundary Update Time
 TIMING SOLV   tsolv  0.0000000000000000E+00 Flow Solver Time
 TIMING REST   trest  0.0000000000000000E+00 Restart File Time
 TIMING OUT    tout   0.0000000000000000E+00 Output file Time
 TIMING TOTAL tottim  0.0000000000000000E+00
  
  processing completed 
  execution normally terminated 
3.7u 0.1s 0:06 62% 0+0k 91+25io 70pf+0w
Finished at   --> Wed Oct  6 16:40:41 EDT 1999