c     -----------------------------------------------------------------

      parameter (n=2100)
      real*8 wl, g, ghigh, loggf, w(n), lgf(n), glo(n)
      real*8 elo(n), gup(n), eup(n), code, wll
      integer i, j, k
      character*4 des
      character*3 family

      open (unit = 100, file='/priv/magus4/bessell/OH/OH_lifbase.dat',
     &     status='old',access='sequential')
      open (unit = 200, file='/priv/magus4/bessell/OH/300-350.sort',
     &     status='old',access='sequential')
      open (unit = 300, file='combine.out')
 
 1000 format(f7.2,2x,f9.6,2x,f4.1,2x,8x,2x,f4.1,2x,
     &             a3,2x,a4,2x,a3) 
 2000 format(2x,f9.4,f7.3,f7.1,f12.3,f5.1,f12.3,f9.2)
 3000 format(2x,f9.4,f7.3,f7.1,f12.3,f5.1,f12.3,f9.2,f8.2,f11.6,
     &             2x,a3,2x,a4)
 4000  format(2x,f9.4,f7.3,f7.1,12x,f5.1,2x,a3,2x,a4)

c*****Read the LifBase data (wavelengths in Angstroms)
 900  read(100,1000,end=5001) wl, loggf, g, ghigh, family,
     &             des, "108"
            wll=wl/10.0 - 0.5
 901        k=0
 902        k=k+1
            read(200,2000,end=5000) w(k),lgf(k),glo(k),elo(k),gup(k),
     &                               eup(k),code
            if(wll.gt.w(k)) then
               go to 902
               else
               if(g.ne.glo(k)) then
                  go to 902
                  else
                     if(gup(k).ne.ghigh) then
                        go to 902
                        else
                     write(300,3000) w(k),lgf(k),glo(k),elo(k),gup(k),
     &                 eup(k),code,wl,loggf,family,des
                  endif
               endif
            endif
       go to 900
 5000  write(300,4000)wl,loggf,glo,gup,family,des
       go to 900
 5001  end