@@ -6,16 +6,16 @@ stop_time = 3.0
66geometry.is_periodic = 0 0
77geometry.coord_sys = 2 # 0 => cart, 1 => RZ 2=>spherical
88geometry.prob_lo = 1.1e6 0
9- geometry.prob_hi = 1.13072e6 3.141592653589793
10- amr.n_cell = 768 2304 #192 1152
9+ geometry.prob_hi = 1.13072e6 0.785398 # pi/4
10+ amr.n_cell = 1536 43200 #This achieves 20cm res in both radial and theta dir
1111
1212# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
1313# 0 = Interior 3 = Symmetry
1414# 1 = Inflow 4 = SlipWall
1515# 2 = Outflow 5 = NoSlipWall
1616# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
1717castro.lo_bc = 3 3 # Inflow in lower R and Symmetry across Theta
18- castro.hi_bc = 2 3 # Outflow in upper boundaries
18+ castro.hi_bc = 2 2 # Outflow in upper boundaries
1919
2020# Allow non-square zones
2121castro.allow_non_unit_aspect_zones = 1
@@ -56,8 +56,8 @@ castro.rotation_include_centrifugal = 0
5656castro.diffuse_temp = 1
5757castro.diffuse_cutoff_density_hi = 5.e4
5858castro.diffuse_cutoff_density = 2.e4
59-
6059castro.diffuse_cond_scale_fac = 1.0
60+ castro.diffuse_use_amrex_mlmg = 0
6161
6262castro.react_rho_min = 1.e2
6363castro.react_rho_max = 1.5e7
@@ -96,15 +96,15 @@ amr.n_error_buf = 2 2 2 2 # number of buffer cells in error est
9696
9797# CHECKPOINT FILES
9898amr.check_file = xrb_spherical_1000Hz_chk # root name of checkpoint file
99- amr.check_int = 250 # number of timesteps between checkpoints
99+ amr.check_int = 50000 # number of timesteps between checkpoints
100100
101101# PLOTFILES
102102amr.plot_file = xrb_spherical_1000Hz_plt # root name of plotfile
103103amr.plot_per = 5.e-3 # number of seconds between plotfiles
104104amr.derive_plot_vars = ALL
105105
106106amr.small_plot_file = xrb_spherical_1000Hz_smallplt # root name of plotfile
107- amr.small_plot_per = 1.e-6 #1.e- 4 # number of seconds between plotfiles
107+ amr.small_plot_per = 1.e-4 # number of seconds between plotfiles
108108amr.small_plot_vars = density Temp
109109amr.derive_small_plot_vars = abar x_velocity y_velocity z_velocity enuc
110110
@@ -121,7 +121,7 @@ problem.T_hi = 2.e8
121121problem.T_lo = 8.e6
122122
123123problem.H_star = 2000.e0
124- problem.atm_delta = 400.0 # 50.0
124+ problem.atm_delta = 50.0
125125
126126problem.fuel1_name = "helium-4"
127127problem.fuel1_frac = 1.0e0
0 commit comments