From 5a13fc5178766ec966562228e07794003d3dae5a Mon Sep 17 00:00:00 2001 From: zhi Date: Tue, 29 Oct 2024 18:33:20 -0400 Subject: [PATCH] add comment about making sliceplots --- Exec/science/xrb_spherical/analysis/slice.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/Exec/science/xrb_spherical/analysis/slice.py b/Exec/science/xrb_spherical/analysis/slice.py index 5ac7905f39..7e55d3848f 100755 --- a/Exec/science/xrb_spherical/analysis/slice.py +++ b/Exec/science/xrb_spherical/analysis/slice.py @@ -7,6 +7,8 @@ import matplotlib.pyplot as plt from yt.frontends.boxlib.api import CastroDataset +from yt.units import cm + """ Given a plot file and field name, it gives slice plots at the top, middle, and bottom of the domain (shell). @@ -54,8 +56,10 @@ def slice(fname:str, field:str, "mid":(r_center*np.sin(theta_center), r_center*np.cos(theta_center)), "bot":(r_center*np.sin(thetar)+0.5*width, r_center*np.cos(thetar))} - sp = yt.SlicePlot(ds, 'phi', field, width=box_widths) + # Note we can also set center during SlicePlot, however then we would enter in [r_center, theta_center, 0] + # rather than the physical R-Z coordinate if we do it via sp.set_center + sp = yt.SlicePlot(ds, 'phi', field, width=box_widths) sp.set_center(centers[loc]) sp.set_cmap(field, "viridis")