trunk/dyn3d/laplacien_rot.f90


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