<div dir="ltr"><div><div><div><div>OH yeah, that is the other funky thing about what I am doing...(I am rolling my eyes at myself for not mentioning this right away...)<br><br></div>I am using another person's enzo run, and he cannot find the parameter file. I made one as best as I could from copying a redshift output into what I thought the format should be, but there could be some problems there. <br><br></div>See attached. Just in case, I am also attaching the redshiftoutput file that I used.<br><br></div>Thank you for looking into this!<br><br></div>Stephanie<br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
<br><div class="gmail_quote">On Thu, May 26, 2016 at 4:20 AM, Britton Smith <span dir="ltr"><<a href="mailto:brittonsmith@gmail.com" target="_blank">brittonsmith@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Stephanie,<div><br></div><div>Could you send along the Enzo parameter file that you're using with this light ray?</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Britton</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, May 26, 2016 at 4:41 AM, Cameron Hummels <span dir="ltr"><<a href="mailto:chummels@gmail.com" target="_blank">chummels@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Stephanie,<div><br></div><div>It looks like there are two problems you're encountering. One, you're trying to include a field "velocity", when there is no such field in your dataset, only 3 scalar fields for "velocity_x", "velocity_y"... These will be included by default when you set use_peculiar_velocity=True as you are doing. So for one, remove the velocity field from your list of fields. The second problem, as you quite rightly identify, is that yt is messing up the order that it's getting the redshifts of the outputs, so it's jumping straight from 0.4 to 0.1. I can't dig into the code right now to see why it is ordering it weirdly, but I might be able to look at it in the next day or two. Britton may also have an idea on what is going on here.</div><div><br></div><div>Also of note is that there are some nice wrapper functions on the LightRay functionality in the new Trident code. In particular, the make_compound_ray function wraps this cleanly. That said, it will not fix this problem, as it is just using the underlying yt LightRay infrastructure: <a href="http://trident.readthedocs.io/en/latest/generated/trident.make_compound_ray.html#trident.make_compound_ray" target="_blank">http://trident.readthedocs.io/en/latest/generated/trident.make_compound_ray.html#trident.make_compound_ray</a></div><div><br></div><div>Cameron</div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Wed, May 25, 2016 at 1:47 PM, Stephanie Tonnesen <span dir="ltr"><<a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div><div><div>Hi all,<br><br></div>I am trying to make a light ray using a few redshift outputs from z=0.4 to z=0.1. Right now I just have two commands in my code taken basically directly from the docs:<br><br>lr = LightRay("CBOX.enzo",simulation_type="Enzo",near_redshift=0.1,<br> far_redshift=0.4,time_data=False,redshift_data=True,find_outputs=True)<br><br>lr.make_light_ray(seed=8934876,fields=['temperature','density','velocity'],use_peculiar_velocity=True)<br><br></div>My current problem is that the first command "LightRay" doesn't seem to be using any of my outputs between 0.4 and 0.1! I am attaching the error file, and you can see that it finds 9 outputs, 7 of which are between and including 0.4 to 0.1. However, then it claims that I am asking it to go straight from 0.399996 to 0.1. <br><br></div><div>I do only have redshift dumps to work with, but I think that those should be used since redshift_data=True<br><br></div><div>Does anyone have any advice for making this work?<br><br></div><div>Thanks,<br></div><div>Stephanie<br></div><div><div><div><div><div><div><div><div><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
</div></div></div></div></div></div></div>
<br></div></div>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><span><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><div>Cameron Hummels<div><span style="font-size:12.8000001907349px">NSF Postdoctoral Fellow</span></div><div><span style="font-size:12.8000001907349px">Department of Astronomy</span></div><div><span style="font-size:12.8000001907349px">California Institute of Technology</span><br></div><div><a href="http://chummels.org" style="font-size:12.8000001907349px" target="_blank">http://chummels.org</a><br></div></div></div></div>
</font></span></div>
<br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org" target="_blank">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/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>
</div></div><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" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>