trunk/dyn3d/laplacien_rotgam.f


! $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	guez	3
2	guez	81	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rotgam.F,v 1.1.1.1
3			! 2004/05/19 12:53:05 lmdzadmin Exp $
4	guez	3
5	guez	81	SUBROUTINE laplacien_rotgam(klevel, rotin, rotout)
6	guez	3
7	guez	81	! 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