anchor.diff:            do 221 j = 1,3
                 strain(i,j) = 0.0d0
     221     continue
  -  222 continue 
  +  222 continue
         do 225 f = 1,nflt
             if (xy.eq.0) then
                 call ddistaz(dist,az,evlo(f),evla(f),stlo,stla)
  
positive.code:       call gcmdln(ffmfile,fltfile,magfile,stafile,haffile,trgfile,
       1            dspfile,stnfile,stsfile,norfile,shrfile,coufile,
       2            gmtfile,flttyp,auto,incr,depth,xy)
        call dbggcmdln(ffmfile,fltfile,magfile,stafile,haffile,
       1                     trgfile,dspfile,stnfile,stsfile,norfile,
       2                     shrfile,coufile,gmtfile,flttyp,auto,incr)
  
  C----
  C Check for input fault source files, and whether output is specified
  
negative.code: C----
  C Verify that an output file is specified.
  C----
        IMPLICIT NONE
        INTEGER kdsp,kstn,ksts,knor,kshr,kcou,kgmt
        COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
        if (kdsp.eq.0.and.kstn.eq.0) then
  