Aquila-Consortium/csiborgtools
0
0
Fork 0
mirror of https://github.com/Richard-Sti/csiborgtools_public.git synced 2025-05-21 01:51:11 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Fix key bug

Browse source
This commit is contained in:
rstiskalek 2023-10-19 13:57:15 +01:00
parent 6eaeff2c7c
commit 637b31cf44
1 changed files with 2 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

4
csiborgtools/read/readsim.py
View file

@ -522,9 +522,9 @@ class QuijoteReader(BaseReader):
vel = fof.GroupVel * (1 + self.read_info(nsnap, nsim)["redshift"]) vel = fof.GroupVel * (1 + self.read_info(nsnap, nsim)["redshift"])
for i, p in enumerate(["x", "y", "z"]): for i, p in enumerate(["x", "y", "z"]):
data[p] = pos[:, i] data[p] = pos[:, i]
data["fof_v" + p] = vel[:, i] data[f"v{p}"] = vel[:, i]
data["group_mass"] = fof.GroupMass * 1e10 data["group_mass"] = fof.GroupMass * 1e10
data["fof_npart"] = fof.GroupLen data["npart"] = fof.GroupLen
# We want to start indexing from 1. Index 0 is reserved for # We want to start indexing from 1. Index 0 is reserved for
# particles unassigned to any FoF group. # particles unassigned to any FoF group.
data["index"] = 1 + numpy.arange(data.size, dtype=numpy.int32) data["index"] = 1 + numpy.arange(data.size, dtype=numpy.int32)