anchor.diff:        PARAMETER (FMAX=150000)
   
   C---- Input files
  -      CHARACTER*40 ffmf      ! faults, finite fault model format
  +      CHARACTER*40 ffmf      ! faults, static out finite fault format
  +      CHARACTER*40 fspf      ! faults, SRCMOD finite fault format
         CHARACTER*40 fltf      ! faults, lon-lat-dep-str-dip-rak-slip-wid-len format
         CHARACTER*40 magf      ! faults, lon-lat-dep-str-dip-rak-mag format
         CHARACTER*40 haff      ! half-space parameters
  
positive.code:       CHARACTER*40 haff,staf,trgf
        LOGICAL ex
        INTEGER nsta,ntrg,ptr,i
        INTEGER kffm,kflt,kmag,ksta,khaf,ktrg,kvol
        INTEGER kdsp,kstn,ksts,knor,kshr,kcou,kgmt
        COMMON /ICHECK/ kffm,kflt,kmag,ksta,khaf,ktrg,kvol
        COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
  C Look for station file, count number of receivers
        if (ksta.eq.0) then
  
negative.code:       LOGICAL ex
        CHARACTER*40 haff
        REAL*8 vp,vs,dens
        INTEGER kffm,kfsp,kflt,kmag,ksta,khaf,ktrg,kvol
        COMMON /ICHECK/ kffm,kfsp,kflt,kmag,ksta,khaf,ktrg,kvol
  
        ! Determine empirical relation used
        p = index(emprel,'p')
  