Initializing CUDA... CUDA initialized with 1 GPU AMReX (20.04-243-gbce44ba71508) 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 InitBaseStateGeometry() 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.391722376 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.0001522888733, 7.31686894e-11 MLMG: Timers: Solve = 1.935680855 Iter = 1.897586089 Bottom = 0.000762493 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.387172672 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.428592444 Iter = 0.392133607 Bottom = 0.000156307 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.000242356683, 5.177443064e-10 MLMG: Timers: Solve = 0.434180726 Iter = 0.394877544 Bottom = 0.000184176 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.309106028e-05, 3.67426241e-11 MLMG: Timers: Solve = 0.622772735 Iter = 0.583281702 Bottom = 0.000255688 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.362207817 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.0002033451456, 8.067998424e-11 MLMG: Timers: Solve = 0.101105325 Iter = 0.091923791 Bottom = 0.007318374 <<< 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 = 2520401.464 MLMG: Initial residual (resid0) = 1493.782876 MLMG: Final Iter. 3 resid, resid/bnorm = 0.002337134676, 9.272866682e-10 MLMG: Timers: Solve = 0.054906807 Iter = 0.047639821 Bottom = 0.003687104 <<< 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 = 663817812.2 MLMG: Initial residual (resid0) = 663817812.2 MLMG: Final Iter. 18 resid, resid/bnorm = 0.06360938677, 9.582356122e-11 MLMG: Timers: Solve = 3.45044934 Iter = 3.409356409 Bottom = 0.001404892 Done calling nodal solver Timestep 0 ends with TIME = 0.000289863906 DT = 0.000289863906 Timing summary: Advection :1.7211055 seconds MAC Proj :0.259071821 seconds Nodal Proj :3.5873013 seconds Reactions :0.17739381 seconds Misc :0.092632773 seconds Base State :0.035897693 seconds Time to advance time step: 5.837991655 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.362779077 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 = 2520420.967 MLMG: Initial residual (resid0) = 2520420.967 MLMG: Final Iter. 6 resid, resid/bnorm = 0.0002864194685, 1.136395357e-10 MLMG: Timers: Solve = 0.100391029 Iter = 0.092839844 Bottom = 0.007365637 <<< 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 = 2520404.547 MLMG: Initial residual (resid0) = 5795.308868 MLMG: Final Iter. 4 resid, resid/bnorm = 0.001145267161, 4.543981489e-10 MLMG: Timers: Solve = 0.070106406 Iter = 0.062518277 Bottom = 0.004888502 <<< 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 = 467308.3034 MLMG: Initial residual (resid0) = 467308.3034 MLMG: Final Iter. 18 resid, resid/bnorm = 3.88992612e-05, 8.324110852e-11 MLMG: Timers: Solve = 3.445455104 Iter = 3.405947007 Bottom = 0.00135688 Done calling nodal solver Timestep 1 ends with TIME = 0.000289863906 DT = 0.000289863906 Timing summary: Advection :1.668581231 seconds MAC Proj :0.272098905 seconds Nodal Proj :3.582868581 seconds Reactions :0.174924992 seconds Misc :0.094112033 seconds Base State :0.036449713 seconds Time to advance time step: 5.793113432 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.365237273 Call to estdt for level 0 gives dt_lev = 2.461551397 Call to estdt for level 1 gives dt_lev = 0.5825185932 Minimum estdt over all levels = 0.5825185932 Call to estdt at beginning of step 2 gives dt =0.5825185932 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 = 2482165.465 MLMG: Initial residual (resid0) = 2482165.465 MLMG: Final Iter. 6 resid, resid/bnorm = 0.0005738855689, 2.31203591e-10 MLMG: Timers: Solve = 0.099344069 Iter = 0.091736372 Bottom = 0.007343132 <<< 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 = 2482179.415 MLMG: Initial residual (resid0) = 2433.708125 MLMG: Final Iter. 4 resid, resid/bnorm = 0.0007517725462, 3.02867932e-10 MLMG: Timers: Solve = 0.069625782 Iter = 0.062310573 Bottom = 0.004928208 <<< 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 = 469051.8379 MLMG: Initial residual (resid0) = 469051.8379 MLMG: Final Iter. 15 resid, resid/bnorm = 3.682986426e-05, 7.851981654e-11 MLMG: Timers: Solve = 2.877570582 Iter = 2.837942511 Bottom = 0.001162484 Done calling nodal solver Timestep 2 ends with TIME = 0.0006087142026 DT = 0.0003188502966 Timing summary: Advection :1.639868261 seconds MAC Proj :0.269428614 seconds Nodal Proj :3.014313201 seconds Reactions :0.17637156 seconds Misc :0.081573839 seconds Base State :0.036829617 seconds Time to advance time step: 5.182025368 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.399377227 Time to regrid: 0.053687252 Call to estdt for level 0 gives dt_lev = 2.549071865 Call to estdt for level 1 gives dt_lev = 0.6639028816 Minimum estdt over all levels = 0.6639028816 Call to estdt at beginning of step 3 gives dt =0.6639028816 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 = 2476179.94 MLMG: Initial residual (resid0) = 2476179.94 MLMG: Final Iter. 5 resid, resid/bnorm = 0.002378463745, 9.605375228e-10 MLMG: Timers: Solve = 0.097126739 Iter = 0.088156655 Bottom = 0.006491522 <<< 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 = 2476171.304 MLMG: Initial residual (resid0) = 2795.767972 MLMG: Final Iter. 3 resid, resid/bnorm = 0.001734971149, 7.006668506e-10 MLMG: Timers: Solve = 0.055291472 Iter = 0.04758689 Bottom = 0.003658829 <<< 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.7616 MLMG: Initial residual (resid0) = 472726.7616 MLMG: Final Iter. 15 resid, resid/bnorm = 4.326957605e-05, 9.153189447e-11 MLMG: Timers: Solve = 2.883479613 Iter = 2.843356862 Bottom = 0.001118181 Done calling nodal solver Timestep 3 ends with TIME = 0.0009594495289 DT = 0.0003507353263 Timing summary: Advection :1.738855575 seconds MAC Proj :0.26218864 seconds Nodal Proj :3.018284959 seconds Reactions :0.186658308 seconds Misc :0.08630828 seconds Base State :0.040446478 seconds Time to advance time step: 5.292816396 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.358571915 Total Time: 32.84893532 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): 9256 [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-243-gbce44ba71508) finalized