libf/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	!
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