Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.12-30-g91a48a24777d) 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.485892879
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.016778938 Iter = 0.974528065 Bottom = 0.000408909
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.454681354
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.7442729427, 7.692883522e-09
MLMG: Timers: Solve = 0.429106364 Iter = 0.389062484 Bottom = 0.00016382
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.01898695052, 3.578901663e-10
MLMG: Timers: Solve = 0.429218089 Iter = 0.388989792 Bottom = 0.000163178
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.001913454252, 5.784313233e-11
MLMG: Timers: Solve = 0.632969772 Iter = 0.592596589 Bottom = 0.000250834
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.44017077
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.0223791264, 1.035810637e-10
MLMG: Timers: Solve = 0.080280407 Iter = 0.072467707 Bottom = 0.007559258
<<< 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.004752771848, 2.199805102e-11
MLMG: Timers: Solve = 0.043056735 Iter = 0.037132902 Bottom = 0.003725676
<<< 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.002826692536, 7.913716888e-11
MLMG: Timers: Solve = 5.82863078 Iter = 5.785807821 Bottom = 0.002343559
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.884985103 seconds
MAC Proj :0.220898392 seconds
Nodal Proj :5.998308167 seconds
Reactions :0.169645225 seconds
Misc :0.107183897 seconds
Base State :0.047035031 seconds
Time to advance time step: 8.381494928
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.440549341
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.02259199321, 1.04566161e-10
MLMG: Timers: Solve = 0.078254043 Iter = 0.072092748 Bottom = 0.007469498
<<< 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.3517590652
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02233632537, 1.033828124e-10
MLMG: Timers: Solve = 0.019746349 Iter = 0.013720237 Bottom = 0.001260053
<<< 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.001655260101, 7.402067436e-11
MLMG: Timers: Solve = 1.91620998 Iter = 1.876114676 Bottom = 0.000829664
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.8612744 seconds
MAC Proj :0.189990942 seconds
Nodal Proj :2.081321644 seconds
Reactions :0.174187738 seconds
Misc :0.12216541 seconds
Base State :0.048609583 seconds
Time to advance time step: 4.42938389
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.449694944
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.02230186388, 1.029863612e-10
MLMG: Timers: Solve = 0.077676691 Iter = 0.071881931 Bottom = 0.007457874
<<< 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.306943228
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01958141199, 9.042375924e-11
MLMG: Timers: Solve = 0.019571974 Iter = 0.013656404 Bottom = 0.001230447
<<< 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.0014101509, 8.655838607e-11
MLMG: Timers: Solve = 1.91780447 Iter = 1.877615935 Bottom = 0.000827604
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060015283 DT = 0.000317429372
Timing summary:
Advection :1.849902612 seconds
MAC Proj :0.189624362 seconds
Nodal Proj :2.086204843 seconds
Reactions :0.171978107 seconds
Misc :0.122078153 seconds
Base State :0.044505531 seconds
Time to advance time step: 4.420246334
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.448764472
Time to regrid: 0.125735542
Call to estdt for level 0 gives dt_lev = 3.523490799
Call to estdt for level 1 gives dt_lev = 0.6591885311
Minimum estdt over all levels = 0.6591885311
Call to estdt at beginning of step 3 gives dt =0.6591885311
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.02218361571, 1.024301027e-10
MLMG: Timers: Solve = 0.078206203 Iter = 0.072202991 Bottom = 0.007462271
<<< 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.2695823647
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01713669087, 7.912655135e-11
MLMG: Timers: Solve = 0.019619088 Iter = 0.013834603 Bottom = 0.001239571
<<< 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.001057911664, 8.060874099e-11
MLMG: Timers: Solve = 2.112930524 Iter = 2.072862157 Bottom = 0.000896918
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551738375 DT = 0.0003491723092
Timing summary:
Advection :1.861960227 seconds
MAC Proj :0.193096088 seconds
Nodal Proj :2.273202201 seconds
Reactions :0.172608749 seconds
Misc :0.124647863 seconds
Base State :0.049082488 seconds
Time to advance time step: 4.625966163
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.446312459
Total Time: 32.3001638
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): 9249
[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.12-30-g91a48a24777d) finalized