From 9cb34607ea9e1011bb65c874d498c98d97685b8f Mon Sep 17 00:00:00 2001
From: sbiri <sbiri@noc.ac.uk>
Date: Mon, 8 Mar 2021 17:01:26 +0000
Subject: [PATCH] Update flux_subs.py
---
flux_subs.py | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/flux_subs.py b/flux_subs.py
index 5042daa..e4cadd1 100755
--- a/flux_subs.py
+++ b/flux_subs.py
@@ -978,8 +978,9 @@ def get_L(L, lat, usr, tsr, qsr, t10n, hin, Ta, sst, qair, qsea, q10n,
monob = 1/np.copy(temp)
elif (L == "ecmwf"):
tsrv = tsr*(1+0.6077*qair)+0.6077*Ta*qsr
- Rb = ((g*hin[1])/np.power(wind, 2))*(2*(Ta-sst)/(Ta+sst-g*hin[1]) +
- 0.6077*(qair-qsea))
+ Rb = g*hin[1]/(wind*wind)*((Ta-sst)/(Ta-0.5*(Ta-sst+g*hin[1] /
+ (1005+1860*qair))) +
+ 0.6077*(qair-qsea))
zo = (0.11*visc_air(Ta)/usr+0.018*np.power(usr, 2)/g)
zot = 0.40*visc_air(Ta)/usr
zol = (Rb*(np.power(np.log((hin[0]+zo)/zo)-psim_calc((hin[0]+zo) /
--
GitLab