Changeset 27

Timestamp:

Jan 19, 2000, 2:35:32 PM (24 years ago)

Author:

lmdz

Message:

modif de physiq et fistrlp.F pour recuperer les flux d'eau precipitante (cf Olivier)

4. Hauglustaine

File:

• LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F (modified) (4 diffs)

LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F

@@ -2 +2 @@
      s                   d_t, d_q, d_ql, rneb, radliq, rain, snow,
      s                   pfrac_impa, pfrac_nucl, pfrac_1nucl,
-     s                   frac_impa, frac_nucl )
+     s                   frac_impa, frac_nucl,
+     s                   prfl, psfl)
 
 c
@@ -32 +33 @@
       REAL rain(klon) ! pluies (mm/s)
       REAL snow(klon) ! neige (mm/s)
+      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)
 cAA
 c Coeffients de fraction lessivee : pour OFF-LINE
@@ -158 +161 @@
 c Initialiser les sorties:
 c
+      DO k = 1, klev+1
+      DO i = 1, klon
+         prfl(i,k) = 0.0
+         psfl(i,k) = 0.0
+      ENDDO
+      ENDDO
+
       DO k = 1, klev
       DO i = 1, klon
@@ -335 +345 @@
          zrfl(i) = zrfl(i)+ MAX(zcond(i)-zoliq(i),0.0)
      .                    * (paprs(i,k)-paprs(i,k+1))/(RG*dtime)
+         IF (zt(i).LT.RTT) THEN
+           psfl(i,k)=zrfl(i)
+         ELSE
+           prfl(i,k)=zrfl(i)
+         ENDIF
       ENDIF
       ENDDO