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 isgrid.set_field_parameter("normal", L)
for grid in pf.h.grids:
grid.set_field_parameter("center", c)
proj = OffAxisProjectionPlot(pf, L, "RadialVelocity", center='max', width=(width, 'pc'),
depth=depth, weight_field="Density")
Cheers,
JohnOn Mon, Dec 16, 2013 at 5:30 PM, Matthew Turk <matthewturk@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@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@gmail.com> wrote:
>>
>> Hi John,
>>
>> Sorry for the delay in replying.
>>
>> On Fri, Dec 13, 2013 at 8:59 AM, John Regan <johnanthonyregan@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@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