--- trunk/Sources/phylmd/ini_histins.f	2017/03/30 14:25:18	217
+++ trunk/phylmd/ini_histins.f	2018/07/20 14:30:23	279
@@ -12,7 +12,7 @@
 
     use clesphys, only: ecrit_ins, ok_instan
     use clesphys2, only: conv_emanuel
-    use dimens_m, only: iim, jjm, llm, nqmx
+    use dimensions, only: iim, jjm, llm, nqmx
     use disvert_m, only: presnivs
     use dynetat0_m, only: day_ref, annee_ref, rlatu, rlonv
     USE histbeg_totreg_m, ONLY : histbeg_totreg
@@ -48,7 +48,7 @@
        CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &
             presnivs/100., nvert)
        
-       CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &
+       CALL histdef(nid_ins, "phis", "surface geopotential", "m2 s-2", &
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &
             "once", zsto, zout)
        CALL histdef(nid_ins, "aire", "Grid area", "-", &
@@ -132,6 +132,8 @@
             "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)
+       CALL histdef(nid_ins, "zxfqcalving", "ice calving", "kg m-2 s-1", &
+            iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
 
        DO nsrf = 1, nbsrf
           call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &
@@ -170,6 +172,14 @@
                "rugosite "//clnsurf(nsrf), "-",   &
                iim, (jjm + 1), nhori, 1, 1, 1, -99, &
                "inst(X)", zsto, zout)
+          call histdef(nid_ins, "u10m_"//clnsurf(nsrf),  &
+               "zonal wind 10 m "//clnsurf(nsrf), "m s-1",   &
+               iim, (jjm + 1), nhori, 1, 1, 1, -99, &
+               "inst(X)", zsto, zout)
+          call histdef(nid_ins, "v10m_"//clnsurf(nsrf),  &
+               "meridional wind 10 m "//clnsurf(nsrf), "m s-1",   &
+               iim, (jjm + 1), nhori, 1, 1, 1, -99, &
+               "inst(X)", zsto, zout)
        END DO
 
        CALL histdef(nid_ins, "rugs", "rugosity", "-", &
@@ -188,8 +198,9 @@
        CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &
             "inst(X)", zsto, zout)
-       CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", &
-            iim, (jjm + 1), nhori, 1, 1, 1, -99, &
+       CALL histdef(nid_ins, "s_capCL", &
+            "Convective available potential energy for atmospheric boundary " &
+            // "layer", "J/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99, &
             "inst(X)", zsto, zout)
        CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &
@@ -200,13 +211,7 @@
        CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &
             "inst(X)", zsto, zout)
-       CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2", &
-            iim, (jjm + 1), nhori, 1, 1, 1, -99, &
-            "inst(X)", zsto, zout)
-       CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2", &
-            iim, (jjm + 1), nhori, 1, 1, 1, -99, &
-            "inst(X)", zsto, zout)
-       CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", &
+       CALL histdef(nid_ins, "qsurf", "Surface Air humidity", "", &
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &
             "inst(X)", zsto, zout)
 
@@ -258,6 +263,7 @@
        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, &