From 8e540c20d5087c6372dd91b25969f7d1ddb601e6 Mon Sep 17 00:00:00 2001
From: sbiri <sbiri@noc.ac.uk>
Date: Mon, 8 Mar 2021 15:05:51 +0000
Subject: [PATCH] Rb in get_L according to eq. 3.25 instead of 3.24
---
flux_subs.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/flux_subs.py b/flux_subs.py
index 24fe5a6..5042daa 100755
--- a/flux_subs.py
+++ b/flux_subs.py
@@ -978,8 +978,8 @@ 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[0]/(Ta*(1+0.61*qair))) *
- ((t10n*(1+0.61*q10n)-sst*(1+0.61*qsea))/np.power(wind, 2)))
+ Rb = ((g*hin[1])/np.power(wind, 2))*(2*(Ta-sst)/(Ta+sst-g*hin[1]) +
+ 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