--- trunk/Sources/phylmd/ini_histins.f	2016/03/11 18:47:26	178
+++ trunk/Sources/phylmd/ini_histins.f	2016/03/29 15:20:23	189
@@ -9,10 +9,12 @@
     ! From phylmd/ini_histins.h, v 1.2 2005/05/25 13:10:09
 
     use clesphys, only: ecrit_ins
+    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
+    use gr_phy_write_m, only: gr_phy_write
     USE histbeg_totreg_m, ONLY : histbeg_totreg
     USE histdef_m, ONLY : histdef
     USE histend_m, ONLY : histend
@@ -40,12 +42,12 @@
 
        CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
 
-       CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlon, zx_lon)
+       zx_lon = gr_phy_write(rlon)
        DO i = 1, iim
           zx_lon(i, 1) = rlon(i+1)
           zx_lon(i, (jjm + 1)) = rlon(i+1)
        ENDDO
-       CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)
+       zx_lat = gr_phy_write(rlat)
        CALL histbeg_totreg("histins", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
             jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)
        write(*, *)'Inst ', itau_phy, zjulian
@@ -260,7 +262,8 @@
             iim, (jjm + 1), nhori, 1, 1, 1, -99, &
             "inst(X)", zsto, zout)
 
-       !IM cf. AM 081204 END
+       if (conv_emanuel) CALL histdef(nid_ins, "ptop", "cloud top pressure", &
+            "Pa", iim, jjm + 1, nhori, 1, 1, 1, -99, "inst(X)", zsto, zout)
 
        ! Champs 3D: