Hi all,<br><br>I would like to plot radiation pressure vectors in the style of p.modify[&quot;velocity&quot;](). Now, my radiation pressure acceleration is in my data as three baryon fields (they appear in pf.h.field_list) : RadAccel1, RadAccel2, RadAccel3 for the three spacial dimensions.<br>
<br>Now I am embarrassed to admit I have not done much yt development in the past so I was hoping someone more familiar with the code could point me in the direction of where I can add code for this? I am hoping it will be as simple as finding all the places code exists for p.modify[&quot;velocity&quot;]() and copy and paste it changing the fields used to the radiation pressure ones. Will it be? Where is the code for p.modify[&quot;velocity&quot;]()? Are there any tricky bits I should be aware of?<br>
<br>Many thanks,<br><br>Libby<br clear="all"><br>-- <br>Elizabeth Harper-Clark MA MSci<br>PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT<br>Sciences and Engineering Coordinator, Teaching Assistants&#39; Training Program, UofT<br>
<br><a href="http://www.astro.utoronto.ca/%7Eh-clark" target="_blank">www.astro.utoronto.ca/~h-clark</a><br><a href="mailto:h-clark@cita.utoronto.ca" target="_blank">h-clark@cita.utoronto.ca</a><br>Astronomy office phone: +1-416-978-5759<br>