Initializing CUDA... CUDA initialized with 1 GPU AMReX (20.06-79-g9b2172017f8a) 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 = model.hse.cool.coulomb model file = model.hse.cool.coulomb reading initial model 640 points found in the initial model 6 variables found in the initial model model file mapping, level: 0 dr of MAESTRO base state = 562500.000000 dr of input file data = 562500.000000 maximum radius (cell-centered) of input model = 359718750.000000 setting r_cutoff to 365 radius at r_cutoff 205593750 Maximum HSE Error = 0.000061 (after putting initial model into base state arrays, and for density < base_cutoff_density) Writing plotfile reacting_bubble-3d_pltInitData after InitData inner sponge: r_sp , r_tp : 186468750, 224718750 Time to write plotfile: 0.59031604 inner sponge: r_sp , r_tp : 186468750, 224718750 Doing initial projection Calling nodal solver MLMG: Initial rhs = 0 MLMG: Initial residual (resid0) = 0 MLMG: No iterations needed MLMG: Timers: Solve = 0.07007434 Iter = 0 Bottom = 0 Done calling nodal solver Writing plotfile reacting_bubble-3d_pltafter_InitProj after InitProj inner sponge: r_sp , r_tp : 186468750, 224718750 Time to write plotfile: 0.513130331 Call to firstdt for level 0 gives dt_lev = 0.0004741795139 Multiplying dt_lev by init_shrink; dt_lev = 0.0004741795139 Minimum firstdt over all levels = 0.0004741795139 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 2054.812665 MLMG: Initial residual (resid0) = 2054.812665 MLMG: Final Iter. 7 resid, resid/bnorm = 3.68354236e-10, 1.792641452e-13 MLMG: Timers: Solve = 3.340050542 Iter = 3.272260959 Bottom = 0.02614741 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 0.1719927642 Minimum estdt over all levels = 0.1719927642 Call to estdt at end of istep_divu_iter = 1 gives dt = 0.1719927642 Multiplying dt by init_shrink; dt = 0.1719927642 Ignoring this new dt since it's larger than the previous dt = 0.0004741795139 Writing plotfile reacting_bubble-3d_pltafter_DivuIter after final DivuIter inner sponge: r_sp , r_tp : 186468750, 224718750 Time to write plotfile: 0.487408333 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 0.0004741795139 Cell Count: Level 0, 1327104 cells inner sponge: r_sp , r_tp : 186468750, 224718750 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 2728631.974 MLMG: Initial residual (resid0) = 2728631.974 MLMG: Final Iter. 8 resid, resid/bnorm = 0.0002633368131, 9.650873243e-11 MLMG: Timers: Solve = 0.216619639 Iter = 0.211377056 Bottom = 0.161159531 <<< 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 = 2728627.802 MLMG: Initial residual (resid0) = 8728.675221 MLMG: Final Iter. 5 resid, resid/bnorm = 0.0002202815958, 8.072980699e-11 MLMG: Timers: Solve = 0.138770069 Iter = 0.136714621 Bottom = 0.105658904 <<< 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 = 1.15755262e+10 MLMG: Initial residual (resid0) = 1.15755262e+10 MLMG: Final Iter. 7 resid, resid/bnorm = 0.002913355827, 2.516823665e-13 MLMG: Timers: Solve = 3.354815979 Iter = 3.278894189 Bottom = 0.03051647 Done calling nodal solver Timestep 0 ends with TIME = 0.0004741795139 DT = 0.0004741795139 Timing summary: Advection :0.199169056 seconds MAC Proj :0.374500926 seconds Nodal Proj :3.465277183 seconds Reactions :0.168766222 seconds Misc :0.038435578 seconds Base State :0.013310348 seconds Time to advance time step: 4.246291772 Writing plotfile 0 after all initialization inner sponge: r_sp , r_tp : 186468750, 224718750 Time to write plotfile: 0.485277213 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 0.0004741795139 Cell Count: Level 0, 1327104 cells inner sponge: r_sp , r_tp : 186468750, 224718750 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1308380.202 MLMG: Initial residual (resid0) = 1308380.202 MLMG: Final Iter. 8 resid, resid/bnorm = 0.0001230611233, 9.405608792e-11 MLMG: Timers: Solve = 0.223120339 Iter = 0.220957967 Bottom = 0.171288402 <<< 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 = 1308380.239 MLMG: Initial residual (resid0) = 7.635504054 MLMG: Final Iter. 5 resid, resid/bnorm = 5.306012864e-05, 4.05540584e-11 MLMG: Timers: Solve = 0.152306565 Iter = 0.150251712 Bottom = 0.119220186 <<< 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 = 5488883.997 MLMG: Initial residual (resid0) = 5488883.997 MLMG: Final Iter. 7 resid, resid/bnorm = 6.844954507e-06, 1.247057601e-12 MLMG: Timers: Solve = 3.364283637 Iter = 3.295959993 Bottom = 0.027566587 Done calling nodal solver Timestep 1 ends with TIME = 0.0004741795139 DT = 0.0004741795139 Timing summary: Advection :0.175768819 seconds MAC Proj :0.38668396 seconds Nodal Proj :3.465399004 seconds Reactions :0.169522631 seconds Misc :0.023023236 seconds Base State :0.014736074 seconds Time to advance time step: 4.220536773 Call to estdt for level 0 gives dt_lev = 0.2012294727 Minimum estdt over all levels = 0.2012294727 Call to estdt at beginning of step 2 gives dt =0.2012294727 dt_growth factor limits the new dt = 0.0005215974653 Timestep 2 starts with TIME = 0.0004741795139 DT = 0.0005215974653 Cell Count: Level 0, 1327104 cells inner sponge: r_sp , r_tp : 186468750, 224718750 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1426181.41 MLMG: Initial residual (resid0) = 1426181.41 MLMG: Final Iter. 9 resid, resid/bnorm = 1.513370907e-05, 1.061134928e-11 MLMG: Timers: Solve = 0.246375457 Iter = 0.244288454 Bottom = 0.188725254 <<< 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 = 1426291.688 MLMG: Initial residual (resid0) = 36463.44575 MLMG: Final Iter. 9 resid, resid/bnorm = 2.261689261e-05, 1.585712993e-11 MLMG: Timers: Solve = 0.241663407 Iter = 0.239597405 Bottom = 0.183958918 <<< 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 = 6063485.841 MLMG: Initial residual (resid0) = 6063485.841 MLMG: Final Iter. 8 resid, resid/bnorm = 5.181147116e-05, 8.544832547e-12 MLMG: Timers: Solve = 3.807559179 Iter = 3.739217538 Bottom = 0.031188915 Done calling nodal solver Timestep 2 ends with TIME = 0.0009957769792 DT = 0.0005215974653 Timing summary: Advection :0.170532019 seconds MAC Proj :0.499509424 seconds Nodal Proj :3.908616444 seconds Reactions :0.168697603 seconds Misc :0.021823675 seconds Base State :0.009830984 seconds Time to advance time step: 4.769317645 Call to estdt for level 0 gives dt_lev = 0.201266633 Minimum estdt over all levels = 0.201266633 Call to estdt at beginning of step 3 gives dt =0.201266633 dt_growth factor limits the new dt = 0.0005737572118 Timestep 3 starts with TIME = 0.0009957769792 DT = 0.0005737572118 Cell Count: Level 0, 1327104 cells inner sponge: r_sp , r_tp : 186468750, 224718750 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 1597824.952 MLMG: Initial residual (resid0) = 1597824.952 MLMG: Final Iter. 9 resid, resid/bnorm = 5.46446845e-05, 3.419941867e-11 MLMG: Timers: Solve = 0.234014615 Iter = 0.231939914 Bottom = 0.176344309 <<< 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 = 1597946.233 MLMG: Initial residual (resid0) = 987353.6977 MLMG: Final Iter. 10 resid, resid/bnorm = 2.140631094e-05, 1.339613968e-11 MLMG: Timers: Solve = 0.261373015 Iter = 0.259316008 Bottom = 0.197561014 <<< 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 = 6669975.86 MLMG: Initial residual (resid0) = 6669975.86 MLMG: Final Iter. 8 resid, resid/bnorm = 2.878462223e-05, 4.315551185e-12 MLMG: Timers: Solve = 3.809889007 Iter = 3.741804148 Bottom = 0.031994884 Done calling nodal solver Timestep 3 ends with TIME = 0.001569534191 DT = 0.0005737572118 Timing summary: Advection :0.17619414 seconds MAC Proj :0.50684175 seconds Nodal Proj :3.913640531 seconds Reactions :0.167652457 seconds Misc :0.174772966 seconds Base State :0.011196896 seconds Time to advance time step: 4.939261556 Writing plotfile 3 inner sponge: r_sp , r_tp : 186468750, 224718750 Time to write plotfile: 0.486701283 Total Time: 27.805238 Unused ParmParse Variables: [TOP]::maestro.v(nvals = 1) :: [1] [TOP]::amr.ref_ratio(nvals = 4) :: [2, 2, 2, 2] [TOP]::amr.check_file(nvals = 1) :: [reacting_bubble-3d_chk] [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] Total GPU global memory (MB): 12066 Free GPU global memory (MB): 7977 [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.06-79-g9b2172017f8a) finalized