Changeset 1746

Timestamp:

Apr 19, 2013, 3:46:37 PM (12 years ago)

Author:

idelkadi

Message:

Correction de bug sur le calcul de l'eau condensee dans fisrtilp.F90

File:

• LMDZ5/trunk/libf/phylmd/fisrtilp.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -441 +441 @@
               zqn(i) = zq(i)
               rneb(i,k) = 1.0                  
-              zcond(i) = MAX(0.0,zqn(i)-zqs(i))
+              zcond(i) = MAX(0.0,zqn(i)-zqs(i))/(1+zdqs(i))
               rhcl(i,k)=1.0
            ELSE
-              zcond(i) = MAX(0.0,zqn(i)-zqs(i))*rneb(i,k)
+              zcond(i) = MAX(0.0,zqn(i)-zqs(i))*rneb(i,k)/(1+zdqs(i))
               rhcl(i,k)=(zqs(i)+zq(i)-zdelq)/2./zqs(i)
            ENDIF