anchor.diff:       1                          + dcos(str-pi/2.0d0)*dsin(rak)*dcos(dip)
         r(3) =                                         dsin(rak)*dsin(dip)
   C Project shear component on slip vector
  -      shear = s(1)*r(1) + s(2)*r(2) + s(3)*r(3)
  +      if (kshr.eq.2) then
  +          shear = dsqrt(s(1)*s(1) + s(2)*s(2) + s(3)*s(3))
  +      else
  +          shear = s(1)*r(1) + s(2)*r(2) + s(3)*r(3)
  +      endif
   C Coulomb stress (recall sign convention: pos = dilation)
         coul = shear + frc*norml
         RETURN
  
positive.code:           write(*,*)
       1 '    onto the rake vector of the target fault.'
            write(*,*)
            write(*,9999)
            write(*,*)
        endif
  
negative.code:       IMPLICIT NONE
        INTEGER kdsp,kstn,ksts,knor,kshr,kcou,kgmt
        COMMON /OCHECK/ kdsp,kstn,ksts,knor,kshr,kcou,kgmt
        if (kdsp.le.0.and.kstn.eq.0) then
            write(*,*) '!! Error: No output file specified'
            call usage('!! Use -disp, -strain, -stress, '//
       1                   '-normal, -shear, -coul, or an -auto option')
  