Geology Reference
In-Depth Information
SUBROUTINE STRIKE1(U1,XI1,XI,X,DIP,R2,R3,Q2,Q3)
C
C This subroutine evaluates the one-component, U1,
C of the double integral for specified limits XI1 and XI.
C
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
DIMENSION X(3)
DOUBLE PRECISION K
C Find distance R from fault to field point and
C distance Q from image fault to field point.
R=(X(1)-XI1)*(X(1)-XI1)+R2*R2+(R3-XI)*(R3-XI)
R=DSQRT(R)
Q=(X(1)-XI1)*(X(1)-XI1)+Q2*Q2+(Q3+XI)*(Q3+XI)
Q=DSQRT(Q)
C Find projection of Q on plane Q3=0.
K=(X(1)-XI1)*(X(1)-XI1)+Q2*Q2
K=DSQRT(K)
C Find fault plane co-ordinates XI2 and XI3.
XI2=XI*DCOS(DIP)
XI3=XI*DSIN(DIP)
C Evaluate U1.
U1=(X(1)-XI1)*(2.D0*R2/(R*(R+R3-XI))
1 -(4.D0*Q2-2.D0*X(3)*DCOS(DIP))/(Q*(Q+Q3+XI))
2 -3.D0*DTAN(DIP)/(Q+X(3)+XI3)+4.D0*Q2*X(3)*DSIN(DIP)/Q**3
3 -4.D0*Q2*Q3*X(3)*DSIN(DIP)*(2.D0*Q+Q3+XI)/(Q**3*(Q+Q3+XI)**2))
4 -6.D0*DTAN(DIP)**2*DATAN(((K-Q2*DCOS(DIP))*(Q-K)
5 +(Q3+XI)*K*DSIN(DIP))/((X(1)-XI1)*(Q3+XI)*DCOS(DIP)))
6 +3.D0*DATAN((X(1)-XI1)*(R3-XI)/(R2*R))
7 -3.D0*DATAN((X(1)-XI1)*(Q3+XI)/(Q2*Q))
RETURN
END
C
SUBROUTINE STRIKE2(U2,XI1,XI,X,DIP,R2,R3,Q2,Q3)
C
C This subroutine evaluates the two-component, U2,
C of the double integral for specified limits XI1 and XI.
C
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
DIMENSION X(3)
C Find fault to field point distance R, and image fault to field point
C distance Q.
R=(X(1)-XI1)*(X(1)-XI1)+R2*R2+(R3-XI)*(R3-XI)
R=DSQRT(R)
Q=(X(1)-XI1)*(X(1)-XI1)+Q2*Q2+(Q3+XI)*(Q3+XI)
Q=DSQRT(Q)
C Find fault plane co-ordinates XI2 and XI3.
XI2=XI*DCOS(DIP)
XI3=XI*DSIN(DIP)
C Evaluate U2.
U2=DSIN(DIP)*(3.D0*DTAN(DIP)*DLOG(Q+X(3)+XI3)/(DCOS(DIP))
1 -DLOG(R+R3-XI)-(1.D0+3.D0*DTAN(DIP)**2)*DLOG(Q+Q3+XI))
2 +2.D0*R2*R2*DSIN(DIP)/(R*(R+R3-XI))+2.D0*R2*DCOS(DIP)/R
3 -2.D0*DSIN(DIP)*(2.D0*X(3)*(Q2*DCOS(DIP)-Q3*DSIN(DIP))
4 +Q2*(Q2+X(2)*DSIN(DIP)))/(Q*(Q+Q3+XI))
5 -3.D0*DTAN(DIP)*(X(2)-XI3)/(Q+X(3)+XI3)
6 +2.D0*(Q2*DCOS(DIP)-Q3*DSIN(DIP)-X(3)*DSIN(DIP)**2)/Q
7 +4.D0*Q2*X(3)*DSIN(DIP)*(X(2)-XI2+Q3*DCOS(DIP))/Q**3
8 -4.D0*Q2*Q2*Q3*X(3)*DSIN(DIP)**2*(2.D0*Q+Q3+XI)/
9 (Q**3*(Q+Q3+XI)**2)
RETURN
END
Search WWH ::
Custom Search