Initializing CUDA... CUDA initialized with 1 GPU AMReX (19.12-36-g87659c492175) initialized Calling Setup() Calling ReadParameters() WARNING: burning_cutoff_density not supplied in the inputs file WARNING: setting burning_cutoff_density = base_cutoff_density reading extern runtime parameters ... ERROR: problem in the namelist Calling VariableSetup() Calling set_method_params() Initializing Helmholtz EOS and using Coulomb corrections. Calling BCSetup() Calling init_base_state_geometry() Calling Init() Calling InitData() ------------------------------------------------------------------------------ reading initial model 320 points found in the initial model file 6 variables found in the initial model file ------------------------------------------------------------------------------ model file mapping (spherical base state) dr of MAESTRO base state = 781250.0000 dr of input file data = 1562500.000 maximum radius (cell-centered) of input model = 499218750.0 setting r_cutoff to 239 radius at r_cutoff 187109375.0000000 Maximum HSE Error = 0.2310976572E-02 (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 : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.801113651 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Doing initial projection Calling nodal solver MLMG: Initial rhs = 188240648.4 MLMG: Initial residual (resid0) = 188240648.4 MLMG: Final Iter. 12 resid, resid/bnorm = 0.01633389294, 8.677133807e-11 MLMG: Timers: Solve = 2.545798195 Iter = 2.509163191 Bottom = 0.000845524 Done calling nodal solver Writing plotfile wdconvect-amr_pltafter_InitProj after InitProj inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.697253133 Call to firstdt for level 0 gives dt_lev = 0.005789678157 Multiplying dt_lev by init_shrink; dt_lev = 0.0005789678157 Call to firstdt for level 1 gives dt_lev = 0.0005789678157 Multiplying dt_lev by init_shrink; dt_lev = 5.789678157e-05 Minimum firstdt over all levels = 5.789678157e-05 Doing initial divu iteration #1 Calling nodal solver MLMG: Initial rhs = 80008358.59 MLMG: Initial residual (resid0) = 80008358.59 MLMG: Final Iter. 2 resid, resid/bnorm = 0.482713986, 6.033294452e-09 MLMG: Timers: Solve = 0.466661618 Iter = 0.431518906 Bottom = 0.000147196 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 5.916319527 Call to estdt for level 1 gives dt_lev = 2.341804478 Minimum estdt over all levels = 2.341804478 Call to estdt at end of istep_divu_iter = 1 gives dt = 2.341804478 Multiplying dt by init_shrink; dt = 0.2341804478 Ignoring this new dt since it's larger than the previous dt = 5.789678157e-05 Doing initial divu iteration #2 Calling nodal solver MLMG: Initial rhs = 44170163.99 MLMG: Initial residual (resid0) = 44170163.99 MLMG: Final Iter. 2 resid, resid/bnorm = 0.02610203723, 5.909427285e-10 MLMG: Timers: Solve = 0.471117426 Iter = 0.431860184 Bottom = 0.000156556 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 5.916319527 Call to estdt for level 1 gives dt_lev = 2.341804478 Minimum estdt over all levels = 2.341804478 Call to estdt at end of istep_divu_iter = 2 gives dt = 2.341804478 Multiplying dt by init_shrink; dt = 0.2341804478 Ignoring this new dt since it's larger than the previous dt = 5.789678157e-05 Doing initial divu iteration #3 Calling nodal solver MLMG: Initial rhs = 26504865.08 MLMG: Initial residual (resid0) = 26504865.08 MLMG: Final Iter. 3 resid, resid/bnorm = 0.001814276911, 6.845071294e-11 MLMG: Timers: Solve = 0.686375077 Iter = 0.647142987 Bottom = 0.000226997 Done calling nodal solver Call to estdt for level 0 gives dt_lev = 5.916319527 Call to estdt for level 1 gives dt_lev = 2.341804478 Minimum estdt over all levels = 2.341804478 Call to estdt at end of istep_divu_iter = 3 gives dt = 2.341804478 Multiplying dt by init_shrink; dt = 0.2341804478 Ignoring this new dt since it's larger than the previous dt = 5.789678157e-05 Writing plotfile wdconvect-amr_pltafter_DivuIter after final DivuIter inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.721706213 Doing initial pressure iteration #1 Timestep 0 starts with TIME = 0 DT = 5.789678157e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 173628411.4 MLMG: Initial residual (resid0) = 173628411.4 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02237717761, 1.288797002e-10 MLMG: Timers: Solve = 0.510276531 Iter = 0.473787605 Bottom = 0.001063723 <<< 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 = 173628411.8 MLMG: Initial residual (resid0) = 2158.123376 MLMG: Final Iter. 2 resid, resid/bnorm = 0.04143755587, 2.386565392e-10 MLMG: Timers: Solve = 0.196153814 Iter = 0.162531845 Bottom = 0.000412809 <<< 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 = 35718515.71 MLMG: Initial residual (resid0) = 35718515.71 MLMG: Final Iter. 35 resid, resid/bnorm = 0.002761640004, 7.731676272e-11 MLMG: Timers: Solve = 7.311943514 Iter = 7.272540264 Bottom = 0.002486704 Done calling nodal solver Timestep 0 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05 Timing summary: Advection :4.692674655 seconds MAC Proj :0.8872916 seconds Nodal Proj :7.544776763 seconds Reactions :0.219515242 seconds Misc :0.114240731 seconds Base State :0.045155249 seconds Time to advance time step: 13.45905378 Writing plotfile 0 after all initialization inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.767565252 Calling Evolve() Timestep 1 starts with TIME = 0 DT = 5.789678157e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 173628606.2 MLMG: Initial residual (resid0) = 173628606.2 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02241813624, 1.291154535e-10 MLMG: Timers: Solve = 0.509998779 Iter = 0.476640177 Bottom = 0.001054222 <<< 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 = 173628606.2 MLMG: Initial residual (resid0) = 0.353521162 MLMG: Final Iter. 1 resid, resid/bnorm = 0.02251002256, 1.296446654e-10 MLMG: Timers: Solve = 0.118225384 Iter = 0.084752723 Bottom = 0.000157371 <<< 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 = 16938159.33 MLMG: Initial residual (resid0) = 16938159.33 MLMG: Final Iter. 10 resid, resid/bnorm = 0.001225439832, 7.234787484e-11 MLMG: Timers: Solve = 2.222047373 Iter = 2.1785341 Bottom = 0.000739102 Done calling nodal solver Timestep 1 ends with TIME = 5.789678157e-05 DT = 5.789678157e-05 Timing summary: Advection :4.627438865 seconds MAC Proj :0.803546637 seconds Nodal Proj :2.449385023 seconds Reactions :0.223487714 seconds Misc :0.133015396 seconds Base State :0.046587063 seconds Time to advance time step: 8.23752734 Writing plotfile 1 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.734728959 Call to estdt for level 0 gives dt_lev = 1.112476815 Call to estdt for level 1 gives dt_lev = 0.2935011253 Minimum estdt over all levels = 0.2935011253 Timestep 2 starts with TIME = 5.789678157e-05 DT = 6.368645973e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 166900886.1 MLMG: Initial residual (resid0) = 166900886.1 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02282487787, 1.367570803e-10 MLMG: Timers: Solve = 0.512430035 Iter = 0.479116923 Bottom = 0.004658936 <<< 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 = 166900886.1 MLMG: Initial residual (resid0) = 0.3012107426 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01909282326, 1.143961767e-10 MLMG: Timers: Solve = 0.118562542 Iter = 0.084965248 Bottom = 0.000209552 <<< 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 = 11608423.25 MLMG: Initial residual (resid0) = 11608423.25 MLMG: Final Iter. 9 resid, resid/bnorm = 0.000867142342, 7.469940776e-11 MLMG: Timers: Solve = 1.998150529 Iter = 1.958601555 Bottom = 0.000649295 Done calling nodal solver Timestep 2 ends with TIME = 0.0001215832413 DT = 6.368645973e-05 Timing summary: Advection :4.615115506 seconds MAC Proj :0.807003306 seconds Nodal Proj :2.226599993 seconds Reactions :0.209171319 seconds Misc :0.129486251 seconds Base State :0.038638979 seconds Time to advance time step: 7.987993464 Writing plotfile 2 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.743602134 Time to regrid: 0.117030122 Call to estdt for level 0 gives dt_lev = 1.539392516 Call to estdt for level 1 gives dt_lev = 0.3670776926 Minimum estdt over all levels = 0.3670776926 Timestep 3 starts with TIME = 0.0001215832413 DT = 7.00551057e-05 Cell Count: Level 0, 262144 cells Level 1, 304128 cells inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 <<< STEP 1 : react state >>> <<< STEP 2 : make w0 >>> <<< STEP 3 : create MAC velocities >>> MLMG: Initial rhs = 166879353.3 MLMG: Initial residual (resid0) = 166879353.3 MLMG: Final Iter. 6 resid, resid/bnorm = 0.02258217894, 1.353203886e-10 MLMG: Timers: Solve = 0.507803136 Iter = 0.474315701 Bottom = 0.001028776 <<< 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 = 166879353.2 MLMG: Initial residual (resid0) = 0.3057069979 MLMG: Final Iter. 1 resid, resid/bnorm = 0.01939887887, 1.162449308e-10 MLMG: Timers: Solve = 0.118232656 Iter = 0.084598494 Bottom = 0.000205069 <<< 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 = 9504180.595 MLMG: Initial residual (resid0) = 9504180.595 MLMG: Final Iter. 9 resid, resid/bnorm = 0.0008270535618, 8.701997543e-11 MLMG: Timers: Solve = 1.972945902 Iter = 1.929136203 Bottom = 0.000693093 Done calling nodal solver Timestep 3 ends with TIME = 0.000191638347 DT = 7.00551057e-05 Timing summary: Advection :4.63735529 seconds MAC Proj :0.804142961 seconds Nodal Proj :2.22552486 seconds Reactions :0.211379785 seconds Misc :0.130101948 seconds Base State :0.040079376 seconds Time to advance time step: 8.009192536 Writing plotfile 3 inner sponge: r_sp , r_tp : 0.151953125000E+09 0.183203125000E+09 outer sponge: r_sp_outer, r_tp_outer: 0.183203125000E+09 0.198828125000E+09 Time to write plotfile: 0.747149963 Total Time: 75.42296518 Unused ParmParse Variables: [TOP]::maestro.v(nvals = 1) :: [1] [TOP]::system.regtest_reduction(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): 9302 [The Arena] space (MB): 9049 [The Device Arena] space (MB): 8 [The Managed Arena] space (MB): 8 [The Pinned Arena] space (MB): 8 AMReX (19.12-36-g87659c492175) finalized