anchor.diff:  
         SUBROUTINE gcmdln(ffmf,fltf,magf,haff,staf,trgf,dspf,stnf,stsf,
        1                  norf,shrf,coulf,gmtf,flttyp,auto,incr,adist,xy,
  -     2                  prog,long)
  +     2                  prog,long,emprel)
         IMPLICIT NONE
         CHARACTER*40 tag
         INTEGER narg,i,ptr
         CHARACTER*40 ffmf,fltf,magf,haff,staf,trgf,
        1             dspf,stnf,stsf,norf,shrf,coulf,
  -     2             gmtf
  +     2             gmtf,emprel
         INTEGER flttyp,auto,xy,prog,long
         REAL*8 incr,adist
         INTEGER kffm,kflt,kmag,ksta,khaf,ktrg
  
positive.code: C----------------------------------------------------------------------C
  
        SUBROUTINE readmag(magf,evlo,evla,evdp,str,dip,rak,dx,dy,
       1                   slip,nflt)
  C----
  C Read shear dislocations from fault file in format:
  C   evlo evla evdp(km) str dip rak mag
  C All units are converted to SI, angles are in degrees.
  C----
        IMPLICIT NONE
        CHARACTER*40 magf
        INTEGER i,f,nflt,FMAX,ln
        PARAMETER (FMAX=1500)
        REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
       1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX),mag
        REAL*8 wid,len,mu
        INTEGER typ
        LOGICAL ex
        mu = 4.3d10 ! shear mod hard-coded, corresponds to autohaf params
        inquire(file='mu.dat',EXIST=ex)
        if (ex) then
  
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,
       2            emprel)
  
  C----
  C Check for input fault source files, and whether output is specified
  