MPI initialized with 16 MPI processes
MPI initialized with thread support level 0
AMReX (22.03-10-gfbdfc840de7c) 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.823268513
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               = 268163993.9
MLMG: Initial residual (resid0) = 268163993.9
MLMG: Final Iter. 5 resid, resid/bnorm = 0.009204054251, 3.432248348e-11
MLMG: Timers: Solve = 0.167332335 Iter = 0.161700486 Bottom = 0.000169916
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.812231637
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               = 100896162.8
MLMG: Initial residual (resid0) = 100896162.8
MLMG: Final Iter. 2 resid, resid/bnorm = 1.087015875, 1.077360967e-08
MLMG: Timers: Solve = 0.07062135 Iter = 0.065805079 Bottom = 6.4917e-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 = 8.768320826
Minimum estdt over all levels = 8.768320826
Call to estdt at end of istep_divu_iter = 1 gives dt = 8.768320826
Multiplying dt by init_shrink; dt = 0.8768320826
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               = 49401163.56
MLMG: Initial residual (resid0) = 49401163.56
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01935658464, 3.918244682e-10
MLMG: Timers: Solve = 0.072089531 Iter = 0.067176299 Bottom = 6.8093e-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 = 8.962174275
Minimum estdt over all levels = 8.962174275
Call to estdt at end of istep_divu_iter = 2 gives dt = 8.962174275
Multiplying dt by init_shrink; dt = 0.8962174275
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               = 29021861.06
MLMG: Initial residual (resid0) = 29021861.06
MLMG: Final Iter. 3 resid, resid/bnorm = 0.001478702134, 5.095132014e-11
MLMG: Timers: Solve = 0.106385124 Iter = 0.101508472 Bottom = 9.5545e-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 = 9.055396423
Minimum estdt over all levels = 9.055396423
Call to estdt at end of istep_divu_iter = 3 gives dt = 9.055396423
Multiplying dt by init_shrink; dt = 0.9055396423
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.807453961
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               = 264654294.1
MLMG: Initial residual (resid0) = 264654294.1
MLMG: Final Iter. 20 resid, resid/bnorm = 0.2198064644, 8.305418399e-10
MLMG: Timers: Solve = 0.20391107 Iter = 0.197929547 Bottom = 0.000881452
<<< 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               = 264654294.5
MLMG: Initial residual (resid0) = 5771.024244
MLMG: Final Iter. 2 resid, resid/bnorm = 0.1596843153, 6.03369447e-10
MLMG: Timers: Solve = 0.025113876 Iter = 0.020737889 Bottom = 9.5524e-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               = 14853703.25
MLMG: Initial residual (resid0) = 14853703.25
MLMG: Final Iter. 14 resid, resid/bnorm = 0.0008678463371, 5.842626062e-11
MLMG: Timers: Solve = 0.473164867 Iter = 0.46817914 Bottom = 0.000449939
Done calling nodal solver

Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection  :0.955921095 seconds
MAC Proj   :0.252961395 seconds
Nodal Proj :0.521421988 seconds
Reactions  :1.050534711 seconds
Misc       :0.220548021 seconds
Base State :0.004782144 seconds
Time to advance time step: 3.027664864

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.80826398
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               = 264654231
MLMG: Initial residual (resid0) = 264654231
MLMG: Final Iter. 20 resid, resid/bnorm = 0.2199404999, 8.310484931e-10
MLMG: Timers: Solve = 0.20205965 Iter = 0.197699126 Bottom = 0.000879746
<<< 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               = 264654230.8
MLMG: Initial residual (resid0) = 8.258030203
MLMG: Final Iter. 2 resid, resid/bnorm = 0.1579990945, 5.9700196e-10
MLMG: Timers: Solve = 0.025215781 Iter = 0.020845297 Bottom = 9.4833e-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               = 19692713.45
MLMG: Initial residual (resid0) = 19692713.45
MLMG: Final Iter. 5 resid, resid/bnorm = 0.001143385301, 5.80613385e-11
MLMG: Timers: Solve = 0.17353025 Iter = 0.168602821 Bottom = 0.0001637
Done calling nodal solver

Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885
Timing summary:
Advection  :0.941951385 seconds
MAC Proj   :0.250234152 seconds
Nodal Proj :0.22204457 seconds
Reactions  :1.052441674 seconds
Misc       :0.223220793 seconds
Base State :0.004666076 seconds
Time to advance time step: 2.720050537

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.823376274
Call to estdt for level 0 gives dt_lev = 3.005544288
Call to estdt for level 1 gives dt_lev = 0.6541866529
Minimum estdt over all levels = 0.6541866529
Call to estdt at beginning of step 2 gives dt =0.6541866529
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               = 263576954
MLMG: Initial residual (resid0) = 263576954
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1794640627, 6.808791889e-10
MLMG: Timers: Solve = 0.055141665 Iter = 0.050819362 Bottom = 0.000232792
<<< 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               = 263576953.8
MLMG: Initial residual (resid0) = 716.4811135
MLMG: Final Iter. 2 resid, resid/bnorm = 0.0134591158, 5.106332556e-11
MLMG: Timers: Solve = 0.024999944 Iter = 0.02072132 Bottom = 8.9229e-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               = 14829531.18
MLMG: Initial residual (resid0) = 14829531.18
MLMG: Final Iter. 6 resid, resid/bnorm = 0.001049363939, 7.07617743e-11
MLMG: Timers: Solve = 0.206032168 Iter = 0.201001499 Bottom = 0.000187059
Done calling nodal solver

Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974
Timing summary:
Advection  :0.972391321 seconds
MAC Proj   :0.103200348 seconds
Nodal Proj :0.254184888 seconds
Reactions  :1.053152201 seconds
Misc       :0.223935888 seconds
Base State :0.006900792 seconds
Time to advance time step: 2.607273324

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.827536726
Time to regrid: 0.131604335
Call to estdt for level 0 gives dt_lev = 4.022408781
Call to estdt for level 1 gives dt_lev = 0.8210389333
Minimum estdt over all levels = 0.8210389333
Call to estdt at beginning of step 3 gives dt =0.8210389333
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               = 263106333.8
MLMG: Initial residual (resid0) = 263106333.8
MLMG: Final Iter. 5 resid, resid/bnorm = 0.1839234261, 6.990459844e-10
MLMG: Timers: Solve = 0.054453286 Iter = 0.050105394 Bottom = 0.000222244
<<< 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               = 263106333.5
MLMG: Initial residual (resid0) = 0.256992593
MLMG: No iterations needed
MLMG: Timers: Solve = 0.004351958 Iter = 0 Bottom = 0
<<< 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               = 12313408.17
MLMG: Initial residual (resid0) = 12313408.17
MLMG: Final Iter. 6 resid, resid/bnorm = 0.000564988474, 4.588400434e-11
MLMG: Timers: Solve = 0.208885242 Iter = 0.203982096 Bottom = 0.000195371
Done calling nodal solver

Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571
Timing summary:
Advection  :0.959021852 seconds
MAC Proj   :0.081829437 seconds
Nodal Proj :0.256995203 seconds
Reactions  :1.051255548 seconds
Misc       :0.224365962 seconds
Base State :0.00463425 seconds
Time to advance time step: 2.573871122

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.829527128

Total Time: 20.20659262
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]

AMReX (22.03-10-gfbdfc840de7c) finalized