!
! $Header$
!
      IF (ok_histNMC(3)) THEN
!c
       ndex3d = 0
       itau_w = itau_phy + itap + start_time * day_step / iphysiq
!ccc
!c  Champs interpolles sur des niveaux de pression du NMC
!c
!c     PARAMETER(nout=3) 
!c nout=1 : in=pdtphys,    out=mth
!c nout=2 : in=pdtphys,    out=day
!c nout=3 : in=pdtphys,    out=hf
!ccc
      CALL histwrite_phy(nid_hfnmc,lNMC,"tnondef",itau_w,                            &
       &tnondef(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"ta",itau_w,                                &
       &twriteSTD3)
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"zg",itau_w,                                &
       &phiwriteSTD3)
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"hus",itau_w,                               &
       &qwriteSTD3)
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"hur",itau_w,                               &
       &rhwriteSTD3)
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"ua",itau_w,                                &
       &uwriteSTD3)
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"va",itau_w,                                &
       &vwriteSTD3)
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"wap",itau_w,                               &
       &wwriteSTD3)
!c
       IF (1.EQ.0) THEN
!c
       DO k=1, nlevSTD
        DO i=1, klon
         IF(tnondef(i,k,3).NE.missing_val) THEN
          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(3)
         ELSE
          zx_tmp_fiNC(i,k) = missing_val
         ENDIF
        ENDDO
       ENDDO !k=1, nlevSTD
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"psbg",itau_w,                              &
       &zx_tmp_fiNC)
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"uv",itau_w,                                &
       &uvsumSTD(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"vq",itau_w,                                &
       &vqsumSTD(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"vT",itau_w,                                &
       &vTsumSTD(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"wq",itau_w,                                &
       &wqsumSTD(:,:,3))
!c
      CALL histwrite_phy(nid_hfnmc,lNMC,"vphi",itau_w,                               &
       &vphisumSTD(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"wT",itau_w,                                &
       &wTsumSTD(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"uxu",itau_w,                               &
       &u2sumSTD(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"vxv",itau_w,                               &
       &v2sumSTD(:,:,3))
!c
       CALL histwrite_phy(nid_hfnmc,lNMC,"TxT",itau_w,                               &
       &T2sumSTD(:,:,3))
!c
       ENDIF !(1.EQ.0) THEN
!c
      if (ok_sync) then
!$OMP MASTER
        call histsync(nid_hfnmc)
!$OMP END MASTER
      endif
!c
      ENDIF !      (ok_histNMC(3)) THEN