<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Hi Tazkera,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">You can always get rays and do the integration yourself:</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">ray = ds.ray(p1,p2)</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">dens = ray['density']</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">where p1 and p2 are tuples of the end points of the ray. If you're using cartesian coordinates then you can also get x, y and z along the ray:</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">x = ray['x']</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">etc.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">then the position is just </div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">r = np.sqrt(x**2 + y**2 + z**2)</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">and you can use a numpy or scipy integration method.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">If you're using data from an AMR run you may need to sort the data:</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">tord = ray['t'])</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">dens = ray['dens'][tord]</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">and similarly for other data from the ray.</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small"><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Regards,</div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif;font-size:small">Jon</div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 1, 2017 at 1:17 AM, <span dir="ltr"><<a href="mailto:yt-users-request@lists.spacepope.org" target="_blank">yt-users-request@lists.spacepope.org</a>></span> wrote:<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Date: Tue, 31 Oct 2017 17:56:39 -0400<br>
From: tazkera haque <<a href="mailto:h.tazkera@gmail.com">h.tazkera@gmail.com</a>><br>
To: Discussion of the yt analysis package<br>
<<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a>><br>
Subject: [yt-users] Integration of a field along line of sight<br>
Message-ID:<br>
<CAFj_7zy_P2FrQSf9DMi0-<wbr>DkXZOOSvdW=<a href="mailto:poWO_JzZ_cgQD%2BaCYQ@mail.gmail.com">poWO_JzZ_cgQD+aCYQ@<wbr>mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear YT people,<br>
<br>
I need to do integration of density field along line of sight. I know the<br>
ProjectionPlot function plots the integration of a field along line of<br>
sight.<br>
<br>
yt.ProjectionPlot(ds, 'z', 'density, center= 'm')<br>
<br>
However, I was wondering if it is possible to get the projection as an<br>
array (or new field) , so that I can make a ProfilePlot of the integrated<br>
field along line of sight and another field ('Mass')<br>
<br>
In other words, I want to do the following:<br>
<br>
sp = ds.sphere("max", (200, "kpc"))<br>
proj = yt.ProfilePlot(sp, 'Integrated_Density_along_<wbr>Line_Of_Sight', 'Mass)<br>
<br>
Thanks<br>
Tazkera<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20171031/038832cb/attachment-0001.html" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>pipermail/yt-users-spacepope.<wbr>org/attachments/20171031/<wbr>038832cb/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://lists.spacepope.org/<wbr>listinfo.cgi/yt-users-<wbr>spacepope.org</a><br>
<br>
<br>
------------------------------<br>
<br>
End of yt-users Digest, Vol 116, Issue 34<br>
******************************<wbr>***********<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">________________________________________________________<br>Jonathan D. Slavin Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a> 60 Garden Street, MS 83<br>phone: (617) 496-7981 Cambridge, MA 02138-1516<br>cell: (781) 363-0035 USA<br>________________________________________________________<br><br></div></div></div></div>
</div></div>