modify - radiation pressure

Hi all, I would like to plot radiation pressure vectors in the style of p.modify["velocity"](). 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. 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["velocity"]() and copy and paste it changing the fields used to the radiation pressure ones. Will it be? Where is the code for p.modify["velocity"]()? Are there any tricky bits I should be aware of? Many thanks, Libby -- Elizabeth Harper-Clark MA MSci PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT www.astro.utoronto.ca/~h-clark <http://www.astro.utoronto.ca/%7Eh-clark> h-clark@cita.utoronto.ca Astronomy office phone: +1-416-978-5759

Hi Libby, If you use QuiverCallback, which would be p.modify["quiver"], you can specify the fields for the x and y axes in the image plane. The call is p.modify["quiver"](field1, field2, factor), where factor controls the scaling of the arrows. If you want to see specifically how p.modify["velocity"] works, check out VelocityCallback in yt/visualization/plot_modifications.py. Britton On Mon, Feb 21, 2011 at 3:12 PM, Elizabeth Harper-Clark < h-clark@astro.utoronto.ca> wrote:
Hi all,
I would like to plot radiation pressure vectors in the style of p.modify["velocity"](). 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.
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["velocity"]() and copy and paste it changing the fields used to the radiation pressure ones. Will it be? Where is the code for p.modify["velocity"]()? Are there any tricky bits I should be aware of?
Many thanks,
Libby
-- Elizabeth Harper-Clark MA MSci PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT
www.astro.utoronto.ca/~h-clark <http://www.astro.utoronto.ca/%7Eh-clark> h-clark@cita.utoronto.ca Astronomy office phone: +1-416-978-5759
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (2)
-
Britton Smith
-
Elizabeth Harper-Clark