diff --git a/AirSeaFluxCode/src/AirSeaFluxCode.py b/AirSeaFluxCode/src/AirSeaFluxCode.py
index 0f8ef48ebcd7a73e8a9464239cfa4e9a8d4a6da2..4a8d1f7df80fef0d10acc4a9b7bbab50fbf916f0 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 ec02db200bdb57cea8fea0abdacd2096e71a3ab9..549bee63d95d104e1de27d6df1c180c49b0644f1 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 665c7a6d1fca6ccb273852886f39696de446e41b..b387d5fbefddb44bceb2ca0aa4f0656e1eb2d0b9 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 cc33b3ee00cda656ca1ac4c408569935afd3b5c2..a72271f6bac67f8cdc2bc3e431ee1ca52d721367 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