subroutine read_nu(Nb,nu_low,nu_hi)
      implicit none
      include 'max.inc'
      include 'formats.inc'
      integer strlen,ios,dim,Nb,i,band
      character*(Nchar_mx) file_nu
      character*(Nchar_mx) strtmp1,strtmp2
      double precision nu_low(1:Nbmx)
      double precision nu_hi(1:Nbmx)

      file_nu='./venus_nu.dat'
      open(10,file=file_nu(1:strlen(file_nu)),
     &     status='old',iostat=ios)
      if (ios.ne.0) then        ! file not found
         write(*,*) 'Error from routine read_nu:'
         write(*,*) 'File could not ne found:'
         write(*,*) file_nu(1:strlen(file_nu))
         stop
      endif
      do i=1,4
         read(10,*)
      enddo
      read(10,33) strtmp1,dim,strtmp2,Nb
      do i=1,3
         read(10,*)
      enddo
      do band=1,Nb
         read(10,*) nu_low(band),nu_hi(band)
      enddo
      close(10)

      return
      end