anchor.diff:       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
  +          open(unit=61,file='mu.dat',status='old')
  +          read(61,*) mu
  +          close(61)
  +      else
  +          mu = 4.3d10
  +      endif
  +C      print *,'SHEAR MOD:',mu
         open (unit=33,file=magf,status='old')
         f = 0
     331 f = f + 1
  
positive.code:       PARAMETER (FMAX=1500)
        REAL*8 evlo(FMAX),evla(FMAX),evdp(FMAX),str(FMAX),dip(FMAX),
       1       rak(FMAX),dx(FMAX),dy(FMAX),slip(FMAX)
        REAL*8 uN,uE,uZ,strain(3,3),stress(3,3),norml,shear,coul
        REAL*8 az,hdsp
        INTEGER kffm,kflt,kmag,ksta,khaf,ktrg
        INTEGER kdsp,kstn,ksts,knor,kshr,kcou,kgmt
  
negative.code:               elseif (noz.eq.1.and.kdsp.eq.2) then
                    az = atan2(uE,uN)*r2d
                    hdsp = dsqrt(uE*uE+uN*uN)
                    write(21,9994) stlo,stla,az,hdsp,uz
                endif
            endif
  C         Compute (3x3) strain matrix
  