--- trunk/Sources/phylmd/ini_histins.f	2016/05/09 19:56:28	191
+++ trunk/Sources/phylmd/ini_histins.f	2016/08/30 12:52:46	206
@@ -155,7 +155,6 @@
             "inst(X)", zsto, zout)
 
        DO nsrf = 1, nbsrf
-
           call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &
                "% "//clnsurf(nsrf), "%",   &
                iim, (jjm + 1), nhori, 1, 1, 1, -99, &
@@ -200,8 +199,8 @@
                "rugosite "//clnsurf(nsrf), "-",   &
                iim, (jjm + 1), nhori, 1, 1, 1, -99, &
                "inst(X)", zsto, zout)
-          !XXX
        END DO
+
        CALL histdef(nid_ins, "rugs", "rugosity", "-", &
             iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
             "inst(X)", zsto, zout)
@@ -251,8 +250,13 @@
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &
             "inst(X)", zsto, zout)
 
-       if (conv_emanuel) CALL histdef(nid_ins, "ptop", "cloud top pressure", &
-            "Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
+       if (conv_emanuel) then
+          CALL histdef(nid_ins, "ptop", "cloud top pressure", &
+               "Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
+          CALL histdef(nid_ins, "dnwd0", "unsaturated downdraft", &
+               "kg/m2/s", iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", &
+               zsto, zout)
+       end if
 
        ! Champs 3D: