Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.08-70-ge33b6cd8b866) initialized
Calling Setup()
Calling ReadParameters()
 WARNING: burning_cutoff_density_lo not supplied in the inputs file
 WARNING: setting burning_cutoff_density_lo = base_cutoff_density
reading extern runtime parameters ...
Calling VariableSetup()
 
 Initializing Helmholtz EOS and using Coulomb corrections.
 
 Calling set_method_params()
Calling BCSetup()
Calling BaseStateGeometry::Init()
 Calling init_base_state_geometry()
Calling Init()
Calling InitData()
initdata model_File = kepler_new_6.25e8.hybrid.hse.320
model file = kepler_new_6.25e8.hybrid.hse.320


reading initial model
320 points found in the initial model
6 variables found in the initial model
model file mapping (spherical base state)
dr of MAESTRO base state =                            781250.000000
dr of input file data =                               1562500.000000

maximum radius (cell-centered) of input model =       499218750.000000
 
setting r_cutoff to 239
radius at r_cutoff 187109375

Maximum HSE Error = 0.002311
(after putting initial model into base state arrays, and
for density < base_cutoff_density)


Writing plotfile wdconvect-amr_pltInitData after InitData
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.427142149
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Doing initial projection
Calling nodal solver
MLMG: Initial rhs               = 223792658.9
MLMG: Initial residual (resid0) = 223792658.9
MLMG: Final Iter. 5 resid, resid/bnorm = 0.01302978396, 5.822257097e-11
MLMG: Timers: Solve = 1.018975889 Iter = 0.977678494 Bottom = 0.000383764
Done calling nodal solver

Writing plotfile wdconvect-amr_pltafter_InitProj after InitProj
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.40451264
Call to firstdt for level 0 gives dt_lev = 0.005789678157
Multiplying dt_lev by init_shrink; dt_lev = 0.0005789678157
Call to firstdt for level 1 gives dt_lev = 0.002885721563
Multiplying dt_lev by init_shrink; dt_lev = 0.0002885721563
Minimum firstdt over all levels = 0.0002885721563
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs               = 96748240.18
MLMG: Initial residual (resid0) = 96748240.18
MLMG: Final Iter. 2 resid, resid/bnorm = 0.7442735396, 7.692889692e-09
MLMG: Timers: Solve = 0.434765253 Iter = 0.394292557 Bottom = 0.00014462
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 5.916319527
Call to estdt for level 1 gives dt_lev = 2.341804478
Minimum estdt over all levels = 2.341804478
Call to estdt at end of istep_divu_iter = 1 gives dt = 2.341804478
Multiplying dt by init_shrink; dt = 0.2341804478
Ignoring this new dt since it's larger than the previous dt = 0.0002885721563
Doing initial divu iteration #2
Calling nodal solver
MLMG: Initial rhs               = 53052451.02
MLMG: Initial residual (resid0) = 53052451.02
MLMG: Final Iter. 2 resid, resid/bnorm = 0.0189867727, 3.578868147e-10
MLMG: Timers: Solve = 0.433354286 Iter = 0.393435442 Bottom = 0.000150627
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 5.916319527
Call to estdt for level 1 gives dt_lev = 2.341804478
Minimum estdt over all levels = 2.341804478
Call to estdt at end of istep_divu_iter = 2 gives dt = 2.341804478
Multiplying dt by init_shrink; dt = 0.2341804478
Ignoring this new dt since it's larger than the previous dt = 0.0002885721563
Doing initial divu iteration #3
Calling nodal solver
MLMG: Initial rhs               = 33080059.37
MLMG: Initial residual (resid0) = 33080059.37
MLMG: Final Iter. 3 resid, resid/bnorm = 0.001913452566, 5.784308136e-11
MLMG: Timers: Solve = 0.637891313 Iter = 0.597967101 Bottom = 0.000218517
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 5.916319527
Call to estdt for level 1 gives dt_lev = 2.341804478
Minimum estdt over all levels = 2.341804478
Call to estdt at end of istep_divu_iter = 3 gives dt = 2.341804478
Multiplying dt by init_shrink; dt = 0.2341804478
Ignoring this new dt since it's larger than the previous dt = 0.0002885721563

Writing plotfile wdconvect-amr_pltafter_DivuIter after final DivuIter
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.385222664
Doing initial pressure iteration #1

Timestep 0 starts with TIME = 0 DT = 0.0002885721563

Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs               = 216054224.6
MLMG: Initial residual (resid0) = 216054224.6
MLMG: Final Iter. 6 resid, resid/bnorm = 0.0223791115, 1.035809947e-10
MLMG: Timers: Solve = 0.104932926 Iter = 0.095445097 Bottom = 0.007428903
<<< STEP 4 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs               = 216054224.3
MLMG: Initial residual (resid0) = 10665.17397
MLMG: Final Iter. 3 resid, resid/bnorm = 0.004752761585, 2.199800351e-11
MLMG: Timers: Solve = 0.056696348 Iter = 0.049096038 Bottom = 0.003714536
<<< STEP 8 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs               = 35718898.93
MLMG: Initial residual (resid0) = 35718898.93
MLMG: Final Iter. 31 resid, resid/bnorm = 0.002826744691, 7.9138629e-11
MLMG: Timers: Solve = 5.908755745 Iter = 5.867579331 Bottom = 0.002164799
Done calling nodal solver

Timestep 0 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection  :1.744925381 seconds
MAC Proj   :0.268178538 seconds
Nodal Proj :6.055365367 seconds
Reactions  :0.175173043 seconds
Misc       :0.073854331 seconds
Base State :0.047063828 seconds
Time to advance time step: 8.317987944

Writing plotfile 0 after all initialization
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.395260948
Calling Evolve()

Timestep 1 starts with TIME = 0 DT = 0.0002885721563

Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs               = 216054534.3
MLMG: Initial residual (resid0) = 216054534.3
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02259200439, 1.045662127e-10
MLMG: Timers: Solve = 0.103816732 Iter = 0.095956311 Bottom = 0.007424479
<<< STEP 4 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs               = 216054534.1
MLMG: Initial residual (resid0) = 0.3517592185
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02233633444, 1.033828544e-10
MLMG: Timers: Solve = 0.026228652 Iter = 0.018570556 Bottom = 0.001244048
<<< STEP 8 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs               = 22362132.14
MLMG: Initial residual (resid0) = 22362132.14
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001655273139, 7.402125742e-11
MLMG: Timers: Solve = 1.942935809 Iter = 1.903137714 Bottom = 0.000720377
Done calling nodal solver

Timestep 1 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection  :1.777736348 seconds
MAC Proj   :0.233318823 seconds
Nodal Proj :2.09782004 seconds
Reactions  :0.187440094 seconds
Misc       :0.082881976 seconds
Base State :0.044919086 seconds
Time to advance time step: 4.379706328

Writing plotfile 1
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.415580165
Call to estdt for level 0 gives dt_lev = 2.560503026
Call to estdt for level 1 gives dt_lev = 0.5423221224
Minimum estdt over all levels = 0.5423221224
Call to estdt at beginning of step 2 gives dt =0.5423221224
dt_growth factor limits the new dt = 0.000317429372

Timestep 2 starts with TIME = 0.0002885721563 DT = 0.000317429372

Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs               = 216551625.1
MLMG: Initial residual (resid0) = 216551625.1
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02230188623, 1.029864644e-10
MLMG: Timers: Solve = 0.105205903 Iter = 0.096258724 Bottom = 0.007506673
<<< STEP 4 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs               = 216551624.9
MLMG: Initial residual (resid0) = 0.3069432486
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01958142804, 9.042383336e-11
MLMG: Timers: Solve = 0.026437167 Iter = 0.018583892 Bottom = 0.001276569
<<< STEP 8 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs               = 16291326.17
MLMG: Initial residual (resid0) = 16291326.17
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001410156488, 8.655872907e-11
MLMG: Timers: Solve = 1.93949286 Iter = 1.899818827 Bottom = 0.00074034
Done calling nodal solver

Timestep 2 ends with TIME = 0.0006060015283 DT = 0.000317429372
Timing summary:
Advection  :1.659178149 seconds
MAC Proj   :0.233323216 seconds
Nodal Proj :2.079937557 seconds
Reactions  :0.172398008 seconds
Misc       :0.073690434 seconds
Base State :0.048763931 seconds
Time to advance time step: 4.218975408

Writing plotfile 2
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.387406478
Time to regrid: 0.057220254
Call to estdt for level 0 gives dt_lev = 3.52349079
Call to estdt for level 1 gives dt_lev = 0.6591885313
Minimum estdt over all levels = 0.6591885313
Call to estdt at beginning of step 3 gives dt =0.6591885313
dt_growth factor limits the new dt = 0.0003491723092

Timestep 3 starts with TIME = 0.0006060015283 DT = 0.0003491723092

Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs               = 216573205.8
MLMG: Initial residual (resid0) = 216573205.8
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02218361199, 1.024300855e-10
MLMG: Timers: Solve = 0.10351307 Iter = 0.095796804 Bottom = 0.007416732
<<< STEP 4 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 4a: thermal conduct >>>
<<< STEP 5 : react state >>>
<<< STEP 6 : make new S and new w0 >>>
<<< STEP 7 : create MAC velocities >>>
MLMG: Initial rhs               = 216573205.6
MLMG: Initial residual (resid0) = 0.2695824104
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01713669203, 7.912655668e-11
MLMG: Timers: Solve = 0.02603994 Iter = 0.0183689 Bottom = 0.001255957
<<< STEP 8 : advect base >>>
            :  density_advance >>>
            :   tracer_advance >>>
            : enthalpy_advance >>>
<<< STEP 8a: thermal conduct >>>
<<< STEP 9 : react state >>>
<<< STEP 10: make new S >>>
<<< STEP 11: update and project new velocity >>>
Calling nodal solver
MLMG: Initial rhs               = 13124031.6
MLMG: Initial residual (resid0) = 13124031.6
MLMG: Final Iter. 11 resid, resid/bnorm = 0.001057909802, 8.060859907e-11
MLMG: Timers: Solve = 2.13885483 Iter = 2.09900251 Bottom = 0.000784759
Done calling nodal solver

Timestep 3 ends with TIME = 0.0009551738375 DT = 0.0003491723092
Timing summary:
Advection  :1.67307331 seconds
MAC Proj   :0.229786817 seconds
Nodal Proj :2.278789877 seconds
Reactions  :0.172228704 seconds
Misc       :0.070765917 seconds
Base State :0.045795915 seconds
Time to advance time step: 4.425123694

Writing plotfile 3
inner sponge: r_sp      , r_tp      : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
Time to write plotfile: 0.387089469

Total Time: 31.5126677
Unused ParmParse Variables:
  [TOP]::maestro.v(nvals = 1)  :: [1]
  [TOP]::amr.check_file(nvals = 1)  :: [wdconvect-amr_chk]
  [TOP]::amr.checkpoint_files_output(nvals = 1)  :: [0]

Total GPU global memory (MB): 12066
Free  GPU global memory (MB): 9396
[The         Arena] space allocated (MB): 9049
[The         Arena] space used      (MB): 0
[The  Device Arena] space allocated (MB): 8
[The  Device Arena] space used      (MB): 0
[The Managed Arena] space allocated (MB): 8
[The Managed Arena] space used      (MB): 0
[The  Pinned Arena] space allocated (MB): 8
[The  Pinned Arena] space used      (MB): 0
AMReX (20.08-70-ge33b6cd8b866) finalized