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