source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/laplacien_gam.f90 @ 5122

Last change on this file since 5122 was 5119, checked in by abarral, 4 months ago

enforce PRIVATE by default in several modules, expose PUBLIC as needed
move eigen.f90 to obsolete/
(lint) aslong the way
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.4 KB
Line 
1! $Header$
2
3SUBROUTINE laplacien_gam(klevel, cuvsga, cvusga, unsaigam, &
4        unsapolnga, unsapolsga, teta, divgra)
5
6  !  P. Le Van
7  !
8  !   ************************************************************
9  !
10  !  ....   calcul de  (div( grad ))   de   teta  .....
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 lmdz_ssum_scopy, ONLY: scopy
16
17  IMPLICIT NONE
18  !
19  INCLUDE "dimensions.h"
20  INCLUDE "paramet.h"
21  INCLUDE "comgeom.h"
22
23  !
24  !    ............     variables  en arguments    ..........
25  !
26  INTEGER :: klevel
27  REAL :: teta(ip1jmp1, klevel), divgra(ip1jmp1, klevel)
28  REAL :: cuvsga(ip1jm), cvusga(ip1jmp1), unsaigam(ip1jmp1), &
29          unsapolnga, unsapolsga
30  !
31  !    ...........    variables  locales    .................
32  !
33  REAL :: ghy(ip1jm, llm), ghx(ip1jmp1, llm)
34  !    ......................................................
35
36  !
37  !
38  !   ...  cvuscugam  = ( cvu/ cu ) ** (- gamdissip )
39  !   ...  cuvscvgam  = ( cuv/ cv ) ** (- gamdissip )  calcules dans inigeom  ..
40  !   ...  unsairegam =  1. /  aire ** (- gamdissip )
41  !
42
43  CALL SCOPY (ip1jmp1 * klevel, teta, 1, divgra, 1)
44  !
45  CALL   grad (klevel, divgra, ghx, ghy)
46  !
47  CALL  diverg_gam (klevel, cuvsga, cvusga, unsaigam, &
48          unsapolnga, unsapolsga, ghx, ghy, divgra)
49
50  !
51
52  RETURN
53END SUBROUTINE laplacien_gam
Note: See TracBrowser for help on using the repository browser.