--- trunk/libf/dyn3d/convflu.f	2008/02/27 13:16:39	3
+++ trunk/libf/dyn3d/convflu.f90	2012/07/26 14:37:37	62
@@ -1,62 +1,51 @@
-!
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/convflu.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
-!
-      SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
-c
-c  P. Le Van
-c
-c
-c    *******************************************************************
-c  ... calcule la (convergence horiz. * aire locale)du flux ayant pour
-c      composantes xflu et yflu ,variables extensives .  ......
-c    *******************************************************************
-c      xflu , yflu et nbniv sont des arguments d'entree pour le s-pg ..
-c      convfl                est  un argument de sortie pour le s-pg .
-c
-c     njxflu  est le nombre de lignes de latitude de xflu, 
-c     ( = jjm ou jjp1 )
-c     nbniv   est le nombre de niveaux vert. de  xflu et de yflu .
-c
-      use dimens_m
-      use paramet_m
-      use comgeom
-      IMPLICIT NONE
-c
-      REAL       xflu,yflu,convfl,convpn,convps
-      INTEGER    l,ij,nbniv
-      DIMENSION  xflu( ip1jmp1,nbniv ),yflu( ip1jm,nbniv ) ,
-     *         convfl( ip1jmp1,nbniv )
-c
-      REAL       SSUM
-c
-c
-c
-      DO 5 l = 1,nbniv
-c
-      DO 2  ij = iip2, ip1jm - 1
-      convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   +
-     *                      yflu(ij +1,l ) - yflu( ij -iim,l )
-   2  CONTINUE
-c
-c
-
-c     ....  correction pour  convfl( 1,j,l)  ......
-c     ....   convfl(1,j,l)= convfl(iip1,j,l) ...
-c
-CDIR$ IVDEP
-      DO 3 ij = iip2,ip1jm,iip1
-      convfl( ij,l ) = convfl( ij + iim,l )
-   3  CONTINUE
-c
-c     ......  calcul aux poles  .......
-c
-      convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
-      convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
-      DO 4 ij = 1,iip1
-      convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
-      convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
-   4  CONTINUE
-c
-   5  CONTINUE
-      RETURN
-      END
+SUBROUTINE convflu( xflu,yflu,nbniv,convfl )
+
+  ! From LMDZ4/libf/dyn3d/convflu.F, version 1.1.1.1 2004/05/19 12:53:05
+
+  !  P. Le Van
+
+  ! Calcule la (convergence horiz. * aire locale) du flux ayant pour
+  ! composantes xflu et yflu ,variables extensives .
+
+  ! nbniv   est le nombre de niveaux vert. de  xflu et de yflu.
+
+  use dimens_m
+  use paramet_m
+  use comgeom
+
+  IMPLICIT NONE
+
+  integer, intent(in):: nbniv
+  REAL, intent(in):: xflu( ip1jmp1,nbniv ), yflu( ip1jm,nbniv )
+  real, intent(out):: convfl( ip1jmp1,nbniv )
+
+  real convpn,convps
+  INTEGER    l,ij
+  REAL       SSUM
+
+  !------------------------------------------------------------------
+
+  DO  l = 1,nbniv
+     DO  ij = iip2, ip1jm - 1
+        convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   + &
+             yflu(ij +1,l ) - yflu( ij -iim,l )
+     end DO
+
+     ! correction pour  convfl( 1,j,l)
+     ! convfl(1,j,l)= convfl(iip1,j,l)
+
+     DO  ij = iip2,ip1jm,iip1
+        convfl( ij,l ) = convfl( ij + iim,l )
+     end DO
+
+     ! calcul aux pôles
+
+     convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
+     convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
+     DO  ij = 1,iip1
+        convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
+        convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
+     end DO
+  end DO
+
+END SUBROUTINE convflu