[yt-users] yt-3.0 devel with RAMSES
Matthew Turk
matthewturk at gmail.com
Sun Oct 20 02:03:12 PDT 2013
Hi Romain,
Did this get fixed with the reboot, like the other issues? My reading of
your command is that it's correct, so this not working is a bug to fix on
our end.
-Matt
On Sun, Oct 20, 2013 at 2:06 AM, Romain Teyssier
<romain.teyssier at gmail.com>wrote:
> I have now installed yt-3.0 properly.
> But I am having issues reading a file with 7 hydro variables.
>
>
> r=load("mike/ramses/output_00040/info_00040.txt",fields=["Density","x-velocity","y-velocity","z-velocity","Pressure","Metal","g"])
> yt : [WARNING ] 2013-10-19 18:59:51,239 RAMSES frontend assumes all
> simulations are cosmological!
> yt : [INFO ] 2013-10-19 18:59:51,240 Parameters: current_time
> = 2.88490981599e+15
> yt : [INFO ] 2013-10-19 18:59:51,240 Parameters: domain_dimensions
> = [128 128 128]
> yt : [INFO ] 2013-10-19 18:59:51,241 Parameters: domain_left_edge
> = [ 0. 0. 0.]
> yt : [INFO ] 2013-10-19 18:59:51,242 Parameters: domain_right_edge
> = [ 1. 1. 1.]
> yt : [INFO ] 2013-10-19 18:59:51,242 Parameters:
> cosmological_simulation = 1
> yt : [INFO ] 2013-10-19 18:59:51,243 Parameters: current_redshift
> = 0.0
> yt : [INFO ] 2013-10-19 18:59:51,243 Parameters: omega_lambda
> = 0.0
> yt : [INFO ] 2013-10-19 18:59:51,243 Parameters: omega_matter
> = 1.0
> yt : [INFO ] 2013-10-19 18:59:51,244 Parameters: hubble_constant
> = 0.01
>
> p=ProjectionPlot(r,'z','Density')
> You are running with the wrong number of fields.
> Please specify these in the load command.
> We are looking for 6 fields.
> The last set of field sizes was: [2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
> 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
> 2, 2, 2, 2, 2, 2, 2, 2, 2, 2]
> ---------------------------------------------------------------------------
> AssertionError Traceback (most recent call last)
> <ipython-input-39-c02d2ecd2249> in <module>()
> ----> 1 p=ProjectionPlot(r,'z','Density')
>
> /Users/teyssier/yt-3.0/yt/visualization/plot_window.py in __init__(self,
> pf, axis, fields, center, width, axes_unit, weight_field, max_level,
> origin, fontsize, field_parameters, data_source)
> 1439 proj = pf.h.proj(fields, axis, weight_field=weight_field,
> 1440 center=center, data_source=data_source,
> -> 1441 field_parameters = field_parameters)
> 1442 PWViewerMPL.__init__(self, proj, bounds, fields=fields,
> origin=origin,
> 1443 fontsize=fontsize)
>
> /Users/teyssier/yt-3.0/yt/data_objects/construction_data_containers.py in
> __init__(self, field, axis, weight_field, center, pf, data_source, style,
> field_parameters)
> 232 self.data_source = data_source
> 233 self.weight_field = weight_field
> --> 234 self.get_data(field)
> 235
> 236 @property
>
> /Users/teyssier/yt-3.0/yt/data_objects/construction_data_containers.py in
> get_data(self, fields)
> 276 # This needs to be parallel_objects-ified
> 277 for chunk in parallel_objects(self.data_source.chunks(
> --> 278 chunk_fields, "io")):
> 279 mylog.debug("Adding chunk (%s) to tree (%0.3e GB
> RAM)", chunk.ires.size,
> 280 get_memory_usage()/1024.)
>
> /Users/teyssier/yt-3.0/yt/utilities/parallel_tools/parallel_analysis_interface.py
> in parallel_objects(objects, njobs, storage, barrier, dynamic)
> 432 # this will prevent intermediate objects from being created.
> 433 oiter = itertools.islice(enumerate(objects), my_new_id, None,
> njobs)
> --> 434 for obj_id, obj in oiter:
> 435 result_id = obj_id * njobs + my_new_id
> 436 if storage is not None:
>
> /Users/teyssier/yt-3.0/yt/data_objects/data_containers.py in chunks(self,
> fields, chunking_style, **kwargs)
> 451 for chunk in self.hierarchy._chunk(self, chunking_style,
> **kwargs):
> 452 with self._chunked_read(chunk):
> --> 453 self.get_data(fields)
> 454 # NOTE: we yield before releasing the context
> 455 yield self
>
> /Users/teyssier/yt-3.0/yt/data_objects/data_containers.py in
> get_data(self, fields)
> 529 # need to be generated.
> 530 read_fluids, gen_fluids =
> self.hierarchy._read_fluid_fields(
> --> 531 fluids, self,
> self._current_chunk)
> 532 self.field_data.update(read_fluids)
> 533
>
> /Users/teyssier/yt-3.0/yt/geometry/geometry_handler.py in
> _read_fluid_fields(self, fields, dobj, chunk)
> 512 selector,
> 513 fields_to_read,
> --> 514 chunk_size)
> 515 for field in fields_to_read:
> 516 ftype, fname = field
>
> /Users/teyssier/yt-3.0/yt/frontends/ramses/io.py in
> _read_fluid_selection(self, chunks, selector, fields, size)
> 39 # and pick off the right vectors
> 40 content = cStringIO.StringIO(f.read())
> ---> 41 rv = subset.fill(content, fields, selector)
> 42 for ft, f in fields:
> 43 d = rv.pop(f)
>
> /Users/teyssier/yt-3.0/yt/frontends/ramses/data_structures.py in
> fill(self, content, fields, selector)
> 302 for field in fields:
> 303 tr[field] = np.zeros(cell_count, 'float64')
> --> 304 for level, offset in enumerate(self.domain.hydro_offset):
> 305 if offset == -1: continue
> 306 content.seek(offset)
>
> /Users/teyssier/yt-3.0/yt/frontends/ramses/data_structures.py in
> hydro_offset(self)
> 100 ('file_ncache', 1, 'I') )
> 101 try:
> --> 102 hvals = fpu.read_attrs(f, header, "=")
> 103 except AssertionError:
> 104 print "You are running with the wrong number
> of fields."
>
> /Users/teyssier/yt-3.0/yt/utilities/fortran_utils.py in read_attrs(f,
> attrs, endian)
> 72 if s1 != s2:
> 73 size = struct.calcsize(endian + "I" + "".join(n*[t]) +
> "I")
> ---> 74 assert(s1 == s2)
> 75 if n == 1: v = v[0]
> 76 if type(a)==tuple:
>
> AssertionError:
>
>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20131020/0ca86919/attachment.htm>
More information about the yt-users
mailing list