From bad9d7d38034074d2a7175147a8c520bc33e3019 Mon Sep 17 00:00:00 2001
From: Volker Hoffmann <volker@cheleb.net>
Date: Sun, 18 Aug 2013 12:43:34 +0200
Subject: halo condition Omega_xy <= 1.0e-6 instead of == 0 [there are no
 perfect zeros due to our sampling and/or the reset]

---
 DiskXYZ.py | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/DiskXYZ.py b/DiskXYZ.py
index e2f070f..b593f2d 100644
--- a/DiskXYZ.py
+++ b/DiskXYZ.py
@@ -491,7 +491,7 @@ class DiskReduce(DiskReduceBase):
 
     def compute_Omega_xy(self):
         self.Omega_xy = self.vcyl_xy["vtheta"] / self.rt["r"]
-        self.Omega_xy[self.Omega_xy == 0.0] = np.nan
+        self.Omega_xy[self.Omega_xy <= 1.0e-6] = np.nan
 
     def compute_h_xy(self):
         self.h_xy = np.sqrt(self.cs2_xy) / self.Omega_xy
-- 
cgit v1.1