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 )
c
c     P. Le Van
c
c   ************************************************************
c   ... calcul de  (rotat x nxgrad)_gam  du rotationnel rotin ..
c   ************************************************************
c     klevel et teta  sont des arguments  d'entree pour le s-prog
c      divgra     est  un argument  de sortie pour le s-prog
c
      use dimens_m
      use paramet_m
      use comgeom
      IMPLICIT NONE
c

c
c    .............   variables  en  arguments    ...........
c
      INTEGER, intent(in):: klevel
      REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel )
c
c   ............     variables   locales     ...............
c
      INTEGER l, ij
      REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
c   ........................................................
c
c

      CALL   nxgrad_gam ( klevel, rotin,   ghx ,   ghy  )
      CALL   rotat_nfil ( klevel, ghx  ,   ghy , rotout )
c
      DO l = 1, klevel
        DO ij = 1, ip1jm
         rotout(ij,l) = rotout(ij,l) * unsairz_gam(ij)
        ENDDO
      ENDDO

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