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

      if (iret .ne. 0) write(stdout,*) 'Error: Open NetCDF file'

      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)

      if (iret .ne. 0) write(stdout,*) 'Error: define dimensions'

      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)

      if (iret .ne. 0) write(stdout,*) 'Error: define axis ids'

      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)

      if (iret .ne. 0) write(stdout,*) 'Error: define variable ids'

      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)

      if (iret .ne. 0) write(stdout,*) 'Error: assign axis attributes'

      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)

      if (iret .ne. 0) write(stdout,*) 'Error: define variable attrib.'

      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)

      if (iret .ne. 0) write(stdout,*) 'Error: write variables'

      call NCCLOS (ncid, iret)