Initializing CUDA...
CUDA initialized with 1 device.
AMReX (23.01-18-gd6f8b2cce0fb) 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.408955903
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.132045632 Iter = 2.084264323 Bottom = 0.000424681
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.382327427
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.201650534, 6.311322694e-08
MLMG: Timers: Solve = 0.471363837 Iter = 0.424857164 Bottom = 0.000119711
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.1141292164, 1.99578864e-09
MLMG: Timers: Solve = 0.468482163 Iter = 0.424008482 Bottom = 0.000108124
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.297881015 Iter = 1.253417862 Bottom = 0.000242909
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.355178983
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.047248769 Iter = 0.041702902 Bottom = 0.004963956
<<< 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.01771812088, 5.891683826e-11
MLMG: Timers: Solve = 0.021960271 Iter = 0.017738703 Bottom = 0.001971883
<<< 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.001556652132, 9.673926727e-11
MLMG: Timers: Solve = 7.380557701 Iter = 7.333778324 Bottom = 0.000995545
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.764174404 seconds
MAC Proj :0.153924509 seconds
Nodal Proj :7.541763217 seconds
Reactions :0.18756516 seconds
Misc :0.079521592 seconds
Base State :0.050160832 seconds
Time to advance time step: 9.727538324
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.356245401
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.284964107, 9.475681659e-10
MLMG: Timers: Solve = 0.04547416 Iter = 0.04106649 Bottom = 0.00500632
<<< 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.014720763 Iter = 0.010373555 Bottom = 0.00086195
<<< 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.001094801351, 3.919260947e-11
MLMG: Timers: Solve = 2.167592017 Iter = 2.120905049 Bottom = 0.000523699
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.783840172 seconds
MAC Proj :0.14257846 seconds
Nodal Proj :2.331383105 seconds
Reactions :0.18555274 seconds
Misc :0.076768701 seconds
Base State :0.051323408 seconds
Time to advance time step: 4.52067538
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.353742461
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.045821977 Iter = 0.041472401 Bottom = 0.005015488
<<< 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.505977073e-11
MLMG: Timers: Solve = 0.014245786 Iter = 0.009572191 Bottom = 0.000867179
<<< 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.886590005e-11
MLMG: Timers: Solve = 2.179465351 Iter = 2.133950133 Bottom = 0.000522484
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974
Timing summary:
Advection :1.783634814 seconds
MAC Proj :0.142528684 seconds
Nodal Proj :2.340545737 seconds
Reactions :0.198182534 seconds
Misc :0.078053676 seconds
Base State :0.05729276 seconds
Time to advance time step: 4.543531682
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.414164486
Time to regrid: 0.060436952
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.047233395 Iter = 0.042877475 Bottom = 0.005063582
<<< 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.014499041 Iter = 0.010302386 Bottom = 0.000852853
<<< 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.131411883 Iter = 2.08610214 Bottom = 0.000444974
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571
Timing summary:
Advection :1.78432083 seconds
MAC Proj :0.147222212 seconds
Nodal Proj :2.291931618 seconds
Reactions :0.190748113 seconds
Misc :0.076759676 seconds
Base State :0.050073435 seconds
Time to advance time step: 4.491552639
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.377130289
Total Time: 32.84957311
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): 10129
[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.01-18-gd6f8b2cce0fb) finalized