trunk/dyn3d/covnat.f

!
! $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
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covnat.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
3			!
4			SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )
5			use dimens_m
6			use paramet_m
7			use comgeom
8			IMPLICIT NONE
9
10			c=======================================================================
11			c
12			c Auteur: F Hourdin Phu LeVan
13			c -------
14			c
15			c Objet:
16			c ------
17			c
18			c *********************************************************************
19			c calcul des compos. naturelles a partir des comp.covariantes
20			c ********************************************************************
21			c
22			c=======================================================================
23
24
25			INTEGER klevel
26			REAL ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel )
27			REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel )
28			INTEGER l,ij
29
30
31			DO l = 1,klevel
32			DO ij = 1, iip1
33			unat (ij,l) =0.
34			END DO
35
36			DO ij = iip2, ip1jm
37			unat( ij,l ) = ucov( ij,l ) / cu(ij)
38			ENDDO
39			DO ij = ip1jm+1, ip1jmp1
40			unat (ij,l) =0.
41			END DO
42
43			DO ij = 1,ip1jm
44			vnat( ij,l ) = vcov( ij,l ) / cv(ij)
45			ENDDO
46
47			ENDDO
48			RETURN
49			END