Geology Reference
In-Depth Information
DO 41 J=1,N1
C(I,J)=D(I,J,M)
41 CONTINUE
40 CONTINUE
C Put radius, density, mean density, ETA and e in
C active locations.
DO 42 I=1,N1
J=NK(M)+I
R(I)=RI(J)
RHO(I)=RHOI(J)
RHOB(I)=RHOBI(J)
ETA(I)=ETAI(J)
E(I)=EI(J)
42 CONTINUE
C Perform final estimate of Radau variable ETA.
C Set exponent in integrand.
EX=4.D0
C Initialize starting value of integral INT.
INT(1,M)=INTI(4,M)
C Form integrand function Q.
DO 43 I=1,N1
C Calculate Froude number or dimensionless parameter m(r).
FR=3.D0*WE*WE/(4.D0*PI*G*RHOB(I))
Q(I)=7.D0*FR*(1.D0+ETA(I))-3.D0*E(I)*(1.D0+ETA(I))**2
1 -4.D0*E(I)
Q(I)=Q(I)*2.D0*(1.D0-RHO(I)/RHOB(I))/35.D0
Q(I)=Q(I)+1.D0+ETA(I)/2.D0-ETA(I)*ETA(I)/10.D0
Q(I)=Q(I)*RHOB(I)/DSQRT(1.D0+ETA(I))
43 CONTINUE
C Perform quadrature.
CALL QUAD(INT,EX,Q,NIS,W,X,N1,C,R,M,M1)
C Store starting value for next region.
INTI(4,MP1)=INT(N1,M)
C Form final estimate of ETA and store in ETA1.
DO 44 I=1,N1
IF(I.EQ.1.AND.M.EQ.1)GO TO 45
ETA(I)=5.D0*INT(I,M)/(RHOB(I)*(R(I)**5))
ETA(I)=ETA(I)*ETA(I)-1.D0
GO TO 46
45 ETA(I)=0.D0
C Store ETA.
46 J=NK(M)+I
ETAI(J)=ETA(I)
44 CONTINUE
C Form final estimate of figure parameter e.
C Set exponent in integral.
EX=-1.D0
C Initialize starting value of integral INT.
INT(1,M)=INTI(5,M)
C Form integrand function Q.
DO 47 I=1,N1
Q(I)=ETA(I)
47 CONTINUE
C Perform quadrature.
CALL QUAD(INT,EX,Q,NIS,W,X,N1,C,R,M,M1)
C Store starting value for next region.
INTI(5,MP1)=INT(N1,M)
C Form final estimate of e(r)/e(0) and store in EI.
DO 48 I=1,N1
IF(I.EQ.1.AND.M.EQ.1)GO TO 49
E(I)=DEXP(INT(I,M))
Search WWH ::
Custom Search