This commit is contained in:
Noe Brucy
2023-05-08 18:06:54 +02:00
parent 125f62a5cc
commit 2d8bc22482

View File

@@ -336,8 +336,8 @@ class Galsec:
for fluid in self.fluids: for fluid in self.fluids:
for i, (delta, name) in enumerate(zip([delta_x, delta_y, delta_z], ["x", "y", "z"])): for i, (delta, name) in enumerate(zip([delta_x, delta_y, delta_z], ["x", "y", "z"])):
pos = self.data[fluid]["position"][:, i] pos = self.data[fluid]["position"][:, i]
minval = np.min(pos) minval = np.floor(np.min(pos))
bin = np.floor((pos + minval) / delta) bin = np.floor((pos - minval) / delta)
# Store the middle value # Store the middle value
self.data[fluid][f"{name}_bin"] = (bin + 0.5) * delta + minval self.data[fluid][f"{name}_bin"] = (bin + 0.5) * delta + minval