anchor.diff:        strain(3,1) = strain(1,3)
         strain(3,2) = strain(2,3)
   
  -      RETURN 
  +      RETURN
         END
   
   C----------------------------------------------------------------------C
   
  +C      SUBROUTINE FNDISP(ux,uy,uz,x,y,stdp,evdp,dipin,rakin,wid,len,
  +C     1                   slip,vp,vs,dens)
         SUBROUTINE o92rect(ux,uy,uz,x,y,stdp,evdp,dipin,rakin,wid,len,
        1                   slip,vp,vs,dens)
  -C      SUBROUTINE DSPFN(UX,UY,UZ,X,Y,STDP,EVDP,DIPIN,RAKIN,WID,LEN,SLIP,VP,VS,DENS)
   C----
   C Static elastic displacement at a location (internal or surface) due
  -C to a finite rectangular source shear dislocation in a uniform, 
  +C to a finite rectangular source shear dislocation in a uniform,
   C isotropic halfspace.
   C
   C INPUTS (ALL IN SI UNITS):
  
positive.code:           f(4,1) =                CA2*3.0d0*x*pq/R5
            f(5,1) =  CA1*s/R3    + CA2*3.0d0*y*pq/R5
            f(6,1) = -CA1*t/R3    + CA2*3.0d0*d*pq/R5
  
            f(1,2) = -3.0d0*xx*q/R5 - CB*I1*sd
            f(2,2) = -3.0d0*xy*q/R5 - CB*I2*sd
  
negative.code:           fx(4,1) =                     CA2*3.0d0*pq*A5/R5
            fx(5,1) = -CA1*3.0d0*x*s/R5 - CA2*15.0d0*xy*pq/R7
            fx(6,1) =  CA1*3.0d0*x*t/R5 - CA2*15.0d0*x*d*pq/R7
  C          fx(7,1) =  CA1*A3/R3        - CA2*3.0d0*q*q*A5/R5
  C          fx(8,1) = -CA1*3.0d0*x*t/R5 + CA2*15.0d0*xy*q*q/R7
  C          fx(9,1) = -CA1*3.0d0*x*s/R5 + CA2*15.0d0*xd*q*q/R7
  C          fx(10,1) = -CA1*A3/R3
  C          fx(11,1) =  CA1*3.0d0*xy/R5
  C          fx(12,1) =  CA1*3.0d0*xd/R5
        endif
  
        RETURN
  