Initializing CUDA... CUDA initialized with 1 GPU AMReX (20.04-222-g0f2e7aae8713) 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.411846971 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.0001522888724, 7.316868897e-11 MLMG: Timers: Solve = 1.906336664 Iter = 1.868815939 Bottom = 0.000796116 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.390620093 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.006086156503, 7.783158865e-09 MLMG: Timers: Solve = 0.422415464 Iter = 0.386359257 Bottom = 0.000159792 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.0002423566884, 5.177443181e-10 MLMG: Timers: Solve = 0.42546141 Iter = 0.38637219 Bottom = 0.000175267 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.309106483e-05, 3.674263686e-11 MLMG: Timers: Solve = 0.610015976 Iter = 0.570987596 Bottom = 0.000259018 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.359621614 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.0002033453202, 8.068005352e-11 MLMG: Timers: Solve = 0.101590135 Iter = 0.092014378 Bottom = 0.007281529 <<< 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 = 2520410.785 MLMG: Initial residual (resid0) = 4105.173899 MLMG: Final Iter. 4 resid, resid/bnorm = 0.001121342182, 4.449045326e-10 MLMG: Timers: Solve = 0.069814011 Iter = 0.06213989 Bottom = 0.004931366 <<< 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 = 663975251.3 MLMG: Initial residual (resid0) = 663975251.3 MLMG: Final Iter. 20 resid, resid/bnorm = 0.06375943781, 9.602682885e-11 MLMG: Timers: Solve = 3.758914586 Iter = 3.71847913 Bottom = 0.001524011 Done calling nodal solver Timestep 0 ends with TIME = 0.000289863906 DT = 0.000289863906 Timing summary: Advection :1.621301287 seconds MAC Proj :0.275626542 seconds Nodal Proj :3.890486268 seconds Reactions :0.172844636 seconds Misc :0.089158977 seconds Base State :0.034247908 seconds Time to advance time step: 6.049916019 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.353783461 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 = 2520411.47 MLMG: Initial residual (resid0) = 2520411.47 MLMG: Final Iter. 6 resid, resid/bnorm = 0.000373349234, 1.48130271e-10 MLMG: Timers: Solve = 0.102273214 Iter = 0.094390625 Bottom = 0.007442071 <<< 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 = 2520430.895 MLMG: Initial residual (resid0) = 6182.069347 MLMG: Final Iter. 4 resid, resid/bnorm = 0.00116662588, 4.628676318e-10 MLMG: Timers: Solve = 0.070281594 Iter = 0.062534892 Bottom = 0.004775682 <<< 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 = 467307.9693 MLMG: Initial residual (resid0) = 467307.9693 MLMG: Final Iter. 18 resid, resid/bnorm = 3.894155816e-05, 8.333168e-11 MLMG: Timers: Solve = 3.383426839 Iter = 3.34472625 Bottom = 0.00134644 Done calling nodal solver Timestep 1 ends with TIME = 0.000289863906 DT = 0.000289863906 Timing summary: Advection :1.588694587 seconds MAC Proj :0.271690004 seconds Nodal Proj :3.512683976 seconds Reactions :0.171787978 seconds Misc :0.08433414 seconds Base State :0.032705308 seconds Time to advance time step: 5.629725897 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.357583129 Call to estdt for level 0 gives dt_lev = 2.461525899 Call to estdt for level 1 gives dt_lev = 0.5825175922 Minimum estdt over all levels = 0.5825175922 Call to estdt at beginning of step 2 gives dt =0.5825175922 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 = 2482166.036 MLMG: Initial residual (resid0) = 2482166.036 MLMG: Final Iter. 6 resid, resid/bnorm = 0.0005592140078, 2.252927482e-10 MLMG: Timers: Solve = 0.100396491 Iter = 0.092633819 Bottom = 0.007362425 <<< 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 = 2482180.367 MLMG: Initial residual (resid0) = 2433.330164 MLMG: Final Iter. 4 resid, resid/bnorm = 0.0007921425859, 3.191317587e-10 MLMG: Timers: Solve = 0.070124209 Iter = 0.062604531 Bottom = 0.004766721 <<< 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 = 469052.6112 MLMG: Initial residual (resid0) = 469052.6112 MLMG: Final Iter. 14 resid, resid/bnorm = 4.605368304e-05, 9.818447214e-11 MLMG: Timers: Solve = 2.644156902 Iter = 2.605514074 Bottom = 0.001094521 Done calling nodal solver Timestep 2 ends with TIME = 0.0006087142026 DT = 0.0003188502966 Timing summary: Advection :1.58161998 seconds MAC Proj :0.268988808 seconds Nodal Proj :2.775335218 seconds Reactions :0.174451479 seconds Misc :0.07732402 seconds Base State :0.033687563 seconds Time to advance time step: 4.878262558 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.377045313 Time to regrid: 0.05293736 Call to estdt for level 0 gives dt_lev = 2.549012516 Call to estdt for level 1 gives dt_lev = 0.6638997325 Minimum estdt over all levels = 0.6638997325 Call to estdt at beginning of step 3 gives dt =0.6638997325 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 = 2476181.06 MLMG: Initial residual (resid0) = 2476181.06 MLMG: Final Iter. 5 resid, resid/bnorm = 0.002411237452, 9.73772674e-10 MLMG: Timers: Solve = 0.085295795 Iter = 0.077623045 Bottom = 0.006148172 <<< 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 = 2476165.294 MLMG: Initial residual (resid0) = 3455.281975 MLMG: Final Iter. 4 resid, resid/bnorm = 0.0005412152968, 2.185699388e-10 MLMG: Timers: Solve = 0.07044147 Iter = 0.062933587 Bottom = 0.004950228 <<< 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 = 472750.3616 MLMG: Initial residual (resid0) = 472750.3616 MLMG: Final Iter. 13 resid, resid/bnorm = 3.50238206e-05, 7.408523281e-11 MLMG: Timers: Solve = 2.459299728 Iter = 2.420716032 Bottom = 0.000989233 Done calling nodal solver Timestep 3 ends with TIME = 0.0009594495289 DT = 0.0003507353263 Timing summary: Advection :1.56892405 seconds MAC Proj :0.255049482 seconds Nodal Proj :2.587404968 seconds Reactions :0.170085474 seconds Misc :0.078714647 seconds Base State :0.033943868 seconds Time to advance time step: 4.660702418 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.394194279 Total Time: 31.90779937 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-222-g0f2e7aae8713) finalized