<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Thanks Britton. I am attaching the lightcone solution file.</div><div><br></div><div>Cheers,</div><div>GB</div><div><br></div><div><div># parameter_filename = Example_Cosmology_Sim.param</div><div><br></div><div># Slice Dataset Redshift depth/box width/degree axis center</div><div>0000 DD0272/DD0272 0.099445 1.008644 0.079617 0 0.990649 0.013154 0.723866</div><div>0001 DD0278/DD0278 0.077086 0.988231 0.062013 2 0.805856 0.062935 0.861305</div><div>0002 DD0284/DD0284 0.055391 0.968416 0.044766 0 0.670546 0.643445 0.723627</div><div>0003 DD0290/DD0290 0.034326 0.949171 0.027864 0 0.611051 0.491222 0.350996</div><div>0004 DD0296/DD0296 0.013859 0.647291 0.011297 2 0.353542 0.185281 0.792768</div></div><br><div><div>On 28 Oct 2014, at 02:04, Britton Smith <<a href="mailto:brittonsmith@gmail.com">brittonsmith@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr">Hi GB,<div><br></div><div>I am forwarding this over to the yt-users list as this is a yt issue and not an enzo issue. I am unable to reproduce the issue you're seeing when I use the light cone cookbook recipe. I'm guessing that you're running into a corner case with the light cone generator that is specific to your simulation data. Could you send along the light cone solution file? It should be named <span style="font-family:arial,sans-serif;font-size:13px">lightcone.dat. There might be some clues in there.</span></div><div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Britton</span></div><div><br></div><div><br><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Gongbo Zhao</b> <span dir="ltr"><<a href="mailto:gong-bo.zhao@port.ac.uk">gong-bo.zhao@port.ac.uk</a>></span><br>Date: Thu, Oct 23, 2014 at 11:35 PM<br>Subject: [enzo-dev] A lightcone problem<br>To: <a href="mailto:brittonsmith@gmail.com">brittonsmith@gmail.com</a><br>Cc: "<a href="mailto:enzo-dev@googlegroups.com">enzo-dev@googlegroups.com</a>" <<a href="mailto:enzo-dev@googlegroups.com">enzo-dev@googlegroups.com</a>><br><br><br>Hi Britton,<br>
<br>
I used the latest yt (version Version = 3.1-dev; Changeset = bf8f9f98b34a (yt) tip) but I am having problems with the lightcone analysis. When I ran it using the following script,<br>
<br>
> from yt.analysis_modules.cosmological_observation.api import \<br>
> LightCone<br>
><br>
> lc = LightCone('Example_Cosmology_Sim.param',<br>
> 'Enzo', 0., 0.1)<br>
> lc.calculate_light_cone_solution(seed=123456789, filename='lightcone.dat')<br>
> field = 'density'<br>
> field_of_view = (600.0, "arcmin")<br>
> resolution = (60.0, "arcsec")<br>
> lc.project_light_cone(field_of_view, resolution,<br>
> field, weight_field=None,<br>
> save_stack=True,<br>
> save_slice_images=True)<br>
<br>
<br>
I got these error messages,<br>
<br>
> yt : [INFO ] 2014-10-24 04:29:36,439 Parameters: domain_dimensions = [128 128 128]<br>
> yt : [INFO ] 2014-10-24 04:29:36,440 Parameters: domain_left_edge = [ 0. 0. 0.]<br>
> yt : [INFO ] 2014-10-24 04:29:36,440 Parameters: domain_right_edge = [ 1. 1. 1.]<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: initial_time = 1.1703706995e+15 s<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: final_time = 4.31424173577e+17 s<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: stop_cycle = 100000<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: cosmological_simulation = 1<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: box_size = 64.0 Mpccm/h<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: omega_lambda = 0.73<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: omega_matter = 0.27<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: hubble_constant = 0.71<br>
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: initial_redshift = 60.0<br>
> yt : [INFO ] 2014-10-24 04:29:36,442 Parameters: final_redshift = 0.0<br>
> yt : [INFO ] 2014-10-24 04:29:36,442 Total datasets: 313.<br>
> yt : [INFO ] 2014-10-24 04:29:49,842 create_cosmology_splice: Used 5 data dumps to get from z = 0.100000 to 0.000000.<br>
> yt : [ERROR ] 2014-10-24 04:29:49,978 Warning: box fraction required to go from z = 0.099445 to 0.077086 is 1.008644<br>
> yt : [ERROR ] 2014-10-24 04:29:49,978 Full box delta z is 0.022168, but it is 0.022359 to the next data dump.<br>
> yt : [INFO ] 2014-10-24 04:29:50,042 Saving light cone solution to lightcone.dat.<br>
> yt : [INFO ] 2014-10-24 04:29:50,217 Parameters: current_time = 272.816550757<br>
> yt : [INFO ] 2014-10-24 04:29:50,217 Parameters: domain_dimensions = [128 128 128]<br>
> yt : [INFO ] 2014-10-24 04:29:50,217 Parameters: domain_left_edge = [ 0. 0. 0.]<br>
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: domain_right_edge = [ 1. 1. 1.]<br>
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: cosmological_simulation = 1<br>
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: current_redshift = 0.0994652856936<br>
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: omega_lambda = 0.73<br>
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: omega_matter = 0.27<br>
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: hubble_constant = 0.71<br>
> yt : [INFO ] 2014-10-24 04:29:50,220 Making projection at z = 0.099445 from DD0272/DD0272.<br>
> Parsing Hierarchy 100% ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Time: 00:00:00<br>
> yt : [INFO ] 2014-10-24 04:29:51,714 Gathering a field list (this may take a moment.)<br>
> yt : [INFO ] 2014-10-24 04:29:52,622 Loading field plugins.<br>
> yt : [INFO ] 2014-10-24 04:29:52,622 Loaded angular_momentum (8 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,622 Loaded astro (14 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,623 Loaded cosmology (20 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,623 Loaded fluid (56 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,624 Loaded fluid_vector (88 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,624 Loaded geometric (102 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,624 Loaded local (102 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,625 Loaded magnetic_field (108 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,625 Loaded my_plugins (108 new fields)<br>
> yt : [INFO ] 2014-10-24 04:29:52,625 Loaded species (110 new fields)<br>
> Traceback (most recent call last):<br>
> File "lightcone.py", line 13, in <module><br>
> save_slice_images=True)<br>
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/analysis_modules/cosmological_observation/light_cone/light_cone.py", line 322, in project_light_cone<br>
> weight_field=weight_field)<br>
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/utilities/parallel_tools/parallel_analysis_interface.py", line 267, in barrierize<br>
> return func(*args, **kwargs)<br>
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/analysis_modules/cosmological_observation/light_cone/light_cone_projection.py", line 84, in _light_cone_projection<br>
> data_source=cut_region)<br>
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 245, in __init__<br>
> self.get_data(field)<br>
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 285, in get_data<br>
> self._initialize_chunk(chunk, tree)<br>
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 369, in _initialize_chunk<br>
> tree.initialize_chunk(i1, i2, ilevel)<br>
> File "QuadTree.pyx", line 330, in yt.utilities.lib.QuadTree.QuadTree.initialize_chunk (yt/utilities/lib/QuadTree.c:4710)<br>
> ValueError: Buffer has wrong number of dimensions (expected 1, got 2)<br>
<br>
It seems to be a buffer dimension problem?<br>
<br>
Thanks,<br>
GB<br>
<span class=""><font color="#888888"><br>
--<br>
You received this message because you are subscribed to the Google Groups "enzo-dev" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:enzo-dev%2Bunsubscribe@googlegroups.com">enzo-dev+unsubscribe@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:enzo-dev@googlegroups.com">enzo-dev@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/enzo-dev" target="_blank">http://groups.google.com/group/enzo-dev</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</font></span></div><br></div></div>
</blockquote></div><br></body></html>