Changeset 5163 for LMDZ6/branches/cirrus/libf/phylmd

Timestamp:

Aug 26, 2024, 11:42:24 AM (4 weeks ago)

Author:

aborella

Message:

Correctif calcul chaleur latente lscp

File:

• LMDZ6/branches/cirrus/libf/phylmd/lmdz_lscp.F90 (modified) (1 diff)

LMDZ6/branches/cirrus/libf/phylmd/lmdz_lscp.F90

@@ -1018 +1018 @@
                         ENDIF
 
-                        qlbef(i)=max(0.,zqn(i)-zqs(i))
+                        IF ( ok_unadjusted_clouds ) THEN
+                          !--AB We relax the saturation adjustment assumption
+                          !-- qvc (grid-mean vapor in cloud) is calculated by the condensation scheme
+                          IF ( rneb(i,k) .GT. eps ) THEN
+                            qlbef(i) = MAX(0., zqn(i) - qvc(i) / rneb(i,k))
+                          ELSE
+                            qlbef(i) = 0.
+                          ENDIF
+                        ELSE
+                          qlbef(i)=max(0.,zqn(i)-zqs(i))
+                        ENDIF
+
                         num = -Tbef(i)+zt(i)+rneb(i,k)*((1-zfice(i))*RLVTT &
                               +zfice(i)*RLSTT)/RCPD/(1.0+RVTMP2*(zq(i)+zmqc(i)))*qlbef(i)