/[lmdze]/trunk/Sources/phylmd/conf_phys.f
ViewVC logotype

Contents of /trunk/Sources/phylmd/conf_phys.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 228 - (show annotations)
Fri Nov 3 12:38:47 2017 UTC (6 years, 6 months ago) by guez
File size: 1698 byte(s)
Bug fix in dynetat0: phisinit to phis.

gcm explodes (stops in hgardfou) in less than one day with iflag_pbl =
7 (Mellor and Yamada 2.0 Fournier) and 11 (corresponding to iflag_pbl
= 31 in LMDZ, call to vdif_kcay). So remove those choices. Not much
used in LMDZ either. Remaining useful choices are iflag = 0, 1, 6, 8,
9.

Remove procedure yamada, which was not used.

1 module conf_phys_m
2
3 implicit none
4
5 integer:: iflag_pbl = 1 ! for the planetary boundary layer
6 REAL:: rad_chau1 = 13., rad_chau2 = 9.
7 real:: epmax = 0.993 ! \'efficacit\'e de pr\'ecipitation
8 integer:: iflag_clw = 0
9
10 contains
11
12 subroutine conf_phys
13
14 ! From phylmd/conf_phys.F90, version 1.7 2005/07/05 07:21:23
15
16 ! Configuration de la "physique" de LMDZ.
17
18 USE clesphys, ONLY: read_clesphys
19 use clesphys2, only: read_clesphys2
20 USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &
21 cld_tau_lsc, coef_eva, ffallv_con, ffallv_lsc, iflag_pdf, reevap_ice
22 use nr_util, only: assert
23 use unit_nml_m, only: unit_nml
24 USE yomcst, ONLY: read_YOMCST
25
26 namelist /conf_phys_nml/ epmax, iflag_clw, cld_lc_lsc, cld_lc_con, &
27 cld_tau_lsc, cld_tau_con, ffallv_lsc, ffallv_con, coef_eva, &
28 reevap_ice, iflag_pdf, iflag_pbl
29
30 namelist /nuagecom/ rad_chau1, rad_chau2
31
32 !-----------------------------------------------------------
33
34 print *, "Call sequence information: conf_phys"
35 call read_clesphys2
36 call read_YOMCST
37
38 cld_lc_lsc = 2.6e-4
39 cld_lc_con = 2.6e-4
40 cld_tau_lsc = 3600.
41 cld_tau_con = 3600.
42 ffallv_lsc = 1.
43 ffallv_con = 1.
44 coef_eva = 2.e-5
45 reevap_ice = .false.
46 iflag_pdf = 0
47
48 print *, "Enter namelist 'conf_phys_nml'."
49 read(unit=*, nml=conf_phys_nml)
50 write(unit_nml, nml=conf_phys_nml)
51
52 call assert(any(iflag_pbl == [0, 1, 6, 8, 9]), &
53 "conf_phys: bad value for iflag_pbl")
54 call read_clesphys
55
56 print *, "Enter namelist 'nuagecom'."
57 read(unit=*, nml=nuagecom)
58 write(unit_nml, nml=nuagecom)
59
60 end subroutine conf_phys
61
62 end module conf_phys_m

  ViewVC Help
Powered by ViewVC 1.1.21