Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.02-404-g433143b8cfa6) 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 InitBaseStateGeometry()
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.52428498
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. 11 resid, resid/bnorm = 0.01346848923, 7.154931382e-11
MLMG: Timers: Solve = 2.166585605 Iter = 2.127761338 Bottom = 0.000888998
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.501112684
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.57
MLMG: Initial residual (resid0) = 80008358.57
MLMG: Final Iter. 2 resid, resid/bnorm = 0.4733840131, 5.916681977e-09
MLMG: Timers: Solve = 0.438697543 Iter = 0.401551609 Bottom = 0.00016781
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 = 44170163.98
MLMG: Initial residual (resid0) = 44170163.98
MLMG: Final Iter. 2 resid, resid/bnorm = 0.03121123224, 7.066134565e-10
MLMG: Timers: Solve = 0.440610767 Iter = 0.40020118 Bottom = 0.000177922
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.08
MLMG: Initial residual (resid0) = 26504865.08
MLMG: Final Iter. 3 resid, resid/bnorm = 0.00104162774, 3.929949226e-11
MLMG: Timers: Solve = 0.633508512 Iter = 0.592634767 Bottom = 0.000254485
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.475200931
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.0223493916, 1.287196687e-10
MLMG: Timers: Solve = 0.101951643 Iter = 0.091491772 Bottom = 0.006777955
<<< 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.124611
MLMG: Final Iter. 2 resid, resid/bnorm = 0.04143720183, 2.386545001e-10
MLMG: Timers: Solve = 0.040127764 Iter = 0.032202908 Bottom = 0.002265865
<<< 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. 33 resid, resid/bnorm = 0.003207540882, 8.980050871e-11
MLMG: Timers: Solve = 6.397483569 Iter = 6.355887382 Bottom = 0.002493367
Done calling nodal solver
Timestep 0 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05
Timing summary:
Advection :1.965245399 seconds
MAC Proj :0.259939085 seconds
Nodal Proj :6.552962916 seconds
Reactions :0.198530311 seconds
Misc :0.128943128 seconds
Base State :0.031585025 seconds
Time to advance time step: 9.106057886
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.474902034
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.02239027107, 1.289549663e-10
MLMG: Timers: Solve = 0.098443276 Iter = 0.090702845 Bottom = 0.006845448
<<< 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.3528474217
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02246488312, 1.293846884e-10
MLMG: Timers: Solve = 0.025899021 Iter = 0.01791166 Bottom = 0.001189479
<<< 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. 8 resid, resid/bnorm = 0.001684285467, 9.943733756e-11
MLMG: Timers: Solve = 1.58992926 Iter = 1.549907292 Bottom = 0.000638147
Done calling nodal solver
Timestep 1 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05
Timing summary:
Advection :1.970955371 seconds
MAC Proj :0.235877804 seconds
Nodal Proj :1.741913344 seconds
Reactions :0.195354804 seconds
Misc :0.137680507 seconds
Base State :0.032673792 seconds
Time to advance time step: 4.282284927
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.478091004
Call to estdt for level 0 gives dt_lev = 1.112476815
Call to estdt for level 1 gives dt_lev = 0.2935011253
Minimum estdt over all levels = 0.2935011253
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.02282220311, 1.367410542e-10
MLMG: Timers: Solve = 0.098671099 Iter = 0.090971331 Bottom = 0.006822451
<<< 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) = 142.3038985
MLMG: Final Iter. 1 resid, resid/bnorm = 0.094277007, 5.648682233e-10
MLMG: Timers: Solve = 0.025460142 Iter = 0.017785833 Bottom = 0.001139566
<<< 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. 7 resid, resid/bnorm = 0.0010750246, 9.26072884e-11
MLMG: Timers: Solve = 1.410169905 Iter = 1.368740052 Bottom = 0.000581241
Done calling nodal solver
Timestep 2 ends with TIME = 0.0001215832413 DT = 6.368645973e-05
Timing summary:
Advection :1.9538228 seconds
MAC Proj :0.239269891 seconds
Nodal Proj :1.564361156 seconds
Reactions :0.196799731 seconds
Misc :0.133447101 seconds
Base State :0.03223432 seconds
Time to advance time step: 4.088154553
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.470875476
Time to regrid: 0.100611775
Call to estdt for level 0 gives dt_lev = 1.539392515
Call to estdt for level 1 gives dt_lev = 0.3670776926
Minimum estdt over all levels = 0.3670776926
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.02259801328, 1.354152736e-10
MLMG: Timers: Solve = 0.099385114 Iter = 0.091654586 Bottom = 0.006889147
<<< 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.3060922084
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01942452794, 1.163986291e-10
MLMG: Timers: Solve = 0.025907817 Iter = 0.018006397 Bottom = 0.001173635
<<< 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.595
MLMG: Initial residual (resid0) = 9504180.595
MLMG: Final Iter. 7 resid, resid/bnorm = 0.0008972726064, 9.440820252e-11
MLMG: Timers: Solve = 1.458166904 Iter = 1.412725235 Bottom = 0.000575588
Done calling nodal solver
Timestep 3 ends with TIME = 0.000191638347 DT = 7.00551057e-05
Timing summary:
Advection :2.126929417 seconds
MAC Proj :0.238659496 seconds
Nodal Proj :1.634334145 seconds
Reactions :0.20237408 seconds
Misc :0.139752765 seconds
Base State :0.033147014 seconds
Time to advance time step: 4.342689417
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.526102277
Total Time: 57.26157345
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): 9350
[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.02-404-g433143b8cfa6) finalized