Geology Reference
In-Depth Information
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 projections of Q on plane X(1)=0 and on plane Q3=0.
H=Q2*Q2+(Q3+XI)*(Q3+XI)
H=DSQRT(H)
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 U2.
U2=DSIN(DIP)*(-DLOG(R+X(1)-XI1)+DLOG(Q+X(1)-XI1)
1 +4.D0*XI3*X(3)/(Q*(Q+X(1)-XI1))+3.D0*(X(1)-XI1)/(Q+X(3)+XI3)
2 +(X(2)-XI2)**2*(2.D0/(R*(R+X(1)-XI1))+4.D0/(Q*(Q+X(1)-XI1))
3 -4.D0*XI3*X(3)*(2.D0*Q+X(1)-XI1)/(Q**3*(Q+X(1)-XI1)**2)))
4 -DCOS(DIP)*((X(2)-XI2)*(2.D0*(X(3)-XI3)/(R*(R+X(1)-XI1))
5 +4.D0*(X(3)-XI3)/(Q*(Q+X(1)-XI1))+4.D0*XI3*X(3)*
6 (X(3)+XI3)*(2.D0*Q+X(1)-XI1)/(Q**3*(Q+X(1)-XI1)**2))
7 +6.D0*DATAN((X(1)-XI1)*(X(2)-XI2)/((H+X(3)+XI3)*(Q+H)))
8 -3.D0*DATAN((X(1)-XI1)*(R3-XI)/(R2*R))
9 +6.D0*DATAN((X(1)-XI1)*(Q3+XI)/(Q2*Q)))
U2=U2+6.D0*DATAN(((K-Q2*DCOS(DIP))*(Q-K)+(Q3+XI)*K*DSIN(DIP))/
1 ((X(1)-XI1)*(Q3+XI)*DCOS(DIP)))/DCOS(DIP)
2 +6.D0*X(3)*(((DSIN(DIP)**2-DCOS(DIP)**2)*(Q3+XI)
3 +2.D0*Q2*DCOS(DIP)*DSIN(DIP))/(Q*(Q+X(1)-XI1))
4 +(X(1)-XI1)*DSIN(DIP)**2/(Q*(Q+Q3+XI)))
RETURN
END
C
SUBROUTINE DIP3(U3,XI1,XI,X,DIP,R2,R3,Q2,Q3)
C
C This subroutine evaluates the three-component, U3,
C of the double integral for specified limits XI1 and XI.
C
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
DIMENSION X(3)
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 X(1)=0.
H=Q2*Q2+(Q3+XI)*(Q3+XI)
H=DSQRT(H)
C Find fault plane co-ordinates XI2 and XI3.
XI2=XI*DCOS(DIP)
XI3=XI*DSIN(DIP)
C Evaluate U3.
U3=DSIN(DIP)*((X(2)-XI2)*(2.D0*(X(3)-XI3)/(R*(R+X(1)-XI3))
1 +4.D0*(X(3)-XI3)/(Q*(Q+X(1)-XI1))-4.D0*XI3*X(3)*(X(3)+XI3)*
2 (2.D0*Q+X(1)-XI1)/(Q**3*(Q+X(1)-XI1)**2))
3 -6.D0*DATAN((X(1)-XI1)*(X(2)-XI2)/((H+X(3)+XI3)*(Q+H)))
4 +3.D0*DATAN((X(1)-XI1)*(R3-XI)/(R2*R))
5 -6.D0*DATAN((X(1)-XI1)*(Q3+XI)/(Q2*Q)))
6 +DCOS(DIP)*(DLOG(R+X(1)-XI1)-DLOG(Q+X(1)-XI1)
7 -2.D0*(X(3)-XI3)**2/(R*(R+X(1)-XI1))
Search WWH ::
Custom Search