Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.01-37-gf349c9281fc4) 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.57965044
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.015492609, 8.230214428e-11
MLMG: Timers: Solve = 0.112722578 Iter = 0.10248314 Bottom = 0.002613071
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.580517503
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.038608601 Iter = 0.030558599 Bottom = 0.000774582
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.2234185169, 5.058131826e-09
MLMG: Timers: Solve = 0.035416446 Iter = 0.030344214 Bottom = 0.000746033
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.002006638679, 7.570831506e-11
MLMG: Timers: Solve = 0.049448668 Iter = 0.044310176 Bottom = 0.001110835
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.564035345
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.02233196562, 1.286193051e-10
MLMG: Timers: Solve = 0.108747004 Iter = 0.096599638 Bottom = 0.007155877
<<< 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.04143702225, 2.386534659e-10
MLMG: Timers: Solve = 0.043652562 Iter = 0.034365241 Bottom = 0.002370477
<<< 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.002904137131, 8.130620975e-11
MLMG: Timers: Solve = 0.615118688 Iter = 0.60988658 Bottom = 0.016101414
Done calling nodal solver
Timestep 0 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05
Timing summary:
Advection :4.215483168 seconds
MAC Proj :0.293332535 seconds
Nodal Proj :0.735566082 seconds
Reactions :0.20920153 seconds
Misc :0.142395527 seconds
Base State :0.035988941 seconds
Time to advance time step: 5.596401195
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.567179413
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.02237275778, 1.288541e-10
MLMG: Timers: Solve = 0.099413461 Iter = 0.091043949 Bottom = 0.006736937
<<< 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.026334938 Iter = 0.017924355 Bottom = 0.001139747
<<< 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.001402715221, 8.281391107e-11
MLMG: Timers: Solve = 0.150479147 Iter = 0.146024691 Bottom = 0.003792177
Done calling nodal solver
Timestep 1 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05
Timing summary:
Advection :3.973372289 seconds
MAC Proj :0.248865702 seconds
Nodal Proj :0.258694806 seconds
Reactions :0.193905004 seconds
Misc :0.145165026 seconds
Base State :0.030972685 seconds
Time to advance time step: 4.820518369
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.543920543
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.02290116064, 1.372141346e-10
MLMG: Timers: Solve = 0.099437233 Iter = 0.09092518 Bottom = 0.006733188
<<< 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.39434843e-10
MLMG: Timers: Solve = 0.040716604 Iter = 0.032454003 Bottom = 0.002241177
<<< 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.124287671 Iter = 0.119715174 Bottom = 0.00310796
Done calling nodal solver
Timestep 2 ends with TIME = 0.0001215832413 DT = 6.368645973e-05
Timing summary:
Advection :3.986454317 seconds
MAC Proj :0.267926477 seconds
Nodal Proj :0.234937227 seconds
Reactions :0.193773499 seconds
Misc :0.139462267 seconds
Base State :0.031878762 seconds
Time to advance time step: 4.822963151
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.574331813
Time to regrid: 0.102849155
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.02257641032, 1.35285821e-10
MLMG: Timers: Solve = 0.099310256 Iter = 0.090959753 Bottom = 0.006743574
<<< 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.026415233 Iter = 0.017917625 Bottom = 0.001140759
<<< 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.0008274372667, 8.706034764e-11
MLMG: Timers: Solve = 0.137584105 Iter = 0.133149584 Bottom = 0.003449781
Done calling nodal solver
Timestep 3 ends with TIME = 0.000191638347 DT = 7.00551057e-05
Timing summary:
Advection :3.991326888 seconds
MAC Proj :0.254211266 seconds
Nodal Proj :0.246855737 seconds
Reactions :0.197560825 seconds
Misc :0.13882174 seconds
Base State :0.031650879 seconds
Time to advance time step: 4.82918089
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.558201109
Total Time: 50.91078649
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): 9340
[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-37-gf349c9281fc4) finalized