From 204083fe5f7b29c22f9b9efc6768116f7c30562b Mon Sep 17 00:00:00 2001 From: Knut Date: Tue, 16 Apr 2024 01:30:35 +0200 Subject: [PATCH] fixed sensible heat flux correction due to rainfall --- src/airsea/fairall.F90 | 2 +- src/airsea/kondo.F90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/airsea/fairall.F90 b/src/airsea/fairall.F90 index cd8e50d0a..d9364eb93 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 83aa11c3f..96a994f9b 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