trunk/dyn3d/laplacien_rotgam.f90


! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rotgam.F,v 1.1.1.1
! 2004/05/19 12:53:05 lmdzadmin Exp $

SUBROUTINE laplacien_rotgam(klevel, rotin, rotout)

  ! P. Le Van

  ! ************************************************************
  ! ... calcul de  (rotat x nxgrad)_gam  du rotationnel rotin ..
  ! ************************************************************
  ! klevel et teta  sont des arguments  d'entree pour le s-prog
  ! divgra     est  un argument  de sortie pour le s-prog

  USE dimens_m
  USE paramet_m
  USE comgeom
  IMPLICIT NONE


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

  INTEGER, INTENT (IN) :: klevel
  REAL rotin(ip1jm, klevel), rotout(ip1jm, klevel)

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

  INTEGER l, ij
  REAL ghy(ip1jm, llm), ghx(ip1jmp1, llm)
  ! ........................................................


  CALL nxgrad_gam(klevel, rotin, ghx, ghy)
  CALL rotat_nfil(klevel, ghx, ghy, rotout)

  DO l = 1, klevel
    DO ij = 1, ip1jm
      rotout(ij, l) = rotout(ij, l)*unsairz_gam(ij)
    END DO
  END DO

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