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)

  ! P. Le Van

  ! ************************************************************
  ! ...  calcul de  ( rotat x nxgrad )  du rotationnel rotin  .
  ! ************************************************************

  ! klevel et rotin  sont des arguments  d'entree pour le s-prog
  ! rotout           est  un argument  de sortie pour le s-prog

  USE dimens_m
  USE paramet_m
  USE comgeom
  USE filtreg_m, ONLY: filtreg
  IMPLICIT NONE


  ! ..........    variables  en  arguments     .............

  INTEGER, INTENT (IN) :: klevel
  REAL rotin(iim + 1, jjm, klevel), rotout(ip1jm, klevel)

  ! ..........    variables   locales       ................

  REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
  ! ........................................................


  CALL filtreg(rotin, jjm, klevel, 2, 1, .FALSE.)

  CALL nxgrad(klevel, rotin, ghx, ghy)
  CALL rotatf(klevel, ghx, ghy, rotout)

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