Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.05-15-ga9d1c7649200) 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()
Initializing Helmholtz EOS and using Coulomb corrections.
Calling set_method_params()
Calling BCSetup()
Calling init_base_state_geometry()
Calling BaseStateGeometry::Init()
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.397130173
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 = 185424520.7
MLMG: Initial residual (resid0) = 185424520.7
MLMG: Final Iter. 10 resid, resid/bnorm = 0.01339600929, 7.224507979e-11
MLMG: Timers: Solve = 1.935538458 Iter = 1.894001822 Bottom = 0.000741424
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.374364305
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 = 72482174.27
MLMG: Initial residual (resid0) = 72482174.27
MLMG: Final Iter. 2 resid, resid/bnorm = 0.5223139888, 7.20610266e-09
MLMG: Timers: Solve = 0.433099185 Iter = 0.392587205 Bottom = 0.00016562
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 = 43334804.49
MLMG: Initial residual (resid0) = 43334804.49
MLMG: Final Iter. 2 resid, resid/bnorm = 0.02298274694, 5.303530778e-10
MLMG: Timers: Solve = 0.436033264 Iter = 0.395371538 Bottom = 0.000176342
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 = 32768245.23
MLMG: Initial residual (resid0) = 32768245.23
MLMG: Final Iter. 3 resid, resid/bnorm = 0.001478955848, 4.513381286e-11
MLMG: Timers: Solve = 0.636469792 Iter = 0.595715518 Bottom = 0.00025317
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.35829767
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 = 224625388.6
MLMG: Initial residual (resid0) = 224625388.6
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02028143406, 9.029003436e-11
MLMG: Timers: Solve = 0.105198648 Iter = 0.095665523 Bottom = 0.00738567
<<< 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 = 224625389.4
MLMG: Initial residual (resid0) = 3940.279649
MLMG: Final Iter. 2 resid, resid/bnorm = 0.04907616749, 2.184800553e-10
MLMG: Timers: Solve = 0.041662663 Iter = 0.033843195 Bottom = 0.002471194
<<< 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 = 35718891.54
MLMG: Initial residual (resid0) = 35718891.54
MLMG: Final Iter. 31 resid, resid/bnorm = 0.003148810356, 8.815532118e-11
MLMG: Timers: Solve = 5.908489913 Iter = 5.865625169 Bottom = 0.002185113
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.823483677 seconds
MAC Proj :0.258893897 seconds
Nodal Proj :6.053884202 seconds
Reactions :0.185649097 seconds
Misc :0.072512418 seconds
Base State :0.04020528 seconds
Time to advance time step: 8.394891111
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.362099449
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 = 224625175.9
MLMG: Initial residual (resid0) = 224625175.9
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02050752938, 9.129666473e-11
MLMG: Timers: Solve = 0.102452884 Iter = 0.0946917 Bottom = 0.007482071
<<< 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 = 224625175.8
MLMG: Initial residual (resid0) = 0.7033640956
MLMG: Final Iter. 1 resid, resid/bnorm = 0.04446259194, 1.979412672e-10
MLMG: Timers: Solve = 0.02608984 Iter = 0.018357372 Bottom = 0.001272311
<<< 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 = 25278219.2
MLMG: Initial residual (resid0) = 25278219.2
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001846274477, 7.303815442e-11
MLMG: Timers: Solve = 1.930804913 Iter = 1.890620335 Bottom = 0.000714226
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.670846205 seconds
MAC Proj :0.23060235 seconds
Nodal Proj :2.067322059 seconds
Reactions :0.173254309 seconds
Misc :0.061997289 seconds
Base State :0.037603861 seconds
Time to advance time step: 4.204455744
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.361326838
Call to estdt for level 0 gives dt_lev = 2.657850076
Call to estdt for level 1 gives dt_lev = 0.6581035052
Minimum estdt over all levels = 0.6581035052
Call to estdt at beginning of step 2 gives dt =0.6581035052
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 = 221100139
MLMG: Initial residual (resid0) = 221100139
MLMG: Final Iter. 6 resid, resid/bnorm = 0.01983714104, 8.972016538e-11
MLMG: Timers: Solve = 0.102696514 Iter = 0.094808725 Bottom = 0.007505861
<<< 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 = 221100138.8
MLMG: Initial residual (resid0) = 0.6525653573
MLMG: Final Iter. 1 resid, resid/bnorm = 0.04140532156, 1.872695412e-10
MLMG: Timers: Solve = 0.025873191 Iter = 0.018332453 Bottom = 0.001250833
<<< 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 = 19841168.99
MLMG: Initial residual (resid0) = 19841168.99
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001567839296, 7.901950214e-11
MLMG: Timers: Solve = 1.929426026 Iter = 1.889134019 Bottom = 0.000758493
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060015283 DT = 0.000317429372
Timing summary:
Advection :1.677001564 seconds
MAC Proj :0.235821024 seconds
Nodal Proj :2.06915251 seconds
Reactions :0.178201632 seconds
Misc :0.061848113 seconds
Base State :0.039005201 seconds
Time to advance time step: 4.222464337
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.358241676
Time to regrid: 0.056087302
Call to estdt for level 0 gives dt_lev = 3.830037618
Call to estdt for level 1 gives dt_lev = 0.8027988744
Minimum estdt over all levels = 0.8027988744
Call to estdt at beginning of step 3 gives dt =0.8027988744
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 = 221012973.3
MLMG: Initial residual (resid0) = 221012973.3
MLMG: Final Iter. 5 resid, resid/bnorm = 0.2166522741, 9.802694878e-10
MLMG: Timers: Solve = 0.087414374 Iter = 0.079719323 Bottom = 0.006259528
<<< 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 = 221012973.1
MLMG: Initial residual (resid0) = 0.6485630063
MLMG: Final Iter. 1 resid, resid/bnorm = 0.0409702386, 1.853748132e-10
MLMG: Timers: Solve = 0.025691215 Iter = 0.017982653 Bottom = 0.001073008
<<< 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 = 15859855.11
MLMG: Initial residual (resid0) = 15859855.11
MLMG: Final Iter. 11 resid, resid/bnorm = 0.001177113503, 7.421968831e-11
MLMG: Timers: Solve = 2.118033328 Iter = 2.077507013 Bottom = 0.000811522
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551738375 DT = 0.0003491723092
Timing summary:
Advection :1.677774986 seconds
MAC Proj :0.216530354 seconds
Nodal Proj :2.260168871 seconds
Reactions :0.175973652 seconds
Misc :0.061864077 seconds
Base State :0.038770961 seconds
Time to advance time step: 4.3927287
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.355930069
Total Time: 31.8580981
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): 9256
[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.05-15-ga9d1c7649200) finalized