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