Index: trunk/LMDZ.MARS/libf/dyn3d/calfis.F
===================================================================
--- trunk/LMDZ.MARS/libf/dyn3d/calfis.F	(revision 1266)
+++ trunk/LMDZ.MARS/libf/dyn3d/calfis.F	(revision 1312)
@@ -43,5 +43,5 @@
 c       pmasse          masse d'air dans chaque maille
 c       pts             surface temperature  (K)
-c       pw              flux vertical (kg m-2)
+c       pw              flux vertical (kg/s)
 c
 c    Output :
@@ -305,28 +305,28 @@
 c On interpole "lineairement" la temperature entre les couches(FF,10/95)
 
-      DO ig=1,ngridmx
-         zvervel(ig,1)=0.
-      END DO
-      DO l=2,llm
-        zvervel(1,l)=(pw(1,1,l)/apoln)
-     &  * r *0.5*(ztfi(1,l)+ztfi(1,l-1)) /zplev(1,l)              
-        ig0=2
-       DO j=2,jjm
-           DO i = 1, iim
-              zvervel(ig0,l) = pw(i,j,l) * unsaire(i,j)
-     &        * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
-              ig0 = ig0 + 1
-           ENDDO
-       ENDDO
-        zvervel(ig0,l)=(pw(1,jjp1,l)/apols)
-     &  * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
-      ENDDO
+!      DO ig=1,ngridmx
+!         zvervel(ig,1)=0.
+!      END DO
+!      DO l=2,llm
+!        zvervel(1,l)=(pw(1,1,l)/apoln)
+!     &  * r *0.5*(ztfi(1,l)+ztfi(1,l-1)) /zplev(1,l)              
+!        ig0=2
+!       DO j=2,jjm
+!           DO i = 1, iim
+!              zvervel(ig0,l) = pw(i,j,l) * unsaire(i,j)
+!     &        * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
+!              ig0 = ig0 + 1
+!           ENDDO
+!       ENDDO
+!        zvervel(ig0,l)=(pw(1,jjp1,l)/apols)
+!     &  * r *0.5*(ztfi(ig0,l)+ztfi(ig0,l-1)) /zplev(ig0,l)              
+!      ENDDO
 
 c    .........  Reindexation : calcul de zvervel au MILIEU des couches
-       DO l=1,llm-1
-	      DO ig=1,ngridmx
-		     zvervel(ig,l) = 0.5*(zvervel(ig,l)+zvervel(ig,l+1))
-          END DO 
-       END DO 
+!       DO l=1,llm-1
+!	      DO ig=1,ngridmx
+!		     zvervel(ig,l) = 0.5*(zvervel(ig,l)+zvervel(ig,l+1))
+!          END DO 
+!       END DO 
 c      (dans la couche llm, on garde la valeur à la limite inférieure llm)
 
@@ -426,5 +426,6 @@
      ,     zplev,zplay,zphi,
      ,     zufi, zvfi,ztfi, zqfi,  
-     ,     zvervel,
+!     ,     zvervel,
+     ,     pw,
 C - sorties
      s     zdufi, zdvfi, zdtfi, zdqfi,zdpsrf,tracer)