cIM    Implemente en modes sequentiel et parallele

       CALL gather(rlat,rlat_glo)
       CALL bcast(rlat_glo)
       CALL gather(rlon,rlon_glo)
       CALL bcast(rlon_glo)

c$OMP MASTER
      if (is_mpi_root) then
c
       zstophy = dtime
       zout = mth_len*un_jour
c
       idayref = day_ref
       CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
c
       CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon_glo,zx_lon)
       if (iim.gt.1) then
       DO i = 1, iim
         zx_lon(i,1) = rlon_glo(i+1)
         zx_lon(i,jjmp1) = rlon_glo(i+1)
       ENDDO
       endif
       CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat_glo,zx_lat)
c
       CALL histbeg("paramLMDZ_phy.nc", 
     .                 np,zx_lon(np:np,1), np,zx_lat(1,np:np),
     .                 1,1,1,1,
     .                 itau_phy, zjulian, dtime,
     .                 nhori, nid_ctesGCM)
c
       CALL histdef(nid_ctesGCM, "R_ecc", 
     .                "Excentricite","-",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "R_peri", 
     .                "Equinoxe","-",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "R_incl", 
     .                "Inclinaison","deg",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "solaire", 
     .                "Constante solaire","W/m2",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "co2_ppm", 
     .                "Concentration du CO2", "ppm",
     .                1,1,nhori, 1,1,1, -99, 32, 
     .                "ave(X)", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "CH4_ppb", 
     .                "Concentration du CH4", "ppb",
     .                1,1,nhori, 1,1,1, -99, 32, 
     .                "ave(X)", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "N2O_ppb",
     .                "Concentration du N2O", "ppb",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "CFC11_ppt",
     .                "Concentration du CFC11", "ppt",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "CFC12_ppt",
     .                "Concentration du CFC12", "ppt",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave(X)", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "bils",
     .                "Surface total heat flux", "W m-2",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "evap",
     .                "Evaporation", "kg m-2 s-1",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "evap_land",
     .                "Land evaporation", "kg m-2 s-1",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "flat",
     .                "Latent heat flux", "W m-2",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "nettop0",
     .                "Clear sky net downward radiatif flux at TOA", 
     .                "W m-2",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "nettop",
     .                "Net downward radiatif flux at TOA", "W m-2",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "precip",
     .                "Total precipitation (liq+sol)", "kg m-2 s-1",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "tsol",
     .                "Surface temperature", "K",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "t2m",
     .                "Temperature at 2m", "K",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout)
c
       CALL histdef(nid_ctesGCM, "prw",
     .                "Precipitable water", "kg m-2",
     .                1,1,nhori, 1,1,1, -99, 32,
     .                "ave", zstophy,zout) 
c=================================================================
c
       CALL histend(nid_ctesGCM)
       
       endif !(is_mpi_root)
c$OMP END MASTER
c=================================================================