Hi Matt, I tried to make some changes following one of your previous suggestions (mixing):
v, c = pf.h.find_max("Density") slc = pf.h.slice(1, c[1], fields=["Density","Temperature","VelocityMagnitude"],center="max")
this is the way how it works... do you have any explanation?
Thanks Stefano
On Mon, Apr 29, 2013 at 3:33 PM, stefano bovino poetaste@gmail.com wrote:
Hi Matt, thanks for the suggestion... the projection works well only if I use c instead of max but the slice still not working, here following the code:
v, c = pf.h.find_max("Density") slc = pf.h.slice(2, c[2], fields=["Density","Temperature","VelocityMagnitude"]) proj = pf.h.proj(2, "Density", weight_field="Density", center=c)
Stefano
On Mon, Apr 29, 2013 at 3:10 PM, Matthew Turk matthewturk@gmail.comwrote:
Hi Stefano,
Changing the argument to ' center = "max" ' for proj should fix the projection, but for slice you will need to change:
v, c = pf.h.find_max("Density")
slc = pf.h.slice(2, c[2], ...
So that it also slices at the point of maximum density. Note also that if you want to slice along x, you would change the 2's to 0's, for y you would change to 1.
-Matt
On Mon, Apr 29, 2013 at 8:59 AM, stefano bovino poetaste@gmail.com wrote:
Hi Kacper, my minimal resolution is less than 1 AU, the problem seems more related
to
the Matt's issue.
The script is the same as in the yt website, here the link:
http://yt-project.org/doc/cookbook/multi_plot_slice_and_proj.py
the output of pf.h.find_max("Density") :
(1.139355201341093e-11, array([ 0.49722083, 0.50645245, 0.49492896]))
Thanks Stefano
---------- Forwarded message ---------- From: Kacper Kowalik xarthisius.kk@gmail.com Date: Mon, Apr 29, 2013 at 1:35 PM Subject: Re: [yt-users] slice plot To: yt-users@lists.spacepope.org
On 04/29/2013 01:22 PM, stefano bovino wrote:
Hi all, when using the following script fromt the yt.project website
multi_plot_slice_and_proj.py<
http://yt-project.org/doc/cookbook/multi_plot_slice_and_proj.py%3EI
have problems if I change the resolution from mpc to pc.
slc_frb = slc.to_frb((1.0, "mpc"), 512) --> slc_frb = slc.to_frb((1.0, "pc"), 512)
I even tried other ways (or scripts) to plot a slice of the density changing the resolution, but the problem is still there. Do you have
any
solution?
I just would like to create a 9 panels plot with slice of density, H2,
and
Temperature with three different width (on the pc scale).
Hi Stefano, what is the size of your computational domain or size of the smallest cell you have? You can obtain that information using: yt stats output_filename
Cheers, Kacper
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org