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

Annotation of /trunk/libf/dyn3d/grad.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 55 - (hide annotations)
Mon Dec 12 13:25:01 2011 UTC (12 years, 6 months ago) by guez
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 guez 3 !
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 guez 55 INTEGER, intent(in):: klevel
20 guez 3 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