[yt-users] Radial velocity Projections
John Regan
johnanthonyregan at gmail.com
Tue Dec 17 04:25:05 PST 2013
Hi Matt,
Thanks for having a look. It seemed to sample the field OK but then fell
over.
Traceback (most recent call last):
File "OffAxisProjection_RadialVelocity_TimeSeries.py", line 46, in
<module>
depth=depth, weight_field="Density")
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/src/yt-hg/yt/visualization/plot_window.py",
line 1149, in __init__
self.set_axes_unit(axes_unit)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/src/yt-hg/yt/visualization/plot_container.py",
line 56, in newfunc
args[0]._setup_plots()
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/src/yt-hg/yt/visualization/plot_window.py",
line 665, in _setup_plots
cax)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/src/yt-hg/yt/visualization/plot_window.py",
line 1416, in __init__
self._init_image(data, cbname, cmap, extent, aspect)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/src/yt-hg/yt/visualization/base_plot_types.py",
line 109, in _init_image
self.cb = self.figure.colorbar(self.image, self.cax)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/figure.py",
line 1451, in colorbar
cb = cbar.colorbar_factory(cax, mappable, **kw)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/colorbar.py",
line 1274, in colorbar_factory
cb = Colorbar(cax, mappable, **kwargs)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/colorbar.py",
line 877, in __init__
ColorbarBase.__init__(self, ax, **kw)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/colorbar.py",
line 317, in __init__
self.draw_all()
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/colorbar.py",
line 338, in draw_all
self._process_values()
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/colorbar.py",
line 652, in _process_values
b = self.norm.inverse(self._uniform_y(self.cmap.N + 1))
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/matplotlib-1.3.0-py2.7-linux-x86_64.egg/matplotlib/colors.py",
line 1003, in inverse
return vmin * ma.power((vmax / vmin), val)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/numpy/ma/core.py",
line 3648, in __mul__
return multiply(self, other)
File
"/homeappl/home/regan/appl_taito/YT/yt-x86_64/lib/python2.7/site-packages/numpy/ma/core.py",
line 941, in __call__
if m:
ValueError: The truth value of an array with more than one element is
ambiguous. Use a.any() or a.all()
The snippet of code that I run is
for grid in pf.h.grids:
grid.set_field_parameter("center", c)
grid.set_field_parameter("normal", L)
proj = OffAxisProjectionPlot(pf, L, "RadialVelocity", center='max',
width=(width, 'pc'),
depth=depth, weight_field="Density")
Cheers,
John
On Mon, Dec 16, 2013 at 5:30 PM, Matthew Turk <matthewturk at gmail.com> wrote:
> Hi John,
>
> Yup, you can. I'm looking at the code right now and I think the issue
> is going to be tricky to maneuver through the code base in a reliable
> way the field parameters. Here's what I'd suggest doing --
>
> for grid in pf.h.grids:
> grid.set_field_parameter("center", c)
>
> just before you call OffAxisProjection. This is something that the
> grid machinery is a bit annoying about the way that things work in
> 2.x, but we're fixing in 3.0. You'll also have to set the "normal"
> field parameter.
>
> Let me know if that works,
>
> -Matt
>
>
> On Mon, Dec 16, 2013 at 9:55 AM, John Regan <johnanthonyregan at gmail.com>
> wrote:
> > Hi Matt,
> >
> > Yes it is the OAP I am interested in. I had originally weighted the
> > projection using the density field.
> >
> > The projected velocity may look funny and it would be interesting to try
> to
> > compare that to the planar. Could I also pass cyl_RadialVelocity to the
> OAP
> > API?
> >
> > Cheers,
> > John
> >
> > On 16 Dec 2013 16:09, "Matthew Turk" <matthewturk at gmail.com> wrote:
> >>
> >> Hi John,
> >>
> >> Sorry for the delay in replying.
> >>
> >> On Fri, Dec 13, 2013 at 8:59 AM, John Regan <johnanthonyregan at gmail.com
> >
> >> wrote:
> >> > Hi All,
> >> >
> >> > I was trying to make an OffAxis radial projection. Naively using
> >> >
> >> > proj = OffAxisProjectionPlot(pf, L, "RadialVelocityKMS",
> center='max',
> >> > width=(width, 'pc'),
> >> > depth=depth)
> >> >
> >> > fails with a message of type
> >> >
> >> > yt.data_objects.field_info_container.NeedsParameter: (['center'])
> >> >
> >> >
> >> > The RadialVelocity probably isn't the field to use here and requires
> >> > some
> >> > manipulation to be a quantity suitable for plotting.
> >> >
> >> > Any suggestions on the best way forward would be much appreciated.
> >> >
> >>
> >> So, a few questions. Do you want to show the way that the radial
> >> velocity is distributed, with respect to an image plane that isn't
> >> axis aligned? If you so, you might have more luck using the
> >> OffAxisSlicePlot, which will take a slice through the domain. The OAP
> >> plot will integrate along the line of sight; you will probably want to
> >> take an average if this is what you want, by supplying the weight
> >> parameter here, as otherwise it'll just sum up the RV values at all
> >> points. I think the issue with it not having the parameter center is
> >> fixable, but we should probably make sure that's what you're looking
> >> for first. It's also going to be doing the full 3D radial velocity,
> >> not the planar velocity, so if your velocity is roughly spherical it
> >> might have odd morphology in the resultant plot.
> >>
> >> If a slice will do what you're looking for, I'd recommend the
> >> OffAxisSlicePlot -- if you indeed do want a radial average, then I'll
> >> take a look at passing the parameter in. (This will also help by
> >> enabling you to pass in the bulk velocity.) And if you want to find
> >> the planar RV, you'll have to use a field like cyl_RadialVelocity.
> >>
> >> -Matt
> >>
> >> > Cheers,
> >> > John
> >> >
> >> > _______________________________________________
> >> > yt-users mailing list
> >> > yt-users at lists.spacepope.org
> >> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >> >
> >> _______________________________________________
> >> yt-users mailing list
> >> yt-users at lists.spacepope.org
> >> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
> >
> > _______________________________________________
> > yt-users mailing list
> > yt-users at lists.spacepope.org
> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> >
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131217/149b57af/attachment.html>
More information about the yt-users
mailing list