Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/clesphys.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/clesphys.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/clesphys.h	(revision 723)
@@ -26,6 +26,6 @@
        LOGICAL ok_isccp, ok_regdyn
        REAL lonmin_ins, lonmax_ins, latmin_ins, latmax_ins
-       INTEGER ecrit_ins, ecrit_hf, ecrit_hf2mth, ecrit_day
-       INTEGER ecrit_mth, ecrit_tra, ecrit_reg 
+       REAL ecrit_ins, ecrit_hf, ecrit_hf2mth, ecrit_day
+       REAL ecrit_mth, ecrit_tra, ecrit_reg 
 
        COMMON/clesphys/cycle_diurne, soil_model, new_oliq,
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_bilKP_ave.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_bilKP_ave.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_bilKP_ave.h	(revision 723)
@@ -5,5 +5,5 @@
 c
          zsto = dtime
-         zout = dtime * FLOAT(ecrit_day)
+         zout = dtime * ecrit_day
          typeval=tave
 c
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histISCCP.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histISCCP.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histISCCP.h	(revision 723)
@@ -47,5 +47,5 @@
      .                 lmaxm1, cldtopres, nvert,'down')
 c
-        IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
+      IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN
 c
 c variables a ecrire
@@ -63,5 +63,5 @@
      .                "ave(X)", zstophy,zout)
 c
-        ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
+      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
 c
          DO k=1, kmaxm1
@@ -81,5 +81,5 @@
      .                "ave(X)", zstophy,zout)
 c
-        ENDIF 
+      ENDIF 
 c
         CALL histend(nid_isccp)
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday.h	(revision 723)
@@ -5,6 +5,6 @@
 c
        zstophy = dtime
-       zstoday= dtime * FLOAT(ecrit_day)
-       zout = dtime * FLOAT(ecrit_day)
+       zstoday= dtime * ecrit_day
+       zout = dtime * ecrit_day
 c
          idayref = day_ref
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday_seri.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday_seri.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histday_seri.h	(revision 723)
@@ -2,8 +2,8 @@
 c $Header$
 c
-      IF (type_run.EQ."AMIP") THEN
+       IF (type_run.EQ."AMIP") THEN
 c
        zstophy = dtime
-       zout = dtime * FLOAT(ecrit_day)
+       zout = dtime * ecrit_day
 c
          idayref = day_ref
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf.h	(revision 723)
@@ -5,8 +5,8 @@
       IF (ok_hf) THEN
 c
-       zstohf = dtime * FLOAT(ecrit_hf)
+       zstohf = dtime * ecrit_hf
        zstophy = dtime
-       zstoday = dtime * FLOAT(ecrit_day)
-       zout = dtime * FLOAT(ecrit_hf)
+       zstoday = dtime * ecrit_day
+       zout = dtime * ecrit_hf
 c
          idayref = day_ref
@@ -190,4 +190,5 @@
       ENDIF !lev_histhf.GE.4
 c
+c
       ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
 c
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf3d.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf3d.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histhf3d.h	(revision 723)
@@ -8,6 +8,6 @@
 cIM     zstohf = dtime * FLOAT(ecrit_hf)
 cIM     zout = dtime * FLOAT(ecrit_hf)
-        zstohf = dtime * FLOAT(ecrit_hf)
-        zout = dtime * FLOAT(ecrit_hf)
+        zstohf = dtime * ecrit_hf
+        zout = dtime * ecrit_hf
 c
 c       PRINT*, 'La frequence de sortie hf3d est de ', ecrit_hf
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmth.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmth.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmth.h	(revision 723)
@@ -5,8 +5,8 @@
 c
        zstophy = dtime
-       zstomth = dtime * FLOAT(ecrit_mth) !valeurs mensuelles "instantannees"
+       zstomth = dtime * ecrit_mth !valeurs mensuelles "instantannees"
        zstorad = dtime * FLOAT(radpas)    !variables stockees toutes les 2h
-       zstohf = dtime *FLOAT(ecrit_hf)   !variables stockees toutes les 6h
-       zout = dtime * FLOAT(ecrit_mth)    !frequence des sorties : 1 val. par mois
+       zstohf = dtime * ecrit_hf   !variables stockees toutes les 6h
+       zout = dtime * ecrit_mth    !frequence des sorties : 1 val. par mois
 c
 cIM    PRINT*,' La zstophy,zstorad,zstohf,zstomth,zout=',
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmthNMC.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmthNMC.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histmthNMC.h	(revision 723)
@@ -5,7 +5,7 @@
 c
        zstophy = dtime
-       zstohf = dtime * FLOAT(ecrit_hf)
-       zstomth = dtime * FLOAT(ecrit_mth)
-       zout = dtime * FLOAT(ecrit_mth)
+       zstohf = dtime * ecrit_hf
+       zstomth = dtime * ecrit_mth
+       zout = dtime * ecrit_mth
 c
          idayref = day_ref
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histrac.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histrac.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/ini_histrac.h	(revision 723)
@@ -39,5 +39,5 @@
 
          zsto = pdtphys
-         zout = pdtphys * FLOAT(ecrit_tra)
+         zout = pdtphys * ecrit_tra
 c
          CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/moy_undefSTD.F
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/moy_undefSTD.F	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/moy_undefSTD.F	(revision 723)
@@ -1,4 +1,4 @@
       SUBROUTINE moy_undefSTD(nlevSTD,itap,
-     $           ecrit_day,ecrit_mth,ecrit_hf2mth,
+     $           dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
      $           tnondef,tsumSTD)
       IMPLICIT none
@@ -28,5 +28,5 @@
       INTEGER nlevSTD, klevSTD, itap
       PARAMETER(klevSTD=17)
-      INTEGER ecrit_day, ecrit_mth, ecrit_hf2mth
+      REAL dtime, ecrit_day, ecrit_mth, ecrit_hf2mth
 c
 c variables locales
@@ -40,10 +40,10 @@
 c calcul 1 fois par jour
 c
-      IF(MOD(itap,ecrit_day).EQ.0) THEN
+      IF(MOD(itap,NINT(ecrit_day/dtime)).EQ.0) THEN
        DO k=1, nlevSTD
         DO i=1, klon
-         IF(tnondef(i,k,1).NE.FLOAT(ecrit_day)) THEN
+         IF(tnondef(i,k,1).NE.(ecrit_day/dtime)) THEN
           tsumSTD(i,k,1)=tsumSTD(i,k,1)/
-     $    (FLOAT(ecrit_day)-tnondef(i,k,1))
+     $    (ecrit_day/dtime-tnondef(i,k,1))
          ELSE
           tsumSTD(i,k,1)=1.E+20
@@ -55,17 +55,17 @@
 c calcul 1 fois par mois
 c
-      IF(MOD(itap,ecrit_mth).EQ.0) THEN
+      IF(MOD(itap,NINT(ecrit_mth/dtime)).EQ.0) THEN
        DO k=1, nlevSTD
         DO i=1, klon
-         IF(tnondef(i,k,2).NE.FLOAT(ecrit_mth)) THEN
+         IF(tnondef(i,k,2).NE.ecrit_mth/dtime) THEN
           tsumSTD(i,k,2)=tsumSTD(i,k,2)/
-     $    (FLOAT(ecrit_mth)-tnondef(i,k,2))
+     $    (ecrit_mth/dtime-tnondef(i,k,2))
          ELSE
           tsumSTD(i,k,2)=1.E+20
          ENDIF !tnondef
 c
-         IF(tnondef(i,k,3).NE.FLOAT(ecrit_hf2mth)) THEN
+         IF(tnondef(i,k,3).NE.NINT(ecrit_hf2mth)) THEN
           tsumSTD(i,k,3)=tsumSTD(i,k,3)/
-     $    (FLOAT(ecrit_hf2mth)-tnondef(i,k,3))
+     $    (ecrit_hf2mth-tnondef(i,k,3))
          ELSE
           tsumSTD(i,k,3)=1.E+20
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/physiq.F
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/physiq.F	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/physiq.F	(revision 723)
@@ -38,8 +38,8 @@
 c   CLEFS CPP POUR LES IO
 c   =====================
+#define histins
 #define histhf
 #define histday
 #define histmth
-#define histins
 #define histISCCP
 #define histREGDYN
@@ -456,13 +456,4 @@
       INTEGER nbapp_isccp,isccppas
 
-#undef histISCCP
-#define histISCCP
-c     data ok_isccp,ecrit_isccp/.true.,0.125/      
-c     data ok_isccp,ecrit_isccp/.true.,1./      
-cIM 190504     data ok_isccp/.true./      
-cIM 190504 #else
-cIM 190504     data ok_isccp/.false./
-cIM 190504 #endif
-
 c sorties statistiques regime dynamique
 c     logical ok_regdyn
@@ -471,14 +462,4 @@
 c     save ok_regdyn, ecrit_regdyn, nid_regdyn
       save nid_regdyn
-
-#undef histREGDYN
-#define histREGDYN
-cIM 190504 #ifdef histREGDYN
-c     data ok_regdyn,ecrit_regdyn/.true.,0.125/
-c     data ok_regdyn,ecrit_regdyn/.true.,1./
-cIM 190504    data ok_regdyn/.true./
-cIM 190504 #else
-cIM 190504   data ok_regdyn/.false./
-cIM 190504 #endif 
 
       REAL zx_tau(kmaxm1), zx_pc(lmaxm1), zx_o500(iwmax)
@@ -1520,6 +1501,4 @@
 #endif
 
-c#undef histmthNMC
-c#define histmthNMC
 #ifdef histmthNMC
 #include "ini_histmthNMC.h"
@@ -2747,6 +2726,5 @@
 c====================================================================
 C
-      call phytrac (iflag_con,
-     I                   rnpb,
+      call phytrac (rnpb,
      I                   itap, 
      I                   julien, 
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/phytrac.F
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/phytrac.F	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/phytrac.F	(revision 723)
@@ -4,6 +4,5 @@
 c
 c
-      SUBROUTINE phytrac (iflag_con,
-     I                    rnpb,
+      SUBROUTINE phytrac (rnpb,
      I                    nstep,
      I                    julien,
@@ -156,5 +155,6 @@
       REAL flxmass_w(klon,klev)
 #endif
-      integer iflag_con
+cIM   integer iflag_con
+#include "clesphys.h"
 
 cAA Rem : nbtr : nombre de vrais traceurs est defini dans dimphy.h
@@ -266,6 +266,6 @@
 C Variables liees a l'ecriture de la bande histoire : phytrac.nc
 c
-      INTEGER ecrit_tra
-      SAVE ecrit_tra   
+cIM   INTEGER ecrit_tra
+cIM   SAVE ecrit_tra   
       logical ok_sync
       parameter (ok_sync = .true.)
@@ -387,6 +387,6 @@
          if (debutphy) then
 
-          ecrit_tra = NINT(86400./pdtphys *ecritphy) 
-          print*,'dans phytrac ',pdtphys,ecritphy,ecrit_tra
+cIM       ecrit_tra = NINT(86400./pdtphys *ecritphy) 
+cIM       print*,'dans phytrac ',pdtphys,ecritphy,ecrit_tra
 
          if(nbtr.lt.nqmax) then
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histmthNMC.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histmthNMC.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histmthNMC.h	(revision 723)
@@ -84,5 +84,5 @@
         DO i=1, klon
          IF(tnondef(i,k,3).NE.1.E+20) THEN
-          zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,3))/FLOAT(ecrit_hf2mth)
+          zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth
          ELSE
           zx_tmp_fi3d(i,k) = 1.E+20
Index: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h
===================================================================
--- LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h	(revision 665)
+++ LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histrac.h	(revision 723)
@@ -6,10 +6,6 @@
       ndex3d = 0
 c
-      zsto = pdtphys
-      zout = pdtphys * FLOAT(ecrit_tra)
       itau_w = itau_phy + nstep
-
       
-      i = NINT(zout/zsto)
       CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
       CALL histwrite(nid_tra,"phis",itau_w,zx_tmp_2d,iim*(jjm+1),ndex2d)