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