Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/grad.f90
- Timestamp:
- Aug 2, 2024, 9:58:25 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/grad.f90
r5105 r5159 3 3 4 4 SUBROUTINE grad(klevel, pg,pgx,pgy ) 5 ! 5 6 6 ! P. Le Van 7 ! 7 8 8 ! ****************************************************************** 9 9 ! .. calcul des composantes covariantes en x et y du gradient de g 10 ! 10 11 11 ! ****************************************************************** 12 12 ! pg est un argument d'entree pour le s-prog 13 13 ! pgx et pgy sont des arguments de sortie pour le s-prog 14 ! 14 15 USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm 16 USE lmdz_paramet 15 17 IMPLICIT NONE 16 18 ! 17 INCLUDE "dimensions.h" 18 INCLUDE "paramet.h" 19 20 19 21 INTEGER :: klevel 20 22 REAL :: pg( ip1jmp1,klevel ) 21 23 REAL :: pgx( ip1jmp1,klevel ) , pgy( ip1jm,klevel ) 22 24 INTEGER :: l,ij 23 ! 24 ! 25 26 25 27 DO l = 1,klevel 26 ! 28 27 29 DO ij = 1, ip1jmp1 - 1 28 30 pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l ) 29 31 END DO 30 ! 32 31 33 ! .... correction pour pgx(ip1,j,l) .... 32 34 ! ... pgx(iip1,j,l)= pgx(1,j,l) .... … … 35 37 pgx( ij,l ) = pgx( ij -iim,l ) 36 38 END DO 37 ! 39 38 40 DO ij = 1,ip1jm 39 41 pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l ) 40 42 END DO 41 ! 43 42 44 END DO 43 45 RETURN
Note: See TracChangeset
for help on using the changeset viewer.