Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.01-220-g45d4129ee9d1) 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 ...
ERROR: problem in the namelist
Calling VariableSetup()
Calling set_method_params()
Initializing Helmholtz EOS and using Coulomb corrections.
Calling BCSetup()
Calling init_base_state_geometry()
Calling Init()
Calling InitData()
------------------------------------------------------------------------------
reading initial model
320 points found in the initial model file
6 variables found in the initial model file
------------------------------------------------------------------------------
model file mapping (spherical base state)
dr of MAESTRO base state = 781250.0000
dr of input file data = 1562500.000
maximum radius (cell-centered) of input model = 499218750.0
setting r_cutoff to 239
radius at r_cutoff 187109375.0000000
Maximum HSE Error = 0.2310976572E-02
(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 : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Time to write plotfile: 0.581728261
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Doing initial projection
Calling nodal solver
MLMG: Initial rhs = 188240648.4
MLMG: Initial residual (resid0) = 188240648.4
MLMG: Final Iter. 7 resid, resid/bnorm = 0.01549260831, 8.230214057e-11
MLMG: Timers: Solve = 0.109664933 Iter = 0.099771111 Bottom = 0.002467841
Done calling nodal solver
Writing plotfile wdconvect-amr_pltafter_InitProj after InitProj
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Time to write plotfile: 0.574123213
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.0005789678157
Multiplying dt_lev by init_shrink; dt_lev = 5.789678157e-05
Minimum firstdt over all levels = 5.789678157e-05
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 80008358.56
MLMG: Initial residual (resid0) = 80008358.56
MLMG: Final Iter. 2 resid, resid/bnorm = 1.348924273, 1.685979186e-08
MLMG: Timers: Solve = 0.037820573 Iter = 0.029512218 Bottom = 0.000715409
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 = 5.789678157e-05
Doing initial divu iteration #2
Calling nodal solver
MLMG: Initial rhs = 44170164.91
MLMG: Initial residual (resid0) = 44170164.91
MLMG: Final Iter. 2 resid, resid/bnorm = 0.2234185168, 5.058131824e-09
MLMG: Timers: Solve = 0.033908572 Iter = 0.029423065 Bottom = 0.000705364
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 = 5.789678157e-05
Doing initial divu iteration #3
Calling nodal solver
MLMG: Initial rhs = 26504865.12
MLMG: Initial residual (resid0) = 26504865.12
MLMG: Final Iter. 3 resid, resid/bnorm = 0.002006638621, 7.570831286e-11
MLMG: Timers: Solve = 0.0476525 Iter = 0.043157229 Bottom = 0.001061197
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 = 5.789678157e-05
Writing plotfile wdconvect-amr_pltafter_DivuIter after final DivuIter
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Time to write plotfile: 0.557427383
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 5.789678157e-05
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 173628411.4
MLMG: Initial residual (resid0) = 173628411.4
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02233195887, 1.286192662e-10
MLMG: Timers: Solve = 0.107573474 Iter = 0.096148402 Bottom = 0.007035059
<<< 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 = 173628411.8
MLMG: Initial residual (resid0) = 2158.124791
MLMG: Final Iter. 2 resid, resid/bnorm = 0.04143702223, 2.386534658e-10
MLMG: Timers: Solve = 0.042406384 Iter = 0.033807704 Bottom = 0.002314365
<<< 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 = 35718515.72
MLMG: Initial residual (resid0) = 35718515.72
MLMG: Final Iter. 43 resid, resid/bnorm = 0.002904127934, 8.130595227e-11
MLMG: Timers: Solve = 0.604466328 Iter = 0.599927996 Bottom = 0.015753299
Done calling nodal solver
Timestep 0 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05
Timing summary:
Advection :2.335868719 seconds
MAC Proj :0.284520927 seconds
Nodal Proj :0.712576929 seconds
Reactions :0.203358194 seconds
Misc :0.133355912 seconds
Base State :0.033506368 seconds
Time to advance time step: 3.670012338
Writing plotfile 0 after all initialization
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Time to write plotfile: 0.544022047
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 5.789678157e-05
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 173628606.2
MLMG: Initial residual (resid0) = 173628606.2
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02237276477, 1.288541402e-10
MLMG: Timers: Solve = 0.104663944 Iter = 0.09596784 Bottom = 0.006973806
<<< 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 = 173628606.2
MLMG: Initial residual (resid0) = 0.3523824454
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02243344788, 1.292036397e-10
MLMG: Timers: Solve = 0.027043724 Iter = 0.018450753 Bottom = 0.001145426
<<< 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 = 16938159.34
MLMG: Initial residual (resid0) = 16938159.34
MLMG: Final Iter. 11 resid, resid/bnorm = 0.001402716152, 8.281396605e-11
MLMG: Timers: Solve = 0.159405696 Iter = 0.154836132 Bottom = 0.004024173
Done calling nodal solver
Timestep 1 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05
Timing summary:
Advection :2.310342348 seconds
MAC Proj :0.260625019 seconds
Nodal Proj :0.27062636 seconds
Reactions :0.203829651 seconds
Misc :0.142397003 seconds
Base State :0.03335732 seconds
Time to advance time step: 3.188179578
Writing plotfile 1
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Time to write plotfile: 0.551958289
Call to estdt for level 0 gives dt_lev = 1.112476819
Call to estdt for level 1 gives dt_lev = 0.2935011252
Minimum estdt over all levels = 0.2935011252
Timestep 2 starts with TIME = 5.789678157e-05 DT = 6.368645973e-05
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 166900886.1
MLMG: Initial residual (resid0) = 166900886.1
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02290113643, 1.372139895e-10
MLMG: Timers: Solve = 0.104997178 Iter = 0.096212251 Bottom = 0.00702865
<<< 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 = 166900886.1
MLMG: Initial residual (resid0) = 8376.291279
MLMG: Final Iter. 2 resid, resid/bnorm = 0.1567925077, 9.394348429e-10
MLMG: Timers: Solve = 0.043037274 Iter = 0.034477589 Bottom = 0.002336517
<<< 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 = 11608423.25
MLMG: Initial residual (resid0) = 11608423.25
MLMG: Final Iter. 9 resid, resid/bnorm = 0.00100406725, 8.649471404e-11
MLMG: Timers: Solve = 0.13151935 Iter = 0.127014712 Bottom = 0.003316729
Done calling nodal solver
Timestep 2 ends with TIME = 0.0001215832413 DT = 6.368645973e-05
Timing summary:
Advection :2.300674403 seconds
MAC Proj :0.276173123 seconds
Nodal Proj :0.244498977 seconds
Reactions :0.204703718 seconds
Misc :0.137636729 seconds
Base State :0.033849321 seconds
Time to advance time step: 3.164045914
Writing plotfile 2
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Time to write plotfile: 0.543851667
Time to regrid: 0.107812609
Call to estdt for level 0 gives dt_lev = 1.539392515
Call to estdt for level 1 gives dt_lev = 0.3670776927
Minimum estdt over all levels = 0.3670776927
Timestep 3 starts with TIME = 0.0001215832413 DT = 7.00551057e-05
Cell Count:
Level 0, 262144 cells
Level 1, 304128 cells
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 166879353.3
MLMG: Initial residual (resid0) = 166879353.3
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02257638052, 1.352856424e-10
MLMG: Timers: Solve = 0.104540928 Iter = 0.095902187 Bottom = 0.006983484
<<< 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 = 166879353.2
MLMG: Initial residual (resid0) = 0.3055612778
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01938926012, 1.161872918e-10
MLMG: Timers: Solve = 0.027503714 Iter = 0.018719298 Bottom = 0.001170329
<<< 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 = 9504180.596
MLMG: Initial residual (resid0) = 9504180.596
MLMG: Final Iter. 10 resid, resid/bnorm = 0.000827436801, 8.706029864e-11
MLMG: Timers: Solve = 0.149685214 Iter = 0.144882848 Bottom = 0.003799954
Done calling nodal solver
Timestep 3 ends with TIME = 0.000191638347 DT = 7.00551057e-05
Timing summary:
Advection :2.311167686 seconds
MAC Proj :0.266909199 seconds
Nodal Proj :0.263362436 seconds
Reactions :0.202545898 seconds
Misc :0.144606549 seconds
Base State :0.03476287 seconds
Time to advance time step: 3.189076936
Writing plotfile 3
inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09
outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09
Time to write plotfile: 0.544236097
Total Time: 44.98424
Unused ParmParse Variables:
[TOP]::maestro.v(nvals = 1) :: [1]
[TOP]::system.regtest_reduction(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): 9341
[The Arena] space (MB): 9049
[The Device Arena] space (MB): 8
[The Managed Arena] space (MB): 8
[The Pinned Arena] space (MB): 8
AMReX (20.01-220-g45d4129ee9d1) finalized