From ee0f2c6f175accd01f9e04062b020019aa54c7a5 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Thu, 1 Jul 2021 10:45:11 -0500
Subject: [PATCH] separate is_day, is_night functions
---
modules/util/util.py | 16 ++++++++--------
1 file changed, 8 insertions(+), 8 deletions(-)
diff --git a/modules/util/util.py b/modules/util/util.py
index e86cd377..7cbed72c 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -185,21 +185,21 @@ def value_to_index(nda, value):
return idx
-# array solzen must be degrees, missing values must NaN. For small 50x50km regions only
-def is_night(solzen, test_angle=80.0, threshold=0.10):
+# array solzen must be degrees, missing values must NaN. For small roughly 50x50km regions only
+def is_day(solzen, test_angle=80.0):
solzen = solzen.flatten()
solzen = solzen[np.invert(np.isnan(solzen))]
- if len(solzen) == 0 or (np.sum(solzen > test_angle) / len(solzen)) > threshold:
- return True
- else:
+ if len(solzen) == 0 or np.sum(solzen <= test_angle) < len(solzen):
return False
+ else:
+ return True
# array solzen must be degrees, missing values must NaN. For small roughly 50x50km regions only
-def is_day(solzen, test_angle=75.0):
+def is_night(solzen, test_angle=100.0):
solzen = solzen.flatten()
solzen = solzen[np.invert(np.isnan(solzen))]
- if len(solzen) == 0 or np.sum(solzen <= test_angle) < len(solzen):
+ if len(solzen) == 0 or np.sum(solzen >= test_angle) < len(solzen):
return False
else:
- return True
+ return True
\ No newline at end of file
--
GitLab