fixed sensible heat flux correction due to rainfall

gotm-model · Apr 15, 2024 · 204083f · 204083f
1 parent 5431037
commit 204083f
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git 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
@@ -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