anchor.diff:                endif
                 i = i + 1
                 call getarg(i,dspf)
  +          elseif (tag(1:4).eq.'-thr') then
  +              i = i + 1
  +              call getarg(i,tag)
  +              read(tag,'(BN,F10.0)') sthr
             elseif (tag(1:3).eq.'-az') then
                 if (kdsp.eq.1) then
                     kdsp = 2
  
positive.code:           write(*,*) '!! Error: no input faults read'
            call usage('!! Check input files for correct format')
        endif
  C     If FFM not used, use unweighted mean coordinates as hylo and hyla
        if (kffm.ne.1) then
            call calcmean(hylo,evlo,nflt)
  
negative.code:       REAL*8 vp,vs,dens
        INTEGER nsta,ntrg,typ(FMAX)
  
        REAL*8 sthr            ! threshold of slip for FFM
  
  C---- Input file markers
        INTEGER kffm           ! 0:no 1:yes
        INTEGER kflt           ! 0:no 1:yes
  