Geology Reference
In-Depth Information
8 -4.D0*((X(3)+XI3)**2-XI3*X(3))/(Q*(Q+X(1)-XI1))-4.D0*XI3*X(3)*
9 (X(3)+XI3)**2*(2.D0*Q+X(1)-XI1)/(Q**3*(Q+X(1)-XI1)**2))
U3=U3+6.D0*X(3)*(DCOS(DIP)*DSIN(DIP)*
1 (2.D0*(Q3+XI)/(Q*(Q+X(1)-XI1))+(X(1)-XI1)/(Q*(Q+Q3+XI)))
2 -Q2*(DSIN(DIP)**2-DCOS(DIP)**2)/(Q*(Q+X(1)-XI1)))
RETURN
END
C
SUBROUTINE STRIKES(U,X,DIP,SLIP,DU,DL,L)
C
C This subroutine calculates the displacement vector U(3),
C at the field point X(3), for a strike-slip fault, dipping at angle DIP,
C with slip magnitude SLIP, upper boundary measured down-dip at DU,
C lower boundary down-dip at DL, fault extending along the strike
C from -L to L.
C
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
DIMENSION U(3),X(3)
DOUBLE PRECISION L
C Set value of pi.
PI=3.141592653589793D0
C Find parameters of field point distances normal and parallel to
C the fault and its image.
R2=X(2)*DSIN(DIP)-X(3)*DCOS(DIP)
R3=X(2)*DCOS(DIP)+X(3)*DSIN(DIP)
Q2=X(2)*DSIN(DIP)+X(3)*DCOS(DIP)
Q3=-X(2)*DCOS(DIP)+X(3)*DSIN(DIP)
C Find U(1).
CALL STRIKE1(U1,L,DL,X,DIP,R2,R3,Q2,Q3)
U(1)=U1
CALL STRIKE1(U1,-L,DU,X,DIP,R2,R3,Q2,Q3)
U(1)=U(1)+U1
CALL STRIKE1(U1,-L,DL,X,DIP,R2,R3,Q2,Q3)
U(1)=U(1)-U1
CALL STRIKE1(U1,L,DU,X,DIP,R2,R3,Q2,Q3)
U(1)=U(1)-U1
C Find U(2).
CALL STRIKE2(U2,L,DL,X,DIP,R2,R3,Q2,Q3)
U(2)=U2
CALL STRIKE2(U2,-L,DU,X,DIP,R2,R3,Q2,Q3)
U(2)=U(2)+U2
CALL STRIKE2(U2,-L,DL,X,DIP,R2,R3,Q2,Q3)
U(2)=U(2)-U2
CALL STRIKE2(U2,L,DU,X,DIP,R2,R3,Q2,Q3)
U(2)=U(2)-U2
C Find U(3).
CALL STRIKE3(U3,L,DL,X,DIP,R2,R3,Q2,Q3)
U(3)=U3
CALL STRIKE3(U3,-L,DU,X,DIP,R2,R3,Q2,Q3)
U(3)=U(3)+U3
CALL STRIKE3(U3,-L,DL,X,DIP,R2,R3,Q2,Q3)
U(3)=U(3)-U3
CALL STRIKE3(U3,L,DU,X,DIP,R2,R3,Q2,Q3)
U(3)=U(3)-U3
C Scale displacements.
U(1)=SLIP*U(1)/(12.D0*PI)
U(2)=SLIP*U(2)/(12.D0*PI)
U(3)=SLIP*U(3)/(12.D0*PI)
RETURN
END
C
Search WWH ::
Custom Search