---------- Forwarded message ----------
From:
Gongbo Zhao <gong-bo.zhao@port.ac.uk>Date: Thu, Oct 23, 2014 at 11:35 PM
Subject: [enzo-dev] A lightcone problem
To:
brittonsmith@gmail.comCc: "
enzo-dev@googlegroups.com" <
enzo-dev@googlegroups.com>
Hi Britton,
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,
> from yt.analysis_modules.cosmological_observation.api import \
> LightCone
>
> lc = LightCone('Example_Cosmology_Sim.param',
> 'Enzo', 0., 0.1)
> lc.calculate_light_cone_solution(seed=123456789, filename='lightcone.dat')
> field = 'density'
> field_of_view = (600.0, "arcmin")
> resolution = (60.0, "arcsec")
> lc.project_light_cone(field_of_view, resolution,
> field, weight_field=None,
> save_stack=True,
> save_slice_images=True)
I got these error messages,
> yt : [INFO ] 2014-10-24 04:29:36,439 Parameters: domain_dimensions = [128 128 128]
> yt : [INFO ] 2014-10-24 04:29:36,440 Parameters: domain_left_edge = [ 0. 0. 0.]
> yt : [INFO ] 2014-10-24 04:29:36,440 Parameters: domain_right_edge = [ 1. 1. 1.]
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: initial_time = 1.1703706995e+15 s
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: final_time = 4.31424173577e+17 s
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: stop_cycle = 100000
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: cosmological_simulation = 1
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: box_size = 64.0 Mpccm/h
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: omega_lambda = 0.73
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: omega_matter = 0.27
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: hubble_constant = 0.71
> yt : [INFO ] 2014-10-24 04:29:36,441 Parameters: initial_redshift = 60.0
> yt : [INFO ] 2014-10-24 04:29:36,442 Parameters: final_redshift = 0.0
> yt : [INFO ] 2014-10-24 04:29:36,442 Total datasets: 313.
> 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.
> 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
> 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.
> yt : [INFO ] 2014-10-24 04:29:50,042 Saving light cone solution to lightcone.dat.
> yt : [INFO ] 2014-10-24 04:29:50,217 Parameters: current_time = 272.816550757
> yt : [INFO ] 2014-10-24 04:29:50,217 Parameters: domain_dimensions = [128 128 128]
> yt : [INFO ] 2014-10-24 04:29:50,217 Parameters: domain_left_edge = [ 0. 0. 0.]
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: domain_right_edge = [ 1. 1. 1.]
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: cosmological_simulation = 1
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: current_redshift = 0.0994652856936
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: omega_lambda = 0.73
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: omega_matter = 0.27
> yt : [INFO ] 2014-10-24 04:29:50,218 Parameters: hubble_constant = 0.71
> yt : [INFO ] 2014-10-24 04:29:50,220 Making projection at z = 0.099445 from DD0272/DD0272.
> Parsing Hierarchy 100% ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Time: 00:00:00
> yt : [INFO ] 2014-10-24 04:29:51,714 Gathering a field list (this may take a moment.)
> yt : [INFO ] 2014-10-24 04:29:52,622 Loading field plugins.
> yt : [INFO ] 2014-10-24 04:29:52,622 Loaded angular_momentum (8 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,622 Loaded astro (14 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,623 Loaded cosmology (20 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,623 Loaded fluid (56 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,624 Loaded fluid_vector (88 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,624 Loaded geometric (102 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,624 Loaded local (102 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,625 Loaded magnetic_field (108 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,625 Loaded my_plugins (108 new fields)
> yt : [INFO ] 2014-10-24 04:29:52,625 Loaded species (110 new fields)
> Traceback (most recent call last):
> File "lightcone.py", line 13, in <module>
> save_slice_images=True)
> 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
> weight_field=weight_field)
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/utilities/parallel_tools/parallel_analysis_interface.py", line 267, in barrierize
> return func(*args, **kwargs)
> 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
> data_source=cut_region)
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 245, in __init__
> self.get_data(field)
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 285, in get_data
> self._initialize_chunk(chunk, tree)
> File "/opt/apps/enzo/new/yt-x86_64/src/yt-hg/yt/data_objects/construction_data_containers.py", line 369, in _initialize_chunk
> tree.initialize_chunk(i1, i2, ilevel)
> File "QuadTree.pyx", line 330, in yt.utilities.lib.QuadTree.QuadTree.initialize_chunk (yt/utilities/lib/QuadTree.c:4710)
> ValueError: Buffer has wrong number of dimensions (expected 1, got 2)
It seems to be a buffer dimension problem?
Thanks,
GB
--
You received this message because you are subscribed to the Google Groups "enzo-dev" group.
To unsubscribe from this group and stop receiving emails from it, send an email to enzo-dev+unsubscribe@googlegroups.com.
To post to this group, send email to enzo-dev@googlegroups.com.
Visit this group at http://groups.google.com/group/enzo-dev.
For more options, visit https://groups.google.com/d/optout.