Index: LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F
===================================================================
--- LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F	(revision 150)
+++ LMDZ.3.3/branches/rel-LF/libf/phylmd/phyetat0.F	(revision 151)
@@ -2,5 +2,5 @@
      .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsol,snow,
      .           albe, evap, rain_fall, snow_fall, solsw, sollw,
-     .           radsol,frugs,agesno,clesphy0,
+     .           fder,radsol,frugs,agesno,clesphy0,
      .           zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,tabcntr0,
      .           t_ancien,q_ancien,ancien_ok)
@@ -35,4 +35,5 @@
       REAL sollw(klon)
       real solsw(klon)
+      real fder(klon)
       REAL frugs(klon,nbsrf)
       REAL agesno(klon)
@@ -769,4 +770,31 @@
 
 c
+c Lecture derive des flux:
+c
+      ierr = NF_INQ_VARID (nid, "fder", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <fder> est absent'
+         PRINT*, 'mis a zero'
+         fder = 0.
+      ELSE
+#ifdef NC_DOUBLE
+        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, fder)
+#else
+        ierr = NF_GET_VAR_REAL(nid, nvarid, fder)
+#endif
+        IF (ierr.NE.NF_NOERR) THEN
+          PRINT*, 'phyetat0: Lecture echouee pour <fder>'
+          CALL abort
+        ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, klon
+         xmin = MIN(fder(i),xmin)
+         xmax = MAX(fder(i),xmax)
+      ENDDO
+      PRINT*,'Derive des flux fder:', xmin, xmax
+
+c
 c Lecture du rayonnement net au sol:
 c
@@ -826,5 +854,5 @@
               xmax = MAX(frugs(i,nsrf),xmax)
            ENDDO
-           PRINT*,'evap du sol RUG**:', nsrf, xmin, xmax
+           PRINT*,'rugosite du sol RUG**:', nsrf, xmin, xmax
          ENDDO
       ELSE
Index: LMDZ.3.3/branches/rel-LF/libf/phylmd/phyredem.F
===================================================================
--- LMDZ.3.3/branches/rel-LF/libf/phylmd/phyredem.F	(revision 150)
+++ LMDZ.3.3/branches/rel-LF/libf/phylmd/phyredem.F	(revision 151)
@@ -2,5 +2,5 @@
      .           rlat,rlon, pctsrf,tsol,tsoil,deltat,qsol,snow,
      .           albedo, evap, rain_fall, snow_fall,
-     .           solsw, sollw,
+     .           solsw, sollw,fder,
      .           radsol,frugs,agesno,
      .           zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,
@@ -37,4 +37,5 @@
       real solsw(klon)
       real sollw(klon)
+      real fder(klon)
       REAL radsol(klon)
       REAL frugs(klon,nbsrf)
@@ -436,4 +437,19 @@
       ierr = NF_REDEF (nid)
 #ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "fder", NF_DOUBLE, 1, idim2,nvarid)
+#else
+      ierr = NF_DEF_VAR (nid, "fder", NF_FLOAT, 1, idim2,nvarid)
+#endif
+      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 14,
+     .                        "Derive de flux")
+      ierr = NF_ENDDEF(nid)
+#ifdef NC_DOUBLE
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,fder)
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,fder)
+#endif
+c
+      ierr = NF_REDEF (nid)
+#ifdef NC_DOUBLE
       ierr = NF_DEF_VAR (nid, "rain_f", NF_DOUBLE, 1, idim2,nvarid)
 #else