Index: LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F
===================================================================
--- LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F	(revision 26)
+++ LMDZ.3.3/trunk/libf/phylmd/fisrtilp_tr.F	(revision 27)
@@ -2,5 +2,6 @@
      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,4 +33,6 @@
       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,4 +161,11 @@
 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,4 +345,9 @@
          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