MPI initialized with           24  MPI processes
 MPI initialized with            1  threads
 PROBINIT max time   =   1.649379730224609E-002
 PROBINIT total time =   1.649379730224609E-002
 
 finished reading eos table
 imax= 211 jmax=  71
 temp(1)   =  1.000E+04 temp(jmax)   =  1.000E+11
 ye*den(1) =  1.000E-10 ye*den(imax) =  1.000E+11
 
                                                            
------------------------------------------------------------------------------
 reading initial model
       14000 points found in the initial model file
           6  variables found in the initial model file
 WARNING: variable not found: iron-56
 WARNING: hydrogen-1 not provided in inputs file
 WARNING: carbon-13 not provided in inputs file
 WARNING: nitrogen-13 not provided in inputs file
 WARNING: nitrogen-14 not provided in inputs file
 WARNING: nitrogen-15 not provided in inputs file
 WARNING: oxygen-14 not provided in inputs file
 WARNING: oxygen-15 not provided in inputs file
 WARNING: oxygen-16 not provided in inputs file
 WARNING: oxygen-17 not provided in inputs file
 WARNING: oxygen-18 not provided in inputs file
 WARNING: fluorine-17 not provided in inputs file
 WARNING: fluorine-18 not provided in inputs file
 WARNING: fluorine-19 not provided in inputs file
 WARNING: neon-18 not provided in inputs file
 WARNING: neon-19 not provided in inputs file
 WARNING: neon-20 not provided in inputs file
 WARNING: magnesium-22 not provided in inputs file
 WARNING: sulfur-30 not provided in inputs file
 WARNING: nickel-56 not provided in inputs file
------------------------------------------------------------------------------
 model file mapping, level:           1
      dr of MAESTRO base state =                              0.5000000000    
      dr of input file data =                                 0.5000000000    
  
      maximum radius (cell-centered) of input model =          6999.750000    
  
 setting r_cutoff to         1689
  
 Maximum HSE Error =    12.0823183940673     
    (after putting initial model into base state arrays, and
     for density < base_cutoff_density)
------------------------------------------------------------------------------
  
 number of MPI processes =           24
 number of threads       =            1
  
 number of dimensions    =            2
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
  
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
------------------------------------------------------------------------------
 physical cutoff densities:
          low density cutoff (for mapping the model) =         10000.00000    
          buoyancy cutoff density                           
              (for zeroing rho - rho_0, centrifugal term) =    50000.00000    
          anelastic cutoff =                                   10000.00000    
          sponge start density =                               250000.0000    
                                                            
 thermodynamics cutoffs:
          EOS temperature floor =                              5000000.000    
          EOS density floor =                                 0.1000000000E-04
------------------------------------------------------------------------------
  
------------------------------------------------------------------------------
 DOING THE INITIAL VELOCITY PROJECTION
------------------------------------------------------------------------------
  
... Level 1 create thermal coeffs:
 PROJ_TYPE IN HGPROJECT:           1
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection  0.0000000000E+00   0.0000000000E+00
... y-velocity before projection  0.0000000000E+00   0.0000000000E+00
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1952243714E-03   0.1952243708E-03
... y-velocity  after projection -0.1652219926E-03   0.5959109948E-03
 
 Call to firstdt for level           1 gives dt_lev =  4.460584682654186E-005
 Multiplying dt_lev by init_shrink; dt_lev =  4.460584682654186E-006
 Minimum firstdt over all levels =  4.460584682654186E-006
 
  
------------------------------------------------------------------------------
 DOING           3 INITIAL DIVU ITERATIONS
------------------------------------------------------------------------------
  
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.580312479911730E-004
 
 PROJ_TYPE IN HGPROJECT:           2
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.1952243714E-03   0.1952243708E-03
... y-velocity before projection -0.1652219926E-03   0.5959109948E-03
 
 doing hgproject with tolerance, eps =   1.000000000000000E-008
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.4712413451E-03   0.4712413450E-03
... y-velocity  after projection -0.4170017939E-03   0.2724799582E-03
 
 Call to estdt at end of istep_divu_iter =           1
 gives dt =  4.460584682654186E-005
 Multiplying dt by init_shrink; dt =  4.460584682654186E-006
 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.580312479911730E-004
 
 PROJ_TYPE IN HGPROJECT:           2
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.4712413451E-03   0.4712413450E-03
... y-velocity before projection -0.4170017939E-03   0.2724799582E-03
 
 doing hgproject with tolerance, eps =   1.000000000000000E-010
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.4712514047E-03   0.4712514046E-03
... y-velocity  after projection -0.4170119195E-03   0.2724906200E-03
 
 Call to estdt at end of istep_divu_iter =           2
 gives dt =  4.460584682654186E-005
 Multiplying dt by init_shrink; dt =  4.460584682654186E-006
 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.580312479911730E-004
 
 PROJ_TYPE IN HGPROJECT:           2
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.4712514047E-03   0.4712514046E-03
... y-velocity before projection -0.4170119195E-03   0.2724906200E-03
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.4712514051E-03   0.4712514051E-03
... y-velocity  after projection -0.4170119194E-03   0.2724906201E-03
 
 Call to estdt at end of istep_divu_iter =           3
 gives dt =  4.460584682654186E-005
 Multiplying dt by init_shrink; dt =  4.460584682654186E-006
 
 
------------------------------------------------------------------------------
 DOING           3 INITIAL PRESSURE ITERATIONS
------------------------------------------------------------------------------
 
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  4.460584682654186E-006
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   6.580312479911730E-004
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.4712514037E-03   0.1120928960E+04
... level / min : MAC vels    1  -0.4712514037E-03  -0.2942382025E+03
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.15559621E+09
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13499340E-01
... level / max : MAC vels    1  -0.2070269500E+03  -0.1687828051E+03
... level / min : MAC vels    1   0.2070269500E+03   0.5885520700E+03
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077798E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602884982E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396856371E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.580532680192355E-004
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.4712514037E-03   0.1120928668E+04
... level / min : MAC vels    1  -0.4712514037E-03  -0.2942384944E+03
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.15559862E+09
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.25592791E-02
... level / max : MAC vels    1  -0.2070266824E+03  -0.1687827611E+03
... level / min : MAC vels    1   0.2070266824E+03   0.5885516571E+03
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077798E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602884988E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396856375E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.4712513931E-03   0.4712513930E-03
... new min/max : y-velocity       -0.2606410225E+05   0.7209535049E+04
 
 PROJ_TYPE IN HGPROJECT:           3
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.4712513931E-03   0.4712513930E-03
... y-velocity before projection -0.2606410225E+05   0.7209535049E+04
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.2859865554E+04   0.2859865554E+04
... y-velocity  after projection -0.4380866017E+04   0.1714660656E+04
 
 Timing summary:
    Advection       :    4.03486871719360       seconds
    MAC   Projection:    1.13203358650208       seconds
    Nodal Projection:    2.95156764984131       seconds
    Thermal         :    2.20339894294739       seconds
    Reactions       :    575.086706161499       seconds
    Misc            :   0.743284940719604       seconds
  
 Time to advance timestep:    586.075458049774       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  4.460584682654186E-006
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   6.580507678635298E-004
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.1429995567E+04   0.1129359756E+05
... level / min : MAC vels    1  -0.1429995567E+04  -0.3143035985E+04
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13288369E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.53079271E+00
... level / max : MAC vels    1  -0.2572750793E+03  -0.2044037791E+03
... level / min : MAC vels    1   0.2572750793E+03   0.6020573981E+03
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077790E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602882509E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396854964E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.580533332964969E-004
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.1429995567E+04   0.1129359757E+05
... level / min : MAC vels    1  -0.1429995567E+04  -0.3143035980E+04
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13288371E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.29673121E+00
... level / max : MAC vels    1  -0.2572815731E+03  -0.2044081626E+03
... level / min : MAC vels    1   0.2572815731E+03   0.6020659379E+03
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077790E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602882506E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396854962E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2859811773E+04   0.2859811773E+04
... new min/max : y-velocity       -0.4774553732E+04   0.3754462557E+04
 
 PROJ_TYPE IN HGPROJECT:           3
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2859811773E+04   0.2859811773E+04
... y-velocity before projection -0.4774553732E+04   0.3754462557E+04
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.2850249451E+04   0.2850249450E+04
... y-velocity  after projection -0.4673732176E+04   0.2546699227E+04
 
 Timing summary:
    Advection       :    4.33721208572388       seconds
    MAC   Projection:    1.03315901756287       seconds
    Nodal Projection:    2.36448311805725       seconds
    Thermal         :    2.85797882080078       seconds
    Reactions       :    574.794075965881       seconds
    Misc            :   0.741092681884766       seconds
  
 Time to advance timestep:    585.971430778503       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  4.460584682654186E-006
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   6.580507716708660E-004
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.1425212886E+04   0.1130964237E+05
... level / min : MAC vels    1  -0.1425212886E+04  -0.3136102834E+04
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13588527E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.54343795E+00
... level / max : MAC vels    1  -0.2875651700E+03  -0.1878958386E+03
... level / min : MAC vels    1   0.2875702931E+03   0.6018607094E+03
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077791E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602882825E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396855144E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.580533336254790E-004
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.1425212886E+04   0.1130964238E+05
... level / min : MAC vels    1  -0.1425212886E+04  -0.3136102829E+04
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13588528E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.30069110E+00
... level / max : MAC vels    1  -0.2875620065E+03  -0.1878996059E+03
... level / min : MAC vels    1   0.2875671296E+03   0.6018692919E+03
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077791E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602882822E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396855142E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2850248101E+04   0.2850248103E+04
... new min/max : y-velocity       -0.4953197750E+04   0.5871232792E+04
 
 PROJ_TYPE IN HGPROJECT:           3
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2850248101E+04   0.2850248103E+04
... y-velocity before projection -0.4953197750E+04   0.5871232792E+04
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.2846388403E+04   0.2846413666E+04
... y-velocity  after projection -0.4755600030E+04   0.6048591180E+04
 
 Timing summary:
    Advection       :    3.82881617546082       seconds
    MAC   Projection:    1.32016110420227       seconds
    Nodal Projection:    2.36971974372864       seconds
    Thermal         :    2.40661478042603       seconds
    Reactions       :    577.819767951965       seconds
    Misc            :   0.739331245422363       seconds
  
 Time to advance timestep:    588.027648925781       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 Writing 1D data to xrb_plt00000
 Writing cell-centered 1D data to xrb_plt00000/model_cc_00000
 Writing edge-centered 1D data to xrb_plt00000/model_ec_00000
 
 
 BEGIN MAIN EVOLUTION LOOP WITH dt =  4.460584682654186E-006
 
 MEMORY STATS AT START OF TIMESTEP            1
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.4712514051E-03   0.2724906201E-03
... level / min : old vels    1  -0.4712514051E-03  -0.4170119194E-03
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  4.460584682654186E-006
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   6.580507720061411E-004
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.1423295321E+04   0.1138680993E+05
... level / min : MAC vels    1  -0.1423282690E+04  -0.3134002359E+04
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13649537E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.54301697E+00
... level / max : MAC vels    1  -0.4791027400E+03  -0.2125226549E+03
... level / min : MAC vels    1   0.4791378601E+03   0.6020304163E+03
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077791E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602882905E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396855189E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.580533325914598E-004
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.1423295321E+04   0.1138680993E+05
... level / min : MAC vels    1  -0.1423282690E+04  -0.3134002355E+04
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13649538E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.30075443E+00
... level / max : MAC vels    1  -0.4791071615E+03  -0.2125266338E+03
... level / min : MAC vels    1   0.4791422808E+03   0.6020389818E+03
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000003E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1342077791E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000031E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000019E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.8602882902E-11
... new min/max : oxygen-17          0.0000000000E+00   0.1700000041E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000000E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800001664E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000169E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900000389E-28
... new min/max : neon-20            0.0000000000E+00   0.1396855188E-20
... new min/max : magnesium-22       0.0000000000E+00   0.2200000003E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968498E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159967E+21   0.1397900334E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2846383441E+04   0.2846408701E+04
... new min/max : y-velocity       -0.9902556257E+04   0.1005107555E+05
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2846383441E+04   0.2846408701E+04
... y-velocity before projection -0.9902556257E+04   0.1005107555E+05
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.3437603196E+04   0.3437684833E+04
... y-velocity  after projection -0.6637825144E+04   0.9503243180E+04
 
 Timing summary:
    Advection       :    4.24980974197388       seconds
    MAC   Projection:    1.22225403785706       seconds
    Nodal Projection:    2.31442403793335       seconds
    Thermal         :    1.92758607864380       seconds
    Reactions       :    580.334980010986       seconds
    Misc            :   0.975484848022461       seconds
  
 Time to advance timestep:    590.789114952087       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            1
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.3437603196E+04   0.3437684833E+04
... min/max : y-velocity       -0.6637825144E+04   0.9503243180E+04
... min/max : gpix              -0.1279097353E+16   0.1279066899E+16
... min/max : gpiy              -0.8575847209E+16   0.2366475815E+16
 
STEP =      1  TIME = 4.4605846827E-06  DT = 4.4605846827E-06
 
 MEMORY STATS AT START OF TIMESTEP            2
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.3437684833E+04   0.9503243180E+04
... level / min : old vels    1  -0.3437603196E+04  -0.6637825144E+04
 
 Call to estdt at beginning of step           2
 gives dt =  1.171453936811436E-005
 dt_growth factor limits the new dt =  4.906643150919605E-006
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  4.906643150919605E-006
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   7.239016924644639E-004
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.5317684596E+04   0.1904649407E+05
... level / min : MAC vels    1  -0.5317558468E+04  -0.2288939360E+05
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.62797945E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.60775799E+00
... level / max : MAC vels    1  -0.1045095584E+05  -0.1172360614E+05
... level / min : MAC vels    1   0.1045115094E+05   0.1643468650E+05
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000011E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1430424662E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000095E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000059E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.2666429140E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000126E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000001E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800005158E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000523E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001206E-28
... new min/max : neon-20            0.0000000000E+00   0.1341580533E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000010E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968671E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159748E+21   0.1397900337E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   7.236506880850782E-004
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.5317684653E+04   0.1904650038E+05
... level / min : MAC vels    1  -0.5317558526E+04  -0.2288938739E+05
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.62798349E+11
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.35025706E+00
... level / max : MAC vels    1  -0.1045113413E+05  -0.1172278554E+05
... level / min : MAC vels    1   0.1045132974E+05   0.1643471402E+05
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000011E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1430424663E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000095E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000059E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.2666429147E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000126E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000001E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800005158E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000523E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001206E-28
... new min/max : neon-20            0.0000000000E+00   0.1341580546E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000010E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968671E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159748E+21   0.1397900337E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.7018321655E+04   0.7018484237E+04
... new min/max : y-velocity       -0.1199618996E+07   0.8883436379E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.7018321655E+04   0.7018484237E+04
... y-velocity before projection -0.1199618996E+07   0.8883436379E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.2828118047E+06   0.2828064330E+06
... y-velocity  after projection -0.1014475220E+07   0.8070222441E+06
 
 Timing summary:
    Advection       :    4.86668109893799       seconds
    MAC   Projection:    1.11878299713135       seconds
    Nodal Projection:    2.77512812614441       seconds
    Thermal         :    3.61002802848816       seconds
    Reactions       :    582.292386293411       seconds
    Misc            :   0.953470706939697       seconds
  
 Time to advance timestep:    595.281323909760       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            2
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.2828118047E+06   0.2828064330E+06
... min/max : y-velocity       -0.1014475220E+07   0.8070222441E+06
... min/max : gpix              -0.9806053563E+17   0.9806238859E+17
... min/max : gpiy              -0.1031709009E+18   0.2309532467E+18
 
STEP =      2  TIME = 9.3672278336E-06  DT = 4.9066431509E-06
 
 MEMORY STATS AT START OF TIMESTEP            3
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.2828064330E+06   0.8070222441E+06
... level / min : old vels    1  -0.2828118047E+06  -0.1014475220E+07
 
 Call to estdt at beginning of step           3
 gives dt =  3.450097061994142E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  3.450097061994142E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   5.088504963135556E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.2907147174E+06   0.8571364573E+06
... level / min : MAC vels    1  -0.2907208333E+06  -0.1019551717E+07
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.18851563E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.28687903E+02
... level / max : MAC vels    1  -0.4438217561E+06  -0.9590409191E+06
... level / min : MAC vels    1   0.4438173777E+06   0.8107744025E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000015E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1479940691E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000131E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000081E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.3678645020E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000174E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800007117E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000722E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001665E-28
... new min/max : neon-20            0.0000000000E+00   0.2553377698E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000014E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968683E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159733E+21   0.1397900305E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   5.095611692916021E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.2907147170E+06   0.8571364617E+06
... level / min : MAC vels    1  -0.2907208329E+06  -0.1019551728E+07
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.18851669E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.50328165E+02
... level / max : MAC vels    1  -0.4438786220E+06  -0.9590526490E+06
... level / min : MAC vels    1   0.4438742420E+06   0.8107871522E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000015E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1479940691E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000131E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000081E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.3678645020E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000174E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800007117E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000722E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001665E-28
... new min/max : neon-20            0.0000000000E+00   0.2553377698E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000014E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968683E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159733E+21   0.1397900305E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2736705684E+06   0.2736639505E+06
... new min/max : y-velocity       -0.6600575247E+06   0.5472389568E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2736705684E+06   0.2736639505E+06
... y-velocity before projection -0.6600575247E+06   0.5472389568E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.2574884764E+06   0.2574807157E+06
... y-velocity  after projection -0.6471051383E+06   0.5399580930E+06
 
 Timing summary:
    Advection       :    4.99487709999084       seconds
    MAC   Projection:   0.951627969741821       seconds
    Nodal Projection:    2.35345816612244       seconds
    Thermal         :    2.64622402191162       seconds
    Reactions       :    535.573900938034       seconds
    Misc            :   0.954963207244873       seconds
  
 Time to advance timestep:    547.048463106155       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            3
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.2574884764E+06   0.2574807157E+06
... min/max : y-velocity       -0.6471051383E+06   0.5399580930E+06
... min/max : gpix              -0.4136330275E+18   0.4136302065E+18
... min/max : gpiy              -0.6662188573E+18   0.6716904277E+18
 
STEP =      3  TIME = 9.7122375398E-06  DT = 3.4500970620E-07
 
 MEMORY STATS AT START OF TIMESTEP            4
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.2574807157E+06   0.5399580930E+06
... level / min : old vels    1  -0.2574884764E+06  -0.6471051383E+06
 
 Call to estdt at beginning of step           4
 gives dt =  5.408795851761535E-007
 dt_growth factor limits the new dt =  3.795106768193556E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  3.795106768193556E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   5.605293977081965E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.2532234880E+06   0.3456253157E+06
... level / min : MAC vels    1  -0.2532331666E+06  -0.4090843340E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.11084767E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.37958518E+02
... level / max : MAC vels    1  -0.2439156876E+06  -0.3259135948E+06
... level / min : MAC vels    1   0.2439188065E+06   0.3539906316E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000016E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1486771803E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000136E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000084E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.3818281169E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000181E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800007387E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000750E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001728E-28
... new min/max : neon-20            0.0000000000E+00   0.2750964690E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000015E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968696E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159716E+21   0.1397900306E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   5.517432226550017E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.2532234599E+06   0.3456257916E+06
... level / min : MAC vels    1  -0.2532331385E+06  -0.4090839786E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.11084439E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.36899112E+02
... level / max : MAC vels    1  -0.2439000500E+06  -0.3259362955E+06
... level / min : MAC vels    1   0.2439031730E+06   0.3539879107E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000016E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1486771803E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000136E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000084E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.3818281176E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000181E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800007387E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000750E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001728E-28
... new min/max : neon-20            0.0000000000E+00   0.2750964707E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000015E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968696E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159717E+21   0.1397900306E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2355773242E+06   0.2355668023E+06
... new min/max : y-velocity       -0.2621497532E+06   0.2748203472E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2355773242E+06   0.2355668023E+06
... y-velocity before projection -0.2621497532E+06   0.2748203472E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1111502977E+06   0.1111441545E+06
... y-velocity  after projection -0.1864074825E+06   0.1885171068E+06
 
 Timing summary:
    Advection       :    4.61441588401794       seconds
    MAC   Projection:   0.915270805358887       seconds
    Nodal Projection:    2.45639586448669       seconds
    Thermal         :    2.40684604644775       seconds
    Reactions       :    536.205208063126       seconds
    Misc            :   0.944188117980957       seconds
  
 Time to advance timestep:    547.089004993439       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            4
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.1111502977E+06   0.1111441545E+06
... min/max : y-velocity       -0.1864074825E+06   0.1885171068E+06
... min/max : gpix              -0.5920012204E+18   0.5920137106E+18
... min/max : gpiy              -0.4896913834E+18   0.4851301605E+18
 
STEP =      4  TIME = 1.0091748217E-05  DT = 3.7951067682E-07
 
 MEMORY STATS AT START OF TIMESTEP            5
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.1111441545E+06   0.1885171068E+06
... level / min : old vels    1  -0.1111502977E+06  -0.1864074825E+06
 
 Call to estdt at beginning of step           5
 gives dt =  6.774348313459897E-007
 dt_growth factor limits the new dt =  4.174617445012912E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  4.174617445012912E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   6.069350248781801E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.1118837412E+06   0.4601873699E+06
... level / min : MAC vels    1  -0.1118897978E+06  -0.5961233403E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.11942092E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.40229554E+02
... level / max : MAC vels    1  -0.1319677968E+06  -0.5696579658E+06
... level / min : MAC vels    1   0.1319675058E+06   0.4072030297E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000017E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1494285800E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000142E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000088E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.3971875332E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000188E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800007684E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000780E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001798E-28
... new min/max : neon-20            0.0000000000E+00   0.2976671356E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000015E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968711E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159698E+21   0.1397900306E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.085892265595775E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.1118837340E+06   0.4601873214E+06
... level / min : MAC vels    1  -0.1118897906E+06  -0.5961234882E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.11941835E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.30023661E+02
... level / max : MAC vels    1  -0.1319813267E+06  -0.5696586553E+06
... level / min : MAC vels    1   0.1319810370E+06   0.4072304951E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000017E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1494285800E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000142E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000088E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.3971875329E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000188E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800007684E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000780E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001798E-28
... new min/max : neon-20            0.0000000000E+00   0.2976671347E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000015E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968711E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159698E+21   0.1397900306E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.1268119698E+06   0.1268196227E+06
... new min/max : y-velocity       -0.7051486836E+06   0.6190282701E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.1268119698E+06   0.1268196227E+06
... y-velocity before projection -0.7051486836E+06   0.6190282701E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1122753510E+06   0.1122790183E+06
... y-velocity  after projection -0.7387124645E+06   0.5700096891E+06
 
 Timing summary:
    Advection       :    3.47820401191711       seconds
    MAC   Projection:   0.921699285507202       seconds
    Nodal Projection:    2.29441189765930       seconds
    Thermal         :    2.46407103538513       seconds
    Reactions       :    531.878621101379       seconds
    Misc            :   0.941529750823975       seconds
  
 Time to advance timestep:    541.772426128387       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            5
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.1122753510E+06   0.1122790183E+06
... min/max : y-velocity       -0.7387124645E+06   0.5700096891E+06
... min/max : gpix              -0.7459885441E+18   0.7459778188E+18
... min/max : gpiy              -0.1030481838E+19   0.8048693688E+18
 
STEP =      5  TIME = 1.0509209961E-05  DT = 4.1746174450E-07
 
 MEMORY STATS AT START OF TIMESTEP            6
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.1122790183E+06   0.5700096891E+06
... level / min : old vels    1  -0.1122753510E+06  -0.7387124645E+06
 
 Call to estdt at beginning of step           6
 gives dt =  4.738041878230316E-007
 dt_growth factor limits the new dt =  4.592079189514203E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  4.592079189514203E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   6.694395811572588E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.1864830220E+06   0.7346373152E+06
... level / min : MAC vels    1  -0.1864870486E+06  -0.8977536226E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.17180966E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.37265713E+02
... level / max : MAC vels    1  -0.3174268056E+06  -0.8807730682E+06
... level / min : MAC vels    1   0.3174220648E+06   0.6650317576E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000017E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1502545271E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000148E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000091E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4140720458E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000196E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800008011E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000813E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001874E-28
... new min/max : neon-20            0.0000000000E+00   0.3234965224E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000016E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968727E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159678E+21   0.1397900306E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   6.732042029354752E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.1864830220E+06   0.7346372750E+06
... level / min : MAC vels    1  -0.1864870486E+06  -0.8977536943E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.17180827E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.47367683E+02
... level / max : MAC vels    1  -0.3174459919E+06  -0.8808412611E+06
... level / min : MAC vels    1   0.3174412280E+06   0.6650059037E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000017E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1502545274E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000148E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000091E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4140720516E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000196E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800008011E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000813E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001874E-28
... new min/max : neon-20            0.0000000000E+00   0.3234965211E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000016E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968727E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159678E+21   0.1397900306E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2874853452E+06   0.2874796369E+06
... new min/max : y-velocity       -0.4389310695E+06   0.4015018393E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2874853452E+06   0.2874796369E+06
... y-velocity before projection -0.4389310695E+06   0.4015018393E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1639866827E+06   0.1639798028E+06
... y-velocity  after projection -0.4906126675E+06   0.4084056028E+06
 
 Timing summary:
    Advection       :    3.72413253784180       seconds
    MAC   Projection:   0.954741954803467       seconds
    Nodal Projection:    3.04112505912781       seconds
    Thermal         :    2.37955498695374       seconds
    Reactions       :    531.701864957809       seconds
    Misc            :   0.937661886215210       seconds
  
 Time to advance timestep:    542.382467031479       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            6
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.1639866827E+06   0.1639798028E+06
... min/max : y-velocity       -0.4906126675E+06   0.4084056028E+06
... min/max : gpix              -0.5139789557E+18   0.5139904392E+18
... min/max : gpiy              -0.1262062031E+19   0.8627404678E+18
 
STEP =      6  TIME = 1.0968417880E-05  DT = 4.5920791895E-07
 
 MEMORY STATS AT START OF TIMESTEP            7
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.1639798028E+06   0.4084056028E+06
... level / min : old vels    1  -0.1639866827E+06  -0.4906126675E+06
 
 Call to estdt at beginning of step           7
 gives dt =  7.134094184993702E-007
 dt_growth factor limits the new dt =  5.051287108465624E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  5.051287108465624E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   7.405279520544997E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.2026872009E+06   0.3518428299E+06
... level / min : MAC vels    1  -0.2026935515E+06  -0.4253066926E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.12947011E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.28450284E+02
... level / max : MAC vels    1  -0.1847290975E+06  -0.4384206268E+06
... level / min : MAC vels    1   0.1847255902E+06   0.3532577386E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000018E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1511602498E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000154E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000095E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4325926624E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000205E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800008371E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000849E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001958E-28
... new min/max : neon-20            0.0000000000E+00   0.3531006187E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000017E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968745E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159655E+21   0.1397900307E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   7.519474621574557E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.2026871789E+06   0.3518424366E+06
... level / min : MAC vels    1  -0.2026935295E+06  -0.4253070856E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.12946948E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.12293819E+02
... level / max : MAC vels    1  -0.1847555794E+06  -0.4384654147E+06
... level / min : MAC vels    1   0.1847520799E+06   0.3532672558E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000018E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1511602508E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000154E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000095E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4325926815E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000205E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800008371E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000849E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900001958E-28
... new min/max : neon-20            0.0000000000E+00   0.3531006155E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000017E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968745E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159655E+21   0.1397900307E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2282078147E+06   0.2282007224E+06
... new min/max : y-velocity       -0.3753947184E+06   0.5389419351E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2282078147E+06   0.2282007224E+06
... y-velocity before projection -0.3753947184E+06   0.5389419351E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1382915421E+06   0.1382822988E+06
... y-velocity  after projection -0.3309010022E+06   0.4072975576E+06
 
 Timing summary:
    Advection       :    4.08501195907593       seconds
    MAC   Projection:   0.944501161575317       seconds
    Nodal Projection:    3.16658806800842       seconds
    Thermal         :    2.43414497375488       seconds
    Reactions       :    534.624751091003       seconds
    Misc            :   0.933863162994385       seconds
  
 Time to advance timestep:    545.631070852280       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            7
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.1382915421E+06   0.1382822988E+06
... min/max : y-velocity       -0.3309010022E+06   0.4072975576E+06
... min/max : gpix              -0.4269110517E+18   0.4269050527E+18
... min/max : gpiy              -0.1086467078E+19   0.6398789565E+18
 
STEP =      7  TIME = 1.1473546591E-05  DT = 5.0512871085E-07
 
 MEMORY STATS AT START OF TIMESTEP            8
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.1382822988E+06   0.4072975576E+06
... level / min : old vels    1  -0.1382915421E+06  -0.3309010022E+06
 
 Call to estdt at beginning of step           8
 gives dt =  7.561128465690018E-007
 dt_growth factor limits the new dt =  5.556415819312187E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  5.556415819312187E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   8.271517539220100E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.1795850847E+06   0.8356799270E+06
... level / min : MAC vels    1  -0.1795926981E+06  -0.6641697099E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13157622E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.45263794E+02
... level / max : MAC vels    1  -0.2342523087E+06  -0.5949307000E+06
... level / min : MAC vels    1   0.2342385093E+06   0.7507201589E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000019E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1521564366E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000162E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000100E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4529525842E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000215E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800008767E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000890E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900002050E-28
... new min/max : neon-20            0.0000000000E+00   0.3871349795E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000017E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968764E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159630E+21   0.1397900307E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   8.280968970963019E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.1795850882E+06   0.8356798823E+06
... level / min : MAC vels    1  -0.1795927016E+06  -0.6641697395E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.13157865E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.47154193E+02
... level / max : MAC vels    1  -0.2342715058E+06  -0.5948967406E+06
... level / min : MAC vels    1   0.2342577132E+06   0.7508566642E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000019E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1521564366E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000162E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000100E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4529525842E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000215E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800008767E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000890E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900002050E-28
... new min/max : neon-20            0.0000000000E+00   0.3871349824E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000017E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968764E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159630E+21   0.1397900307E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.1789520374E+06   0.1789406618E+06
... new min/max : y-velocity       -0.4360865508E+06   0.5586961769E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.1789520374E+06   0.1789406618E+06
... y-velocity before projection -0.4360865508E+06   0.5586961769E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1236483948E+06   0.1236518820E+06
... y-velocity  after projection -0.3419063242E+06   0.4591243752E+06
 
 Timing summary:
    Advection       :    3.48886275291443       seconds
    MAC   Projection:    1.21695303916931       seconds
    Nodal Projection:    2.98878002166748       seconds
    Thermal         :    1.93809103965759       seconds
    Reactions       :    527.425149202347       seconds
    Misc            :   0.937372207641602       seconds
  
 Time to advance timestep:    537.918967962265       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            8
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.1236483948E+06   0.1236518820E+06
... min/max : y-velocity       -0.3419063242E+06   0.4591243752E+06
... min/max : gpix              -0.4403923141E+18   0.4403700656E+18
... min/max : gpiy              -0.5376835741E+18   0.5611187024E+18
 
STEP =      8  TIME = 1.2029188173E-05  DT = 5.5564158193E-07
 
 MEMORY STATS AT START OF TIMESTEP            9
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.1236518820E+06   0.4591243752E+06
... level / min : old vels    1  -0.1236483948E+06  -0.3419063242E+06
 
 Call to estdt at beginning of step           9
 gives dt =  7.622972108436919E-007
 dt_growth factor limits the new dt =  6.112057401243406E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  6.112057401243406E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   9.109161640455507E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.2025598359E+06   0.5139966032E+06
... level / min : MAC vels    1  -0.2025529146E+06  -0.2761248282E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.14706169E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.73449199E+02
... level / max : MAC vels    1  -0.2032236805E+06  -0.2264208698E+06
... level / min : MAC vels    1   0.2032218895E+06   0.5162165580E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000020E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1532523018E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000170E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000105E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4753514521E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000225E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800009202E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000934E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900002152E-28
... new min/max : neon-20            0.0000000000E+00   0.4263473616E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000018E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968786E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159602E+21   0.1397900308E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   8.970755219858907E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.2025598438E+06   0.5139969959E+06
... level / min : MAC vels    1  -0.2025529225E+06  -0.2761246833E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.14706471E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.10379127E+03
... level / max : MAC vels    1  -0.2032081396E+06  -0.2263938104E+06
... level / min : MAC vels    1   0.2032063473E+06   0.5161714638E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000020E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1532523006E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000170E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000105E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4753514308E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000225E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800009202E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000934E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900002152E-28
... new min/max : neon-20            0.0000000000E+00   0.4263473447E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000018E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968786E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159602E+21   0.1397900308E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2438163417E+06   0.2438295527E+06
... new min/max : y-velocity       -0.4930187588E+06   0.3772779276E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2438163417E+06   0.2438295527E+06
... y-velocity before projection -0.4930187588E+06   0.3772779276E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1170112039E+06   0.1169894519E+06
... y-velocity  after projection -0.3083044236E+06   0.4294971437E+06
 
 Timing summary:
    Advection       :    3.66676592826843       seconds
    MAC   Projection:   0.936353206634521       seconds
    Nodal Projection:    2.92580294609070       seconds
    Thermal         :    2.02882480621338       seconds
    Reactions       :    534.097352266312       seconds
    Misc            :   0.940403938293457       seconds
  
 Time to advance timestep:    544.479822874069       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP            9
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.1170112039E+06   0.1169894519E+06
... min/max : y-velocity       -0.3083044236E+06   0.4294971437E+06
... min/max : gpix              -0.3290616963E+18   0.3290738998E+18
... min/max : gpiy              -0.6277974652E+18   0.2582828122E+18
 
STEP =      9  TIME = 1.2640393913E-05  DT = 6.1120574012E-07
 
 MEMORY STATS AT START OF TIMESTEP           10
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
... level / max : old vels    1   0.1169894519E+06   0.4294971437E+06
... level / min : old vels    1  -0.1170112039E+06  -0.3083044236E+06
 
 Call to estdt at beginning of step          10
 gives dt =  8.149013173895082E-007
 dt_growth factor limits the new dt =  6.723263141367747E-007
 
inner sponge: topsponge_lo_r      , topsponge_hi_r      :   0.645625000000E+04    0.662075000000E+04
 level:            1
    number of boxes =           64
    maximum zones   =          512        2048
 Number of valid cells at level                   1     1048576
 Number of valid + ghost cells at level           1     1098304
 
 <<< STEP  1 : react state     
 <<< CALLING advance_timestep with dt =  6.723263141367747E-007
 <<< STEP  2 : make w0 >>> 
 ... max CFL of w0:   9.867853480570769E-005
 
 <<< STEP  3 : create MAC velocities>>> 
 
... level / max : MAC vels    1   0.1690632993E+06   0.5293526686E+06
... level / min : MAC vels    1  -0.1690870401E+06  -0.7182897202E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.12492635E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.93874877E+02
... level / max : MAC vels    1  -0.1337387820E+06  -0.6296524279E+06
... level / min : MAC vels    1   0.1337265578E+06   0.5687757862E+06
 
 <<< STEP  4 : advect base        
... Level 1 create thermal coeffs:
             :  density_advance >>> 
             :   tracer_advance >>> 
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000021E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1544441140E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000178E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000110E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4997825224E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000237E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800009681E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000982E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900002264E-28
... new min/max : neon-20            0.0000000000E+00   0.4718739234E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000019E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968809E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>> 
... Level 1 update:
... new min/max : rho * H            0.1786159572E+21   0.1397900308E+25
 
 <<< STEP  4a: thermal conduct >>>
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  5 : react state     
 <<< STEP  6 : make new S and new w0 >>> 
... Level 1 create thermal coeffs:
 ... max CFL of w0:   9.832779707320983E-005
 
 <<< STEP  7 : create MAC velocities >>> 
 
... level / max : MAC vels    1   0.1690632995E+06   0.5293526927E+06
... level / min : MAC vels    1  -0.1690870403E+06  -0.7182896722E+06
 
... before mac_projection: max of [div (coeff * UMAC) - RHS)] 0.12492491E+13
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
...  after mac_projection: max of [div (coeff * UMAC) - RHS)] 0.62478883E+02
... level / max : MAC vels    1  -0.1336800684E+06  -0.6293903887E+06
... level / min : MAC vels    1   0.1336678435E+06   0.5685366373E+06
 
 <<< STEP  8 : advect base   
             :  density_advance >>>
             :   tracer_advance >>>
... Level 1 update:
... new min/max : hydrogen-1         0.0000000000E+00   0.1000000000E-29
... new min/max : helium-4           0.5000000000E+00   0.9999999999E+00
... new min/max : carbon-12          0.9999999999E-10   0.5000000000E+00
... new min/max : carbon-13          0.0000000000E+00   0.1300000021E-28
... new min/max : nitrogen-13        0.0000000000E+00   0.1544441141E-28
... new min/max : nitrogen-14        0.0000000000E+00   0.1400000178E-28
... new min/max : nitrogen-15        0.0000000000E+00   0.1500000110E-28
... new min/max : oxygen-14          0.0000000000E+00   0.1400000000E-28
... new min/max : oxygen-15          0.0000000000E+00   0.1500000000E-28
... new min/max : oxygen-16          0.0000000000E+00   0.4997825237E-10
... new min/max : oxygen-17          0.0000000000E+00   0.1700000237E-28
... new min/max : oxygen-18          0.0000000000E+00   0.1800000002E-28
... new min/max : fluorine-17        0.0000000000E+00   0.1700000000E-28
... new min/max : fluorine-18        0.0000000000E+00   0.1800009681E-28
... new min/max : fluorine-19        0.0000000000E+00   0.1900000982E-28
... new min/max : neon-18            0.0000000000E+00   0.1800000000E-28
... new min/max : neon-19            0.0000000000E+00   0.1900002264E-28
... new min/max : neon-20            0.0000000000E+00   0.4718739245E-19
... new min/max : magnesium-22       0.0000000000E+00   0.2200000019E-28
... new min/max : sulfur-30          0.0000000000E+00   0.3000000000E-28
... new min/max : nickel-56          0.0000000000E+00   0.5600000000E-28
... new min/max : density            0.9954009627E+04   0.6965968809E+07
... new min/max : tracer             0.0000000000E+00   0.0000000000E+00
             : enthalpy_advance >>>
... Level 1 update:
... new min/max : rho * H            0.1786159572E+21   0.1397900308E+25
 
 <<< STEP  8a: thermal conduct >>>
... Level 1 create thermal coeffs:
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
 <<< STEP  9 : react state 
 <<< STEP 10 : make new S >>>
... Level 1 create thermal coeffs:
 <<< STEP 11 : update and project new velocity >>>
... Level 1 update:
... new min/max : x-velocity       -0.2020556357E+06   0.2020294644E+06
... new min/max : y-velocity       -0.3963540665E+06   0.2290486025E+06
 
 PROJ_TYPE IN HGPROJECT:           4
 PROJ_TYPE IN HGPROJECT: ST_DENSE
... x-velocity before projection -0.2020556357E+06   0.2020294644E+06
... y-velocity before projection -0.3963540665E+06   0.2290486025E+06
 
 doing hgproject with tolerance, eps =   1.000000000000000E-012
 F90mg: Coarse problem domain for bottom_solver = 4: 
BOX[{0, 0}, {7, 31}]
    ... Original boxes           64
    ... New      boxes            4
 # cells on each side             8
... x-velocity  after projection -0.1009080420E+06   0.1009026719E+06
... y-velocity  after projection -0.2971897181E+06   0.2143177760E+06
 
 Timing summary:
    Advection       :    3.91132712364197       seconds
    MAC   Projection:   0.955620765686035       seconds
    Nodal Projection:    2.90385723114014       seconds
    Thermal         :    2.37325096130371       seconds
    Reactions       :    533.840269088745       seconds
    Misc            :   0.944639682769775       seconds
  
 Time to advance timestep:    544.477264165878       seconds
 
FAB byte spread across CPUs: [  85993824 ...  128990736]
 
 MEMORY STATS AT END OF TIMESTEP           10
  
    multifab:        95841664       95841664             15             15
         fab:         2995052        4492578             30             45
    boxarray:             988            988             22             22
      layout:               0              0              1              1
    boxassoc:            1006           1840             35             35
     fgassoc:               0              0              0              0
   syncassoc:              60            108              6              6
   copyassoc:              15            178              7              7
   fluxassoc:               0              0              0              0
 
At level   1
... min/max : x-velocity       -0.1009080420E+06   0.1009026719E+06
... min/max : y-velocity       -0.2971897181E+06   0.2143177760E+06
... min/max : gpix              -0.2158392963E+18   0.2158108541E+18
... min/max : gpiy              -0.8131900109E+18   0.7626335563E+18
 
STEP =     10  TIME = 1.3312720227E-05  DT = 6.7232631414E-07
 
 MEMORY STATS AT END OF PROGRAM
  
    multifab:               0              0              0              0
         fab:               0              0              0              0
    boxarray:               0              0              0              0
      layout:               0              0              0              0
    boxassoc:               0              0              0              0
     fgassoc:               0              0              0              0
   syncassoc:               0              0              0              0
   copyassoc:               0              0              0              0
   fluxassoc:               0              0              0              0
  
 Run Time =    8026.53070092201