subroutine sread(ib, wav, str, wid, eng, nlines,stot,wav0
     & ,wmin,wmax)

      implicit double precision (o-z, a-h)

      dimension wav(1), str(1), wid(1), eng(1)

c         Opens and reads line list for band

      if(ib .eq. 1) then 
         open(unit=20, file='/home/asgard/yelle/lines/ch4.dat') 
      else if(ib .eq. 2) then
         open(unit=20, file='/home/asgard/yelle/lines/ch4a.dat')
      else if(ib .eq. 3) then
         open(unit=20, file='/home/asgard/yelle/lines/ch4b.dat')
      else if(ib .eq. 4) then
         open(unit=20, file='/home/asgard/yelle/lines/ch4c.dat')
      else if(ib .eq. 5) then
         open(unit=20, file='/home/asgard/yelle/lines/c2h2.dat')
      else if(ib .eq. 6) then 
         open(unit=20, file='/home/asgard/yelle/lines/c2h2a.dat')
      else if(ib .eq. 7) then
         open(unit=20, file='/home/asgard/yelle/lines/c2h6.dat')
      else if(ib .eq. 8) then
         open(unit=20, file='/home/asgard/yelle/lines/hcnb.dat')
      else if(ib .eq. 9) then
         open(unit=20, file='/home/asgard/yelle/lines/hcna.dat')
      end if
         read(unit=20, fmt=*) nlines
         stot = 0.d0
         wav0 = 0.d0
         do n = 1, nlines
            if(ib .ge. 7) then
            read(unit=20, fmt=910) wav(n), str(n), wid(n), eng(n)
            else
            read(unit=20, fmt=*) wav(n), str(n),  wid(n),  eng(n)
            end if
            stot = stot + str(n)
            wav0 = wav0 + str(n)*wav(n)
c            if(ib .eq. 2) then
c              write(unit=*,fmt='(i4,1p2e13.5)')
c     &         n,wav(n),str(n)
c            end if
         end do
      close(unit=20) 

 910  format(f10.6,e10.3,f5.3,f10.3)

      wav0 = wav0/stot 
      wmin = wav(1)
      wmax = wav(nlines)

      return 
      end