--- trunk/Sources/phylmd/ini_histins.f	2017/02/27 15:44:55	213
+++ trunk/Sources/phylmd/ini_histins.f	2017/03/30 14:25:18	217
@@ -6,7 +6,7 @@
 
 contains
 
-  subroutine ini_histins(dtime)
+  subroutine ini_histins(dtime, ok_newmicro)
 
     ! From phylmd/ini_histins.h, version 1.2, 2005/05/25 13:10:09
 
@@ -26,7 +26,8 @@
     USE ymds2ju_m, only: ymds2ju
 
     REAL, intent(in):: dtime ! pas temporel de la physique (s)
-
+    logical, intent(in):: ok_newmicro
+    
     ! Local:
     real zjulian, zsto, zout
     integer nhori, nvert, nsrf, iq, it
@@ -35,7 +36,7 @@
 
     print *, 'Call sequence information: ini_histins'
 
-    IF (ok_instan) THEN
+    test_ok_instan: IF (ok_instan) THEN
        zsto = dtime * ecrit_ins
        zout = dtime * ecrit_ins
        CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
@@ -129,6 +130,8 @@
        CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", &
             iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
             "inst(X)", zsto, zout)
+       CALL histdef(nid_ins, "msnow", "surface snow amount", "kg/m2", &
+            iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
 
        DO nsrf = 1, nbsrf
           call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &
@@ -255,6 +258,14 @@
        CALL histdef(nid_ins, "dtlw0", "CS LW radiation dT", &
             "K/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
             zsto, zout)
+       if (ok_newmicro) then
+          CALL histdef(nid_ins, "re", "cloud droplet effective radius", &
+               "micrometer", iim, jjm + 1, nhori, llm, 1, llm, nvert, &
+               "inst(X)", zsto, zout)
+          CALL histdef(nid_ins, "fl", &
+               "denominator of Cloud droplet effective radius", "", iim, &
+               jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
+       end if
 
        DO it = 1, nqmx - 2
           ! champ 2D
@@ -279,7 +290,7 @@
        ENDDO
 
        CALL histend(nid_ins)
-    ENDIF
+    ENDIF test_ok_instan
 
   end subroutine ini_histins