--- trunk/Sources/dyn3d/gcm.f	2015/09/09 10:41:47	168
+++ trunk/Sources/dyn3d/gcm.f	2016/12/07 16:44:53	208
@@ -9,13 +9,13 @@
   ! l'advection de "q", en modifiant "iadv" dans "traceur.def".
 
   use comconst, only: daysec, dtvr, iniconst
+  use comdissnew, only: read_comdissnew
   use comgeom, only:  aire_2d, inigeom
   use comgeomphy, only: airephy
   use conf_gcm_m, only: day_step, iperiod, iecri, iphysiq, nday, periodav, &
        conf_gcm, iflag_phys
   use conf_guide_m, only: conf_guide
   use dimens_m, only: iim, jjm, llm, nqmx
-  use dimphy, only: klon
   USE disvert_m, ONLY : disvert
   use dynetat0_m, only: dynetat0, day_ini
   use dynredem0_m, only: dynredem0
@@ -58,6 +58,7 @@
   open(unit_nml, file="used_namelists.txt", status="replace", action="write")
 
   CALL conf_gcm
+  call read_comdissnew
 
   print *, "Enter namelist 'main_nml'."
   read (unit=*, nml=main_nml)
@@ -80,7 +81,7 @@
   call init_dyn_phy
 
   ! Initialisation de la physique :
-  IF (iflag_phys == 1) THEN
+  IF (iflag_phys) THEN
      airephy = pack(aire_2d, dyn_phy)
      CALL suphec
      call yoethf