!
! $Header: /home/cvsroot/LMDZ4/libf/phylmd/ini_histrac.h,v 1.2 2004/10/18 15:42:20 lmdzadmin Exp $
!
         zsto = pdtphys
	 z4julian = zjulian
         CALL ymds2ju(annee_ref, 1, idayref, zero, z4julian)
c
         CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlon,zx_lon)
         DO i = 1, iim
            zx_lon(i,1) = xlon(i+1)
            zx_lon(i,jjm+1) = xlon(i+1)
         ENDDO
         CALL gr_fi_ecrit(1,klon,iim,jjm+1,xlat,zx_lat)
	 z4_lon = zx_lon
	 z4_lat = zx_lat
	 zpresnivs=presnivs
         CALL histbeg("histrac", iim,z4_lon(:,1), jjm+1,z4_lat(1,:),
     .                 1,iim,1,jjm+1, itau_phy, z4julian, zsto,
     .                 nhori, nid_tra)
         CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",
     .                 klev, zpresnivs, nvert)

         zout = pdtphys * FLOAT(ecrit_tra)
c
         CALL histdef(nid_tra, "phis", "Surface geop. height", "-",
     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     .                "once",  zsto,zout)
c
         CALL histdef(nid_tra, "aire", "Grid area", "-",
     .                iim,jjm+1,nhori, 1,1,1, -99, 32,
     .                "once",  zsto,zout)
c
         CALL histend(nid_tra)
         ndex2d = 0
         ndex3d = 0
         ndex = 0