<div dir="ltr">Hi Britton,<div><br></div><div>Thank you so much for your detailed reply! I've read documentation about Light Cone Generator, Light Ray Generator and Absorption Spectrum. I also downloaded the dataset and ran the example scripts with it. There might be two typos here:</div>
<div>(1) In Section "Adding Features to the Spectrum" on <a href="http://yt-project.org/doc/analysis_modules/absorption_spectrum.html">http://yt-project.org/doc/analysis_modules/absorption_spectrum.html</a> gamme should be gamma.</div>
<div>(2) When I ran make_light_ray.py I got an error</div><div>"OSError: [Errno 2] No such file or directory: 'halo_analysis/enzo_tiny_cosmology/RD0000'"</div><div>I need to do "cp -r enzo_tiny_cosmology halo_analysis/" to make it work.</div>
<div><br></div><div>Back to my enzo parameter file, which is attached. I correct that line and know that I don't have enough dataset, so I ran plan_cosmology_splice function to it but I got the following error:</div><div>
<br></div><div><div>Traceback (most recent call last):</div><div> File "plan_cosmology_splice.py", line 2, in <module></div><div> my_splice = CosmologySplice('enzo_run/32cube.enzo', 'Enzo')</div>
<div> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/analysis_modules/cosmological_observation/cosmology_splice.py", line 44, in __init__</div><div> find_outputs=find_outputs)</div><div> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/convenience.py", line 131, in simulation</div>
<div> find_outputs=find_outputs)</div><div> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/frontends/enzo/simulation_handling.py", line 81, in __init__</div><div> find_outputs=find_outputs)</div><div> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/data_objects/time_series.py", line 350, in __init__</div>
<div> self._get_all_outputs(find_outputs=find_outputs)</div><div> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/frontends/enzo/simulation_handling.py", line 434, in _get_all_outputs</div><div> self.all_outputs = self.all_time_outputs + self.all_redshift_outputs</div>
<div>AttributeError: 'EnzoSimulation' object has no attribute 'all_time_outputs'</div></div><div><br></div><div>Since the Enzo I need to use is a modified version of Enzo 1.5, I don't know if that's the problem and how I could fix it. Could you please give me some advice?</div>
<div><br></div><div>Also I am a little bit confused about how light ray and absorption spectrum works and have a few questions:</div><div>(1) What's the initial spectrum before any absorption? </div><div>(2) Where does the light ray originate from? Is it shoot from a random position in the first dataset? Can we set it manually like on the surface of the box? Does it go out of the box and back into it until the next output redshift?</div>
<div>(3) In the output file lightraysolution.txt, what do the parameters(dl/box, Start x/y/z, End x/y/z) mean?</div><div><br></div><div>Sorry that I have so many questions... I wonder if there is some existing documents or paper that have a more detailed description. </div>
<div><br></div><div>Thank you again for your help! Have a nice weekend!</div><div><br></div><div>Pengfei</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/9/6 Britton Smith <span dir="ltr"><<a href="mailto:brittonsmith@gmail.com" target="_blank">brittonsmith@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi Pengfei,<br><br></div>The near_redshift and far_redshift are free parameters for you to set, and just control the redshift interval over which to make the light cone projection. You might also want to look at the narrative documentation on the light cone generator, which you can find here:<br>
<a href="http://yt-project.org/doc/analysis_modules/light_cone_generator.html" target="_blank">http://yt-project.org/doc/analysis_modules/light_cone_generator.html</a><br></div>Note, that you need to use this with a simulation that made enough data outputs to connect the redshift interval you specify. The "enzo_tiny_cosmology" dataset available at <a href="http://yt-project.org/data/" target="_blank">http://yt-project.org/data/</a> was configured to work with the light cone recipe in the cookbook, so I would suggest downloading that.<br>
<br></div>As for the error, this is happening because there is a line in your enzo parameter file that contains more than one "=" sign. I have just submitted a pull request that fixes it, so this should be fixed soon. In the mean time, you can simply comment out that line in your parameter file and everything should work fine.<br>
<br></div>Britton<br></div><div class="gmail_extra"><br><br><div class="gmail_quote"><div><div class="h5">On Fri, Sep 6, 2013 at 3:02 AM, Pengfei Chen <span dir="ltr"><<a href="mailto:madcpf@gmail.com" target="_blank">madcpf@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 class="h5"><div dir="ltr">Hi all,<br><br>I downloaded the example script of light cone projectoin(<a href="http://yt-project.org/doc/cookbook/cosmological_analysis.html#light-cone-projection" target="_blank">http://yt-project.org/doc/cookbook/cosmological_analysis.html#light-cone-projection</a>) to my laptop and changed the Enzo simulation parameter file and redshifts. But when I ran it the following error occurred:<br>
<br> File "light_cone_projection.py", line 19, in <module><br> find_outputs=False)<br> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/analysis_modules/cosmological_observation/light_cone/light_cone.py", line 158, in __init__<br>
find_outputs=find_outputs)<br> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/analysis_modules/cosmological_observation/cosmology_splice.py", line 44, in __init__<br> find_outputs=find_outputs)<br> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/convenience.py", line 131, in simulation<br>
find_outputs=find_outputs)<br> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/frontends/enzo/simulation_handling.py", line 81, in __init__<br> find_outputs=find_outputs)<br> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/data_objects/time_series.py", line 343, in __init__<br>
self._parse_parameter_file()<br> File "/Users/madcpf/Documents/dev_yt/src/yt-hg/yt/frontends/enzo/simulation_handling.py", line 273, in _parse_parameter_file<br> param, vals = (i.strip() for i in line.split("="))<br>
ValueError: need more than 1 value to unpack<br><br>Does the parameters "near_redshift" and "far_redshift" in LightCone need to be the same with "CosmologyFinalRedshift" and "CosmologyInitialRedshift" in enzo parameter file?<br>
Does anyone know how to solve this?<br><br>Thank you all!<span><font color="#888888"><br><br>Pengfei<br></font></span></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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></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" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>