!
! calcul moyennes globales
!
       zx_tmp_fi2d=bils*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gbils)
       zx_tmp_fi2d=evap*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gevap)
       zx_tmp_fi2d(:)=fevap(:, is_ter)*cell_area(:)
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gevapt)
       zx_tmp_fi2d=zxfluxlat*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,glat)
       zx_tmp_fi2d=(topsw0-toplw0)*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gnet0)
       zx_tmp_fi2d=(topsw-toplw)*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gnet)
       zx_tmp_fi2d=(rain_fall+snow_fall)*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,grain)
       zx_tmp_fi2d=zxtsol*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gtsol)
       zx_tmp_fi2d=zt2m*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gt2m)
       zx_tmp_fi2d=prw*cell_area
       CALL global_mean(zx_tmp_fi2d,cell_area,.TRUE.,gprw)
!
!$OMP MASTER
      if (is_mpi_root) then      
!
      ndex2d = 0
      itau_w = itau_phy + itap + int(start_time * day_step_phy)
!
! Variables globales
!
      zx_tmp_0d=R_ecc
      CALL histwrite(nid_ctesGCM,"R_ecc",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=R_peri
      CALL histwrite(nid_ctesGCM,"R_peri",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=R_incl
      CALL histwrite(nid_ctesGCM,"R_incl",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=solaire
      CALL histwrite(nid_ctesGCM,"solaire",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=co2_ppm
      CALL histwrite(nid_ctesGCM,"co2_ppm",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=CH4_ppb
      CALL histwrite(nid_ctesGCM,"CH4_ppb",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=N2O_ppb
      CALL histwrite(nid_ctesGCM,"N2O_ppb",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=CFC11_ppt
      CALL histwrite(nid_ctesGCM,"CFC11_ppt",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
      zx_tmp_0d=CFC12_ppt
      CALL histwrite(nid_ctesGCM,"CFC12_ppt",itau_w, &
                     zx_tmp_0d,np,ndex2d)
!
!=================================================================
! moyennes globales
!
      CALL histwrite(nid_ctesGCM,"bils",itau_w, &
                     gbils,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"evap",itau_w, &
                     gevap,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"evap_land",itau_w, &
                     gevapt,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"flat",itau_w, &
                     glat,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"nettop0",itau_w, &
                     gnet0,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"nettop",itau_w, &
                     gnet,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"precip",itau_w, &
                     grain,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"tsol",itau_w, &
                     gtsol,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"t2m",itau_w, &
                     gt2m,np,ndex2d)
      CALL histwrite(nid_ctesGCM,"prw",itau_w, &
                     gprw,np,ndex2d)
!=================================================================
!
      if (ok_sync) then
        call histsync(nid_ctesGCM)
      endif
!
      endif !(is_mpi_root) then      
!$OMP END MASTER