/[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 56 - (hide annotations)
Tue Jan 10 19:02:02 2012 UTC (12 years, 4 months ago) by guez
File size: 1147 byte(s) 
Imported "writehist.f" from LMDZ.

Moved module variable "histaveid" from "com_io_dyn" to "initdynav_m".

In "inithist", access directly module variables from "com_io_dyn"
instead of going through the arguments. Copying from LMDZ, write "u"
and scalar variables to separate files. Create a new variable for the
new file in "com_io_dyn". Copying from LMDZ, change the vertical axes
of the three files.

Removed some useless initializations in "dissip".

In "bilan_dyn", removed useless variable "time". Avoiding the
approximate test on "dt_cum" being a multiple of "dt_app", just
compute "ncum" from known usage of "bilan_dyn" and compute "dt_cum"
from "ncum". Change "periodav" from real to integer in
"conf_gcm_m". Since "day_step" is required to be a multiple of
"iperiod", so is "ncum".
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 guez 56 ! temps de dissipation des plus petites longueurs d'ondes pour u, v (grad div)
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