From 5f1c433f1181c640aeb302ed2b232757ff304760 Mon Sep 17 00:00:00 2001 From: joshi994 Date: Sat, 29 Jun 2024 14:04:06 -0700 Subject: [PATCH] Heat flux in slab module are updated --- phys/module_sf_slab.F | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/phys/module_sf_slab.F b/phys/module_sf_slab.F index 7304b5af44..71527aba0e 100644 --- a/phys/module_sf_slab.F +++ b/phys/module_sf_slab.F @@ -407,7 +407,7 @@ SUBROUTINE SLAB1D(J,T1D,QV1D,P1D,FLHC,FLQC, & ESG=SVP1*EXP(SVP2*(TG0(I)-SVPT0)/(TG0(I)-SVP3)) QSG=EP2*ESG/(PSFC(I)-ESG) THG=TSK(I)*(100./PSFC(I))**ROVCP - HFX(I)=FLHC(I)*(THG-THX(I)) + HFX(I)=FLHC(I)*CAPG(I)*(THG-THX(I)) QFX(I)=FLQC(I)*(QSG-QX(I)) LH(I)=QFX(I)*XLV QS(I)=HFX(I)+QFX(I)*XLV @@ -440,7 +440,7 @@ SUBROUTINE SLAB1D(J,T1D,QV1D,P1D,FLHC,FLQC, & -SVP3)) QSG=EP2*ESG/(PS-ESG) ! UPDATE FLUXES FOR NEW GROUND TEMPERATURE - HFXT=FLHC(I)*(THG-THX(I)) + HFXT=FLHC(I)*CAPG(I)*(THG-THX(I)) QFXT=FLQC(I)*(QSG-QX(I)) QS(I)=HFXT+QFXT*XLV ! SUM HFX AND QFX OVER SOIL TIMESTEPS