!
! $Id: write_histmthNMC.h 1665 2012-10-09 13:35:26Z fairhead $
!
      IF (ok_histNMC(1)) 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_mthnmc,lNMC,"tnondef",itau_w,                           &
       &tnondef(:,:,1))
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"ta",itau_w,                                &
       &twriteSTD(:,:,1))
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"zg",itau_w,                                &
       &phiwriteSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"hus",itau_w,                              &
       &qwriteSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"hur",itau_w,                              &
       &rhwriteSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"ua",itau_w,                               &
       &uwriteSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"va",itau_w,                               &
       &vwriteSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"wap",itau_w,                              &
       &wwriteSTD(:,:,1))
!c
       DO k=1, nlevSTD
        DO i=1, klon
         IF(tnondef(i,k,1).NE.missing_val) THEN
          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,1))/freq_moyNMC(1)
         ELSE
          zx_tmp_fiNC(i,k) = missing_val
         ENDIF
        ENDDO
       ENDDO !k=1, nlevSTD
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"psbg",itau_w,                              &
       &zx_tmp_fiNC)
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"uv",itau_w,                                &
       &uvsumSTD(:,:,1))
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"vq",itau_w,                                &
       &vqsumSTD(:,:,1))
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"vT",itau_w,                                &
       &vTsumSTD(:,:,1))
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"wq",itau_w,                                &
       &wqsumSTD(:,:,1))
!c
      CALL histwrite_phy(nid_mthnmc,lNMC,"vphi",itau_w,                              &
       &vphisumSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"wT",itau_w,                               &
       &wTsumSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"uxu",itau_w,                              &
       &u2sumSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"vxv",itau_w,                              &
       &v2sumSTD(:,:,1))
!c
       CALL histwrite_phy(nid_mthnmc,lNMC,"TxT",itau_w,                              &
       &T2sumSTD(:,:,1))
!c
       DO k=1, nlevSTD
        DO i=1, klon
         IF(O3sumSTD(i,k,1).NE.missing_val) THEN
          zx_tmp_fiNC(i,k) = O3sumSTD(i,k,1) * 1.e+9
         ELSE
          zx_tmp_fiNC(i,k) = missing_val
         ENDIF
        ENDDO
       ENDDO !k=1, nlevSTD
       CALL histwrite_phy(nid_mthnmc,lNMC,"tro3",itau_w,                             &
       & zx_tmp_fiNC)
!c
       if (read_climoz == 2) THEN
       DO k=1, nlevSTD
        DO i=1, klon
         IF(O3daysumSTD(i,k,1).NE.missing_val) THEN
          zx_tmp_fiNC(i,k) = O3daysumSTD(i,k,1) * 1.e+9
         ELSE
          zx_tmp_fiNC(i,k) = missing_val
         ENDIF
        ENDDO
       ENDDO !k=1, nlevSTD
!c
        CALL histwrite_phy(nid_mthnmc,lNMC,"tro3_daylight",                          &
       &itau_w, zx_tmp_fiNC)
       endif 
!c
      if (ok_sync) then
!$OMP MASTER
        call histsync(nid_mthnmc)
!$OMP END MASTER
      endif
!c
      ENDIF !(ok_histNMC(1)) THEN