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

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

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
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 guez 3 module comdissnew
2    
3     implicit none
4    
5 guez 54 LOGICAL:: lstardis = .TRUE. ! choix de l'opérateur de dissipation star
6 guez 3
7 guez 54 INTEGER:: nitergdiv = 1
8     ! nombre d'itérations de l'opérateur de dissipation gradiv
9 guez 3
10 guez 54 INTEGER:: nitergrot = 2
11     ! nombre d'itérations de l'opérateur de dissipation nxgradrot
12 guez 3
13 guez 54 INTEGER:: niterh = 2
14     ! nombre d'itérations de l'opérateur de dissipation divgrad
15 guez 3
16 guez 54 REAL:: tetagdiv = 7200. ! in s
17     ! temps de dissipation des plus petites longueurs d'ondes pour u, v (gradiv)
18 guez 3
19 guez 54 REAL:: tetagrot = 7200. ! in s
20     ! temps de dissipation des plus petites longueurs d'ondes pour u, v
21 guez 27 ! (nxgradrot)
22 guez 3
23 guez 54 REAL:: tetatemp = 7200. ! in s
24 guez 27 ! temps de dissipation des plus petites longueurs d'ondes pour h (divgrad)
25 guez 3
26 guez 54 REAL:: coefdis = 0.
27 guez 3
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