Index: LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90	(revision 2978)
+++ LMDZ5/trunk/libf/phylmd/phyaqua_mod.F90	(revision 2979)
@@ -349,4 +349,5 @@
     wake_deltaq = 0.
     wake_s = 0.
+    wake_dens = 0. 
     wake_cstar = 0.
     wake_pe = 0.
@@ -358,5 +359,5 @@
     ale_bl_trig =0. 
     alp_bl =0.
-
+    treedrg(:,:,:)=0.
 
 
Index: LMDZ5/trunk/libf/phylmd/phyetat0.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/phyetat0.F90	(revision 2978)
+++ LMDZ5/trunk/libf/phylmd/phyetat0.F90	(revision 2979)
@@ -323,18 +323,20 @@
   ENDIF
 !FC
-     IF (ifl_pbltree>0) then
+  IF (ifl_pbltree>0) THEN
 !CALL get_field("FTER", pctsrf(:, is_ter), found)
-     treedrg(:,1:klev,1:nbsrf)= 0.0
-   CALL get_field("treedrg_ter", drg_ter(:,:), found)
+    treedrg(:,1:klev,1:nbsrf)= 0.0
+    CALL get_field("treedrg_ter", drg_ter(:,:), found)
 !  found=phyetat0_srf(1,treedrg,"treedrg","drag from vegetation" , 0.)
-  !lecture du profile de freinage des arbres
-   IF (.not. found ) then
-     treedrg(:,1:klev,1:nbsrf)= 0.0
-   else
-     treedrg(:,1:klev,is_ter)= drg_ter(:,:)
+    !lecture du profile de freinage des arbres
+    IF (.not. found ) THEN
+      treedrg(:,1:klev,1:nbsrf)= 0.0
+    ELSE
+      treedrg(:,1:klev,is_ter)= drg_ter(:,:)
 !     found=phyetat0_srf(klev,treedrg,"treedrg","freinage arbres",0.)
-   ENDIF
-      endif
-
+    ENDIF
+  ELSE
+    ! initialize treedrg(), because it will be written in restartphy.nc
+    treedrg(:,:,:) = 0.0 
+  ENDIF
 
   ! Lecture de l'age de la neige:
Index: LMDZ5/trunk/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2978)
+++ LMDZ5/trunk/libf/phylmd/physiq_mod.F90	(revision 2979)
@@ -2082,5 +2082,7 @@
        CALL zenang_an(iflag_cycle_diurne.GE.1,jH_cur, &
             latitude_deg,longitude_deg,rmu0,fract)
-       JrNt = 1.0
+       swradcorr(:) = 1.0
+       JrNt(:) = 1.0
+       zrmu0(:) = rmu0(:)
     ELSE
        ! recode par Olivier Boucher en sept 2015