[yt-users] yt-3.0 devel with RAMSES

[Romain Teyssier]

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: