/[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 56 - (show 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 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 (grad div)
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