snapshot...

modules/util/util.py

@@ -151,3 +151,13 @@ def value_to_index(nda, value):
     diff = np.abs(nda - value)
     idx = np.argmin(diff)
     return idx
+
+
+# array solzen must be degrees, missing values must NaN
+def is_night(solzen, threshold):
+    solzen = solzen.flatten()
+    solzen = solzen[np.invert(np.isnan(solzen))]
+    if np.max(solzen) > threshold:
+        return True
+    else:
+        return False