[yt-users] Error in reading the hydrodynamics dataset in Ramses

Matthew Turk matthewturk at gmail.com
Wed Nov 20 13:51:10 PST 2013


Hi Junhwan,

How many variables *does* your RAMSES dataset have?

-Matt

On Wed, Nov 20, 2013 at 4:50 PM, Junhwan Choi (최준환)
<choi.junhwan at gmail.com> wrote:
> Thank you, Kearn.
>
> I tried your suggestion that I used
> ds = load("/lustre/widow2/proj/ast031/test11_rhd_fesc=0.05/output_00048/info_00048.txt",
> fields = ["Density","x-velocity", "y-velocity",
> "z-velocity","Pressure"]).
> But, it results in the identical error message.
> I still do not know what I need to do.
>
> Regard,
> Junhwan
>
>
>
> On Wed, Nov 20, 2013 at 10:37 AM,  <k.grisdale at surrey.ac.uk> wrote:
>> Hi Junhwan
>>
>> I think the problem you might the having is you need to specify the fields you read into yt. To do this just modify you load command to this:
>>
>> ds = load("/lustre/widow2/proj/ast031/test11_rhd_fesc=0.05/output_00048/info_00048.txt”, fields = ["Density","x-velocity", "y-velocity", "z-velocity", "Pressure”])
>>
>> (Assuming that you are running a 3d Simulation with NVAR set as 5)
>>
>> Hope that helps
>>
>> Kearn
>>
>>
>> On 20 Nov 2013, at 16:26, Junhwan Choi (최준환) <choi.junhwan at gmail.com> wrote:
>>
>>> Hi yt users,
>>>
>>> I try to make a basic visualization for the Ramses simulation using
>>> following script:
>>> ==========================
>>> from yt.mods import *
>>>
>>> ds = load("/lustre/widow2/proj/ast031/test11_rhd_fesc=0.05/output_00048/info_00048.txt")
>>> center = [0.5, 0.5, 0.5]
>>>
>>> pw = ProjectionPlot(ds, "x", ("deposit", "all_density"))
>>> pw.zoom(1.01)
>>> pw.save("allview")
>>>
>>> pw = ProjectionPlot(ds, "x", "Density")
>>> pw.zoom(1.01)
>>> pw.save("allviewGas")
>>> ============================
>>> It WORKS for the particle projection (for allview) but it does not
>>> work for gas projection.
>>>
>>> When the yt try to deal with gas dataset, I got the following error
>>> message and the script crash.
>>> 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]
>>> Traceback (most recent call last):
>>>  File "DenVis_DM_Gas.py", line 10, in <module>
>>>    pw = ProjectionPlot(ds, "x", ("gas", "Density"))
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/visualization/plot_window.py",
>>> line 1441, in __init__
>>>    field_parameters = field_parameters)
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/data_objects/construction_data_containers.py",
>>> line 234, in __init__
>>>    self.get_data(field)
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/data_objects/construction_data_containers.py",
>>> line 278, in get_data
>>>    chunk_fields, "io")):
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/utilities/parallel_tools/parallel_analysis_interface.py",
>>> line 434, in parallel_objects
>>>    for obj_id, obj in oiter:
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/data_objects/data_containers.py",
>>> line 453, in chunks
>>>    self.get_data(fields)
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/data_objects/data_containers.py",
>>> line 531, in get_data
>>>    fluids, self, self._current_chunk)
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/geometry/geometry_handler.py",
>>> line 520, in _read_fluid_fields
>>>    chunk_size)
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/frontends/ramses/io.py",
>>> line 41, in _read_fluid_selection
>>>    rv = subset.fill(content, fields, selector)
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/frontends/ramses/data_structures.py",
>>> line 304, in fill
>>>    for level, offset in enumerate(self.domain.hydro_offset):
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/frontends/ramses/data_structures.py",
>>> line 102, in hydro_offset
>>>    hvals = fpu.read_attrs(f, header, "=")
>>>  File "/autofs/na3_home1/jhchoi/common/src/yt-hg/yt/utilities/fortran_utils.py",
>>> line 74, in read_attrs
>>>    assert(s1 == s2)
>>> AssertionError
>>>
>>> I do not understand what is the problem since it is such a basic visualization.
>>> Since I am new for Ramses data, I may take a simple mistake.
>>> Can anyone explain what the error message means and how to resolve the problem?
>>>
>>> Thank you in advance,
>>> Junhwan
>>> _______________________________________________
>>> 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
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org



More information about the yt-users mailing list