source: LMDZ5/trunk/libf/dyn3d_common/grad.F @ 2433

Last change on this file since 2433 was 1945, checked in by lguez, 11 years ago

Duplicated files moved to dyn3d_common.

  • 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
RevLine 
[524]1!
2! $Header$
3!
4      SUBROUTINE  grad(klevel, pg,pgx,pgy )
5c
6c      P. Le Van
7c
8c    ******************************************************************
9c     .. calcul des composantes covariantes en x et y du gradient de g
10c
11c    ******************************************************************
12c             pg        est un   argument  d'entree pour le s-prog
13c       pgx  et  pgy    sont des arguments de sortie pour le s-prog
14c
15      IMPLICIT NONE
16c
17#include "dimensions.h"
18#include "paramet.h"
19      INTEGER klevel
20      REAL  pg( ip1jmp1,klevel )
21      REAL pgx( ip1jmp1,klevel ) , pgy( ip1jm,klevel )
22      INTEGER  l,ij
23c
24c
25      DO 6 l = 1,klevel
26c
27      DO 2  ij = 1, ip1jmp1 - 1
28      pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
29   2  CONTINUE
30c
31c    .... correction pour  pgx(ip1,j,l)  ....
32c    ...    pgx(iip1,j,l)= pgx(1,j,l)  ....
33CDIR$ IVDEP
34      DO 3  ij = iip1, ip1jmp1, iip1
35      pgx( ij,l ) = pgx( ij -iim,l )
36   3  CONTINUE
37c
38      DO 4 ij = 1,ip1jm
39      pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
40   4  CONTINUE
41c
42   6  CONTINUE
43      RETURN
44      END
Note: See TracBrowser for help on using the repository browser.