From 988fe5b8ded65db94b0147391b4c140f6000b18e Mon Sep 17 00:00:00 2001
From: sbiri <sbiri@noc.ac.uk>
Date: Mon, 23 Aug 2021 08:20:04 +0100
Subject: [PATCH] added flag "l" to catch values where z/L>1000
---
AirSeaFluxCode.py | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/AirSeaFluxCode.py b/AirSeaFluxCode.py
index 8f1afec..afb2ef4 100644
--- a/AirSeaFluxCode.py
+++ b/AirSeaFluxCode.py
@@ -130,7 +130,8 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10,
41. lv latent heat of vaporization (Jkg−1)
40. flag ("n": normal, "o": out of nominal range,
"u": u10n<0, "q":q10n<0
- "m": missing, "l": Rib<-0.5 or Rib>0.2,
+ "m": missing,
+ "l": Rib<-0.5 or Rib>0.2 or z/L>1000,
"r" : rh>100%,
"i": convergence fail at n)
@@ -471,8 +472,10 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10,
flag = np.where((q10n < 0) & (flag == "n"), "q",
np.where((q10n < 0) & (flag != "n"), flag+[","]+["q"],
flag))
- flag = np.where(((Rb < -0.5) | (Rb > 0.2)) & (flag == "n"), "l",
- np.where(((Rb < -0.5) | (Rb > 0.2)) &
+ flag = np.where(((Rb < -0.5) | (Rb > 0.2) | ((hin[0]/monob) > 1000)) &
+ (flag == "n"), "l",
+ np.where(((Rb < -0.5) | (Rb > 0.2) |
+ ((hin[0]/monob) > 1000)) &
(flag != "n"), flag+[","]+["l"], flag))
if (out == 1):
flag = np.where((itera == -1) & (flag == "n"), "i",
--
GitLab