<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto"><div>Everything works now like a charm.</div><div>Thanks !<br><br>Romain </div><div><br>Le 20 oct. 2013 à 11:03, Matthew Turk <<a href="mailto:matthewturk@gmail.com">matthewturk@gmail.com</a>> a écrit :<br><br></div><blockquote type="cite"><div><div dir="ltr">Hi Romain,<div><br></div><div>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.</div><div><br></div>
<div>-Matt</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Oct 20, 2013 at 2:06 AM, Romain Teyssier <span dir="ltr"><<a href="mailto:romain.teyssier@gmail.com" target="_blank">romain.teyssier@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I have now installed yt-3.0 properly.<br>
But I am having issues reading a file with 7 hydro variables.<br>
<br>
r=load("mike/ramses/output_00040/info_00040.txt",fields=["Density","x-velocity","y-velocity","z-velocity","Pressure","Metal","g"])<br>
yt : [WARNING ] 2013-10-19 18:59:51,239 RAMSES frontend assumes all simulations are cosmological!<br>
yt : [INFO ] 2013-10-19 18:59:51,240 Parameters: current_time = 2.88490981599e+15<br>
yt : [INFO ] 2013-10-19 18:59:51,240 Parameters: domain_dimensions = [128 128 128]<br>
yt : [INFO ] 2013-10-19 18:59:51,241 Parameters: domain_left_edge = [ 0. 0. 0.]<br>
yt : [INFO ] 2013-10-19 18:59:51,242 Parameters: domain_right_edge = [ 1. 1. 1.]<br>
yt : [INFO ] 2013-10-19 18:59:51,242 Parameters: cosmological_simulation = 1<br>
yt : [INFO ] 2013-10-19 18:59:51,243 Parameters: current_redshift = 0.0<br>
yt : [INFO ] 2013-10-19 18:59:51,243 Parameters: omega_lambda = 0.0<br>
yt : [INFO ] 2013-10-19 18:59:51,243 Parameters: omega_matter = 1.0<br>
yt : [INFO ] 2013-10-19 18:59:51,244 Parameters: hubble_constant = 0.01<br>
<br>
p=ProjectionPlot(r,'z','Density')<br>
You are running with the wrong number of fields.<br>
Please specify these in the load command.<br>
We are looking for 6 fields.<br>
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]<br>
---------------------------------------------------------------------------<br>
AssertionError Traceback (most recent call last)<br>
<ipython-input-39-c02d2ecd2249> in <module>()<br>
----> 1 p=ProjectionPlot(r,'z','Density')<br>
<br>
/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)<br>
1439 proj = pf.h.proj(fields, axis, weight_field=weight_field,<br>
1440 center=center, data_source=data_source,<br>
-> 1441 field_parameters = field_parameters)<br>
1442 PWViewerMPL.__init__(self, proj, bounds, fields=fields, origin=origin,<br>
1443 fontsize=fontsize)<br>
<br>
/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)<br>
232 self.data_source = data_source<br>
233 self.weight_field = weight_field<br>
--> 234 self.get_data(field)<br>
235<br>
236 @property<br>
<br>
/Users/teyssier/yt-3.0/yt/data_objects/construction_data_containers.py in get_data(self, fields)<br>
276 # This needs to be parallel_objects-ified<br>
277 for chunk in parallel_objects(self.data_source.chunks(<br>
--> 278 chunk_fields, "io")):<br>
279 mylog.debug("Adding chunk (%s) to tree (%0.3e GB RAM)", chunk.ires.size,<br>
280 get_memory_usage()/1024.)<br>
<br>
/Users/teyssier/yt-3.0/yt/utilities/parallel_tools/parallel_analysis_interface.py in parallel_objects(objects, njobs, storage, barrier, dynamic)<br>
432 # this will prevent intermediate objects from being created.<br>
433 oiter = itertools.islice(enumerate(objects), my_new_id, None, njobs)<br>
--> 434 for obj_id, obj in oiter:<br>
435 result_id = obj_id * njobs + my_new_id<br>
436 if storage is not None:<br>
<br>
/Users/teyssier/yt-3.0/yt/data_objects/data_containers.py in chunks(self, fields, chunking_style, **kwargs)<br>
451 for chunk in self.hierarchy._chunk(self, chunking_style, **kwargs):<br>
452 with self._chunked_read(chunk):<br>
--> 453 self.get_data(fields)<br>
454 # NOTE: we yield before releasing the context<br>
455 yield self<br>
<br>
/Users/teyssier/yt-3.0/yt/data_objects/data_containers.py in get_data(self, fields)<br>
529 # need to be generated.<br>
530 read_fluids, gen_fluids = self.hierarchy._read_fluid_fields(<br>
--> 531 fluids, self, self._current_chunk)<br>
532 self.field_data.update(read_fluids)<br>
533<br>
<br>
/Users/teyssier/yt-3.0/yt/geometry/geometry_handler.py in _read_fluid_fields(self, fields, dobj, chunk)<br>
512 selector,<br>
513 fields_to_read,<br>
--> 514 chunk_size)<br>
515 for field in fields_to_read:<br>
516 ftype, fname = field<br>
<br>
/Users/teyssier/yt-3.0/yt/frontends/ramses/io.py in _read_fluid_selection(self, chunks, selector, fields, size)<br>
39 # and pick off the right vectors<br>
40 content = cStringIO.StringIO(f.read())<br>
---> 41 rv = subset.fill(content, fields, selector)<br>
42 for ft, f in fields:<br>
43 d = rv.pop(f)<br>
<br>
/Users/teyssier/yt-3.0/yt/frontends/ramses/data_structures.py in fill(self, content, fields, selector)<br>
302 for field in fields:<br>
303 tr[field] = np.zeros(cell_count, 'float64')<br>
--> 304 for level, offset in enumerate(self.domain.hydro_offset):<br>
305 if offset == -1: continue<br>
306 content.seek(offset)<br>
<br>
/Users/teyssier/yt-3.0/yt/frontends/ramses/data_structures.py in hydro_offset(self)<br>
100 ('file_ncache', 1, 'I') )<br>
101 try:<br>
--> 102 hvals = fpu.read_attrs(f, header, "=")<br>
103 except AssertionError:<br>
104 print "You are running with the wrong number of fields."<br>
<br>
/Users/teyssier/yt-3.0/yt/utilities/fortran_utils.py in read_attrs(f, attrs, endian)<br>
72 if s1 != s2:<br>
73 size = struct.calcsize(endian + "I" + "".join(n*[t]) + "I")<br>
---> 74 assert(s1 == s2)<br>
75 if n == 1: v = v[0]<br>
76 if type(a)==tuple:<br>
<br>
AssertionError:<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
</blockquote></div><br></div>
</div></blockquote><blockquote type="cite"><div><span>_______________________________________________</span><br><span>yt-users mailing list</span><br><span><a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a></span><br><span><a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a></span><br></div></blockquote></body></html>