anchor.diff:       1 '  -h/-?    Online help (this screen)'
         write(*,*)
         write(*,*)
  -     1 'polyfit numerically integrates a time series using trapezoidal',
  +     1 'numint numerically integrates a time series using trapezoidal',
        2 ' integration'
         write(*,*)
         write(*,*)
  
positive.code:       fint(1) = 0.0d0
        do 103 i = 2,ntot
            dx = x(i) - x(i-1)
            print *,fint(i-1),dx,f(i),f(i-1)
            fint(i) = fint(i-1) + dx*0.5d0*(f(i)+f(i-1))
    103 continue
  
  
negative.code:       call gcmdln(ifile,ofile)
        inquire(file=ifile,EXIST=ex)
        if (.not.ex) then
            write(*,*) 'No input file:',ifile
            call usage()
        endif
  
  