Changeset 517 for LMDZ.3.3/branches/rel-LF/libf/phylmd/fisrtilp_tr.F

Timestamp:

Apr 16, 2004, 5:43:38 PM (20 years ago)

Author:

lmdzadmin

Message:

Inclusion des modifications de O. Boucher et de J. Quaas pour le calcul des
premiers effets directs et indirects dus aux aerosols
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/fisrtilp_tr.F (modified) (3 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/fisrtilp_tr.F

@@ -5 +5 @@
      s                   pfrac_impa, pfrac_nucl, pfrac_1nucl,
      s                   frac_impa, frac_nucl,
-     s                   prfl, psfl)
+     s                   prfl, psfl,
+     s                   RHcl) ! relative humidity in clear sky (needed for aer optical properties; aeropt.F)
 
 c
@@ -37 +38 @@
       REAL prfl(klon,klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s)
       REAL psfl(klon,klev+1) ! flux d'eau precipitante aux interfaces (kg/m2/s)
+      
+Cjq   For aerosol opt properties needed (see aeropt.F)
+      REAL RHcl(klon,klev)
+      
 cAA
 c Coeffients de fraction lessivee : pour OFF-LINE
@@ -290 +295 @@
             rneb(i,k) = MAX(0.0,MIN(1.0,rneb(i,k)))
             zcond(i) = MAX(0.0,zqn(i)-zqs(i))*rneb(i,k)/(1.+zdqs(i))
+            
+c--Olivier
+            RHcl(i,k)=(zqs(i)+zq(i)-zdelq)/2./zqs(i)
+            IF (rneb(i,k) .LE. 0.0) RHcl(i,k)=zq(i)/zqs(i)
+            IF (rneb(i,k) .GE. 1.0) RHcl(i,k)=1.0
+c--fin
+            
          ENDDO
       ELSE