--- trunk/libf/dyn3d/covnat.f	2008/02/27 13:16:39	3
+++ trunk/dyn3d/covnat.f90	2013/11/15 18:45:49	76
@@ -1,49 +1,48 @@
-!
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covnat.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
-!
-      SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )
-      use dimens_m
-      use paramet_m
-      use comgeom
-      IMPLICIT NONE
-
-c=======================================================================
-c
-c   Auteur:  F Hourdin Phu LeVan
-c   -------
-c
-c   Objet:
-c   ------
-c
-c  *********************************************************************
-c    calcul des compos. naturelles a partir des comp.covariantes
-c  ********************************************************************
-c
-c=======================================================================
-
-
-      INTEGER klevel
-      REAL ucov( ip1jmp1,klevel ),  vcov( ip1jm,klevel )
-      REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel )
-      INTEGER   l,ij
-
-
-      DO l = 1,klevel
-         DO ij = 1, iip1
-            unat (ij,l) =0.
-         END DO
-
-         DO ij = iip2, ip1jm
-            unat( ij,l ) = ucov( ij,l ) / cu(ij)
-         ENDDO
-         DO ij = ip1jm+1, ip1jmp1  
-            unat (ij,l) =0.
-         END DO
-
-         DO ij = 1,ip1jm
-            vnat( ij,l ) = vcov( ij,l ) / cv(ij)
-         ENDDO
-
-      ENDDO
-      RETURN
-      END
+module covnat_m
+
+  IMPLICIT NONE
+
+contains
+
+  SUBROUTINE covnat(klevel, ucov, vcov, unat, vnat)
+
+    ! From LMDZ4/libf/dyn3d/covnat.F, version 1.1.1.1 2004/05/19 12:53:07
+
+    USE comgeom, ONLY: cu, cv
+    USE paramet_m, ONLY: iip1, iip2, ip1jm, ip1jmp1
+
+    ! Authors: F. Hourdin, Phu Le Van
+
+    ! Objet : calcul des composantes naturelles à partir des
+    ! composantes covariantes.
+
+    INTEGER, intent(in):: klevel
+    REAL, intent(in):: ucov(ip1jmp1, klevel), vcov(ip1jm, klevel)
+    REAL, intent(out):: unat(ip1jmp1, klevel), vnat(ip1jm, klevel)
+
+    ! Local:
+    INTEGER l, ij
+
+    !------------------------------------------------------------------
+
+    DO l = 1, klevel
+       DO ij = 1, iip1
+          unat(ij, l) =0.
+       END DO
+
+       DO ij = iip2, ip1jm
+          unat(ij, l) = ucov(ij, l) / cu(ij)
+       ENDDO
+
+       DO ij = ip1jm+1, ip1jmp1 
+          unat(ij, l) =0.
+       END DO
+
+       DO ij = 1, ip1jm
+          vnat(ij, l) = vcov(ij, l) / cv(ij)
+       ENDDO
+    ENDDO
+
+  END SUBROUTINE covnat
+
+end module covnat_m