diff --git a/src/airsea/fairall.F90 b/src/airsea/fairall.F90
index cd8e50d0..d9364eb9 100644
--- a/src/airsea/fairall.F90
+++ b/src/airsea/fairall.F90
@@ -249,7 +249,7 @@ subroutine fairall(sst,airt,u10,v10,precip,evap,taux,tauy,qe,qh)
                x3 = qa * L /(rgas * ta_K * ta_K)
                cd_rain = 1.0/(1.0+const06*(x3*L*x1)/(cpa*x2))
                cd_rain = cd_rain*cpw*((tw-ta) + (qs-qa)*L/cpa)
-               qe = qe - rainfall * cd_rain
+               qh = qh - rainfall * cd_rain
             end if
 
 !           Compute turbulent latent heat flux (W/m2), qe.
diff --git a/src/airsea/kondo.F90 b/src/airsea/kondo.F90
index 83aa11c3..96a994f9 100644
--- a/src/airsea/kondo.F90
+++ b/src/airsea/kondo.F90
@@ -198,7 +198,7 @@ subroutine kondo(sst,airt,u10,v10,precip,evap,taux,tauy,qe,qh)
       x3 = qa * L /(rgas * ta_K * ta_K)
       cd_rain = 1.0/(1.0+const06*(x3*L*x1)/(cpa*x2))
       cd_rain = cd_rain*cpw*((tw-ta) + (qs-qa)*L/cpa)
-      qe = qe - rainfall * cd_rain
+      qh = qh - rainfall * cd_rain
    end if
 
 !  calculation of evaporation/condensation in m/s