<p dir="ltr">Hi Matt,</p>
<p dir="ltr">Yes it is the OAP I am interested in. I had originally weighted the projection using the density field. </p>
<p dir="ltr">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?</p>
<p dir="ltr">Cheers,<br>
John</p>
<div class="gmail_quote">On 16 Dec 2013 16:09, "Matthew Turk" <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi John,<br>
<br>
Sorry for the delay in replying.<br>
<br>
On Fri, Dec 13, 2013 at 8:59 AM, John Regan <<a href="mailto:johnanthonyregan@gmail.com">johnanthonyregan@gmail.com</a>> wrote:<br>
> Hi All,<br>
><br>
> I was trying to make an OffAxis radial projection. Naively using<br>
><br>
>  proj = OffAxisProjectionPlot(pf, L, "RadialVelocityKMS", center='max',<br>
> width=(width, 'pc'),<br>
>                                      depth=depth)<br>
><br>
> fails with a message of type<br>
><br>
> yt.data_objects.field_info_container.NeedsParameter: (['center'])<br>
><br>
><br>
> The RadialVelocity probably isn't the field to use here and requires some<br>
> manipulation to be a quantity suitable for plotting.<br>
><br>
> Any suggestions on the best way forward would be much appreciated.<br>
><br>
<br>
So, a few questions.  Do you want to show the way that the radial<br>
velocity is distributed, with respect to an image plane that isn't<br>
axis aligned?  If you so, you might have more luck using the<br>
OffAxisSlicePlot, which will take a slice through the domain.  The OAP<br>
plot will integrate along the line of sight; you will probably want to<br>
take an average if this is what you want, by supplying the weight<br>
parameter here, as otherwise it'll just sum up the RV values at all<br>
points.  I think the issue with it not having the parameter center is<br>
fixable, but we should probably make sure that's what you're looking<br>
for first.  It's also going to be doing the full 3D radial velocity,<br>
not the planar velocity, so if your velocity is roughly spherical it<br>
might have odd morphology in the resultant plot.<br>
<br>
If a slice will do what you're looking for, I'd recommend the<br>
OffAxisSlicePlot -- if you indeed do want a radial average, then I'll<br>
take a look at passing the parameter in.  (This will also help by<br>
enabling you to pass in the bulk velocity.)  And if you want to find<br>
the planar RV, you'll have to use a field like cyl_RadialVelocity.<br>
<br>
-Matt<br>
<br>
> Cheers,<br>
> John<br>
><br>
> _______________________________________________<br>
> yt-users mailing list<br>
> <a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
> <a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
><br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div>