/[lmdze]/trunk/libf/dyn3d/Dissipation/comdissnew.f90
ViewVC logotype

Contents of /trunk/libf/dyn3d/Dissipation/comdissnew.f90

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
File size: 1145 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 module comdissnew
2
3 implicit none
4
5 LOGICAL:: lstardis = .TRUE. ! choix de l'opérateur de dissipation star
6
7 INTEGER:: nitergdiv = 1
8 ! nombre d'itérations de l'opérateur de dissipation gradiv
9
10 INTEGER:: nitergrot = 2
11 ! nombre d'itérations de l'opérateur de dissipation nxgradrot
12
13 INTEGER:: niterh = 2
14 ! nombre d'itérations de l'opérateur de dissipation divgrad
15
16 REAL:: tetagdiv = 7200. ! in s
17 ! temps de dissipation des plus petites longueurs d'ondes pour u, v (gradiv)
18
19 REAL:: tetagrot = 7200. ! in s
20 ! temps de dissipation des plus petites longueurs d'ondes pour u, v
21 ! (nxgradrot)
22
23 REAL:: tetatemp = 7200. ! in s
24 ! temps de dissipation des plus petites longueurs d'ondes pour h (divgrad)
25
26 REAL:: coefdis = 0.
27
28 contains
29
30 subroutine read_comdissnew
31
32 namelist /comdissnew_nml/lstardis, nitergdiv, nitergrot, niterh, &
33 tetagdiv, tetagrot, tetatemp, coefdis
34
35 !-------------------------------------------------
36
37 print *, "Enter namelist 'comdissnew_nml'."
38 read(unit=*, nml=comdissnew_nml)
39 write(unit=*, nml=comdissnew_nml)
40
41 end subroutine read_comdissnew
42
43 end module comdissnew
  ViewVC Help
Powered by ViewVC 1.1.21