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.com> 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 <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>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