trunk/dyn3d/covcont.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covcont.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
!
      SUBROUTINE covcont (klevel,ucov, vcov, ucont, vcont )
      use dimens_m
      use paramet_m
      use comgeom
      IMPLICIT NONE

c=======================================================================
c
c   Auteur:  P. Le Van
c   -------
c
c   Objet:
c   ------
c
c  *********************************************************************
c    calcul des compos. contravariantes a partir des comp.covariantes
c  ********************************************************************
c
c=======================================================================


      INTEGER klevel
      REAL, intent(in):: ucov( ip1jmp1,klevel ),  vcov( ip1jm,klevel )
      REAL ucont( ip1jmp1,klevel ), vcont( ip1jm,klevel )
      INTEGER   l,ij


      DO 10 l = 1,klevel

      DO 2  ij = iip2, ip1jm
      ucont( ij,l ) = ucov( ij,l ) * unscu2( ij )
   2  CONTINUE

      DO 4 ij = 1,ip1jm
      vcont( ij,l ) = vcov( ij,l ) * unscv2( ij )
   4  CONTINUE

  10  CONTINUE
      RETURN
      END
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/covcont.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
3			!
4			SUBROUTINE covcont (klevel,ucov, vcov, ucont, vcont )
5			use dimens_m
6			use paramet_m
7			use comgeom
8			IMPLICIT NONE
9
10			c=======================================================================
11			c
12			c Auteur: P. Le Van
13			c -------
14			c
15			c Objet:
16			c ------
17			c
18			c *********************************************************************
19			c calcul des compos. contravariantes a partir des comp.covariantes
20			c ********************************************************************
21			c
22			c=======================================================================
23
24
25			INTEGER klevel
26			REAL, intent(in):: ucov( ip1jmp1,klevel ), vcov( ip1jm,klevel )
27			REAL ucont( ip1jmp1,klevel ), vcont( ip1jm,klevel )
28			INTEGER l,ij
29
30
31			DO 10 l = 1,klevel
32
33			DO 2 ij = iip2, ip1jm
34			ucont( ij,l ) = ucov( ij,l ) * unscu2( ij )
35			2 CONTINUE
36
37			DO 4 ij = 1,ip1jm
38			vcont( ij,l ) = vcov( ij,l ) * unscv2( ij )
39			4 CONTINUE
40
41			10 CONTINUE
42			RETURN
43			END