c     -----------------------------------------------------------------
      program matchCH
      parameter (n=13000)
      real*8 wl(n), Jl(n), Ju(n), loggf(n), w(n), lgf(n), jlo(n)
      real*8 elo(n), jup(n), eup(n), code, wll, cod, dwve, dlgf
      real*8 wlx, wll
      integer i, k, kp
      character*4 des(n)
      character*3 family(n)
      character*5 G(n),U(n)

      open (unit = 100, file='/priv/magus4/bessell/OH/CHBX_lifbase.dat',
     &     status='old',access='sequential')
      open (unit = 200, file='/priv/magus4/bessell/OH/kurucz_CHall.dat',
     &     status='old',access='sequential')
      open (unit = 300, file='combine_CHBX.dat')

 1000 format(2x,f9.4,f7.3,f5.1,12x,f5.1,12x,f9.2,a3,2x,a4)
c 2000 format(2x,f9.4,f7.3,f5.1,f12.3,f5.1,f12.3,F9.2)
 2000 format(1x,f9.4,f7.3,f5.1,f10.3,f6.1,f11.3,1x,f3.0,A5,3x,A5)
c 3000 format(2x,f9.4,f7.3,f7.1,f12.3,f5.1,f12.3,f9.2,f8.2,f11.6,
c     &             2x,a3,2x,a4,2f7.3)
 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,2f7.3,1x,a5,a5)
 4000  format(2x,f9.4,f7.3,f5.1,12x,f5.1,2x,a3,2x,a4)

c*****Read the LifBase data (wavelengths in Angstroms)
       i=0
       kp=0
 900  i=1+1
c
c  lifbase data
c
      read(100,1000,end=5001) wl(i), loggf(i), Jl(i), Ju(i),cod, 
     &                           family(i), des(i)
c            wll=wl(i)/10.0 - 0.22
c            wlx=wl(i)/10.0 + 0.1
            wll=wl(i)/10.0 - 0.02
            wlx=wl(i)/10.0 + 0.04
 901        k=kp
 902        if(w(k).gt.wlx) then
               go to 5000
               else
                  k=k+1
            endif
c
c     Kurucz data
c
          read(200,2000,end=5000) w(k),lgf(k),jlo(k),elo(k),jup(k),
     &                               eup(k),code,G(k),U(k)
c            type *,wl(i),10.0*w(k)
c            type *, Jl(i),jlo(k)
c            type *, Ju(i),jup(k)
            if(wll.gt.w(k)) then
               go to 902
               else
               if(Jl(i).ne.jlo(k)) then
                  go to 902
                  else
                     if(Ju(i).ne.jup(k)) then
                        go to 902
                        else
                     dlgf=loggf(i)-lgf(k)
                     dwve=wl(i) - 10.0*w(k)
                     write(300,3000) w(k),lgf(k),jlo(k),elo(k),jup(k),
     &  eup(k),code,wl(i),loggf(i),family(i),des(i),dwve,dlgf,G(k),U(k)
                     write(6,3000) w(k),lgf(k),jlo(k),elo(k),jup(k),
     &  eup(k),code,wl(i),loggf(i),family(i),des(i),dwve,dlgf,G(k),U(k)
c                   go to 900
                  endif
               endif
            endif
       rewind 200
       go to 900
 5000  write(300,4000)wl(i),loggf(i),Jl(i),Ju(i),family(i),des(i)
       rewind 200
       go to 900
 5001  end