anchor.diff:        RETURN
         END
   
  -
   C----------------------------------------------------------------------C
  -C Source size (moment) subroutines                                     C
  +C------------------- Source Moment Subroutines ------------------------C
   C----------------------------------------------------------------------C
   
         SUBROUTINE moment0(Mss,Mds,slip,rakin,area)
   C----
  -C Divide slip into strike slip and dip slip, and calculate moment of
  +C Divide slip into strike-slip and dip-slip, and calculate moment of
   C each component (POINT SOURCE).
  +C
  +C Note: These values are not strictly seismic moments.
   C----
         IMPLICIT none
         REAL*8 pi,d2r
  
positive.code:       R7 = R5*R2
  
        Rd = R + d
  
        A3 = 1.0d0 - 3.0d0*x*x/R2
        A5 = 1.0d0 - 5.0d0*x*x/R2
        A7 = 1.0d0 - 7.0d0*x*x/R2
        B3 = 1.0d0 - 3.0d0*y*y/R2
        B5 = 1.0d0 - 5.0d0*y*y/R2 
        B7 = 1.0d0 - 7.0d0*y*y/R2
        C3 = 1.0d0 - 3.0d0*d*d/R2
        C5 = 1.0d0 - 5.0d0*d*d/R2
        C7 = 1.0d0 - 7.0d0*d*d/R2
  
  C----
  C p, q, s, t = distances in rotated reference frame.
  C p = fault dip-parallel distance
  C q = fault dip-perpendicular distance
  C----
        p = y*cd + d*sd
        q = y*sd - d*cd
        s = p*sd + q*cd
  
negative.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
  C Note: These values are not strictly seismic moments.
  C----
        IMPLICIT none
        REAL*8 pi,d2r
  