5 |
implicit none |
implicit none |
6 |
|
|
7 |
LOGICAL:: cycle_diurne = .TRUE. |
LOGICAL:: cycle_diurne = .TRUE. |
8 |
! Cette option permet d'éteindre le cycle diurne. |
! Cette option permet d'\'eteindre le cycle diurne. Peut \^etre utile |
9 |
! Peut être utile pour accélérer le code. |
! pour acc\'el\'erer le code. |
10 |
|
|
11 |
LOGICAL:: soil_model = .TRUE. |
LOGICAL:: soil_model = .TRUE. |
12 |
! Choix du modele de sol (Thermique ?) |
! Choix du modele de sol (Thermique ?) |
|
! Option qui pourait un string afin de pouvoir |
|
|
! plus de choix ! Ou meme une liste d'options |
|
13 |
|
|
14 |
LOGICAL:: new_oliq = .TRUE. |
LOGICAL:: new_oliq = .TRUE. |
15 |
! Permet de mettre en route la nouvelle parametrisation de l'eau liquide |
! Permet de mettre en route la nouvelle parametrisation de l'eau liquide |
16 |
|
|
17 |
LOGICAL:: ok_orodr = .TRUE. ! pour l'orographie |
! Pour l'orographie: |
18 |
LOGICAL:: ok_orolf = .TRUE. ! pour l'orographie |
LOGICAL:: ok_orodr = .TRUE. |
19 |
|
LOGICAL:: ok_orolf = .TRUE. |
20 |
|
|
21 |
LOGICAL:: ok_limitvrai = .FALSE. |
LOGICAL:: ok_limitvrai = .FALSE. |
22 |
! On peut forcer le modele a lire le fichier SST de la bonne |
! On peut forcer le modele a lire le fichier SST de la bonne |
23 |
! annee. C'est une tres bonne idee, pourquoi ne pas mettre toujours |
! annee. |
|
! a y ? |
|
24 |
|
|
25 |
INTEGER:: nbapp_rad = 12 |
INTEGER:: nbapp_rad = 12 |
26 |
! (nombre d'appels des routines de rayonnements par jour) |
! nombre d'appels des routines de rayonnements par jour |
27 |
|
|
28 |
INTEGER:: iflag_con = 2 |
logical:: conv_emanuel = .true. ! convection scheme of Emanuel, else Tiedtke |
|
! Flag pour la convection : |
|
|
! 1 LMD, |
|
|
! 2 Tiedtke, |
|
|
! 3 CCM(NCAR) |
|
|
! 3 KE |
|
|
! 4 KE vect |
|
29 |
|
|
30 |
contains |
contains |
31 |
|
|
32 |
subroutine read_clesphys2 |
subroutine read_clesphys2 |
33 |
|
|
34 |
namelist /clesphys2_nml/cycle_diurne, soil_model, new_oliq, & |
use unit_nml_m, only: unit_nml |
35 |
ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad, iflag_con |
use nr_util, only: assert |
36 |
|
use conf_gcm_m, only: day_step, iphysiq |
37 |
|
|
38 |
|
namelist /clesphys2_nml/cycle_diurne, soil_model, new_oliq, ok_orodr, & |
39 |
|
ok_orolf, ok_limitvrai, nbapp_rad, conv_emanuel |
40 |
|
|
41 |
!------------------------------------ |
!------------------------------------ |
42 |
|
|
43 |
print *, "Enter namelist 'clesphys2_nml'." |
print *, "Enter namelist 'clesphys2_nml'." |
44 |
read(unit=*, nml=clesphys2_nml) |
read(unit=*, nml=clesphys2_nml) |
45 |
write(unit=*, nml=clesphys2_nml) |
write(unit_nml, nml=clesphys2_nml) |
46 |
|
call assert(mod(day_step / iphysiq, nbapp_rad) == 0, & |
47 |
select case (iflag_con) |
"read_clesphys2 nbapp_rad") |
48 |
case (1) |
call assert(nbapp_rad >= 4 .or. .not. cycle_diurne, & |
49 |
PRINT *, 'Schéma convection LMD' |
"read_clesphys2: minimum 4 calls to radiative transfer per day if " & |
50 |
case (2) |
// "cycle diurne") |
|
PRINT *, 'Schéma convection Tiedtke' |
|
|
case (3) |
|
|
PRINT *, 'Schéma convection CCM' |
|
|
END select |
|
51 |
|
|
52 |
end subroutine read_clesphys2 |
end subroutine read_clesphys2 |
53 |
|
|