ADAS Subroutine eiqip
SUBROUTINE EIQIP(EI,EIJ,EM,Z,PHI,SC,WI,WJ,R,EIQ,FLAG)
C
IMPLICIT REAL*8 (A-H,O-Z)
C
C-----------------------------------------------------------------------
C
C ****************** FORTRAN77 SUBROUTINE: EIQIP **********************
C
C PURPOSE UNKNOWN
C
C NOTES: THIS ROUTINE IS NOT YET PROPERLY ANNOTATED
C
C UNIX-IDL PORT:
C
C VERSION: 1.1 DATE: 16-1-96
C MODIFIED: TIM HAMMOND (TESSELLA SUPPORT SERVICES PLC)
C - FIRST VERSION
C
C-----------------------------------------------------------------------
C
FLAG=1.0
EF=EI-EIJ
T=EF/EI
C=1.0
IF(T-0.2)1,4,4
1 IF(Z-0.01)2,3,3
2 C=DSQRT(4.0*EF/EIJ)
GO TO 21
3 IF(EM-1.5)21,21,20
20 C=DABS(EIJ)/(EI*DSQRT(EF)+EF*DSQRT(EI))
C=Z*DSQRT(EM)*(1.1056/DSQRT(EIJ)-C)
C=DEXP(3.142*C)
21 EI=1.25*EIJ
EF=0.25*EIJ
GO TO 8
4 IF(T-5.0)8,8,5
5 IF(Z-0.01)6,6,7
6 C=DSQRT(-4.0*EI/EIJ)
GO TO 23
7 IF(EM-1.5)23,23,22
22 C=DABS(EIJ)/(EI*DSQRT(EF)+EF*DSQRT(EI))
C=Z*DSQRT(EM)*(1.1056/DSQRT(-EIJ)-C)
C=DEXP(3.142*C)
23 EI=-0.25*EIJ
EF=-1.25*EIJ
8 TI=DSQRT(EI)
TF=DSQRT(EF)
TIF=DABS(EIJ)/(TI+TF)
XI=Z*TIF/(TI*TF)
D=TIF*R
E=TI*TF
IF(EM-1.5)25,24,24
24 EM2=DSQRT(EM)
TIF=EM2*TIF
XI=-EM2*XI
D=EM2*D
25 CONTINUE
T=Z+E*R
T1=XIP(XI,D)/(T*T)
T2=4.0*PHI*E*SC*EM
IF(WI-WJ)10,10,9
9 T2=T2*WI/WJ
10 P=T2*T1
EIQW=8.0*C*PHI*YIP(XI,D)*EM
IF(P-0.5)15,15,11
11 R1=R
12 A=R1
FLAG=0.0
VA=0.5-P
R1=R1+R1
T=Z+E*R1
D=TIF*R1
P=T2*XIP(XI,D)/(T*T)
IF(P-0.5)13,14,12
13 B=R1
VB=0.5-P
CALL ZERO1(A,B,VA,VB,0.01D0,R1,XI,Z,E,TIF,T2)
14 T=Z+E*R1
D=TIF*R1
T1=XIP(XI,D)/(T*T)
15 EIQ=8.0*C*PHI*(YIP(XI,D)+0.5*(T*T-Z*Z)*T1)*EM
EIQ=DMIN1(EIQ,EIQW)
RETURN
END
REAL*8 EI, EIJ, EIQ, EM
REAL*8 FLAG, PHI, R, SC
REAL*8 WI, WJ, Z