MPI initialized with 4 MPI processes
MPI initialized with thread support level 0
OMP initialized with 4 OMP threads
AMReX (23.10-13-gecaa14456e3f) initialized
Calling Setup()
Calling ReadParameters()
reading extern runtime parameters ...
Calling VariableSetup()
Calling BCSetup()
Calling BaseStateGeometry::Init()
Calling Init()
Calling InitData()
initdata model_File = model.hse.cool.coulomb
model file = model.hse.cool.coulomb
reading initial model
640 points found in the initial model
6 variables found in the initial model
model file mapping, level: 0
dr of MAESTRO base state = 562500.000000
dr of input file data = 562500.000000
maximum radius (cell-centered) of input model = 359718750.000000
setting r_cutoff to 365
radius at r_cutoff 205593750
Maximum HSE Error = 0.000061
(after putting initial model into base state arrays, and
for density < base_cutoff_density)
Writing plotfile reacting_bubble-3d_pltInitData after InitData
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 1.61143551
inner sponge: r_sp , r_tp : 186468750, 224718750
Doing initial projection
Calling nodal solver
MLMG: Initial rhs = 0
MLMG: Initial residual (resid0) = 0
MLMG: No iterations needed
MLMG: Timers: Solve = 0.019145208 Iter = 0 Bottom = 0
Done calling nodal solver
Writing plotfile reacting_bubble-3d_pltafter_InitProj after InitProj
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 1.481854138
Call to firstdt for level 0 gives dt_lev = 0.0004741795139
Multiplying dt_lev by init_shrink; dt_lev = 0.0004741795139
Minimum firstdt over all levels = 0.0004741795139
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 2054.815795
MLMG: Initial residual (resid0) = 2054.815795
MLMG: Final Iter. 7 resid, resid/bnorm = 5.60984148e-10, 2.730094587e-13
MLMG: Timers: Solve = 1.024214189 Iter = 1.002528156 Bottom = 0.03503439
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 0.1719927642
Minimum estdt over all levels = 0.1719927642
Call to estdt at end of istep_divu_iter = 1 gives dt = 0.1719927642
Multiplying dt by init_shrink; dt = 0.1719927642
Ignoring this new dt since it's larger than the previous dt = 0.0004741795139
Writing plotfile reacting_bubble-3d_pltafter_DivuIter after final DivuIter
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 1.513356044
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 0.0004741795139
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 2728631.973
MLMG: Initial residual (resid0) = 2728631.973
MLMG: Final Iter. 9 resid, resid/bnorm = 1.566065475e-05, 5.739379626e-12
MLMG: Timers: Solve = 0.32604496 Iter = 0.302369991 Bottom = 0.053432602
<<< 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 = 2728627.792
MLMG: Initial residual (resid0) = 8729.745178
MLMG: Final Iter. 5 resid, resid/bnorm = 0.0002251297119, 8.250656708e-11
MLMG: Timers: Solve = 0.136947091 Iter = 0.117131862 Bottom = 0.021517005
<<< 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 = 1.157552615e+10
MLMG: Initial residual (resid0) = 1.157552615e+10
MLMG: Final Iter. 7 resid, resid/bnorm = 0.002693429606, 2.326831256e-13
MLMG: Timers: Solve = 1.010887902 Iter = 0.991900082 Bottom = 0.026768097
Done calling nodal solver
Timestep 0 ends with TIME = 0.0004741795139 DT = 0.0004741795139
Timing summary:
Advection :2.136988317 seconds
MAC Proj :0.552795736 seconds
Nodal Proj :1.129390475 seconds
Reactions :2.10071173 seconds
Misc :0.475024771 seconds
Base State :0.004436585 seconds
Time to advance time step: 6.419443814
Writing plotfile 0 after all initialization
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 1.467941612
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 0.0004741795139
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 1308380.202
MLMG: Initial residual (resid0) = 1308380.202
MLMG: Final Iter. 8 resid, resid/bnorm = 0.0001225989981, 9.370288387e-11
MLMG: Timers: Solve = 0.225291502 Iter = 0.206392474 Bottom = 0.036767114
<<< 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 = 1308380.236
MLMG: Initial residual (resid0) = 0.9822512265
MLMG: Final Iter. 4 resid, resid/bnorm = 3.542963896e-05, 2.707900806e-11
MLMG: Timers: Solve = 0.123082471 Iter = 0.103930499 Bottom = 0.017753892
<<< 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 = 5488883.985
MLMG: Initial residual (resid0) = 5488883.985
MLMG: Final Iter. 7 resid, resid/bnorm = 1.254607923e-06, 2.285724979e-13
MLMG: Timers: Solve = 1.028983652 Iter = 1.00679711 Bottom = 0.026075145
Done calling nodal solver
Timestep 1 ends with TIME = 0.0004741795139 DT = 0.0004741795139
Timing summary:
Advection :2.054342658 seconds
MAC Proj :0.430767398 seconds
Nodal Proj :1.146422052 seconds
Reactions :2.094659315 seconds
Misc :0.499498267 seconds
Base State :0.002510461 seconds
Time to advance time step: 6.253564654
Call to estdt for level 0 gives dt_lev = 0.2012294718
Minimum estdt over all levels = 0.2012294718
Call to estdt at beginning of step 2 gives dt =0.2012294718
dt_growth factor limits the new dt = 0.0005215974653
Timestep 2 starts with TIME = 0.0004741795139 DT = 0.0005215974653
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 1426181.405
MLMG: Initial residual (resid0) = 1426181.405
MLMG: Final Iter. 8 resid, resid/bnorm = 0.000137500756, 9.641182772e-11
MLMG: Timers: Solve = 0.224191102 Iter = 0.205359083 Bottom = 0.036536642
<<< 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 = 1426184.06
MLMG: Initial residual (resid0) = 5132.459686
MLMG: Final Iter. 6 resid, resid/bnorm = 3.080249007e-05, 2.159783644e-11
MLMG: Timers: Solve = 0.179085911 Iter = 0.160305202 Bottom = 0.028258375
<<< 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 = 6063789.798
MLMG: Initial residual (resid0) = 6063789.798
MLMG: Final Iter. 7 resid, resid/bnorm = 1.490382713e-06, 2.457840332e-13
MLMG: Timers: Solve = 1.011440192 Iter = 0.992235728 Bottom = 0.025979248
Done calling nodal solver
Timestep 2 ends with TIME = 0.0009957769792 DT = 0.0005215974653
Timing summary:
Advection :2.042228512 seconds
MAC Proj :0.484311425 seconds
Nodal Proj :1.117146481 seconds
Reactions :2.109191535 seconds
Misc :0.526548357 seconds
Base State :0.004989934 seconds
Time to advance time step: 6.301060765
Call to estdt for level 0 gives dt_lev = 0.20128497
Minimum estdt over all levels = 0.20128497
Call to estdt at beginning of step 3 gives dt =0.20128497
dt_growth factor limits the new dt = 0.0005737572118
Timestep 3 starts with TIME = 0.0009957769792 DT = 0.0005737572118
Cell Count:
Level 0, 1327104 cells
inner sponge: r_sp , r_tp : 186468750, 224718750
<<< STEP 1 : react state >>>
<<< STEP 2 : make w0 >>>
<<< STEP 3 : create MAC velocities >>>
MLMG: Initial rhs = 1597942.993
MLMG: Initial residual (resid0) = 1597942.993
MLMG: Final Iter. 8 resid, resid/bnorm = 0.0001546049898, 9.675250647e-11
MLMG: Timers: Solve = 0.220563848 Iter = 0.201559252 Bottom = 0.032779548
<<< 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 = 1597945.914
MLMG: Initial residual (resid0) = 5863.170738
MLMG: Final Iter. 6 resid, resid/bnorm = 3.355263198e-05, 2.099735147e-11
MLMG: Timers: Solve = 0.173019759 Iter = 0.154286004 Bottom = 0.02762951
<<< 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 = 6670310.291
MLMG: Initial residual (resid0) = 6670310.291
MLMG: Final Iter. 7 resid, resid/bnorm = 4.630721957e-06, 6.942288672e-13
MLMG: Timers: Solve = 1.060366841 Iter = 1.040539788 Bottom = 0.032190362
Done calling nodal solver
Timestep 3 ends with TIME = 0.001569534191 DT = 0.0005737572118
Timing summary:
Advection :2.032796195 seconds
MAC Proj :0.47450644 seconds
Nodal Proj :1.173903898 seconds
Reactions :2.060059584 seconds
Misc :0.488040711 seconds
Base State :0.002621639 seconds
Time to advance time step: 6.255840215
Writing plotfile 3
inner sponge: r_sp , r_tp : 186468750, 224718750
Time to write plotfile: 1.640067969
Total Time: 36.75246318
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) :: [reacting_bubble-3d_chk]
[TOP]::amr.checkpoint_files_output(nvals = 1) :: [0]
AMReX (23.10-13-gecaa14456e3f) finalized