Commit ab9085f1 authored by sbiri's avatar sbiri
Browse files

Update AirSeaFluxCode.py

parent 978547b9
...@@ -128,7 +128,7 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10, ...@@ -128,7 +128,7 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10,
40. flag ("n": normal, "o": out of nominal range, 40. flag ("n": normal, "o": out of nominal range,
"u": u10n<0, "q":q10n<0 "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,
"rh" : rh>100%, "r" : rh>100%,
"i": convergence fail at n) "i": convergence fail at n)
2021 / Author S. Biri 2021 / Author S. Biri
...@@ -187,7 +187,7 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10, ...@@ -187,7 +187,7 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10,
flag = np.empty(spd.shape, dtype="object") flag = np.empty(spd.shape, dtype="object")
flag[:] = "n" flag[:] = "n"
flag = np.where(np.isnan(spd+T+SST+hum[1]+P+Rs+Rl), "m", flag) flag = np.where(np.isnan(spd+T+SST+hum[1]+P+Rs+Rl), "m", flag)
flag = np.where(rh > 100, "rh", flag) flag = np.where(rh > 100, "r", flag)
dt = Ta - sst dt = Ta - sst
dq = qair - qsea dq = qair - qsea
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment