Initializing AMReX (25.04-16-gae7024c40673)...
MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (25.04-16-gae7024c40673) 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.876536135
inner sponge: r_sp , r_tp : 151953125, 183203125
outer sponge: r_sp_outer, r_tp_outer: 183203125, 198828125
[1m[32mDoing initial projection[0m
Calling nodal solver
MLMG: Initial rhs = 214086148.6
MLMG: Initial residual (resid0) = 214086148.6
MLMG: Final Iter. 8 resid, resid/bnorm = 0.01809367537, 8.451586192e-11
MLMG: Timers: Solve = 0.235211526 Iter = 0.22856584 Bottom = 0.000165907
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.832115479
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
[1m[32mDoing initial divu iteration #1[0m
Calling nodal solver
MLMG: Initial rhs = 98262295.13
MLMG: Initial residual (resid0) = 98262295.13
MLMG: Final Iter. 2 resid, resid/bnorm = 0.3533489313, 3.595976777e-09
MLMG: Timers: Solve = 0.06378394 Iter = 0.057729753 Bottom = 4.0748e-05
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.10291158
Minimum estdt over all levels = 10.10291158
Call to estdt at end of istep_divu_iter = 1 gives dt = 10.10291158
Multiplying dt by init_shrink; dt = 1.010291158
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
[1m[32mDoing initial divu iteration #2[0m
Calling nodal solver
MLMG: Initial rhs = 57185020
MLMG: Initial residual (resid0) = 57185020
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01771998301, 3.098710644e-10
MLMG: Timers: Solve = 0.06460272 Iter = 0.058103771 Bottom = 4.2482e-05
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.47579658
Minimum estdt over all levels = 10.47579658
Call to estdt at end of istep_divu_iter = 2 gives dt = 10.47579658
Multiplying dt by init_shrink; dt = 1.047579658
Ignoring this new dt since it's larger than the previous dt = 0.0002885750885
[1m[32mDoing initial divu iteration #3[0m
Calling nodal solver
MLMG: Initial rhs = 36948393.74
MLMG: Initial residual (resid0) = 36948393.74
MLMG: Final Iter. 3 resid, resid/bnorm = 0.002141558866, 5.796081099e-11
MLMG: Timers: Solve = 0.094804089 Iter = 0.088666832 Bottom = 6.0272e-05
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: 1.590064762
[1m[32mDoing initial pressure iteration #1[0m
[1m[32m
Timestep 0 starts with TIME = 0 DT = 0.0002885750885[0m
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.2849493027, 9.47522151e-10
MLMG: Timers: Solve = 0.058503624 Iter = 0.051657537 Bottom = 0.000150571
<<< 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.175235
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01771806707, 5.891665934e-11
MLMG: Timers: Solve = 0.025759497 Iter = 0.021530605 Bottom = 5.9504e-05
<<< 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.78
MLMG: Initial residual (resid0) = 16091212.78
MLMG: Final Iter. 31 resid, resid/bnorm = 0.001162230037, 7.222762218e-11
MLMG: Timers: Solve = 0.846254502 Iter = 0.840954356 Bottom = 0.000636104
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.055036363 seconds
MAC Proj :0.113020146 seconds
Nodal Proj :0.901112257 seconds
Reactions :1.085834428 seconds
Misc :0.25984108 seconds
Base State :0.008404518 seconds
Time to advance time step: 3.445042149
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.803331853
[1m[32mBeginning main evolution[0m
[1m[32m
Timestep 1 starts with TIME = 0 DT = 0.0002885750885[0m
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.2849341929, 9.47468695e-10
MLMG: Timers: Solve = 0.057243912 Iter = 0.05259255 Bottom = 0.00014405
<<< 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.4609831659
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02871837036, 9.54948812e-11
MLMG: Timers: Solve = 0.015889174 Iter = 0.011430803 Bottom = 2.7961e-05
<<< 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. 7 resid, resid/bnorm = 0.0019977591, 7.151744206e-11
MLMG: Timers: Solve = 0.201478062 Iter = 0.196096811 Bottom = 0.000141217
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection :1.000253893 seconds
MAC Proj :0.099260333 seconds
Nodal Proj :0.257454914 seconds
Reactions :1.118108305 seconds
Misc :0.230535372 seconds
Base State :0.00539958 seconds
Time to advance time step: 2.730102922
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.825834034
Call to estdt for level 0 gives dt_lev = 2.502711294
Call to estdt for level 1 gives dt_lev = 0.5259217833
Minimum estdt over all levels = 0.5259217833
Call to estdt at beginning of step 2 gives dt =0.5259217833
dt_growth factor limits the new dt = 0.0003174325974
[1m[32m
Timestep 2 starts with TIME = 0.0002885750885 DT = 0.0003174325974[0m
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.2824965268, 9.583373463e-10
MLMG: Timers: Solve = 0.062178827 Iter = 0.05769771 Bottom = 0.000151801
<<< 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.4040347729
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02510970458, 8.518181776e-11
MLMG: Timers: Solve = 0.016184144 Iter = 0.011793533 Bottom = 2.8386e-05
<<< 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.001258467135, 5.858120261e-11
MLMG: Timers: Solve = 0.151131172 Iter = 0.144849679 Bottom = 0.000104163
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974
Timing summary:
Advection :1.006106505 seconds
MAC Proj :0.104352858 seconds
Nodal Proj :0.207414286 seconds
Reactions :1.106911336 seconds
Misc :0.233005375 seconds
Base State :0.005474355 seconds
Time to advance time step: 2.658381994
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.815468792
Time to regrid: 0.124174441
Call to estdt for level 0 gives dt_lev = 3.21006914
Call to estdt for level 1 gives dt_lev = 0.6510223114
Minimum estdt over all levels = 0.6510223114
Call to estdt at beginning of step 3 gives dt =0.6510223114
dt_growth factor limits the new dt = 0.0003491758571
[1m[32m
Timestep 3 starts with TIME = 0.0006060076859 DT = 0.0003491758571[0m
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.2731440663, 9.331719007e-10
MLMG: Timers: Solve = 0.058682172 Iter = 0.053906111 Bottom = 0.000142374
<<< 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.3805782406
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02359678427, 8.06162709e-11
MLMG: Timers: Solve = 0.016513459 Iter = 0.011912404 Bottom = 2.7239e-05
<<< 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.001384173593, 8.238923472e-11
MLMG: Timers: Solve = 0.151252588 Iter = 0.145953984 Bottom = 0.000101717
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571
Timing summary:
Advection :1.007465848 seconds
MAC Proj :0.10144255 seconds
Nodal Proj :0.209583255 seconds
Reactions :1.119589006 seconds
Misc :0.228942308 seconds
Base State :0.006437196 seconds
Time to advance time step: 2.667572962
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.824493058
Total Time: 21.89715164
Unused ParmParse Variables:
[TOP]::amr.check_file(nvals = 1) :: [wdconvect-amr_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
AMReX (25.04-16-gae7024c40673) finalized