From d23f00c3c7972accc08deaa9eb32b77b34f98718 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sat, 10 Apr 2021 15:20:30 -0500
Subject: [PATCH] initial commit for aircraft icing
---
modules/util/util.py | 10 ++++++++++
1 file changed, 10 insertions(+)
diff --git a/modules/util/util.py b/modules/util/util.py
index d21fc755..fd405d33 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -161,3 +161,13 @@ def is_night(solzen, test_angle=80.0, threshold=0.10):
return True
else:
return False
+
+
+# array solzen must be degrees, missing values must NaN. For small roughly 50x50km regions only
+def is_day(solzen, test_angle=75.0):
+ solzen = solzen.flatten()
+ solzen = solzen[np.invert(np.isnan(solzen))]
+ if len(solzen) == 0 or np.sum(solzen <= test_angle) < len(solzen):
+ return False
+ else:
+ return True
--
GitLab