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

Annotation of /trunk/dyn3d/laplacien_gam.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 55 - (hide annotations)
Mon Dec 12 13:25:01 2011 UTC (12 years, 5 months ago) by guez
Original Path: trunk/libf/dyn3d/laplacien_gam.f
File size: 1352 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/laplacien_gam.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3     !
4     SUBROUTINE laplacien_gam ( klevel, cuvsga, cvusga, unsaigam ,
5     * unsapolnga, unsapolsga, teta, divgra )
6    
7     c P. Le Van
8     c
9     c ************************************************************
10     c
11     c .... calcul de (div( grad )) de teta .....
12     c ************************************************************
13     c klevel et teta sont des arguments d'entree pour le s-prog
14     c divgra est un argument de sortie pour le s-prog
15     c
16     use dimens_m
17     use paramet_m
18     use comgeom
19     IMPLICIT NONE
20     c
21    
22     c
23     c ............ variables en arguments ..........
24     c
25 guez 55 INTEGER, intent(in):: klevel
26 guez 3 REAL teta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
27     REAL cuvsga(ip1jm) , cvusga( ip1jmp1 ),unsaigam(ip1jmp1),
28     * unsapolnga, unsapolsga
29     c
30     c ........... variables locales .................
31     c
32     REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
33     c ......................................................
34    
35     CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 )
36     c
37     CALL grad ( klevel, divgra, ghx, ghy )
38     c
39     CALL diverg_gam ( klevel, cuvsga, cvusga, unsaigam ,
40     * unsapolnga, unsapolsga, ghx , ghy , divgra )
41    
42     c
43    
44    
45     RETURN
46     END

  ViewVC Help
Powered by ViewVC 1.1.21