anchor.diff:  
   C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -C
   
  -      SUBROUTINE spaceconst(vp,vs,dens)
  -C----
  -C Calculate half-space constants.
  -C Reminder: in Poisson halfspace, vs = vp/sqrt(3), lamda = mu, a = 2/3
  +      SUBROUTINE hafspc(vp,vs,dens)
  +C----
  +C Calculate half-space constants (Lame parameters, etc.)
  +C
  +C Note: in Poisson half-space, vs = vp/sqrt(3), lamda = mu, a = 2/3
   C----
         IMPLICIT none
         REAL*8 vp,vs,dens,mu,lamda
  
positive.code: 
        SUBROUTINE moment1(Mss,Mds,slip,rakin)
  C----
  C Divide slip into strike slip and dip slip, and calculate moment of
  C each component (FINITE SOURCE).
  C----
        IMPLICIT none
        REAL*8 pi,d2r
  
negative.code:       R7 = R5*R2
  
        Rd = R + d
        Rd2 = Rd*Rd
        Rd3 = Rd2*Rd
        Rd4 = Rd3*Rd
  
  C p, q, s, t = distances in rotated reference frame.
  C p = fault dip-parallel distance
  C q = fault dip-perpendicular distance
        p = y*cd + d*sd
        q = y*sd - d*cd
        s = p*sd + q*cd
  