Geology Reference
In-Depth Information
CALL DIP1(U1,L,DU,X,DIP,R2,R3,Q2,Q3)
U(1)=U(1)-U1
C Find U(2).
CALL DIP2(U2,L,DL,X,DIP,R2,R3,Q2,Q3)
U(2)=U2
CALL DIP2(U2,-L,DU,X,DIP,R2,R3,Q2,Q3)
U(2)=U(2)+U2
CALL DIP2(U2,-L,DL,X,DIP,R2,R3,Q2,Q3)
U(2)=U(2)-U2
CALL DIP2(U2,L,DU,X,DIP,R2,R3,Q2,Q3)
U(2)=U(2)-U2
C Find U(3).
CALL DIP3(U3,L,DL,X,DIP,R2,R3,Q2,Q3)
U(3)=U3
CALL DIP3(U3,-L,DU,X,DIP,R2,R3,Q2,Q3)
U(3)=U(3)+U3
CALL DIP3(U3,-L,DL,X,DIP,R2,R3,Q2,Q3)
U(3)=U(3)-U3
CALL DIP3(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
SUBROUTINE DIP1(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)
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 fault plane co-ordinates XI2 and XI3.
XI2=XI*DCOS(DIP)
XI3=XI*DSIN(DIP)
C Evaluate U1.
U1=(X(2)-XI2)*DSIN(DIP)*(2.D0/R+4.D0/Q-4.D0*XI3*X(3)/Q**3
1 -3.D0/(Q+X(3)+XI3))-DCOS(DIP)*(3.D0*DLOG(Q+X(3)+XI3)
2 +2.D0*(X(3)-XI3)/R+4.D0*(X(3)-XI3)/Q
3 +4.D0*XI3*X(3)*(X(3)+XI3)/Q**3)
4 +3.D0*(DLOG(Q+X(3)+XI3)-DSIN(DIP)*DLOG(Q+Q3+XI))/DCOS(DIP)
5 +6.D0*X(3)*(DCOS(DIP)/Q-Q2*DSIN(DIP)/(Q*(Q+Q3+XI)))
RETURN
END
C
SUBROUTINE DIP2(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)
DOUBLE PRECISION K
Search WWH ::
Custom Search