program gcm
      character*80 s80
      character*2 fn(3)
      data fn/'AB','AB','CD'/
      write(6,*)'extracts rotationally-resolved spectrum from gaussian'
      if(iargc().ne.1)then
       write(6,*)'usage: gcm <output>'
      else
       call getarg(1,s80)
       open(9,file=s80,status='old')
1      read(9,800,end=99,err=99)s80
800    format(a80)
       it=0
       if(s80(2:26).eq.'Axis Y = Molar absorption')it=1
       if(s80(2:26).eq.'Intensity: Molar absorpti')it=2
       if(s80(2:26).eq.'Intensity: Diff. in molar')it=3
       if(it.gt.0)then
        if(it.eq.1)read(9,*)
        read(9,*)
        open(91,file=fn(it)//'cm.prn')
        open(92,file=fn(it)//'nm.prn')
        i=0
2       read(9,800,end=99,err=99)s80
        if(s80(2:3).eq.'--')then
         close(9)
         close(91)
         close(92)
         write(6,*)i,' points to '//fn(it)//'nm.prn'
         stop
        endif
        read(s80,*)x,y
        write(91,910)x,y
        write(92,910)1.0e7/x,y
910     format(f17.4,e20.5)
        i=i+1
        goto 2 
       endif
       goto 1
99     close(9)
      endif
      end