diff --git a/Exec/science/xrb_spherical/analysis/slice.py b/Exec/science/xrb_spherical/analysis/slice.py index 42e5b12a14..e459771a8e 100755 --- a/Exec/science/xrb_spherical/analysis/slice.py +++ b/Exec/science/xrb_spherical/analysis/slice.py @@ -5,7 +5,7 @@ import yt import numpy as np import matplotlib.pyplot as plt -from yt.units import cm +from yt.frontends.boxlib.api import CastroDataset """ Given a plot file and field name, it gives slice plots at the top, @@ -23,7 +23,7 @@ def slice(fname:str, field:str, loc:str="top") -> None: loc: location on the domain. {top, mid, bot} """ - ds = yt.load(fname, hint="castro") + ds = CastroDataset(fname) currentTime = ds.current_time.in_units("s") print(f"Current time of this plot file is {currentTime} s") @@ -39,7 +39,7 @@ def slice(fname:str, field:str, loc:str="top") -> None: theta_center = 0.5 * (thetar + thetal) # Domain width of the slice plot - width = (2.0*dr, 2.0*dr) + width = (3.0*dr, 3.0*dr) loc = loc.lower() loc_options = ["top", "mid", "bot"] @@ -48,9 +48,9 @@ def slice(fname:str, field:str, loc:str="top") -> None: raise Exception("loc parameter must be top, mid or bot") # Centers for the Top, Mid and Bot panels - centers = {"top":(r_center*np.sin(thetal)+0.8*dr, r_center*np.cos(thetal)), + centers = {"top":(r_center*np.sin(thetal)+1.48*dr, r_center*np.cos(thetal)), "mid":(r_center*np.sin(theta_center), r_center*np.cos(theta_center)), - "bot":(r_center*np.sin(thetar)+0.8*dr, r_center*np.cos(thetar))} + "bot":(r_center*np.sin(thetar)+1.48*dr, r_center*np.cos(thetar))} sp = yt.SlicePlot(ds, 'phi', field, width=width) diff --git a/Exec/science/xrb_spherical/inputs.He.1000Hz b/Exec/science/xrb_spherical/inputs.He.1000Hz index d55d64e14d..db5c4ab4a5 100644 --- a/Exec/science/xrb_spherical/inputs.He.1000Hz +++ b/Exec/science/xrb_spherical/inputs.He.1000Hz @@ -5,8 +5,8 @@ stop_time = 3.0 # PROBLEM SIZE & GEOMETRY geometry.is_periodic = 0 0 geometry.coord_sys = 2 # 0 => cart, 1 => RZ 2=>spherical -geometry.prob_lo = 1.1e6 0.392699 #0.0 -geometry.prob_hi = 1.13072e6 2.748893 #3.14159 +geometry.prob_lo = 1.1e6 0 +geometry.prob_hi = 1.13072e6 3.1415926 amr.n_cell = 192 1152 # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<