From 360352b9b8a69db81b7888d8bce30a30acf5e61e Mon Sep 17 00:00:00 2001
From: sbiri <sbiri@noc.ac.uk>
Date: Thu, 11 Mar 2021 07:17:01 +0000
Subject: [PATCH] Update AirSeaFluxCode.py, flux_subs.py files
---
AirSeaFluxCode.py | 3 ---
flux_subs.py | 2 +-
2 files changed, 1 insertion(+), 4 deletions(-)
diff --git a/AirSeaFluxCode.py b/AirSeaFluxCode.py
index 58394a4..835127c 100644
--- a/AirSeaFluxCode.py
+++ b/AirSeaFluxCode.py
@@ -406,7 +406,6 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10,
zo = ref_ht/np.exp(kappa/np.sqrt(cd10n))
# adjust neutral cdn at any output height
cdn = np.power(kappa/np.log(hout/zo), 2)
- # cd = np.power(kappa/(np.log(hout/zo)-psim), 2)
cd = cd_calc(cdn, h_out[0], h_out[0], psim)
# solve for zot, zoq from ct10n, cq10n
zot = ref_ht/(np.exp(kappa**2/(ct10n*np.log(ref_ht/zo))))
@@ -414,8 +413,6 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10,
# adjust neutral ctn, cqn at any output height
ctn =np.power(kappa, 2)/(np.log(h_out[0]/zo)*np.log(h_out[1]/zot))
cqn =np.power(kappa, 2)/(np.log(h_out[0]/zo)*np.log(h_out[2]/zoq))
- # ct = np.power(kappa, 2)/((np.log(h_out[0]/zo)-psim)*(np.log(h_out[1]/zot)-psit))
- # cq = np.power(kappa, 2)/((np.log(h_out[0]/zo)-psim)*(np.log(h_out[2]/zoq)-psiq))
ct, cq = ctcq_calc(cdn, cd, ctn, cqn, h_out[1], h_out[2], h_out[1],
psit, psiq)
uref = (spd-usr/kappa*(np.log(h_in[0]/h_out[0])-psim +
diff --git a/flux_subs.py b/flux_subs.py
index f3e137f..bede98c 100755
--- a/flux_subs.py
+++ b/flux_subs.py
@@ -123,7 +123,7 @@ def cd_calc(cdn, hin, hout, psim):
cdn : float
neutral drag coefficient
hin : float
- wind speed height [m]
+ wind speed height [m]
hout : float
reference height [m]
psim : float
--
GitLab