Initializing CUDA... CUDA initialized with 1 GPU AMReX (20.04-249-gdccaf5b22f14) 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 ... ERROR: problem in the namelist Calling VariableSetup() Initializing Helmholtz EOS and using Coulomb corrections. Calling set_method_params() Calling BCSetup() Calling init_base_state_geometry() 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 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.392821901 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 = 2081339.362 MLMG: Initial residual (resid0) = 2081339.362 MLMG: Final Iter. 10 resid, resid/bnorm = 0.0001522888713, 7.316868842e-11 MLMG: Timers: Solve = 1.930391027 Iter = 1.892509855 Bottom = 0.00072644 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.396582268 Call to firstdt for level 0 gives dt_lev = 0.005815502985 Multiplying dt_lev by init_shrink; dt_lev = 0.0005815502985 Call to firstdt for level 1 gives dt_lev = 0.00289863906 Multiplying dt_lev by init_shrink; dt_lev = 0.000289863906 Minimum firstdt over all levels = 0.000289863906 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 781964.8306 MLMG: Initial residual (resid0) = 781964.8306 MLMG: Final Iter. 2 resid, resid/bnorm = 0.006086156514, 7.783158879e-09 MLMG: Timers: Solve = 0.428608197 Iter = 0.392309829 Bottom = 0.000166313 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 1.19799377 Call to estdt for level 1 gives dt_lev = 0.4128889473 Minimum estdt over all levels = 0.4128889473 Call to estdt at end of istep_divu_iter = 1 gives dt = 0.4128889473 Multiplying dt by init_shrink; dt = 0.04128889473 Ignoring this new dt since it's larger than the previous dt = 0.000289863906 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 468101.107 MLMG: Initial residual (resid0) = 468101.107 MLMG: Final Iter. 2 resid, resid/bnorm = 0.0002423566957, 5.177443336e-10 MLMG: Timers: Solve = 0.435606118 Iter = 0.396229832 Bottom = 0.000183287 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 1.19799377 Call to estdt for level 1 gives dt_lev = 0.4128889473 Minimum estdt over all levels = 0.4128889473 Call to estdt at end of istep_divu_iter = 2 gives dt = 0.4128889473 Multiplying dt by init_shrink; dt = 0.04128889473 Ignoring this new dt since it's larger than the previous dt = 0.000289863906 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 356290.8367 MLMG: Initial residual (resid0) = 356290.8367 MLMG: Final Iter. 3 resid, resid/bnorm = 1.309105392e-05, 3.674260623e-11 MLMG: Timers: Solve = 0.639658314 Iter = 0.60039684 Bottom = 0.000253156 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 1.19799377 Call to estdt for level 1 gives dt_lev = 0.4128889473 Minimum estdt over all levels = 0.4128889473 Call to estdt at end of istep_divu_iter = 3 gives dt = 0.4128889473 Multiplying dt by init_shrink; dt = 0.04128889473 Ignoring this new dt since it's larger than the previous dt = 0.000289863906 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.372456808 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 0.000289863906 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 = 2520391.488 MLMG: Initial residual (resid0) = 2520391.488 MLMG: Final Iter. 6 resid, resid/bnorm = 0.0002033452038, 8.068000733e-11 MLMG: Timers: Solve = 0.102346433 Iter = 0.092876104 Bottom = 0.007352433 <<< 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 = 2520391.394 MLMG: Initial residual (resid0) = 3.099700325 MLMG: Final Iter. 2 resid, resid/bnorm = 4.180637188e-05, 1.658725386e-11 MLMG: Timers: Solve = 0.04055721 Iter = 0.033065991 Bottom = 0.002464884 <<< 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 = 663704574.4 MLMG: Initial residual (resid0) = 663704574.4 MLMG: Final Iter. 11 resid, resid/bnorm = 0.06578297998, 9.911485098e-11 MLMG: Timers: Solve = 2.123266958 Iter = 2.08284544 Bottom = 0.000857583 Done calling nodal solver Timestep 0 ends with TIME = 0.000289863906 DT = 0.000289863906 Timing summary: Advection :1.693340235 seconds MAC Proj :0.249158228 seconds Nodal Proj :2.259371575 seconds Reactions :0.175240255 seconds Misc :0.093526856 seconds Base State :0.039702461 seconds Time to advance time step: 4.471154497 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.418432678 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 0.000289863906 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 = 2520429.043 MLMG: Initial residual (resid0) = 2520429.043 MLMG: Final Iter. 6 resid, resid/bnorm = 0.000212295563, 8.422993045e-11 MLMG: Timers: Solve = 0.100515352 Iter = 0.09289242 Bottom = 0.007354168 <<< 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 = 2520428.948 MLMG: Initial residual (resid0) = 3.108290578 MLMG: Final Iter. 2 resid, resid/bnorm = 4.589231685e-05, 1.82081375e-11 MLMG: Timers: Solve = 0.040300396 Iter = 0.032861181 Bottom = 0.002431811 <<< 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 = 467226.6822 MLMG: Initial residual (resid0) = 467226.6822 MLMG: Final Iter. 9 resid, resid/bnorm = 3.914281115e-05, 8.377691737e-11 MLMG: Timers: Solve = 1.741776705 Iter = 1.7024174 Bottom = 0.000702646 Done calling nodal solver Timestep 1 ends with TIME = 0.000289863906 DT = 0.000289863906 Timing summary: Advection :1.603915032 seconds MAC Proj :0.241877372 seconds Nodal Proj :1.87330774 seconds Reactions :0.173313788 seconds Misc :0.09539482 seconds Base State :0.036745568 seconds Time to advance time step: 3.988291122 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.351529216 Call to estdt for level 0 gives dt_lev = 2.45872793 Call to estdt for level 1 gives dt_lev = 0.5822724056 Minimum estdt over all levels = 0.5822724056 Call to estdt at beginning of step 2 gives dt =0.5822724056 dt_growth factor limits the new dt = 0.0003188502966 Timestep 2 starts with TIME = 0.000289863906 DT = 0.0003188502966 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 = 2482204.383 MLMG: Initial residual (resid0) = 2482204.383 MLMG: Final Iter. 5 resid, resid/bnorm = 0.002290578326, 9.228000491e-10 MLMG: Timers: Solve = 0.08539487 Iter = 0.077818774 Bottom = 0.006111079 <<< 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 = 112216276.3 MLMG: Initial residual (resid0) = 4947703770 MLMG: Final Iter. 5 resid, resid/resid0 = 4.8988773, 9.901314889e-10 MLMG: Timers: Solve = 0.08931147 Iter = 0.081873353 Bottom = 0.006111776 <<< 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 = 469119.7189 MLMG: Initial residual (resid0) = 469119.7189 MLMG: Final Iter. 20 resid, resid/bnorm = 3.768790884e-05, 8.03375073e-11 MLMG: Timers: Solve = 3.803894762 Iter = 3.764774681 Bottom = 0.001436764 Done calling nodal solver Timestep 2 ends with TIME = 0.0006087142026 DT = 0.0003188502966 Timing summary: Advection :1.556832161 seconds MAC Proj :0.274580955 seconds Nodal Proj :3.94128035 seconds Reactions :0.178385192 seconds Misc :0.082179024 seconds Base State :0.036494975 seconds Time to advance time step: 6.033806312 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.356734601 Time to regrid: 0.05609651 Call to estdt for level 0 gives dt_lev = 2.543837552 Call to estdt for level 1 gives dt_lev = 0.6636224091 Minimum estdt over all levels = 0.6636224091 Call to estdt at beginning of step 3 gives dt =0.6636224091 dt_growth factor limits the new dt = 0.0003507353263 Timestep 3 starts with TIME = 0.0006087142026 DT = 0.0003507353263 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 = 2476270.349 MLMG: Initial residual (resid0) = 2476270.349 MLMG: Final Iter. 6 resid, resid/bnorm = 0.0009233292658, 3.728709453e-10 MLMG: Timers: Solve = 0.100517652 Iter = 0.092902646 Bottom = 0.007145072 <<< 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 = 2476270.424 MLMG: Initial residual (resid0) = 3.663909508 MLMG: Final Iter. 2 resid, resid/bnorm = 6.453310902e-05, 2.606060646e-11 MLMG: Timers: Solve = 0.040008814 Iter = 0.032609058 Bottom = 0.002276828 <<< 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 = 472726.3531 MLMG: Initial residual (resid0) = 472726.3531 MLMG: Final Iter. 20 resid, resid/bnorm = 3.542462036e-05, 7.49368427e-11 MLMG: Timers: Solve = 3.80596085 Iter = 3.765883661 Bottom = 0.001482871 Done calling nodal solver Timestep 3 ends with TIME = 0.0009594495289 DT = 0.0003507353263 Timing summary: Advection :1.655231499 seconds MAC Proj :0.242133533 seconds Nodal Proj :3.942076541 seconds Reactions :0.171972135 seconds Misc :0.084967401 seconds Base State :0.037924406 seconds Time to advance time step: 6.096967636 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.353565658 Total Time: 31.39718765 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): 9257 [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.04-249-gdccaf5b22f14) finalized