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 wrote:

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 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 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>I

...

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

Hi Stefano, I think I can help explain why that line you last emailed is the correct one. The c[1] and center='max' do two different things. The c[1] tells it to slice the y dimension at y=c[1]. The center='max' tells it where to center the image in the x-z plane. If you don't have the center='max', then you'll end up at the center of the domain I think, where you might have poor resolution. Since the projection doesn't need a specified c[1] like value, it correctly centers simply using the center='max' or center=c command. Hopefully that clears it up a bit. Best, Sam On Mon, Apr 29, 2013 at 8:17 AM, stefano bovino wrote:

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 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 wrote:

...

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 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 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>I

...

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

_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org