Subroutine STG (AL1,B2,E,G,TN,PO,TO,RT,CP,D1K,RK,BB1,BB2,FDU,DEL,TM,POI,TOI,HO,UL1,UL2,RC,ROI,PN,C2,AL2,DL,DLL)
PI=3.141593
POL1=1.
FI=0.96
PSI=0.940
AM1=0.97
AM2=0.950
RC=RK
2 UKCF=FI*COS(AL1)/(2.*SQRT(1.-RK))
CF=PI*D1K*TN/UKCF
HO=CF**2/2
T1T=TO-1(1.-RC)*HO/CP
TYPE 3,PO,TO,TM,RT
3 FORMAT (4F14.4)
CALL GAS(PO,TO,TM,RT,T1T,P1,V1T)
TYPE 3,P1,V1T,RC,G
A=G*V1T/(AM1*CF*PI*SIN(AL1)*SQRT(1.-RC))
TYPE 3,A,D1K,AL1,AL1
UL1=(-D1K+SQRT(D1K**2+4.*A))/2
RCR=1.-(1.-RK)*(((D1K+UL1)/D1K)**(-1.8))
TYPE 4,UL1,RCR
4 FORMAT (2F10.4)
FI=0.980-0.009*BB1/UL1
AM1=0.982-0.005*BB1/UL1
RR=RCR-RC
RC=RCR
TYPE 4,RR,RK
IF (ABS(RR)-0.0005) 1,1,2
CONTINUE
F1=PI*(D1K+UL1)*UL1*SIN(AL1)
C1=FI*CF*SQRT(1.-RC)
U1=PI*(D1K+UL1)*TN
TYPE 3,F1,C1,U1,AL1
CALL TRK(C1,U1,AL1,C1U,W1,B1)
TYPE 3,C1U,W1,B1,B1
W2T=SQRT(W1**2+RC*CF**2)
T2T=TO-HO/CP
TYPE 4,W2T,T2T
CALL GAS (PO,TO,TM,RT,T2T,P2,V2T)
TYPE 4,P2,V2T
B=G*V2T/(AM2*PI*PSI*SIN(B2)*W2T)
TYPE 4,B,B
D2K=D1K-0.001
UL2=(-D2K+SQRT(D2K**2+4.*B))/2.
AM2=0.965-0.015*BB2/UL2
PSI=0.957-0.011*BB2/UL2
TYPE 3,UL2,AM2,PSI,D2K
W2=PSI*W2T
U2=PI*(D2K+UL2)*TN
CALL TRK (W2,U2,B2,W2U,C2,AL2)
TYPE 3,W2,U2,B2,W2U
TYPE 4,C2,AL2
POL=U2*(C1U+W2U-U2)/HO
TYPE 3,POL,POL1
PP=POL1-POL
POL1=POL
TYPE 4,PP,PP
IF (ABS(PP)-0.001) 5,5,1
5 U2CF=UKCF*(D2K+UL2)/D2K
ATP=0.0007*(D2K+UL2)*(U2CF**3)/(PI*UL2*SIN(B2))
ADU=FDU*POL/(AM1*F1)
RP=1.-(1.-RK)*((1.+2*UL2/D2K)**(-1.8))
TYPE 3,U2CF,ATP,ADU,RP
APU=PI*(D2K+2.*UL2)*DEL*POL*SQRT(RP/(1.-RC))/(AM1*F1)
TYPE 4,APU,PP
PAUSE
ROI=POL-ATP-ADU-APU
PN=G*HO*ROI/1000.
TOJ=TO-HO*ROI/CP
TOOJ=TOJ+(C2**2)/(2.*CP)
CALL GAS (P2,TOJ,TM,RT,TOOJ,POOJ,V)
POI=POOJ
TOI=TOOJ
DL=UL2-UL1
DLL=(D2K+UL2)/UL2
RETURN
END