subroutine read_nuk(kfile,nu,k,nl)
      implicit none
      include 'max.inc'
      include 'formats.inc'
      integer nl,strlen,i
      double precision tmp
      double precision nu(1:Npmx)
      double precision k(1:Npmx)
      character*100 kfile,sfile

      sfile='./nlines'
      call get_nlines(kfile,sfile,nl)

      if (nl.gt.Npmx) then
         write(*,*) 'Error from routine read_nuk:'
         write(*,*) 'nl=',nl
         stop
      endif

      write(*,*) 'reading file:',kfile(1:strlen(kfile))
      open(10,file=kfile(1:strlen(kfile)))
      do i=1,nl
         read(10,50) nu(i),tmp,k(i)
      enddo
      close(10)
      write(*,*) '...done'

      return
      end