trunk/phylmd/clesphys2.f

module clesphys2

  ! From version 1.3 2005/06/06 13:16:33

  implicit none

  LOGICAL:: new_oliq = .TRUE.
  ! Permet de mettre en route la nouvelle parametrisation de l'eau liquide

  ! Pour l'orographie:
  LOGICAL:: ok_orodr = .TRUE.
  LOGICAL:: ok_orolf = .TRUE.

  LOGICAL:: ok_limitvrai = .FALSE.
  ! On peut forcer le modele a lire le fichier SST de la bonne
  ! annee.

  INTEGER:: nbapp_rad = 12
  ! nombre d'appels des routines de rayonnements par jour

  logical:: conv_emanuel = .true. ! convection scheme of Emanuel, else Tiedtke

contains

  subroutine read_clesphys2

    use unit_nml_m, only: unit_nml
    use nr_util, only: assert
    use conf_gcm_m, only: day_step, iphysiq

    namelist /clesphys2_nml/new_oliq, ok_orodr, ok_orolf, ok_limitvrai, &
         nbapp_rad, conv_emanuel

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

    print *, "Enter namelist 'clesphys2_nml'."
    read(unit=*, nml=clesphys2_nml)
    write(unit_nml, nml=clesphys2_nml)
    call assert(mod(day_step / iphysiq, nbapp_rad) == 0, &
         "read_clesphys2 nbapp_rad")
    call assert(nbapp_rad >= 4, &
         "read_clesphys2: minimum 4 calls to radiative transfer per day")

  end subroutine read_clesphys2

end module clesphys2
1	module clesphys2
2
3	! From version 1.3 2005/06/06 13:16:33
4
5	implicit none
6
7	LOGICAL:: new_oliq = .TRUE.
8	! Permet de mettre en route la nouvelle parametrisation de l'eau liquide
9
10	! Pour l'orographie:
11	LOGICAL:: ok_orodr = .TRUE.
12	LOGICAL:: ok_orolf = .TRUE.
13
14	LOGICAL:: ok_limitvrai = .FALSE.
15	! On peut forcer le modele a lire le fichier SST de la bonne
16	! annee.
17
18	INTEGER:: nbapp_rad = 12
19	! nombre d'appels des routines de rayonnements par jour
20
21	logical:: conv_emanuel = .true. ! convection scheme of Emanuel, else Tiedtke
22
23	contains
24
25	subroutine read_clesphys2
26
27	use unit_nml_m, only: unit_nml
28	use nr_util, only: assert
29	use conf_gcm_m, only: day_step, iphysiq
30
31	namelist /clesphys2_nml/new_oliq, ok_orodr, ok_orolf, ok_limitvrai, &
32	nbapp_rad, conv_emanuel
33
34	!------------------------------------
35
36	print *, "Enter namelist 'clesphys2_nml'."
37	read(unit=*, nml=clesphys2_nml)
38	write(unit_nml, nml=clesphys2_nml)
39	call assert(mod(day_step / iphysiq, nbapp_rad) == 0, &
40	"read_clesphys2 nbapp_rad")
41	call assert(nbapp_rad >= 4, &
42	"read_clesphys2: minimum 4 calls to radiative transfer per day")
43
44	end subroutine read_clesphys2
45
46	end module clesphys2