source: LMDZ6/trunk/libf/dyn3d_common/laplacien_rotgam.F @ 5150

Last change on this file since 5150 was 4593, checked in by yann meurdesoif, 17 months ago

Replace #include (c preprocessor) by INCLUDE (fortran keyword)

in phylmd (except rrtm and ecrad) filtrez, dy3dmem and dyn3dcommon

Other directories will follow
YM
  • 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.1 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE laplacien_rotgam ( klevel, rotin, rotout )
5c
6c     P. Le Van
7c
8c   ************************************************************
9c   ... calcul de  (rotat x nxgrad)_gam  du rotationnel rotin ..
10c   ************************************************************
11c     klevel et teta  sont des arguments  d'entree pour le s-prog
12c      divgra     est  un argument  de sortie pour le s-prog
13c
14      IMPLICIT NONE
15c
16      INCLUDE "dimensions.h"
17      INCLUDE "paramet.h"
18      INCLUDE "comgeom.h"
19
20c
21c    .............   variables  en  arguments    ...........
22c
23      INTEGER klevel
24      REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel )
25c
26c   ............     variables   locales     ...............
27c
28      INTEGER l, ij
29      REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
30c   ........................................................
31c
32c
33
34      CALL   nxgrad_gam ( klevel, rotin,   ghx ,   ghy  )
35      CALL   rotat_nfil ( klevel, ghx  ,   ghy , rotout )
36c
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
Note: See TracBrowser for help on using the repository browser.