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

Contents of /trunk/dyn3d/grad.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 55 - (show annotations)
Mon Dec 12 13:25:01 2011 UTC (12 years, 5 months ago) by guez
Original Path: trunk/libf/dyn3d/grad.f
File size: 1170 byte(s) 
-- In procedure "bilan_dyn", replaced average of "zvq" by integral of
"zvq", following a comment of Francis Codron :

Le calcul actuel donne des unités peu pratiques : transports de
chaleur en K m / s par exemple. C'est bien pour les sorties à 2
dimensions, latitude et pression, car alors le transport ne dépend pas
de l'espacement des niveaux, mieux pour comparer ou tracer en latitude
et pression. Par contre, quand on somme sur la verticale, on
préfèrerait avoir des transports d'énergie en watts, ou au moins an K
kg / s (à multiplier par "Cp" ou "L"). On doit pouvoir recalculer le
transport intégré à partir des fichiers de sortie, mais c'est embêtant
(calcul de "cv").

-- Gathered files in directory Dissipation.
1 !
2 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/grad.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3 !
4 SUBROUTINE grad(klevel, pg,pgx,pgy )
5 c
6 c P. Le Van
7 c
8 c ******************************************************************
9 c .. calcul des composantes covariantes en x et y du gradient de g
10 c
11 c ******************************************************************
12 c pg est un argument d'entree pour le s-prog
13 c pgx et pgy sont des arguments de sortie pour le s-prog
14 c
15 use dimens_m
16 use paramet_m
17 IMPLICIT NONE
18 c
19 INTEGER, intent(in):: klevel
20 REAL pg( ip1jmp1,klevel )
21 REAL pgx( ip1jmp1,klevel ) , pgy( ip1jm,klevel )
22 INTEGER l,ij
23 c
24 c
25 DO 6 l = 1,klevel
26 c
27 DO 2 ij = 1, ip1jmp1 - 1
28 pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
29 2 CONTINUE
30 c
31 c .... correction pour pgx(ip1,j,l) ....
32 c ... pgx(iip1,j,l)= pgx(1,j,l) ....
33 CDIR$ IVDEP
34 DO 3 ij = iip1, ip1jmp1, iip1
35 pgx( ij,l ) = pgx( ij -iim,l )
36 3 CONTINUE
37 c
38 DO 4 ij = 1,ip1jm
39 pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
40 4 CONTINUE
41 c
42 6 CONTINUE
43 RETURN
44 END
  ViewVC Help
Powered by ViewVC 1.1.21