Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.06-1-gb83904febeb2) 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.452208244
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               = 185424520.7
MLMG: Initial residual (resid0) = 185424520.7
MLMG: Final Iter. 10 resid, resid/bnorm = 0.01339600951, 7.224508097e-11
MLMG: Timers: Solve = 1.919709922 Iter = 1.878650075 Bottom = 0.000746524
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.411180973
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               = 72482174.27
MLMG: Initial residual (resid0) = 72482174.27
MLMG: Final Iter. 2 resid, resid/bnorm = 0.5223139888, 7.20610266e-09
MLMG: Timers: Solve = 0.426172249 Iter = 0.387419151 Bottom = 0.000169421
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               = 43334804.49
MLMG: Initial residual (resid0) = 43334804.49
MLMG: Final Iter. 2 resid, resid/bnorm = 0.02298274693, 5.303530777e-10
MLMG: Timers: Solve = 0.426134022 Iter = 0.387397346 Bottom = 0.000168862
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               = 32768245.23
MLMG: Initial residual (resid0) = 32768245.23
MLMG: Final Iter. 3 resid, resid/bnorm = 0.001478955848, 4.513381286e-11
MLMG: Timers: Solve = 0.612525917 Iter = 0.573657957 Bottom = 0.000252974
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.394627162
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               = 224625388.6
MLMG: Initial residual (resid0) = 224625388.6
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02028144896, 9.02901007e-11
MLMG: Timers: Solve = 0.113231654 Iter = 0.102708212 Bottom = 0.007769666
<<< 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               = 224625389.4
MLMG: Initial residual (resid0) = 3940.279649
MLMG: Final Iter. 2 resid, resid/bnorm = 0.04907616747, 2.184800552e-10
MLMG: Timers: Solve = 0.044639172 Iter = 0.036281386 Bottom = 0.002591466
<<< 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               = 35718891.54
MLMG: Initial residual (resid0) = 35718891.54
MLMG: Final Iter. 31 resid, resid/bnorm = 0.003148811782, 8.815536111e-11
MLMG: Timers: Solve = 5.8155375 Iter = 5.77452942 Bottom = 0.00227466
Done calling nodal solver

Timestep 0 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection  :1.700868481 seconds
MAC Proj   :0.268517844 seconds
Nodal Proj :5.959825999 seconds
Reactions  :0.177251447 seconds
Misc       :0.074537687 seconds
Base State :0.042712608 seconds
Time to advance time step: 8.181448533

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.406859293
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               = 224625175.9
MLMG: Initial residual (resid0) = 224625175.9
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02050751448, 9.129659839e-11
MLMG: Timers: Solve = 0.110480054 Iter = 0.102359709 Bottom = 0.007775315
<<< 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               = 224625175.8
MLMG: Initial residual (resid0) = 0.7033640956
MLMG: Final Iter. 1 resid, resid/bnorm = 0.04446259194, 1.979412672e-10
MLMG: Timers: Solve = 0.028451294 Iter = 0.019980716 Bottom = 0.001303813
<<< 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               = 25278219.2
MLMG: Initial residual (resid0) = 25278219.2
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001846272498, 7.303807613e-11
MLMG: Timers: Solve = 1.910490754 Iter = 1.87170066 Bottom = 0.000757619
Done calling nodal solver

Timestep 1 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection  :1.779800885 seconds
MAC Proj   :0.246384446 seconds
Nodal Proj :2.05663589 seconds
Reactions  :0.176926718 seconds
Misc       :0.079359391 seconds
Base State :0.046562457 seconds
Time to advance time step: 4.339558877

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.403343436
Call to estdt for level 0 gives dt_lev = 2.657850077
Call to estdt for level 1 gives dt_lev = 0.6581035052
Minimum estdt over all levels = 0.6581035052
Call to estdt at beginning of step 2 gives dt =0.6581035052
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               = 221100139
MLMG: Initial residual (resid0) = 221100139
MLMG: Final Iter. 6 resid, resid/bnorm = 0.01983711123, 8.972003059e-11
MLMG: Timers: Solve = 0.111480135 Iter = 0.103143176 Bottom = 0.007804702
<<< 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               = 221100138.8
MLMG: Initial residual (resid0) = 0.6525653412
MLMG: Final Iter. 1 resid, resid/bnorm = 0.04140532438, 1.87269554e-10
MLMG: Timers: Solve = 0.028177741 Iter = 0.019868415 Bottom = 0.001303893
<<< 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               = 19841168.99
MLMG: Initial residual (resid0) = 19841168.99
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001567842555, 7.901966643e-11
MLMG: Timers: Solve = 1.990877923 Iter = 1.947107567 Bottom = 0.000749076
Done calling nodal solver

Timestep 2 ends with TIME = 0.0006060015283 DT = 0.000317429372
Timing summary:
Advection  :1.752489523 seconds
MAC Proj   :0.244204642 seconds
Nodal Proj :2.143015755 seconds
Reactions  :0.179662256 seconds
Misc       :0.073623552 seconds
Base State :0.04206025 seconds
Time to advance time step: 4.393614843

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.418507088
Time to regrid: 0.066148796
Call to estdt for level 0 gives dt_lev = 3.830037616
Call to estdt for level 1 gives dt_lev = 0.8027988743
Minimum estdt over all levels = 0.8027988743
Call to estdt at beginning of step 3 gives dt =0.8027988743
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               = 221012973.3
MLMG: Initial residual (resid0) = 221012973.3
MLMG: Final Iter. 5 resid, resid/bnorm = 0.216652289, 9.802695552e-10
MLMG: Timers: Solve = 0.089674764 Iter = 0.081664221 Bottom = 0.006173823
<<< 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               = 221012973.1
MLMG: Initial residual (resid0) = 0.6485630108
MLMG: Final Iter. 1 resid, resid/bnorm = 0.04097023881, 1.853748141e-10
MLMG: Timers: Solve = 0.026413586 Iter = 0.018515195 Bottom = 0.001073093
<<< 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               = 15859855.11
MLMG: Initial residual (resid0) = 15859855.11
MLMG: Final Iter. 11 resid, resid/bnorm = 0.001177113154, 7.421966629e-11
MLMG: Timers: Solve = 2.18304822 Iter = 2.139428915 Bottom = 0.000801244
Done calling nodal solver

Timestep 3 ends with TIME = 0.0009551738375 DT = 0.0003491723092
Timing summary:
Advection  :1.838974476 seconds
MAC Proj   :0.223641495 seconds
Nodal Proj :2.33565426 seconds
Reactions  :0.182858745 seconds
Misc       :0.07575526 seconds
Base State :0.051423406 seconds
Time to advance time step: 4.657498731

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.453392475

Total Time: 32.59733462
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): 9256
[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.06-1-gb83904febeb2) finalized