aboutsummaryrefslogtreecommitdiffstats
path: root/ReduceH5
diff options
context:
space:
mode:
authorVolker Hoffmann <volker@cheleb.net>2014-12-10 17:36:04 +0100
committerVolker Hoffmann <volker@cheleb.net>2014-12-10 17:36:04 +0100
commita855d0ce832bf2a9c77cb63e755582e51b4e4422 (patch)
tree35f06ac34b1333865adec7f86298ab7caa965987 /ReduceH5
parent166f4e5547a25784078d41135f4a1a8800e65f2b (diff)
add close encounter counts
Diffstat (limited to 'ReduceH5')
-rw-r--r--ReduceH5/reduce_genga.py2
1 files changed, 2 insertions, 0 deletions
diff --git a/ReduceH5/reduce_genga.py b/ReduceH5/reduce_genga.py
index 61a3f1c..fc32168 100644
--- a/ReduceH5/reduce_genga.py
+++ b/ReduceH5/reduce_genga.py
@@ -77,6 +77,7 @@ for istep, nstep in enumerate(nsteps):
vx = genga[:,7].astype("float64")
vy = genga[:,8].astype("float64")
vz = genga[:,9].astype("float64")
+ ce_count = genga[:,19].astype("float64")
if args.barycentric:
x, vx = kh.helio2bary(x, vx, m)
@@ -103,6 +104,7 @@ for istep, nstep in enumerate(nsteps):
f5.create_dataset("particles/Omega", data=Omega)
f5.create_dataset("particles/omega", data=omega)
f5.create_dataset("particles/M0", data=M0)
+ f5.create_dataset("particles/ce_count", data=ce_count)
f5["particles/x"].attrs["units"] = "AU"
f5["particles/y"].attrs["units"] = "AU"
f5["particles/z"].attrs["units"] = "AU"