anchor.diff:        CHARACTER*200 clip,line
         REAL*8 x1,x2,y1,y2,dx,dy,dz,dd,ref(5),dist,az,lo,la
         REAL*8 x,y,z,xx(1000),yy(1000),zsave
  -      INTEGER i,j,nx,ny,p,xsec,xz,nn,io,datain
  +      INTEGER i,j,k,nx,ny,p,xsec,xz,nn,io,datain,expflg
   
   C----
   C Get parameters from the command line
  
positive.code:       write(*,*)
       1 '-o OFILE               Output to file (default prints to ',
       2              'standard output)'
        write (*,*)
       1 '-p                     Print results to standard output ',
       2                               '(overrides -o)'
  
negative.code: C     clip:  
  C----
        call gcmdln(x1,x2,nx,dx,y1,y2,ny,dy,z,ofile,p,ref,xsec,xz,clip,
       1            datain,expflg)
        if (datain.eq.1) goto 9001
  C Check that limits and increments are properly defined
        if (x2.lt.x1) call usage('!! Error: X2 must be greater than X1')
  