Initializing CUDA...
CUDA initialized with 1 GPU
AMReX (20.10-30-gc0710a2bc472) 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.480059369
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.013487254 Iter = 0.97268434 Bottom = 0.000480695
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.438087206
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.432871383 Iter = 0.393333488 Bottom = 0.000227684
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.432305682 Iter = 0.392765183 Bottom = 0.000226167
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.636294698 Iter = 0.596805842 Bottom = 0.000339619
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.460746872
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.02237912267, 1.035810464e-10
MLMG: Timers: Solve = 0.078688905 Iter = 0.071039401 Bottom = 0.0072965
<<< 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.004752751237, 2.199795562e-11
MLMG: Timers: Solve = 0.04218823 Iter = 0.036434916 Bottom = 0.003650852
<<< 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.002826735377, 7.91383683e-11
MLMG: Timers: Solve = 5.88393646 Iter = 5.84175113 Bottom = 0.002320383
Done calling nodal solver
Timestep 0 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.720613007 seconds
MAC Proj :0.214362749 seconds
Nodal Proj :6.04387416 seconds
Reactions :0.166547695 seconds
Misc :0.103419479 seconds
Base State :0.04370258 seconds
Time to advance time step: 8.249289789
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.430046127
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.02259197459, 1.045660748e-10
MLMG: Timers: Solve = 0.076612816 Iter = 0.070770933 Bottom = 0.007280015
<<< 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.01930922 Iter = 0.013562587 Bottom = 0.001214839
<<< 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.001655275002, 7.402134071e-11
MLMG: Timers: Solve = 1.932845258 Iter = 1.893257685 Bottom = 0.000854098
Done calling nodal solver
Timestep 1 ends with TIME = 0.0002885721563 DT = 0.0002885721563
Timing summary:
Advection :1.702006205 seconds
MAC Proj :0.18414577 seconds
Nodal Proj :2.09018417 seconds
Reactions :0.165332936 seconds
Misc :0.119879047 seconds
Base State :0.042573226 seconds
Time to advance time step: 4.262035013
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.438576647
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.0223018676, 1.029863784e-10
MLMG: Timers: Solve = 0.076473577 Iter = 0.070690081 Bottom = 0.007261447
<<< 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.019202766 Iter = 0.013469657 Bottom = 0.00121061
<<< 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.001410153694, 8.655855757e-11
MLMG: Timers: Solve = 1.933906696 Iter = 1.894319137 Bottom = 0.000807819
Done calling nodal solver
Timestep 2 ends with TIME = 0.0006060015283 DT = 0.000317429372
Timing summary:
Advection :1.709811141 seconds
MAC Proj :0.183368107 seconds
Nodal Proj :2.08743955 seconds
Reactions :0.165726446 seconds
Misc :0.120203193 seconds
Base State :0.048645654 seconds
Time to advance time step: 4.26702747
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.437832591
Time to regrid: 0.118078087
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.02218363434, 1.024301887e-10
MLMG: Timers: Solve = 0.07688294 Iter = 0.07098333 Bottom = 0.007307449
<<< 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.019269309 Iter = 0.013523363 Bottom = 0.001210279
<<< 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.001057907008, 8.060838618e-11
MLMG: Timers: Solve = 2.145441239 Iter = 2.106162785 Bottom = 0.000867557
Done calling nodal solver
Timestep 3 ends with TIME = 0.0009551738375 DT = 0.0003491723092
Timing summary:
Advection :1.700901906 seconds
MAC Proj :0.184247208 seconds
Nodal Proj :2.299256272 seconds
Reactions :0.163078659 seconds
Misc :0.118318931 seconds
Base State :0.047337509 seconds
Time to advance time step: 4.466255333
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.432355871
Total Time: 31.74316621
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): 9249
[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.10-30-gc0710a2bc472) finalized