Make median calculation more robust

actris-cloudnet · Apr 24, 2024 · 351c3e5 · 351c3e5
1 parent 854ac9d
commit 351c3e5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/mwrpy/atmos.py b/mwrpy/atmos.py
@@ -179,7 +179,7 @@ def find_lwcl_free(lev1: dict) -> tuple[np.ndarray, np.ndarray]:
         tb_rat = tb_rat.rolling(
             pd.tseries.frequencies.to_offset("20min"), center=True, min_periods=100
         ).max()
-        index[tb_mx["Tb"] < np.median(tb_rat["Tb"]) * 0.1] = 0
+        index[tb_mx["Tb"] < np.nanmedian(tb_rat["Tb"]) * 0.1] = 0
 
         df = pd.DataFrame({"index": index}, index=ind)
         df = df.bfill(limit=120)