--- trunk/Sources/phylmd/conf_phys.f	2017/02/27 15:44:55	213
+++ trunk/Sources/phylmd/conf_phys.f	2017/11/06 17:20:45	229
@@ -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
@@ -19,6 +22,7 @@
     use clesphys2, only: read_clesphys2
     USE comfisrtilp, ONLY: cld_lc_con, cld_lc_lsc, cld_tau_con, &
          cld_tau_lsc, coef_eva, ffallv_con, ffallv_lsc, iflag_pdf, reevap_ice
+    use nr_util, only: assert
     use unit_nml_m, only: unit_nml
     USE yomcst, ONLY: read_YOMCST
 
@@ -48,6 +52,8 @@
     read(unit=*, nml=conf_phys_nml)
     write(unit_nml, nml=conf_phys_nml)
 
+    call assert(any(iflag_pbl == [0, 1, 6, 8, 9]), &
+         "conf_phys: bad value for iflag_pbl")
     call read_clesphys
 
     print *, "Enter namelist 'nuagecom'."