anchor.diff:  
   C----
   C Compare number of lines in stations.txt and targetparams.txt
  -C and use line count of stations.txt to initialize progress indicator
   C----
         stact = 0
     112 read (22,*,end=113)
  
positive.code: 
  C======================================================================C
  
        SUBROUTINE gcmdln(flttyp,verbose)
        IMPLICIT none
        CHARACTER*5 tag
        INTEGER narg,i,flttyp
        LOGICAL verbose
  
        flttyp = 1
        verbose = .false.
  
        narg = iargc()
        i = 0
     11 i = i + 1
        if (i.gt.narg) goto 12
            call getarg(i,tag)
            if (tag(1:2).eq.'-p') then
                flttyp = 0
            elseif (tag(1:2).eq.'-f') then
                flttyp = 1
            elseif (tag(1:2).eq.'-V') then
                verbose = .true.
            elseif (tag(1:2).eq.'-h'.or.tag(1:2).eq.'-?') then
                call usage() 
            endif
        goto 11
     12 continue
  
        RETURN
        END
  
  C----------------------------------------------------------------------C
  
        SUBROUTINE readstaticout(nflt,evlo,evla,evdp,str,dip,rak,dx,dy,
       1                         area,slip,maxflt)
        IMPLICIT none
        CHARACTER*1 dum
        INTEGER ct,seg,nseg,nx,ny,i,maxflt,nflt,ptr
        CHARACTER*10 dxc,dyc
        REAL*8 evlo(maxflt),evla(maxflt),evdp(maxflt),str(maxflt),
       1       dip(maxflt),rak(maxflt),dx(maxflt),dy(maxflt),area(maxflt),
       2       slip(maxflt)
  C----
  C Read header with number of fault segments
  C----
        ct = 0
        read (21,*) dum,dum,dum,dum,nseg
  C----
  
negative.code:       REAL*8 pi,d2r
        PARAMETER (pi=4.0d0*datan(1.0d0),d2r=pi/1.8d2)
  
        INTEGER f,nflt,maxflt,trgmatchsta
        PARAMETER (maxflt=1000)
        REAL*8 evlo(maxflt),evla(maxflt),evdp(maxflt),str(maxflt),
       1       dip(maxflt),rak(maxflt),slip(maxflt),dx(maxflt),
  