initial commit for aircraft icing

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