trunk/dyn3d/laplacien_gam.f


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

SUBROUTINE laplacien_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, &
    unsapolsga, teta, divgra)

  ! P. Le Van

  ! ************************************************************

  ! ....   calcul de  (div( grad ))   de   teta  .....
  ! ************************************************************
  ! klevel et teta  sont des arguments  d'entree pour le s-prog
  ! divgra     est  un argument  de sortie pour le s-prog

  USE grad_m, ONLY: grad
  USE dimens_m
  USE paramet_m
  USE comgeom
  IMPLICIT NONE


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

  INTEGER, INTENT (IN) :: klevel
  REAL teta(ip1jmp1, klevel), divgra(ip1jmp1, klevel)
  REAL cuvsga(ip1jm), cvusga(ip1jmp1), unsaigam(ip1jmp1), unsapolnga, &
    unsapolsga

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

  REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
  ! ......................................................

  CALL scopy(ip1jmp1*klevel, teta, 1, divgra, 1)

  CALL grad(klevel, divgra, ghx, ghy)

  CALL diverg_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, unsapolsga, &
    ghx, ghy, divgra)


  RETURN
END SUBROUTINE laplacien_gam
1
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_gam.F,v 1.1.1.1
3	! 2004/05/19 12:53:06 lmdzadmin Exp $
4
5	SUBROUTINE laplacien_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, &
6	unsapolsga, teta, divgra)
7
8	! P. Le Van
9
10	! ************************************************************
11
12	! .... calcul de (div( grad )) de teta .....
13	! ************************************************************
14	! klevel et teta sont des arguments d'entree pour le s-prog
15	! divgra est un argument de sortie pour le s-prog
16
17	USE grad_m, ONLY: grad
18	USE dimens_m
19	USE paramet_m
20	USE comgeom
21	IMPLICIT NONE
22
23
24
25	! ............ variables en arguments ..........
26
27	INTEGER, INTENT (IN) :: klevel
28	REAL teta(ip1jmp1, klevel), divgra(ip1jmp1, klevel)
29	REAL cuvsga(ip1jm), cvusga(ip1jmp1), unsaigam(ip1jmp1), unsapolnga, &
30	unsapolsga
31
32	! ........... variables locales .................
33
34	REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel)
35	! ......................................................
36
37	CALL scopy(ip1jmp1*klevel, teta, 1, divgra, 1)
38
39	CALL grad(klevel, divgra, ghx, ghy)
40
41	CALL diverg_gam(klevel, cuvsga, cvusga, unsaigam, unsapolnga, unsapolsga, &
42	ghx, ghy, divgra)
43
44
45
46
47	RETURN
48	END SUBROUTINE laplacien_gam