Hey all,
I was hoping that someone might have a suggestion or answer for an
issue I'm having with generating simple radial profiles from
simulation data. The profiles automatically generate axis bounds and
while I can set the bounds on the x axis using the x_bounds parameter,
I don't know how, or if it's even possible, to set a similar
restriction for the y axis. (I have tried simply setting a y_bounds
parameter, but the script breaks stating that y_bounds is unrecognized).
Here is an example of the code:
from yt.mods import * # set up our namespace
for n in range(1,242):
pf = load("/work/01198/earln/RSW_nC_1-0/DD%04i/DD%04i" % (n,n))
pc = PlotCollection(pf, center=[0.5, 0.5, 0.5])
pc.add_profile_sphere(10.0, "mpc",
["Radius", "RadialVelocity"], weight="CellMassMsun",
x_bounds=(1.0e+18,1.0e+20), y_bounds=(1e-3, 1.e+8))
pc.save("/work/01198/earln/RSW_nC_1-0/graph-data_radvel/
profiles/%s" % pf)
Any suggestions or help would be greatly appreciated,
Nicholas Earl