c     IF (ok_mensuelNMC) THEN
      IF (ok_mensuel) THEN
c
      ndex2d = 0
      ndex3d = 0
c
c Champs 2D:
c

        zsto = dtime
        zout = dtime * ecrit_mth
        itau_w = itau_phy + itap

      i = NINT(zout/zsto)
c
c  Champs interpolles sur des niveaux de pression du NMC
c110304 BEG
      DO k=1, nlevSTD
c
         bb=clevSTD(k)
c
         IF(k.GE.2) THEN
          aa=clevSTD(k)
          bb=aa(1:lnblnk1(aa))
         ENDIF
c
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlevSTD(:,k),zx_tmp_2d)
      CALL histwrite(nid_nmc,"t"//bb,itau_w,zx_tmp_2d,
     $       iim*jjmp1,ndex2d)
c
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, philevSTD(:,k),zx_tmp_2d)
      CALL histwrite(nid_nmc,"phi"//bb,itau_w,zx_tmp_2d,
     $       iim*jjmp1,ndex2d)
c
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlevSTD(:,k),zx_tmp_2d)
      CALL histwrite(nid_nmc,"q"//bb,itau_w,zx_tmp_2d,
     $       iim*jjmp1,ndex2d)
c
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlevSTD(:,k),zx_tmp_2d)
      CALL histwrite(nid_nmc,"rh"//bb,itau_w,zx_tmp_2d,
     $       iim*jjmp1,ndex2d)
c
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulevSTD(:,k),zx_tmp_2d)
      CALL histwrite(nid_nmc,"u"//bb,itau_w,zx_tmp_2d,
     $       iim*jjmp1,ndex2d)
c
      CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlevSTD(:,k),zx_tmp_2d)
      CALL histwrite(nid_nmc,"v"//bb,itau_w,zx_tmp_2d,
     $       iim*jjmp1,ndex2d)
c
      ENDDO !k=1, nlevSTD
ccc
c110304 END
c
      if (ok_sync) then
        call histsync(nid_nmc)
      endif

      ENDIF