anchor.diff:       1                   hylo,hyla,nflt)
   C----
   C Read shear dislocations from FFM in standard subfault format.
  +C All units are converted to SI, angles are in degrees.
   C----
         IMPLICIT NONE
         CHARACTER*40 ffmf,du,dxc,dyc
  
positive.code:               prog = 1
                if (ksta.eq.0) staf = 'autosta.dat'
                ksta = 1
                i = i + 1
                call getarg(i,tag)
                read(tag,'(BN,I3)') auto
                i = i + 1
                call getarg(i,tag)
                read(tag,'(BN,F10.0)') val
                i = i + 1
                call getarg(i,tag)
                read(tag,'(BN,F10.0)') incr
            elseif (tag(1:4).eq.'-gmt') then
                kgmt = 1
                i = i + 1
  
negative.code:       INTEGER kdsp,kstn,ksts,knor,kshr,kcou,kgmt
        COMMON /ICHECK/ kffm,kflt,kmag,ksta,khaf,ktrg
        COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
  C Look for station file, count number of receivers
        if (ksta.eq.0) then
            write(*,*) '!! Error: No station file specified'
            call usage('!! Use -sta STAFILE or an -auto option')
  