icount      = INT(tmp(5))

      print*, 'timesteps= ',icount

      allocate (time(icount))

      allocate (xpart(npart,icount))

      allocate (ypart(npart,icount))

      allocate (kpart(npart,icount))

      allocate (temp(npart,Nr,icount))

         allocate (uvel(npart,Nr,icount))

         allocate (vvel(npart,Nr,icount))

         allocate (salt(npart,Nr,icount))

      allocate (zeta(npart,icount))

         do n=1,icount

      do m=2,icount

      print*,icount,' x ',npart,' = ',icount*npart,' records expected,',

      print*,'==> ',icount*npart-itotalrecord,' float records missing'

      edges(1) = icount

      edges(2) = icount

      edges(3) = icount

       icount      = INT(tmp(5))

      print*, 'timesteps= ',icount

         icount=icount-1

         print*, 'preliminary --> use timesteps= ',icount

      allocate (time(icount))

      allocate (xpart(npart,icount))

      allocate (ypart(npart,icount))

      allocate (kpart(npart,icount))

      allocate (uvel(npart,icount))

      allocate (vvel(npart,icount))

      allocate (temp(npart,icount))

      allocate (salt(npart,icount))

      allocate (zeta(npart,icount))

         do n=1,icount

      do m=2,icount

               if (preflag .and. (np .gt. icount .or. np .lt. 1))

      print*,icount,' x ',npart,' = ',icount*npart,' records expected,',

      print*,'==> ',icount*npart-itotalrecord,' float records missing'

      edges(1) = icount

      edges(2) = icount