diff --git a/AirSeaFluxCode.py b/AirSeaFluxCode.py
index 9997d78eed38be6c4dce85b6359029a705058008..cd131a8c18aa83d4f78043d05a49819177a1f1ed 100644
--- a/AirSeaFluxCode.py
+++ b/AirSeaFluxCode.py
@@ -407,6 +407,7 @@ def AirSeaFluxCode(spd, T, SST, lat=None, hum=None, P=None, hin=18, hout=10,
         utmp = np.copy(u10n)
         utmp = np.where(utmp < 0, np.nan, utmp)
         itera[ind] = np.ones(1)*it
+        rho = (0.34838*P)/(tv10n)
         sensible = -rho*cp*usr*tsr
         latent = -rho*lv*usr*qsr
         if (gust[0] == 1):