PROGRAM reform
*************************************************************************
* reformat from inverse wavenumber to meV.                              *
* the standard input file must have the format x y (z)                  *
* The outputfile comes as: x y                                          *
*************************************************************************
      REAL X(100000),y(100000),xmn,xmx
      INTEGER I,mm
      character*40 flin
      read(*,'(a40)') flin
      read(*,*) xmn,xmx 
      open(23,file=flin)   
      mm=100000 
      do 10 i=1,mm                                                  
       READ(23,*,END=11) X(i),y(i)
10    continue
      close(23)   
11    mm=i-1
      do 25 i=1,mm
       if ((x(i).gt.xmn).and.(x(i).lt.xmx)) 
     *   write(*,*) x(i)/8.0657,y(i)
25    continue
      END