--- trunk/Sources/phylmd/conf_phys.f	2017/11/02 15:47:03	227
+++ trunk/phylmd/conf_phys.f	2018/02/05 10:39:38	254
@@ -2,7 +2,10 @@
 
   implicit none
 
-  integer:: iflag_pbl = 1 ! for the planetary boundary layer
+  integer, protected:: iflag_pbl = 1 ! for the planetary boundary layer
+  ! 6 : Mellor and Yamada 2.0
+  ! 8 : Mellor and Yamada 2.5
+
   REAL:: rad_chau1 = 13., rad_chau2 = 9.
   real:: epmax = 0.993 ! \'efficacit\'e de pr\'ecipitation
   integer:: iflag_clw = 0
@@ -49,7 +52,7 @@
     read(unit=*, nml=conf_phys_nml)
     write(unit_nml, nml=conf_phys_nml)
 
-    call assert(iflag_pbl <=2 .or. iflag_pbl >= 6, &
+    call assert(any(iflag_pbl == [0, 1, 6, 8, 9]), &
          "conf_phys: bad value for iflag_pbl")
     call read_clesphys