dyn3d/Dissipation/comdissnew.f

module comdissnew

  implicit none

  INTEGER:: nitergdiv = 1
  ! nombre d'itérations de l'opérateur de dissipation gradiv

  INTEGER:: nitergrot = 2
  ! nombre d'itérations de l'opérateur de dissipation nxgradrot

  INTEGER:: niterh = 2
  ! nombre d'itérations de l'opérateur de dissipation divgrad

  REAL:: tetagdiv = 7200. ! in s
  ! temps de dissipation des plus petites longueurs d'ondes pour u, v (grad div)

  REAL:: tetagrot = 7200. ! in s
  ! temps de dissipation des plus petites longueurs d'ondes pour u, v
  ! (nxgradrot)

  REAL:: tetatemp = 7200. ! in s
  ! temps de dissipation des plus petites longueurs d'ondes pour h (divgrad) 

  REAL:: coefdis = 0.

contains

  subroutine read_comdissnew

    use unit_nml_m, only: unit_nml
    
    namelist /comdissnew_nml/nitergdiv, nitergrot, niterh, tetagdiv, &
         tetagrot, tetatemp, coefdis

    !-------------------------------------------------

    print *, "Enter namelist 'comdissnew_nml'."
    read(unit=*, nml=comdissnew_nml)
    write(unit_nml, nml=comdissnew_nml)

  end subroutine read_comdissnew

end module comdissnew
1	guez	3	module comdissnew
2
3			implicit none
4
5	guez	54	INTEGER:: nitergdiv = 1
6			! nombre d'itérations de l'opérateur de dissipation gradiv
7	guez	3
8	guez	54	INTEGER:: nitergrot = 2
9			! nombre d'itérations de l'opérateur de dissipation nxgradrot
10	guez	3
11	guez	54	INTEGER:: niterh = 2
12			! nombre d'itérations de l'opérateur de dissipation divgrad
13	guez	3
14	guez	54	REAL:: tetagdiv = 7200. ! in s
15	guez	56	! temps de dissipation des plus petites longueurs d'ondes pour u, v (grad div)
16	guez	3
17	guez	54	REAL:: tetagrot = 7200. ! in s
18			! temps de dissipation des plus petites longueurs d'ondes pour u, v
19	guez	27	! (nxgradrot)
20	guez	3
21	guez	54	REAL:: tetatemp = 7200. ! in s
22	guez	27	! temps de dissipation des plus petites longueurs d'ondes pour h (divgrad)
23	guez	3
24	guez	54	REAL:: coefdis = 0.
25	guez	3
26			contains
27
28			subroutine read_comdissnew
29
30	guez	57	use unit_nml_m, only: unit_nml
31
32	guez	64	namelist /comdissnew_nml/nitergdiv, nitergrot, niterh, tetagdiv, &
33			tetagrot, tetatemp, coefdis
34	guez	3
35			!-------------------------------------------------
36
37			print *, "Enter namelist 'comdissnew_nml'."
38			read(unit=*, nml=comdissnew_nml)
39	guez	57	write(unit_nml, nml=comdissnew_nml)
40	guez	3
41			end subroutine read_comdissnew
42
43			end module comdissnew