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
  use rotatf_m, only: rotatf

  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, direct = .true., intensive = .true.)

  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	use rotatf_m, only: rotatf
21
22	IMPLICIT NONE
23
24
25
26	! .......... variables en arguments .............
27
28	INTEGER, INTENT (IN) :: klevel
29	REAL rotin(iim + 1, jjm, klevel), rotout(ip1jm, klevel)
30
31	! .......... variables locales ................
32
33	REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
34	! ........................................................
35
36
37	CALL filtreg(rotin, direct = .true., intensive = .true.)
38
39	CALL nxgrad(klevel, rotin, ghx, ghy)
40	CALL rotatf(klevel, ghx, ghy, rotout)
41
42	RETURN
43	END SUBROUTINE laplacien_rot