[yt-users] yt-3.0 devel with RAMSES
Romain Teyssier
romain.teyssier at gmail.com
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.htm>
More information about the yt-users
mailing list