Re: [yt-users] Radial velocity Projections

17 Dec 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  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
...
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"  wrote:
...
Hi John,
Sorry for the delay in replying.
On Fri, Dec 13, 2013 at 8:59 AM, John Regan 
...
wrote:
...
Hi All,
I was trying to make an OffAxis radial projection. Naively using
proj = OffAxisProjectionPlot(pf, L, "RadialVelocityKMS",
On Mon, Dec 16, 2013 at 9:55 AM, John Regan 
wrote:
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@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
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org