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	guez	3
2	guez	81	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rot.F,v 1.1.1.1
3			! 2004/05/19 12:53:07 lmdzadmin Exp $
4	guez	3
5	guez	81	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	guez	106	REAL rotin(iim + 1, jjm, klevel), rotout(ip1jm, klevel)
28	guez	81
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