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
      use filtreg_m, only: filtreg
      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.)

      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	use filtreg_m, only: filtreg
19	IMPLICIT NONE
20	c
21
22	c
23	c .......... variables en arguments .............
24	c
25	INTEGER klevel
26	REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel )
27	c
28	c .......... variables locales ................
29	c
30	REAL ghy(ip1jm,klevel), ghx(ip1jmp1,klevel)
31	c ........................................................
32	c
33	c
34	CALL filtreg ( rotin , jjm, klevel, 2, 1, .FALSE.)
35
36	CALL nxgrad ( klevel, rotin, ghx , ghy )
37	CALL rotatf ( klevel, ghx , ghy , rotout )
38	c
39	RETURN
40	END