MPI initialized with 2 MPI processes
MPI initialized with thread support level 0
OMP initialized with 4 OMP threads
AMReX (22.01-6-g5726371bbb1a) 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 ...
Calling VariableSetup()
Calling BCSetup()
Calling BaseStateGeometry::Init()
Calling Init()
Calling InitData()
initdata model_File =
cutoff densities:
low density cutoff (for mapping the model) = 1e-10
buoyancy cutoff density
(for zeroing rho - rho_0, centrifugal term) = 5e-10
anelastic cutoff = 1e-10
Maximum HSE Error = 4.547473509e-13
(after putting initial model into base state arrays, and
for density < base_cutoff_density)
Writing plotfile RT_pltInitData after InitData
Time to write plotfile: 0.06851243973
Doing initial projection
Calling nodal solver
MLMG: Initial rhs = 0
MLMG: Initial residual (resid0) = 0
MLMG: No iterations needed
MLMG: Timers: Solve = 0.002491389 Iter = 0 Bottom = 0
Done calling nodal solver
Writing plotfile RT_pltafter_InitProj after InitProj
Time to write plotfile: 0.06131410599
Call to firstdt for level 0 gives dt_lev = 0.0006089835859
Multiplying dt_lev by init_shrink; dt_lev = 6.089835859e-05
Minimum firstdt over all levels = 6.089835859e-05
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 0
MLMG: Initial residual (resid0) = 0
MLMG: No iterations needed
MLMG: Timers: Solve = 0.002404138 Iter = 0 Bottom = 0
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 0.08660444773
Minimum estdt over all levels = 0.08660444773
Call to estdt at end of istep_divu_iter = 1 gives dt = 0.08660444773
Multiplying dt by init_shrink; dt = 0.008660444773
Ignoring this new dt since it's larger than the previous dt = 6.089835859e-05
Writing plotfile RT_pltafter_DivuIter after final DivuIter
Time to write plotfile: 0.06060862541
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 6.089835859e-05
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 0.002421164472
MLMG: Initial residual (resid0) = 0.002421164472
MLMG: Final Iter. 9 resid, resid/bnorm = 1.751907767e-13, 7.235806522e-11
MLMG: Timers: Solve = 0.026833606 Iter = 0.025556974 Bottom = 0.003031849
<<< 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 = 0.002421205691
MLMG: Initial residual (resid0) = 2.151449931e-06
MLMG: Final Iter. 7 resid, resid/bnorm = 9.927909435e-14, 4.100399017e-11
MLMG: Timers: Solve = 0.020637072 Iter = 0.019340899 Bottom = 0.002154427
<<< 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 = 58.12657035
MLMG: Initial residual (resid0) = 58.12657035
MLMG: Final Iter. 8 resid, resid/bnorm = 2.15578666e-11, 3.708780077e-13
MLMG: Timers: Solve = 0.069372673 Iter = 0.068180167 Bottom = 0.004471901
Done calling nodal solver
Timestep 0 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05
Timing summary:
Advection :0.07403898239 seconds
MAC Proj :0.05349493027 seconds
Nodal Proj :0.07539820671 seconds
Reactions :0.007260084152 seconds
Misc :0.02550411224 seconds
Base State :0.0003921985626 seconds
Time to advance time step: 0.2395207882
Writing plotfile 0 after all initialization
Time to write plotfile: 0.1504199505
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 6.089835859e-05
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 2.394238647e-05
MLMG: Initial residual (resid0) = 2.394238647e-05
MLMG: Final Iter. 9 resid, resid/bnorm = 3.181792119e-16, 1.328936914e-11
MLMG: Timers: Solve = 0.037866167 Iter = 0.03577853 Bottom = 0.006008658
<<< 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 = 2.394218844e-05
MLMG: Initial residual (resid0) = 1.928423574e-06
MLMG: Final Iter. 9 resid, resid/bnorm = 6.036515824e-16, 2.521288244e-11
MLMG: Timers: Solve = 0.036787828 Iter = 0.034747024 Bottom = 0.005286864
<<< 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 = 0.003539812449
MLMG: Initial residual (resid0) = 0.003539812449
MLMG: Final Iter. 8 resid, resid/bnorm = 1.336170757e-15, 3.77469365e-13
MLMG: Timers: Solve = 0.07185152 Iter = 0.070109136 Bottom = 0.005706657
Done calling nodal solver
Timestep 1 ends with TIME = 6.089835859e-05 DT = 6.089835859e-05
Timing summary:
Advection :0.07461643219 seconds
MAC Proj :0.07995247841 seconds
Nodal Proj :0.07813286781 seconds
Reactions :0.006771564484 seconds
Misc :0.01695632935 seconds
Base State :0.0004260540009 seconds
Time to advance time step: 0.2650647163
Call to estdt for level 0 gives dt_lev = 0.3128878797
Minimum estdt over all levels = 0.3128878797
Call to estdt at beginning of step 2 gives dt =0.3128878797
dt_growth factor limits the new dt = 6.698819445e-05
Timestep 2 starts with TIME = 6.089835859e-05 DT = 6.698819445e-05
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 3.050043585e-05
MLMG: Initial residual (resid0) = 3.050043585e-05
MLMG: Final Iter. 9 resid, resid/bnorm = 2.284813777e-15, 7.491085662e-11
MLMG: Timers: Solve = 0.039276383 Iter = 0.037525748 Bottom = 0.005746822
<<< 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 = 3.059063412e-05
MLMG: Initial residual (resid0) = 5.66083759e-06
MLMG: Final Iter. 10 resid, resid/bnorm = 4.892619004e-16, 1.59938463e-11
MLMG: Timers: Solve = 0.045905586 Iter = 0.044016242 Bottom = 0.005840238
<<< 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 = 0.003894013591
MLMG: Initial residual (resid0) = 0.003894013591
MLMG: Final Iter. 8 resid, resid/bnorm = 1.456300358e-15, 3.739844056e-13
MLMG: Timers: Solve = 0.071610186 Iter = 0.069906551 Bottom = 0.005596405
Done calling nodal solver
Timestep 2 ends with TIME = 0.000127886553 DT = 6.698819445e-05
Timing summary:
Advection :0.08489513397 seconds
MAC Proj :0.09297704697 seconds
Nodal Proj :0.07922697067 seconds
Reactions :0.008228778839 seconds
Misc :0.01756882668 seconds
Base State :0.0004982948303 seconds
Time to advance time step: 0.2834761143
Call to estdt for level 0 gives dt_lev = 0.3130286362
Minimum estdt over all levels = 0.3130286362
Call to estdt at beginning of step 3 gives dt =0.3130286362
dt_growth factor limits the new dt = 7.36870139e-05
Timestep 3 starts with TIME = 0.000127886553 DT = 7.36870139e-05
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 3.168800511e-05
MLMG: Initial residual (resid0) = 3.168800511e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 5.650717727e-16, 1.783235552e-11
MLMG: Timers: Solve = 0.041473312 Iter = 0.039832886 Bottom = 0.006227368
<<< 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 = 3.178722337e-05
MLMG: Initial residual (resid0) = 1.012935749e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 6.83363312e-16, 2.149804983e-11
MLMG: Timers: Solve = 0.041607521 Iter = 0.040043554 Bottom = 0.005930906
<<< 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 = 0.004283414244
MLMG: Initial residual (resid0) = 0.004283414244
MLMG: Final Iter. 8 resid, resid/bnorm = 1.625435897e-15, 3.794720297e-13
MLMG: Timers: Solve = 0.073033194 Iter = 0.071370851 Bottom = 0.00562424
Done calling nodal solver
Timestep 3 ends with TIME = 0.0002015735669 DT = 7.36870139e-05
Timing summary:
Advection :0.07200503349 seconds
MAC Proj :0.08906912804 seconds
Nodal Proj :0.0795788765 seconds
Reactions :0.006958961487 seconds
Misc :0.02866125107 seconds
Base State :0.0007336139679 seconds
Time to advance time step: 0.2766151428
Call to estdt for level 0 gives dt_lev = 0.3130285358
Minimum estdt over all levels = 0.3130285358
Call to estdt at beginning of step 4 gives dt =0.3130285358
dt_growth factor limits the new dt = 8.105571529e-05
Timestep 4 starts with TIME = 0.0002015735669 DT = 8.105571529e-05
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 3.570538736e-05
MLMG: Initial residual (resid0) = 3.570538736e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 7.571110825e-16, 2.120439347e-11
MLMG: Timers: Solve = 0.041202768 Iter = 0.039627343 Bottom = 0.006280951
<<< 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 = 3.581452755e-05
MLMG: Initial residual (resid0) = 1.504387741e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 8.870129024e-16, 2.476684639e-11
MLMG: Timers: Solve = 0.041082392 Iter = 0.039524342 Bottom = 0.005938324
<<< 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 = 0.004711755906
MLMG: Initial residual (resid0) = 0.004711755906
MLMG: Final Iter. 8 resid, resid/bnorm = 1.773321073e-15, 3.763609806e-13
MLMG: Timers: Solve = 0.067739623 Iter = 0.066157739 Bottom = 0.005637406
Done calling nodal solver
Timestep 4 ends with TIME = 0.0002826292822 DT = 8.105571529e-05
Timing summary:
Advection :0.07154297829 seconds
MAC Proj :0.08817362785 seconds
Nodal Proj :0.07386922836 seconds
Reactions :0.006035804749 seconds
Misc :0.01495885849 seconds
Base State :0.0004301071167 seconds
Time to advance time step: 0.2549021244
Call to estdt for level 0 gives dt_lev = 0.313028505
Minimum estdt over all levels = 0.313028505
Call to estdt at beginning of step 5 gives dt =0.313028505
dt_growth factor limits the new dt = 8.916128682e-05
Timestep 5 starts with TIME = 0.0002826292822 DT = 8.916128682e-05
Cell Count:
Level 0, 131072 cells
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 4.013180326e-05
MLMG: Initial residual (resid0) = 4.013180326e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 9.681417181e-16, 2.412405224e-11
MLMG: Timers: Solve = 0.041146476 Iter = 0.03959826 Bottom = 0.00612412
<<< 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 = 4.025185748e-05
MLMG: Initial residual (resid0) = 2.044985444e-05
MLMG: Final Iter. 10 resid, resid/bnorm = 1.111034482e-15, 2.760206737e-11
MLMG: Timers: Solve = 0.041092726 Iter = 0.039455509 Bottom = 0.005885782
<<< 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 = 0.005182931901
MLMG: Initial residual (resid0) = 0.005182931901
MLMG: Final Iter. 8 resid, resid/bnorm = 1.960671209e-15, 3.782938395e-13
MLMG: Timers: Solve = 0.068076875 Iter = 0.066378532 Bottom = 0.005604447
Done calling nodal solver
Timestep 5 ends with TIME = 0.000371790569 DT = 8.916128682e-05
Timing summary:
Advection :0.07080745697 seconds
MAC Proj :0.08815407753 seconds
Nodal Proj :0.07432770729 seconds
Reactions :0.006016016006 seconds
Misc :0.01533389091 seconds
Base State :0.000440120697 seconds
Time to advance time step: 0.2549612522
Writing plotfile 5
Time to write plotfile: 0.05202865601
Total Time: 2.023777485
Unused ParmParse Variables:
[TOP]::maestro.v(nvals = 1) :: [1]
[TOP]::amr.ref_ratio(nvals = 4) :: [2, 2, 2, 2]
[TOP]::amr.check_file(nvals = 1) :: [RT_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
AMReX (22.01-6-g5726371bbb1a) finalized