trunk/dyn3d/laplacien_rot.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rot.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
!
      SUBROUTINE laplacien_rot ( klevel, rotin, rotout,ghx,ghy )
c
c    P. Le Van
c
c   ************************************************************
c    ...  calcul de  ( rotat x nxgrad )  du rotationnel rotin  .
c   ************************************************************
c
c     klevel et rotin  sont des arguments  d'entree pour le s-prog
c      rotout           est  un argument  de sortie pour le s-prog
c
      use dimens_m
      use paramet_m
      use comgeom
      IMPLICIT NONE
c

c 
c   ..........    variables  en  arguments     .............
c
      INTEGER klevel
      REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel )
c
c   ..........    variables   locales       ................
c
      REAL ghy(ip1jm,klevel), ghx(ip1jmp1,klevel)
c   ........................................................
c
c
      CALL  filtreg ( rotin ,   jjm, klevel,   2, 1, .FALSE., 1 )

      CALL   nxgrad ( klevel, rotin,   ghx ,  ghy               )
      CALL   rotatf  ( klevel, ghx  ,   ghy , rotout             )
c
      RETURN
      END
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rot.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
3			!
4			SUBROUTINE laplacien_rot ( klevel, rotin, rotout,ghx,ghy )
5			c
6			c P. Le Van
7			c
8			c ************************************************************
9			c ... calcul de ( rotat x nxgrad ) du rotationnel rotin .
10			c ************************************************************
11			c
12			c klevel et rotin sont des arguments d'entree pour le s-prog
13			c rotout est un argument de sortie pour le s-prog
14			c
15			use dimens_m
16			use paramet_m
17			use comgeom
18			IMPLICIT NONE
19			c
20
21			c
22			c .......... variables en arguments .............
23			c
24			INTEGER klevel
25			REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel )
26			c
27			c .......... variables locales ................
28			c
29			REAL ghy(ip1jm,klevel), ghx(ip1jmp1,klevel)
30			c ........................................................
31			c
32			c
33			CALL filtreg ( rotin , jjm, klevel, 2, 1, .FALSE., 1 )
34
35			CALL nxgrad ( klevel, rotin, ghx , ghy )
36			CALL rotatf ( klevel, ghx , ghy , rotout )
37			c
38			RETURN
39			END