Hi,
I am working on data from ramses and I have the following output files from
my simulations:
amr
grav
hydro
part and
info.txt
However, when I try to make projection plots, I am obtaining the error:
FutureWarning: Conversion of the second argument of issubdtype from `float`
> to `np.floating` is deprecated. In future, it will be treated as
> `np.float64 == np.dtype(float).type`.
>
> from ._conv import register_converters as _register_converters
>
> yt : [INFO ] 2018-01-19 16:31:46,982 Parameters: current_time
> = 0.254337758876396
>
> yt : [INFO ] 2018-01-19 16:31:46,982 Parameters: domain_dimensions
> = [32 32 32]
>
> yt : [INFO ] 2018-01-19 16:31:46,983 Parameters: domain_left_edge
> = [0. 0. 0.]
>
> yt : [INFO ] 2018-01-19 16:31:46,983 Parameters: domain_right_edge
> = [1. 1. 1.]
>
> yt : [INFO ] 2018-01-19 16:31:46,983 Parameters:
> cosmological_simulation = 0
>
> Traceback (most recent call last):
>
> File "pyt.py", line 11, in <module>
>
> p = yt.ProjectionPlot(ds, "x", "density")
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/visualization/plot_window.py",
> line 1441, in __init__
>
> validate_mesh_fields(test_data_source, fields)
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/visualization/plot_window.py",
> line 132, in validate_mesh_fields
>
> canonical_fields = data_source._determine_fields(fields)
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/data_objects/data_containers.py",
> line 1127, in _determine_fields
>
> finfo = self.ds._get_field_info("unknown", fname)
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/data_objects/static_output.py",
> line 711, in _get_field_info
>
> self.index
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/data_objects/static_output.py",
> line 504, in index
>
> self, dataset_type=self.dataset_type)
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/frontends/ramses/data_structures.py",
> line 357, in __init__
>
> super(RAMSESIndex, self).__init__(ds, dataset_type)
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/geometry/geometry_handler.py",
> line 50, in __init__
>
> self._setup_geometry()
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/geometry/oct_geometry_handler.py",
> line 25, in _setup_geometry
>
> self._initialize_oct_handler()
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/frontends/ramses/data_structures.py",
> line 361, in _initialize_oct_handler
>
> for i in range(self.dataset['ncpu'])]
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/frontends/ramses/data_structures.py",
> line 361, in <listcomp>
>
> for i in range(self.dataset['ncpu'])]
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/frontends/ramses/data_structures.py",
> line 67, in __init__
>
> self._read_amr_header()
>
> File
> "/Users/xxxx/Documents/yt-conda/lib/python3.6/site-packages/yt/frontends/ramses/data_structures.py",
> line 191, in _read_amr_header
>
> hvals.update(fpu.read_attrs(f, header))
>
> File
> "/Users/samhitha/Documents/yt-conda/lib/python3.6/site-packages/yt/utilities/fortran_utils.py",
> line 86, in read_attrs
>
> assert(s1 == s2)
>
> AssertionError
>
I tried specifying the fields externally using:
ds =
> yt.load("/Users/xxxx/SIM001/output_00002/info_00002.txt",fields=["density",
> "velocity_x", "velocity_y", "velocity_z", "thermal_pressure", "var6",
> "var7"])
>
p = yt.ProjectionPlot(ds, "x", "density")
PS: Rebooting the system didn't work in my case. Any help is greatly
appreciated.