anchor.diff:        do 201 i = 1,nsta
             read(12,*) stlo,stla,stdp
             stdp = stdp*1.0d3
  -C         If NTRG=1, read target fault parameters here
  +C         If NTRG=NSTA, read target fault parameters here
             if (ntrg.eq.nsta.and.ktrg.eq.1.and.nsta.ne.1) then
                 read(13,*) trgstr,trgdip,trgrak,frict
             endif
  
positive.code:       REAL*8 pi,r2d,d2r
        PARAMETER (pi=4.0d0*atan(1.0d0),r2d=180.0d0/pi,d2r=pi/180.0d0)
        CHARACTER*40 staf
        REAL*8 xmin,xmax,ymin,ymax,incr,x0,y0,str0,dx,x1,y1,z
        INTEGER xy
        open(unit=96,file=staf,status='unknown')
        z = ymin
    961 if (z.le.ymax) then
            dx = xmin
    962     if (dx.le.xmax) then
                if (xy.eq.0) then
                    call dlola(x1,y1,x0,y0,dx,str0)
                    x1 = x1*r2d
                    y1 = y1*r2d
                else
                    x1 = x0 + dx*dsin(str0*d2r)
                    y1 = y0 + dx*dcos(str0*d2r)
                endif
                write(96,*) x1,y1,z
                dx = dx + incr
                goto 962
            endif
            z = z + incr
            goto 961
        endif
        close(96)
        RETURN
  
negative.code: C Parse command line
  C----
        call gcmdln(ffmf,fltf,magf,haff,staf,trgf,dspf,stnf,stsf,norf,
       1            shrf,coulf,gmtf,flttyp,auto,incr,adist,xy,prog,long)
  
  C----
  C Check for input fault source files, and whether output is specified
  