From cc425860c1fd27424c483940ecf580bac6d6a68a Mon Sep 17 00:00:00 2001
From: sbiri <sbiri@noc.ac.uk>
Date: Tue, 5 Jul 2022 05:37:51 +0000
Subject: [PATCH] Update AirSeaFluxCode/src/AirSeaFluxCode.py,
 AirSeaFluxCode/src/AirSeaFluxCode_dev.py, AirSeaFluxCode/src/flux_subs.py,
 AirSeaFluxCode/src/flux_subs_dev.py files

---
 AirSeaFluxCode/src/AirSeaFluxCode.py     | 4 ++--
 AirSeaFluxCode/src/AirSeaFluxCode_dev.py | 4 ++--
 AirSeaFluxCode/src/flux_subs.py          | 6 +++---
 AirSeaFluxCode/src/flux_subs_dev.py      | 6 +++---
 4 files changed, 10 insertions(+), 10 deletions(-)

diff --git a/AirSeaFluxCode/src/AirSeaFluxCode.py b/AirSeaFluxCode/src/AirSeaFluxCode.py
index 0f8ef48..4a8d1f7 100644
--- a/AirSeaFluxCode/src/AirSeaFluxCode.py
+++ b/AirSeaFluxCode/src/AirSeaFluxCode.py
@@ -692,8 +692,8 @@ def AirSeaFluxCode(spd, T, SST, SST_fl, meth, lat=None, hum=None, P=None,
         gust : int
             4x1 [x, beta, zi, ustb] x=0 gustiness is OFF, x=1-5 gustiness is ON
             and use gustiness factor: 1. Fairall et al. 2003, 2. GF is removed
-            from TSFs u10n, uref, 3. GF=1, 4. following Zeng et al. 1998 or
-            Brodeau et al. 2006, 5. following C35 matlab code;
+            from TSFs u10n, uref, 3. GF=1, 4. following ECMWF, 
+            4. following Zeng et al. 1998, 6. following C35 matlab code;
             beta gustiness parameter, default is 1.2,
             zi PBL height (m) default is 600,
             min is the value for gust speed in stable conditions,
diff --git a/AirSeaFluxCode/src/AirSeaFluxCode_dev.py b/AirSeaFluxCode/src/AirSeaFluxCode_dev.py
index ec02db2..549bee6 100644
--- a/AirSeaFluxCode/src/AirSeaFluxCode_dev.py
+++ b/AirSeaFluxCode/src/AirSeaFluxCode_dev.py
@@ -709,8 +709,8 @@ def AirSeaFluxCode_dev(spd, T, SST, SST_fl, meth, lat=None, hum=None, P=None,
         gust : int
             4x1 [x, beta, zi, ustb] x=0 gustiness is OFF, x=1-5 gustiness is ON
             and use gustiness factor: 1. Fairall et al. 2003, 2. GF is removed
-            from TSFs u10n, uref, 3. GF=1, 4. following Zeng et al. 1998 or
-            Brodeau et al. 2006, 5. following C35 matlab code;
+            from TSFs u10n, uref, 3. GF=1, 4. following ECMWF, 
+            4. following Zeng et al. 1998, 6. following C35 matlab code;
             beta gustiness parameter, default is 1.2,
             zi PBL height (m) default is 600,
             min is the value for gust speed in stable conditions,
diff --git a/AirSeaFluxCode/src/flux_subs.py b/AirSeaFluxCode/src/flux_subs.py
index 665c7a6..b387d5f 100755
--- a/AirSeaFluxCode/src/flux_subs.py
+++ b/AirSeaFluxCode/src/flux_subs.py
@@ -610,9 +610,9 @@ def apply_GF(gust, spd, wind, step):
         1: gustiness is switched ON following Fairall et al.
         2: gustiness is switched ON and GF is removed from TSFs u10n, uref
         3: gustiness is switched ON and GF=1
-        4: gustiness is switched ON following Zeng et al. 1998 or
-           Brodeau et al. 2006
-        5: gustiness is switched ON following C35 matlab code
+        4: gustiness is switched ON following ECMWF 
+        5: gustiness is switched ON following Zeng et al. (1998) 
+        6: gustiness is switched ON following C35 matlab code
     spd : float
         wind speed                      [ms^{-1}]
     wind : float
diff --git a/AirSeaFluxCode/src/flux_subs_dev.py b/AirSeaFluxCode/src/flux_subs_dev.py
index cc33b3e..a72271f 100644
--- a/AirSeaFluxCode/src/flux_subs_dev.py
+++ b/AirSeaFluxCode/src/flux_subs_dev.py
@@ -681,9 +681,9 @@ def apply_GF(gust, spd, wind, step):
         1: gustiness is switched ON following Fairall et al.
         2: gustiness is switched ON and GF is removed from TSFs u10n, uref
         3: gustiness is switched ON and GF=1
-        4: gustiness is switched ON following Zeng et al. 1998 or
-           Brodeau et al. 2006
-        5: gustiness is switched ON following C35 matlab code
+        4: gustiness is switched ON following ECMWF 
+        5: gustiness is switched ON following Zeng et al. (1998) 
+        6: gustiness is switched ON following C35 matlab code
     spd : float
         wind speed                      [ms^{-1}]
     wind : float
-- 
GitLab