Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.10-38-g6f07fe4e84a4) 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.474810219
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.01302976906, 5.822250438e-11
MLMG: Timers: Solve = 1.014711627 Iter = 0.974148803 Bottom = 0.000533133
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.447599425
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.7442735399, 7.692889695e-09
MLMG: Timers: Solve = 0.433082034 Iter = 0.393660253 Bottom = 0.000227343
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.01898677269, 3.578868144e-10
MLMG: Timers: Solve = 0.433210936 Iter = 0.393752219 Bottom = 0.000233342
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.63701777 Iter = 0.597587582 Bottom = 0.000352202
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.432137167
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.02237911895, 1.035810292e-10
MLMG: Timers: Solve = 0.079654698 Iter = 0.071966368 Bottom = 0.007397376
<<< 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.17408
MLMG: Final Iter. 3 resid, resid/bnorm = 0.004752751176, 2.199795534e-11
MLMG: Timers: Solve = 0.042553464 Iter = 0.036724505 Bottom = 0.003652248
<<< 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.91
MLMG: Initial residual (resid0) = 35718898.91
MLMG: Final Iter. 31 resid, resid/bnorm = 0.002826683223, 7.913690818e-11
MLMG: Timers: Solve = 5.882937128 Iter = 5.842013664 Bottom = 0.002193071
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.798905536 seconds
MAC Proj :0.218108534 seconds
Nodal Proj :6.044961935 seconds
Reactions :0.167775378 seconds
Misc :0.103281757 seconds
Base State :0.045104615 seconds
Time to advance time step: 8.33350392
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.431963298
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.02259196714, 1.045660403e-10
MLMG: Timers: Solve = 0.077282725 Iter = 0.071477578 Bottom = 0.007373305
<<< 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.35175923
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02233633684, 1.033828655e-10
MLMG: Timers: Solve = 0.019213013 Iter = 0.013472818 Bottom = 0.001231337
<<< 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.001655275002, 7.402134071e-11
MLMG: Timers: Solve = 1.938326085 Iter = 1.898887232 Bottom = 0.000875756
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.757027451 seconds
MAC Proj :0.185915754 seconds
Nodal Proj :2.097704116 seconds
Reactions :0.165420836 seconds
Misc :0.118257307 seconds
Base State :0.044283478 seconds
Time to advance time step: 4.324815569
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.480095709
Call to estdt for level 0 gives dt_lev = 2.560503025
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.02230187505, 1.029864128e-10
MLMG: Timers: Solve = 0.077200954 Iter = 0.071373993 Bottom = 0.007377629
<<< 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.3069432653
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01958142943, 9.042383976e-11
MLMG: Timers: Solve = 0.019289176 Iter = 0.013556805 Bottom = 0.001227273
<<< 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.001410154626, 8.655861473e-11
MLMG: Timers: Solve = 1.937162581 Iter = 1.897768508 Bottom = 0.000889283
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060015283 DT = 0.000317429372
Timing summary:
Advection :1.794965009 seconds
MAC Proj :0.187770053 seconds
Nodal Proj :2.099444102 seconds
Reactions :0.16731584 seconds
Misc :0.119142156 seconds
Base State :0.042047574 seconds
Time to advance time step: 4.369133005
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.43817257
Time to regrid: 0.121417582
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.02218363807, 1.024302059e-10
MLMG: Timers: Solve = 0.077381356 Iter = 0.071431472 Bottom = 0.007357
<<< 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.2695824229
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01713669178, 7.912655553e-11
MLMG: Timers: Solve = 0.01935583 Iter = 0.013589357 Bottom = 0.001227446
<<< 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.001057907939, 8.060845714e-11
MLMG: Timers: Solve = 2.124509712 Iter = 2.085092379 Bottom = 0.000937597
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551738375 DT = 0.0003491723092
Timing summary:
Advection :1.80220171 seconds
MAC Proj :0.185910217 seconds
Nodal Proj :2.288458039 seconds
Reactions :0.168048065 seconds
Misc :0.119550367 seconds
Base State :0.045386384 seconds
Time to advance time step: 4.564620331
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.438614671
Total Time: 31.97098105
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): 9248
[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.10-38-g6f07fe4e84a4) finalized