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
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
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
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
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 <poetaste@gmail.com> 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 <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
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (3)
-
Matthew Turk -
Sam Skillman -
stefano bovino