/[lmdze]/trunk/dyn3d/grad.f
ViewVC logotype

Diff of /trunk/dyn3d/grad.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/libf/dyn3d/grad.f revision 55 by guez, Mon Dec 12 13:25:01 2011 UTC trunk/dyn3d/grad.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC
# Line 1  Line 1 
1  !  module grad_m
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/grad.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $  
3  !    IMPLICIT NONE
4        SUBROUTINE  grad(klevel, pg,pgx,pgy )  
5  c  contains
6  c      P. Le Van  
7  c    SUBROUTINE grad(klevel, g, gx, gy)
8  c    ******************************************************************  
9  c     .. calcul des composantes covariantes en x et y du gradient de g      ! From LMDZ4/libf/dyn3d/grad.F, version 1.1.1.1 2004/05/19 12:53:05
10  c      ! P. Le Van
11  c    ******************************************************************  
12  c             pg        est un   argument  d'entree pour le s-prog      ! Calcul des composantes covariantes en x et y du gradient de g.
13  c       pgx  et  pgy    sont des arguments de sortie pour le s-prog  
14  c      USE dimens_m, ONLY: iim
15        use dimens_m      USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1
16        use paramet_m  
17        IMPLICIT NONE      INTEGER, intent(in):: klevel
18  c      REAL, intent(in):: g(ip1jmp1, klevel)
19        INTEGER, intent(in):: klevel      REAL, intent(out):: gx(ip1jmp1, klevel) , gy(ip1jm, klevel)
20        REAL  pg( ip1jmp1,klevel )  
21        REAL pgx( ip1jmp1,klevel ) , pgy( ip1jm,klevel )      ! Local:
22        INTEGER  l,ij      INTEGER l, ij
23  c  
24  c      !----------------------------------------------------------------
25        DO 6 l = 1,klevel  
26  c      DO l = 1, klevel
27        DO 2  ij = 1, ip1jmp1 - 1         DO ij = 1, ip1jmp1 - 1
28        pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )            gx(ij, l) = g(ij +1, l) - g(ij, l)
29     2  CONTINUE         end DO
30  c  
31  c    .... correction pour  pgx(ip1,j,l)  ....         ! correction pour gx(ip1, j, l)
32  c    ...    pgx(iip1,j,l)= pgx(1,j,l)  ....         ! gx(iip1, j, l)= gx(1, j, l)
33  CDIR$ IVDEP         DO ij = iip1, ip1jmp1, iip1
34        DO 3  ij = iip1, ip1jmp1, iip1            gx(ij, l) = gx(ij -iim, l)
35        pgx( ij,l ) = pgx( ij -iim,l )         end DO
36     3  CONTINUE  
37  c         DO ij = 1, ip1jm
38        DO 4 ij = 1,ip1jm            gy(ij, l) = g(ij, l) - g(ij +iip1, l)
39        pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )         end DO
40     4  CONTINUE      end DO
41  c  
42     6  CONTINUE    END SUBROUTINE grad
43        RETURN  
44        END  end module grad_m
Legend:
Removed from v.55  
changed lines
  Added in v.76
  ViewVC Help
Powered by ViewVC 1.1.21