Initializing CUDA... CUDA initialized with 1 GPU AMReX (21.12-30-ga38ccbc90843) 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() Calling BCSetup() 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 WARNING: variable not found: ash WARNING: magnesium-24 not provided in inputs file 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.414855451 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 = 223792109.4 MLMG: Initial residual (resid0) = 223792109.4 MLMG: Final Iter. 4 resid, resid/bnorm = 0.01715009583, 7.66340506e-11 MLMG: Timers: Solve = 1.313118348 Iter = 1.251465389 Bottom = 0.000257876 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.34831205 Call to firstdt for level 0 gives dt_lev = 0.005789736987 Multiplying dt_lev by init_shrink; dt_lev = 0.0005789736987 Call to firstdt for level 1 gives dt_lev = 0.002885750885 Multiplying dt_lev by init_shrink; dt_lev = 0.0002885750885 Minimum firstdt over all levels = 0.0002885750885 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 96748113.54 MLMG: Initial residual (resid0) = 96748113.54 MLMG: Final Iter. 1 resid, resid/bnorm = 10.35489981, 1.070294751e-07 MLMG: Timers: Solve = 0.384495245 Iter = 0.324262961 Bottom = 6.2001e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 11.21229293 Minimum estdt over all levels = 11.21229293 Call to estdt at end of istep_divu_iter = 1 gives dt = 11.21229293 Multiplying dt by init_shrink; dt = 1.121229293 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 53052389.3 MLMG: Initial residual (resid0) = 53052389.3 MLMG: Final Iter. 1 resid, resid/bnorm = 0.2186066424, 4.12058053e-09 MLMG: Timers: Solve = 0.383991909 Iter = 0.325428134 Bottom = 6.3459e-05 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 11.22745362 Minimum estdt over all levels = 11.22745362 Call to estdt at end of istep_divu_iter = 2 gives dt = 11.22745362 Multiplying dt by init_shrink; dt = 1.122745362 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 33080017.02 MLMG: Initial residual (resid0) = 33080017.02 MLMG: Final Iter. 3 resid, resid/bnorm = 0.002709050849, 8.189387713e-11 MLMG: Timers: Solve = 1.012892308 Iter = 0.953381486 Bottom = 0.000185751 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 21.60061365 Call to estdt for level 1 gives dt_lev = 11.20744185 Minimum estdt over all levels = 11.20744185 Call to estdt at end of istep_divu_iter = 3 gives dt = 11.20744185 Multiplying dt by init_shrink; dt = 1.120744185 Ignoring this new dt since it's larger than the previous dt = 0.0002885750885 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.365865353 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 0.0002885750885 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 = 216054180.2 MLMG: Initial residual (resid0) = 216054180.2 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02244902775, 1.039046212e-10 MLMG: Timers: Solve = 0.058443691 Iter = 0.053251538 Bottom = 0.007302247 <<< 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 = 216054180 MLMG: Initial residual (resid0) = 4063.444174 MLMG: Final Iter. 2 resid, resid/bnorm = 0.04702966002, 2.176753073e-10 MLMG: Timers: Solve = 0.023228376 Iter = 0.01872256 Bottom = 0.002436055 <<< 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 = 14621038.7 MLMG: Initial residual (resid0) = 14621038.7 MLMG: Final Iter. 19 resid, resid/bnorm = 0.0008469979512, 5.793008065e-11 MLMG: Timers: Solve = 5.954711522 Iter = 5.894147278 Bottom = 0.001956721 Done calling nodal solver Timestep 0 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :2.158972838 seconds MAC Proj :0.186291549 seconds Nodal Proj :6.162947898 seconds Reactions :0.242224809 seconds Misc :0.093318546 seconds Base State :0.062483378 seconds Time to advance time step: 8.844322308 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.335203562 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 0.0002885750885 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 = 216054489.9 MLMG: Initial residual (resid0) = 216054489.9 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02244241908, 1.038738843e-10 MLMG: Timers: Solve = 0.05772077 Iter = 0.053002578 Bottom = 0.00734454 <<< 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 = 216054489.7 MLMG: Initial residual (resid0) = 0.3527048053 MLMG: Final Iter. 1 resid, resid/bnorm = 0.02240421663, 1.036970658e-10 MLMG: Timers: Solve = 0.015535709 Iter = 0.010631474 Bottom = 0.001232465 <<< 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 = 22362103.76 MLMG: Initial residual (resid0) = 22362103.76 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001080892519, 4.833590481e-11 MLMG: Timers: Solve = 1.632422074 Iter = 1.572267083 Bottom = 0.000628337 Done calling nodal solver Timestep 1 ends with TIME = 0.0002885750885 DT = 0.0002885750885 Timing summary: Advection :2.200189185 seconds MAC Proj :0.171069216 seconds Nodal Proj :1.84154408 seconds Reactions :0.240677718 seconds Misc :0.087275555 seconds Base State :0.062582838 seconds Time to advance time step: 4.541299673 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.356324052 Call to estdt for level 0 gives dt_lev = 2.560516355 Call to estdt for level 1 gives dt_lev = 0.5423247709 Minimum estdt over all levels = 0.5423247709 Call to estdt at beginning of step 2 gives dt =0.5423247709 dt_growth factor limits the new dt = 0.0003174325974 Timestep 2 starts with TIME = 0.0002885750885 DT = 0.0003174325974 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 = 216551218.7 MLMG: Initial residual (resid0) = 216551218.7 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02222341299, 1.026242804e-10 MLMG: Timers: Solve = 0.056831266 Iter = 0.052270783 Bottom = 0.007320958 <<< 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 = 216551218.6 MLMG: Initial residual (resid0) = 0.4579974225 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01965397086, 9.075899453e-11 MLMG: Timers: Solve = 0.014963831 Iter = 0.010506974 Bottom = 0.001205048 <<< 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 = 16291309.39 MLMG: Initial residual (resid0) = 16291309.39 MLMG: Final Iter. 5 resid, resid/bnorm = 0.001264760614, 7.763406759e-11 MLMG: Timers: Solve = 1.639893073 Iter = 1.580369668 Bottom = 0.000611378 Done calling nodal solver Timestep 2 ends with TIME = 0.0006060076859 DT = 0.0003174325974 Timing summary: Advection :2.177720981 seconds MAC Proj :0.174429318 seconds Nodal Proj :1.845592268 seconds Reactions :0.23528598 seconds Misc :0.083501075 seconds Base State :0.063365967 seconds Time to advance time step: 4.517106376 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.32507609 Time to regrid: 0.073055187 Call to estdt for level 0 gives dt_lev = 3.523509294 Call to estdt for level 1 gives dt_lev = 0.6591917221 Minimum estdt over all levels = 0.6591917221 Call to estdt at beginning of step 3 gives dt =0.6591917221 dt_growth factor limits the new dt = 0.0003491758571 Timestep 3 starts with TIME = 0.0006060076859 DT = 0.0003491758571 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 = 216572797.4 MLMG: Initial residual (resid0) = 216572797.4 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02191495895, 1.011898042e-10 MLMG: Timers: Solve = 0.057684729 Iter = 0.053094079 Bottom = 0.007290582 <<< 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 = 216572797.2 MLMG: Initial residual (resid0) = 0.2694439256 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01715105083, 7.919300601e-11 MLMG: Timers: Solve = 0.014657038 Iter = 0.010243138 Bottom = 0.00125234 <<< 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 = 13124019.97 MLMG: Initial residual (resid0) = 13124019.97 MLMG: Final Iter. 6 resid, resid/bnorm = 0.0005545876167, 4.225744993e-11 MLMG: Timers: Solve = 1.958279251 Iter = 1.898738596 Bottom = 0.000670087 Done calling nodal solver Timestep 3 ends with TIME = 0.0009551835431 DT = 0.0003491758571 Timing summary: Advection :2.204333664 seconds MAC Proj :0.175121279 seconds Nodal Proj :2.166679253 seconds Reactions :0.240991052 seconds Misc :0.088304437 seconds Base State :0.06174471 seconds Time to advance time step: 4.876014605 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.344616905 Total Time: 31.64015365 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): 10130 [The Arena] space allocated (MB): 9049 [The Arena] space used (MB): 0 [The Pinned Arena] space allocated (MB): 33 [The Pinned Arena] space used (MB): 0 AMReX (21.12-30-ga38ccbc90843) finalized