anchor.diff:  C     clip:  
   C----
         call gcmdln(x1,x2,nx,dx,y1,y2,ny,dy,z,ofile,p,ref,xsec,xz,clip,
  -     1            datain)
  +     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')
  
positive.code:       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
  
  C----
  C Get parameters from the command line
  
negative.code:           i = i + 1
            call getarg(i,tag)
            read (tag,'(BN,F12.0)') dy
        elseif (tag(1:4).eq.'-exp') then
            i = i + 1
            call getarg(i,tag)
            if(tag.eq.'x') then
                expflg = expflg + 1
            elseif (tag.eq.'y') then
                expflg = expflg + 10
            else
                call usage('!! Error: no -exp option'//trim(tag))
            endif
        elseif (tag(1:4).eq.'-dip') then
            do 102 j = 1,5
                call getarg(i+j,tag)
  