--- trunk/dyn3d/dudv1.f	2014/03/05 14:57:53	82
+++ trunk/dyn3d/dudv1.f	2014/03/11 15:09:02	88
@@ -1,52 +1,52 @@
+module dudv1_m
 
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/dudv1.F,v 1.1.1.1 2004/05/19
-! 12:53:06 lmdzadmin Exp $
-
-SUBROUTINE dudv1(vorpot, pbaru, pbarv, du, dv)
-  USE dimens_m
-  USE paramet_m
   IMPLICIT NONE
 
-  ! -----------------------------------------------------------------------
+contains
 
-  ! Auteur:   P. Le Van
-  ! -------
+  SUBROUTINE dudv1(vorpot, pbaru, pbarv, du, dv)
 
-  ! Objet:
-  ! ------
-  ! calcul du terme de  rotation
-  ! ce terme est ajoute a  d(ucov)/dt et a d(vcov)/dt  ..
-  ! vorpot, pbaru et pbarv sont des arguments d'entree  pour le s-pg ..
-  ! du  et dv              sont des arguments de sortie pour le s-pg ..
+    ! From LMDZ4/libf/dyn3d/dudv1.F, version 1.1.1.1, 2004/05/19 12:53:06 
 
-  ! -----------------------------------------------------------------------
+    ! Author: P. Le Van
 
+    ! Objet: calcul du terme de rotation. Ce terme est ajouté à
+    ! d(ucov)/dt et à d(vcov)/dt.
 
-  REAL vorpot(ip1jm, llm), pbaru(ip1jmp1, llm), pbarv(ip1jm, llm), &
-    du(ip1jmp1, llm), dv(ip1jm, llm)
-  INTEGER l, ij
+    USE dimens_m, ONLY: iim, llm
+    USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1
 
+    REAL, intent(in):: vorpot(ip1jm, llm)
+    REAL, intent(in):: pbaru(ip1jmp1, llm), pbarv(ip1jm, llm)
 
-  DO l = 1, llm
+    real du(ip1jmp1, llm), dv(ip1jm, llm)
+    ! du et dv sont des arguments de sortie pour le s-pg 
 
-    DO ij = iip2, ip1jm - 1
-      du(ij, l) = 0.125*(vorpot(ij-iip1,l)+vorpot(ij,l))* &
-        (pbarv(ij-iip1,l)+pbarv(ij-iim,l)+pbarv(ij,l)+pbarv(ij+1,l))
-    END DO
+    ! Local:
+    INTEGER l, ij
 
-    DO ij = 1, ip1jm - 1
-      dv(ij+1, l) = -0.125*(vorpot(ij,l)+vorpot(ij+1,l))*(pbaru(ij,l)+pbaru( &
-        ij+1,l)+pbaru(ij+iip1,l)+pbaru(ij+iip2,l))
-    END DO
+    !----------------------------------------------------------------------
 
-    ! .... correction  pour  dv( 1,j,l )  .....
-    ! ....   dv(1,j,l)= dv(iip1,j,l) ....
+    DO l = 1, llm
+       DO ij = iip2, ip1jm - 1
+          du(ij, l) = 0.125 * (vorpot(ij - iip1, l) + vorpot(ij, l)) &
+               * (pbarv(ij - iip1, l) + pbarv(ij - iim, l) + pbarv(ij, l) &
+               + pbarv(ij + 1, l))
+       END DO
 
-    ! DIR$ IVDEP
-    DO ij = 1, ip1jm, iip1
-      dv(ij, l) = dv(ij+iim, l)
+       DO ij = 1, ip1jm - 1
+          dv(ij + 1, l) = - 0.125 * (vorpot(ij, l) + vorpot(ij + 1, l)) &
+               * (pbaru(ij, l) + pbaru(ij + 1, l) + pbaru(ij + iip1, l) &
+               + pbaru(ij + iip2, l))
+       END DO
+
+       ! correction pour dv(1, j, l) 
+       ! dv(1, j, l) = dv(iip1, j, l) 
+       DO ij = 1, ip1jm, iip1
+          dv(ij, l) = dv(ij + iim, l)
+       END DO
     END DO
 
-  END DO
-  RETURN
-END SUBROUTINE dudv1
+  END SUBROUTINE dudv1
+
+end module dudv1_m