Initializing CUDA...
CUDA initialized with 1 device.
AMReX (23.02-42-gd8e6401642ff) initialized
Calling Setup()
Calling ReadParameters()
reading extern runtime parameters ...
Calling VariableSetup()
Calling BCSetup()
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
WARNING: variable not found: ash
WARNING: magnesium-24 not provided in inputs file
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.404756409
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 = 214086148.6
MLMG: Initial residual (resid0) = 214086148.6
MLMG: Final Iter. 5 resid, resid/bnorm = 0.005460712517, 2.550707999e-11
MLMG: Timers: Solve = 2.139281456 Iter = 2.091320714 Bottom = 0.000354042
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.373000628
Call to firstdt for level 0 gives dt_lev = 0.005789736987
Multiplying dt_lev by init_shrink; dt_lev = 0.0005789736987
Call to firstdt for level 1 gives dt_lev = 0.002885750885
Multiplying dt_lev by init_shrink; dt_lev = 0.0002885750885
Minimum firstdt over all levels = 0.0002885750885
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 98262295.15
MLMG: Initial residual (resid0) = 98262295.15
MLMG: Final Iter. 1 resid, resid/bnorm = 6.201650538, 6.311322698e-08
MLMG: Timers: Solve = 0.473858598 Iter = 0.427010067 Bottom = 0.000121726
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 21.60061365
Call to estdt for level 1 gives dt_lev = 10.10291125
Minimum estdt over all levels = 10.10291125
Call to estdt at end of istep_divu_iter = 1 gives dt = 10.10291125
Multiplying dt by init_shrink; dt = 1.010291125
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
Doing initial divu iteration #2
Calling nodal solver
MLMG: Initial rhs = 57185021.56
MLMG: Initial residual (resid0) = 57185021.56
MLMG: Final Iter. 1 resid, resid/bnorm = 0.1141292159, 1.995788632e-09
MLMG: Timers: Solve = 0.469690515 Iter = 0.425033713 Bottom = 0.000110625
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 21.60061365
Call to estdt for level 1 gives dt_lev = 10.47579657
Minimum estdt over all levels = 10.47579657
Call to estdt at end of istep_divu_iter = 2 gives dt = 10.47579657
Multiplying dt by init_shrink; dt = 1.047579657
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
Doing initial divu iteration #3
Calling nodal solver
MLMG: Initial rhs = 36948393.66
MLMG: Initial residual (resid0) = 36948393.66
MLMG: Final Iter. 3 resid, resid/bnorm = 0.003112264979, 8.423275469e-11
MLMG: Timers: Solve = 1.295777799 Iter = 1.251097591 Bottom = 0.000241481
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 21.60061365
Call to estdt for level 1 gives dt_lev = 10.62248182
Minimum estdt over all levels = 10.62248182
Call to estdt at end of istep_divu_iter = 3 gives dt = 10.62248182
Multiplying dt by init_shrink; dt = 1.062248182
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
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.350148746
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 0.0002885750885
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 = 300731019.7
MLMG: Initial residual (resid0) = 300731019.7
MLMG: Final Iter. 5 resid, resid/bnorm = 0.284978807, 9.476202596e-10
MLMG: Timers: Solve = 0.047165289 Iter = 0.04238751 Bottom = 0.005101272
<<< 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 = 300731020.2
MLMG: Initial residual (resid0) = 1282.174763
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01771812089, 5.891683828e-11
MLMG: Timers: Solve = 0.021173272 Iter = 0.017070845 Bottom = 0.00199276
<<< 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 = 16091212.77
MLMG: Initial residual (resid0) = 16091212.77
MLMG: Final Iter. 18 resid, resid/bnorm = 0.001556653995, 9.673938302e-11
MLMG: Timers: Solve = 7.39478229 Iter = 7.346495256 Bottom = 0.001027526
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.709780476 seconds
MAC Proj :0.152622301 seconds
Nodal Proj :7.551797409 seconds
Reactions :0.185951014 seconds
Misc :0.078662276 seconds
Base State :0.047008735 seconds
Time to advance time step: 9.679316488
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.352234308
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 0.0002885750885
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 = 300732039.4
MLMG: Initial residual (resid0) = 300732039.4
MLMG: Final Iter. 5 resid, resid/bnorm = 0.2849641219, 9.475682154e-10
MLMG: Timers: Solve = 0.045971405 Iter = 0.041715159 Bottom = 0.005056432
<<< 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 = 300732039.2
MLMG: Initial residual (resid0) = 0.4613403517
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02874185051, 9.557295786e-11
MLMG: Timers: Solve = 0.014050725 Iter = 0.009667985 Bottom = 0.000907625
<<< 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 = 27933872.39
MLMG: Initial residual (resid0) = 27933872.39
MLMG: Final Iter. 5 resid, resid/bnorm = 0.001094801817, 3.919262614e-11
MLMG: Timers: Solve = 2.152242006 Iter = 2.106759772 Bottom = 0.000542678
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.73546847 seconds
MAC Proj :0.142840285 seconds
Nodal Proj :2.308010394 seconds
Reactions :0.184466253 seconds
Misc :0.074958285 seconds
Base State :0.051325931 seconds
Time to advance time step: 4.446271663
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.352808588
Call to estdt for level 0 gives dt_lev = 2.502711311
Call to estdt for level 1 gives dt_lev = 0.5259217834
Minimum estdt over all levels = 0.5259217834
Call to estdt at beginning of step 2 gives dt =0.5259217834
dt_growth factor limits the new dt = 0.0003174325974
Timestep 2 starts with TIME = 0.0002885750885 DT = 0.0003174325974
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 = 294777750.2
MLMG: Initial residual (resid0) = 294777750.2
MLMG: Final Iter. 5 resid, resid/bnorm = 0.2824536264, 9.581918114e-10
MLMG: Timers: Solve = 0.045502781 Iter = 0.041099938 Bottom = 0.005001093
<<< 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 = 294777750
MLMG: Initial residual (resid0) = 0.4034976006
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02507372783, 8.505977072e-11
MLMG: Timers: Solve = 0.013803394 Iter = 0.009479024 Bottom = 0.000884879
<<< 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 = 21482439.4
MLMG: Initial residual (resid0) = 21482439.4
MLMG: Final Iter. 5 resid, resid/bnorm = 0.00126458313, 5.886590003e-11
MLMG: Timers: Solve = 2.129901394 Iter = 2.084934253 Bottom = 0.000512068
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974
Timing summary:
Advection :1.735032663 seconds
MAC Proj :0.141833173 seconds
Nodal Proj :2.287667055 seconds
Reactions :0.185480572 seconds
Misc :0.07591345 seconds
Base State :0.047622622 seconds
Time to advance time step: 4.426408447
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.361732264
Time to regrid: 0.057489622
Call to estdt for level 0 gives dt_lev = 3.210069119
Call to estdt for level 1 gives dt_lev = 0.6510223113
Minimum estdt over all levels = 0.6510223113
Call to estdt at beginning of step 3 gives dt =0.6510223113
dt_growth factor limits the new dt = 0.0003491758571
Timestep 3 starts with TIME = 0.0006060076859 DT = 0.0003491758571
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 = 292704984.1
MLMG: Initial residual (resid0) = 292704984.1
MLMG: Final Iter. 5 resid, resid/bnorm = 0.2731596529, 9.33225151e-10
MLMG: Timers: Solve = 0.044840436 Iter = 0.040625751 Bottom = 0.005020868
<<< 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 = 292704983.8
MLMG: Initial residual (resid0) = 0.3807753602
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02360978279, 8.066067917e-11
MLMG: Timers: Solve = 0.014323992 Iter = 0.009960773 Bottom = 0.00092479
<<< 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 = 16800418.13
MLMG: Initial residual (resid0) = 16800418.13
MLMG: Final Iter. 5 resid, resid/bnorm = 0.001389111143, 8.268312921e-11
MLMG: Timers: Solve = 2.130738948 Iter = 2.08564447 Bottom = 0.000506733
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571
Timing summary:
Advection :1.737400724 seconds
MAC Proj :0.143879867 seconds
Nodal Proj :2.285084522 seconds
Reactions :0.186833483 seconds
Misc :0.074973146 seconds
Base State :0.050657865 seconds
Time to advance time step: 4.428685625
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.355607373
Total Time: 32.39430928
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): 12056
Free GPU global memory (MB): 10127
[The Arena] space allocated (MB): 9042
[The Arena] space used (MB): 0
[The Device Arena] space allocated (MB): 8
[The Device Arena] space used (MB): 0
[The Pinned Arena] space allocated (MB): 33
[The Pinned Arena] space used (MB): 0
AMReX (23.02-42-gd8e6401642ff) finalized