Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.08-92-gd48f68f1a641) 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()
Initializing Helmholtz EOS and using Coulomb corrections.
Calling set_method_params()
Calling BCSetup()
Calling BaseStateGeometry::Init()
Calling init_base_state_geometry()
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
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.473135536
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 = 223792658.9
MLMG: Initial residual (resid0) = 223792658.9
MLMG: Final Iter. 5 resid, resid/bnorm = 0.01302976906, 5.822250438e-11
MLMG: Timers: Solve = 1.000394441 Iter = 0.959212924 Bottom = 0.000404507
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.436317278
Call to firstdt for level 0 gives dt_lev = 0.005789678157
Multiplying dt_lev by init_shrink; dt_lev = 0.0005789678157
Call to firstdt for level 1 gives dt_lev = 0.002885721563
Multiplying dt_lev by init_shrink; dt_lev = 0.0002885721563
Minimum firstdt over all levels = 0.0002885721563
Doing initial divu iteration #1
Calling nodal solver
MLMG: Initial rhs = 96748240.18
MLMG: Initial residual (resid0) = 96748240.18
MLMG: Final Iter. 2 resid, resid/bnorm = 0.7442735399, 7.692889695e-09
MLMG: Timers: Solve = 0.425131635 Iter = 0.38663921 Bottom = 0.000168108
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 5.916319527
Call to estdt for level 1 gives dt_lev = 2.341804478
Minimum estdt over all levels = 2.341804478
Call to estdt at end of istep_divu_iter = 1 gives dt = 2.341804478
Multiplying dt by init_shrink; dt = 0.2341804478
Ignoring this new dt since it's larger than the previous dt = 0.0002885721563
Doing initial divu iteration #2
Calling nodal solver
MLMG: Initial rhs = 53052451.02
MLMG: Initial residual (resid0) = 53052451.02
MLMG: Final Iter. 2 resid, resid/bnorm = 0.01898677269, 3.578868144e-10
MLMG: Timers: Solve = 0.425690375 Iter = 0.387199442 Bottom = 0.000170391
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 5.916319527
Call to estdt for level 1 gives dt_lev = 2.341804478
Minimum estdt over all levels = 2.341804478
Call to estdt at end of istep_divu_iter = 2 gives dt = 2.341804478
Multiplying dt by init_shrink; dt = 0.2341804478
Ignoring this new dt since it's larger than the previous dt = 0.0002885721563
Doing initial divu iteration #3
Calling nodal solver
MLMG: Initial rhs = 33080059.37
MLMG: Initial residual (resid0) = 33080059.37
MLMG: Final Iter. 3 resid, resid/bnorm = 0.001913452566, 5.784308136e-11
MLMG: Timers: Solve = 0.626452951 Iter = 0.587952785 Bottom = 0.000255958
Done calling nodal solver
Call to estdt for level 0 gives dt_lev = 5.916319527
Call to estdt for level 1 gives dt_lev = 2.341804478
Minimum estdt over all levels = 2.341804478
Call to estdt at end of istep_divu_iter = 3 gives dt = 2.341804478
Multiplying dt by init_shrink; dt = 0.2341804478
Ignoring this new dt since it's larger than the previous dt = 0.0002885721563
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.43972457
Doing initial pressure iteration #1
Timestep 0 starts with TIME = 0 DT = 0.0002885721563
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 = 216054224.6
MLMG: Initial residual (resid0) = 216054224.6
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02237911895, 1.035810292e-10
MLMG: Timers: Solve = 0.10267438 Iter = 0.092983896 Bottom = 0.007412299
<<< 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 = 216054224.3
MLMG: Initial residual (resid0) = 10665.17408
MLMG: Final Iter. 3 resid, resid/bnorm = 0.004752751117, 2.199795506e-11
MLMG: Timers: Solve = 0.0554889 Iter = 0.047889919 Bottom = 0.003733679
<<< 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 = 35718898.91
MLMG: Initial residual (resid0) = 35718898.91
MLMG: Final Iter. 31 resid, resid/bnorm = 0.002826690674, 7.913711676e-11
MLMG: Timers: Solve = 5.815004401 Iter = 5.774991544 Bottom = 0.002278558
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.94762168 seconds
MAC Proj :0.265082129 seconds
Nodal Proj :5.981274843 seconds
Reactions :0.175258157 seconds
Misc :0.104587882 seconds
Base State :0.042265706 seconds
Time to advance time step: 8.474310194
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.440138118
Calling Evolve()
Timestep 1 starts with TIME = 0 DT = 0.0002885721563
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 = 216054534.3
MLMG: Initial residual (resid0) = 216054534.3
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02259196341, 1.045660231e-10
MLMG: Timers: Solve = 0.100549954 Iter = 0.092934234 Bottom = 0.007380401
<<< 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 = 216054534.1
MLMG: Initial residual (resid0) = 0.35175923
MLMG: Final Iter. 1 resid, resid/bnorm = 0.02233633684, 1.033828655e-10
MLMG: Timers: Solve = 0.025353541 Iter = 0.017841762 Bottom = 0.001236125
<<< 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 = 22362132.14
MLMG: Initial residual (resid0) = 22362132.14
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001655276865, 7.402142401e-11
MLMG: Timers: Solve = 1.909288996 Iter = 1.870870097 Bottom = 0.000787772
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.931932613 seconds
MAC Proj :0.226621311 seconds
Nodal Proj :2.074231153 seconds
Reactions :0.173415859 seconds
Misc :0.122296022 seconds
Base State :0.04362206 seconds
Time to advance time step: 4.52899378
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.445010369
Call to estdt for level 0 gives dt_lev = 2.560503025
Call to estdt for level 1 gives dt_lev = 0.5423221224
Minimum estdt over all levels = 0.5423221224
Call to estdt at beginning of step 2 gives dt =0.5423221224
dt_growth factor limits the new dt = 0.000317429372
Timestep 2 starts with TIME = 0.0002885721563 DT = 0.000317429372
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 = 216551625.1
MLMG: Initial residual (resid0) = 216551625.1
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02230187505, 1.029864128e-10
MLMG: Timers: Solve = 0.100468971 Iter = 0.092832839 Bottom = 0.00737869
<<< 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 = 216551624.9
MLMG: Initial residual (resid0) = 0.3069432653
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01958142943, 9.042383976e-11
MLMG: Timers: Solve = 0.025575245 Iter = 0.01787282 Bottom = 0.001248883
<<< 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 = 16291326.17
MLMG: Initial residual (resid0) = 16291326.17
MLMG: Final Iter. 10 resid, resid/bnorm = 0.001410155557, 8.65586719e-11
MLMG: Timers: Solve = 1.908249169 Iter = 1.869658944 Bottom = 0.000769205
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060015283 DT = 0.000317429372
Timing summary:
Advection :1.923620336 seconds
MAC Proj :0.226864238 seconds
Nodal Proj :2.072449029 seconds
Reactions :0.17441979 seconds
Misc :0.122979599 seconds
Base State :0.046546886 seconds
Time to advance time step: 4.520822791
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.450060137
Time to regrid: 0.125634767
Call to estdt for level 0 gives dt_lev = 3.52349079
Call to estdt for level 1 gives dt_lev = 0.6591885313
Minimum estdt over all levels = 0.6591885313
Call to estdt at beginning of step 3 gives dt =0.6591885313
dt_growth factor limits the new dt = 0.0003491723092
Timestep 3 starts with TIME = 0.0006060015283 DT = 0.0003491723092
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 = 216573205.8
MLMG: Initial residual (resid0) = 216573205.8
MLMG: Final Iter. 6 resid, resid/bnorm = 0.02218363807, 1.024302059e-10
MLMG: Timers: Solve = 0.100123804 Iter = 0.09264518 Bottom = 0.007428352
<<< 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 = 216573205.6
MLMG: Initial residual (resid0) = 0.2695824229
MLMG: Final Iter. 1 resid, resid/bnorm = 0.01713669178, 7.912655553e-11
MLMG: Timers: Solve = 0.02544758 Iter = 0.017782645 Bottom = 0.00123752
<<< 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 = 13124031.6
MLMG: Initial residual (resid0) = 13124031.6
MLMG: Final Iter. 11 resid, resid/bnorm = 0.001057907939, 8.060845714e-11
MLMG: Timers: Solve = 2.093031145 Iter = 2.054470414 Bottom = 0.000849017
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551738375 DT = 0.0003491723092
Timing summary:
Advection :1.932108458 seconds
MAC Proj :0.228562924 seconds
Nodal Proj :2.258572977 seconds
Reactions :0.175053552 seconds
Misc :0.122422569 seconds
Base State :0.043059688 seconds
Time to advance time step: 4.717205213
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.496268305
Total Time: 32.72186055
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): 12066
Free GPU global memory (MB): 9386
[The Arena] space allocated (MB): 9049
[The Arena] space used (MB): 0
[The Device Arena] space allocated (MB): 8
[The Device Arena] space used (MB): 0
[The Managed Arena] space allocated (MB): 8
[The Managed Arena] space used (MB): 0
[The Pinned Arena] space allocated (MB): 8
[The Pinned Arena] space used (MB): 0
AMReX (20.08-92-gd48f68f1a641) finalized