[yt-users] yt-3.0 devel with RAMSES

Sun Oct 20 03:03:18 PDT 2013

Everything works now like a charm.
Thanks !

Romain 

> Le 20 oct. 2013 à 11:03, Matthew Turk <matthewturk at gmail.com> a écrit :
> 
> 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
> 
> _______________________________________________
> 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/23e61afa/attachment.html>