--- trunk/Sources/phylmd/ini_histins.f	2016/04/14 15:15:56	190
+++ trunk/Sources/phylmd/ini_histins.f	2017/01/12 12:31:31	212
@@ -2,16 +2,17 @@
 
   implicit none
 
+  integer, save:: nid_ins
+
 contains
 
-  subroutine ini_histins(dtime, ok_instan, nid_ins, itau_phy)
+  subroutine ini_histins(dtime)
 
     ! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09
 
-    use clesphys, only: ecrit_ins
+    use clesphys, only: ecrit_ins, ok_instan
     use clesphys2, only: conv_emanuel
     use dimens_m, only: iim, jjm, llm, nqmx
-    use dimphy, only: klon
     use disvert_m, only: presnivs
     use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv
     USE histbeg_totreg_m, ONLY : histbeg_totreg
@@ -21,12 +22,10 @@
     use indicesol, only: nbsrf, clnsurf
     use iniadvtrac_m, only: tname, ttext
     use nr_util, only: pi
+    use phyetat0_m, only: itau_phy
     USE ymds2ju_m, only: ymds2ju
 
     REAL, intent(in):: dtime ! pas temporel de la physique (s)
-    logical, intent(in):: ok_instan
-    integer, intent(out):: nid_ins
-    integer, intent(in):: itau_phy
 
     ! Local:
     real zjulian, zsto, zout
@@ -155,7 +154,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 +198,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,11 +249,19 @@
             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:
 
+       CALL histdef(nid_ins, "tro3", "ozone mole fraction", "-", &
+            iim, jjm + 1, nhori, llm, 1, llm, nvert, "inst(X)", zsto, zout)
+
        CALL histdef(nid_ins, "temp", "Temperature", "K", &
             iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
             "inst(X)", zsto, zout)