diff --git a/AirSeaFluxCode/src/AirSeaFluxCode.py b/AirSeaFluxCode/src/AirSeaFluxCode.py
index 4a8d1f7df80fef0d10acc4a9b7bbab50fbf916f0..a1b5f916c34f123850076c4ff323acdef3129bc6 100644
--- a/AirSeaFluxCode/src/AirSeaFluxCode.py
+++ b/AirSeaFluxCode/src/AirSeaFluxCode.py
@@ -416,7 +416,6 @@ class S88:
 
         self.GustFact = self.wind/self.spd
         self.usr_gust = np.copy(self.usr)
-        self.usr_nogust = self.usr/self.GFo
         # include lapse rate adjustment as theta is well-mixed
         self.tref = self.theta-self.tlapse*self.h_out[1]-self.tsr/kappa * \
             (np.log(self.h_in[1]/self.h_out[1])-self.psit +
@@ -719,15 +718,13 @@ def AirSeaFluxCode(spd, T, SST, SST_fl, meth, lat=None, hum=None, P=None,
             set 1 to keep points
         out_var : str
            optional. user can define pandas array of variables to be output.
-           the default full pandas array is :
+           the default full pandas array, with cskin=0 gust=0, is :
                out_var = ("tau", "sensible", "latent", "monob", "cd", "cd10n",
                            "ct", "ct10n", "cq", "cq10n", "tsrv", "tsr", "qsr",
                            "usr", "psim", "psit", "psiq", "psim_ref", "psit_ref",
                            "psiq_ref", "u10n", "t10n", "q10n", "zo", "zot", "zoq",
-                           "uref", "tref", "qref", "dter", "dqer", "dtwl", "tkt",
-                           "qair", "qsea", "Rl", "Rs", "Rnl", "ug", "usrGF",
-                           "GustFact", "Rb", "rh", "rho", "cp", "lv", "theta",
-                           "itera")
+                           "uref", "tref", "qref", "qair", "qsea", "Rb", "rh",
+                           "rho", "cp", "lv", "theta", "itera")
             the "limited" pandas array is:
                 out_var = ("tau", "sensible", "latent", "uref", "tref", "qref")
             the user can define a custom pandas array of variables to  output.
@@ -778,7 +775,7 @@ def AirSeaFluxCode(spd, T, SST, SST_fl, meth, lat=None, hum=None, P=None,
                        37. downward shortwave radiation (Rs)
                        38. downward net longwave radiation (Rnl)
                        39. gust wind speed (ug)
-                       40. star wind speed/GustFact (usrGF)
+                       40. star wind speed with gust (usr_gust)
                        41. Gustiness Factor (GustFact)
                        42. Bulk Richardson number (Rb)
                        43. relative humidity (rh)
diff --git a/AirSeaFluxCode/src/AirSeaFluxCode_dev.py b/AirSeaFluxCode/src/AirSeaFluxCode_dev.py
index 549bee63d95d104e1de27d6df1c180c49b0644f1..72d93c318312b94c348c04489892c0f203f87afe 100644
--- a/AirSeaFluxCode/src/AirSeaFluxCode_dev.py
+++ b/AirSeaFluxCode/src/AirSeaFluxCode_dev.py
@@ -426,7 +426,6 @@ class S88:
 
         self.GustFact = self.wind/self.spd
         self.usr_gust = np.copy(self.usr)
-        self.usr_nogust = self.usr/self.GFo
         # include lapse rate adjustment as theta is well-mixed
         self.tref = self.theta-self.tlapse*self.h_out[1]-self.tsr/kappa * \
             (np.log(self.h_in[1]/self.h_out[1])-self.psit +
@@ -795,7 +794,7 @@ def AirSeaFluxCode_dev(spd, T, SST, SST_fl, meth, lat=None, hum=None, P=None,
                        37. downward shortwave radiation (Rs)
                        38. downward net longwave radiation (Rnl)
                        39. gust wind speed (ug)
-                       40. star wind speed/GustFact (usrGF)
+                       40. star wind with gust (usr_gust)
                        41. Gustiness Factor (GustFact)
                        42. Bulk Richardson number (Rb)
                        43. relative humidity (rh)
diff --git a/AirSeaFluxCode/src/util_subs.py b/AirSeaFluxCode/src/util_subs.py
index dbe44557f5297b36d4c016679b440b081b826fd1..ddf4f2e9c01bdf342ad614702ed422b8dbb55e29 100644
--- a/AirSeaFluxCode/src/util_subs.py
+++ b/AirSeaFluxCode/src/util_subs.py
@@ -191,7 +191,7 @@ def get_outvars(out_var, cskin, gust):
         elif cskin == 0 and gust[0] != 0:  # skin OFF and gust ON
             res_vars = ("tau", "sensible", "latent", "monob", "cd", "cd10n",
                         "ct", "ct10n", "cq", "cq10n", "tsrv", "tsr", "qsr",
-                        "usr", "usr_gust", "usr_nogust","ug", "GustFact",
+                        "usr_gust", "ug", "GustFact",
                         "psim", "psit", "psiq", "psim_ref", "psit_ref",
                         "psiq_ref", "u10n", "t10n", "q10n", "zo", "zot", "zoq",
                         "uref", "tref", "qref", "qair", "qsea",  "Rb", "rh",
@@ -206,7 +206,7 @@ def get_outvars(out_var, cskin, gust):
         else:
             res_vars = ("tau", "sensible", "latent", "monob", "cd", "cd10n",
                         "ct", "ct10n", "cq", "cq10n", "tsrv", "tsr", "qsr",
-                        "usr", "usr_gust", "usr_nogust","ug", "GustFact",
+                        "usr_gust", "ug", "GustFact",
                         "psim", "psit", "psiq", "psim_ref", "psit_ref",
                         "psiq_ref", "u10n", "t10n", "q10n", "zo", "zot", "zoq",
                         "uref", "tref", "qref", "dter", "dqer", "dtwl", "tkt",