ncid = nccre(outname2(1:IL+4), NCCLOB, iret)

      partdim = ncddef(ncid, 'Particles', npart, iret)

      Dep_tdim  = ncddef(ncid, 'Depth_t',     Nr, iret)

      Dep_wm1dim= ncddef(ncid, 'Depth_wm1',   Nr, iret)

      Dep_wdim  = ncddef(ncid, 'Depth_w',   Nr+1, iret)

      Timedim = ncddef(ncid, 'Time', NCUNLIM, iret)

      partid  = ncvdef (ncid,'Particles',NCFLOAT,1,dims,iret)

      Dep_tid  = ncvdef (ncid,'Depth_t',    NCFLOAT,1,dims,iret)

      Dep_wid  = ncvdef (ncid,'Depth_w',    NCFLOAT,1,dims,iret)

      Dep_wm1id= ncvdef (ncid,'Depth_wm1',  NCFLOAT,1,dims,iret)

      Timeid   = ncvdef (ncid,'Time',   NCFLOAT,1,dims,iret)

      xpartid = ncvdef (ncid,'xpart', NCFLOAT,2,dims,iret)

      ypartid = ncvdef (ncid,'ypart', NCFLOAT,2,dims,iret)

      kpartid = ncvdef (ncid,'kpart', NCFLOAT,2,dims,iret)

      zetaid  = ncvdef (ncid,'zeta',  NCFLOAT,2,dims,iret)

      uvelid  = ncvdef (ncid,'uvel',  NCFLOAT,3,dims,iret)

      vvelid  = ncvdef (ncid,'vvel',  NCFLOAT,3,dims,iret)

      tempid  = ncvdef (ncid,'temp',  NCFLOAT,3,dims,iret)

      saltid  = ncvdef (ncid,'salt',  NCFLOAT,3,dims,iret)

      call NCAPTC(ncid, partid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, partid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPTC(ncid, Timeid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, Timeid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPTC(ncid, Timeid,'time_origin',NCCHAR, 20,startDate, iret)

      call NCAPTC(ncid, Dep_tid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, Dep_tid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPTC(ncid, Dep_tid, 'positive',  NCCHAR, 4, 'down',iret)     

     & call NCAPTC(ncid, Dep_tid, 'edges',NCCHAR, 7,'Depth_w',iret)     

      call NCAPTC(ncid, Dep_wm1id, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, Dep_wm1id, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPTC(ncid, Dep_wm1id, 'positive',  NCCHAR, 4, 'down',iret)     

      call NCAPTC(ncid, Dep_wid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, Dep_wid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPTC(ncid, Dep_wid, 'positive',  NCCHAR, 4, 'down',iret)     

      integer  iret, ncid, VARid

      call NCAPTC(ncid, xpartid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, xpartid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, xpartid,'missing_value',NCFLOAT,1,floatval,iret)

      call NCAPT (ncid, xpartid,'_FillValue', NCFLOAT, 1,floatval, iret)

      call NCAPTC(ncid, ypartid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, ypartid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, ypartid,'missing_value',NCFLOAT,1,floatval,iret)

      call NCAPT (ncid, ypartid,'_FillValue', NCFLOAT, 1,floatval, iret)

      call NCAPTC(ncid, kpartid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, kpartid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, kpartid,'missing_value',NCFLOAT,1,floatval,iret)

      call NCAPT (ncid, kpartid,'_FillValue', NCFLOAT, 1,floatval, iret)

      call NCAPTC(ncid, tempid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, tempid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, tempid, 'missing_value',NCFLOAT,1,floatval,iret)

      call NCAPT (ncid, tempid, '_FillValue', NCFLOAT, 1,floatval, iret)

      call NCAPTC(ncid, saltid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, saltid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, saltid, 'missing_value',NCFLOAT,1,floatval,iret)

      call NCAPT (ncid, saltid, '_FillValue', NCFLOAT, 1,floatval, iret)

      call NCAPTC(ncid, uvelid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, uvelid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, uvelid, 'missing_value',NCFLOAT,1,floatval,iret)

      call NCAPT (ncid, uvelid, '_FillValue', NCFLOAT, 1,floatval, iret)

      call NCAPTC(ncid, vvelid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, vvelid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, vvelid, 'missing_value',NCFLOAT,1,floatval,iret)

      call NCAPT (ncid, vvelid, '_FillValue', NCFLOAT, 1,floatval, iret)

      call NCAPTC(ncid, zetaid, 'long_name', NCCHAR, 24, name, iret) 

      call NCAPTC(ncid, zetaid, 'units',     NCCHAR, 16, unit, iret) 

      call NCAPT (ncid, zetaid,'missing_value',NCFLOAT, 1,floatval,iret)

      call NCAPT (ncid, zetaid,'_FillValue', NCFLOAT, 1, floatval, iret)

      call NCAPTC(ncid, NCGLOBAL, 'title',   NCCHAR, 60, expnam, iret)

      call NCAPTC(ncid, NCGLOBAL, 'history', NCCHAR, 32, stamp, iret)

      call NCENDF(ncid, iret)

      call NCVPT(ncid, partid, corner, edges, pnum, iret)

      call NCVPT(ncid, Timeid, corner, edges, Time, iret)

      call NCVPT(ncid, Dep_tid, corner, edges, Dep_t, iret)

      call NCVPT(ncid, Dep_wid, corner, edges, Dep_w, iret)

      call NCVPT(ncid, Dep_wm1id, corner, edges, Dep_wm1, iret)

      call NCVPT(ncid, VARid, corner, edges, xpart, iret)

      call NCVPT(ncid, VARid, corner, edges, ypart, iret)

      call NCVPT(ncid, VARid, corner, edges, kpart, iret)

      call NCVPT(ncid, VARid, corner, edges, zeta, iret)

      call NCVPT(ncid, VARid, corner, edges, temp, iret)

         call NCVPT(ncid, VARid, corner, edges, temp, iret)

         call NCVPT(ncid, VARid, corner, edges, temp, iret)

         call NCVPT(ncid, VARid, corner, edges, temp, iret)

         call NCVPT(ncid, VARid, corner, edges, salt, iret)

         call NCVPT(ncid, VARid, corner, edges, uvel, iret)

         call NCVPT(ncid, VARid, corner, edges, vvel, iret)

      call NCCLOS (ncid, iret)