COMPUTER PROGRAM USER'S MANUAL FOR ... - CAFE Foundation

FOHTRAN IV G LEVEL 20.1 MAIN
0 135 NT S= 10
0136 TI PSDC4 1) =700.
0137
0138
TI PSPO=700.
GO TO 320 .
0139
0140
310 TlPSPD=TS(fC)-DTS(IC)
NTS=NDTS I IC 1
0141 320 DO 600 ITS=l*NTS
0 142 TIPSPD=TlPSPD+DTSlIC)
OATE = 72034 10/08/04
C MACH NUMBER CALCULATION AND AUVANCE RATIO J
0143
0 144
LHS~l~=.001512*VKTASlIC~*FC~IC~
ZMS(ZI=TIPSPD*fC( IC)/1120.
0 145 ZMl=ZMSI 1)
0146
0147
0148
ZJI=5.309*VKTASllC)/TlPSPD
340
IFlZJI.EQ.0.) ZHL=ZMS(Z)
IF ISTALITIIC).LE..SO.AND.ZJI.LE.5.0) GO TO 342
0149 IF(STALITLIC).GT..50.AND.ZJI.L€.3.0) GO TO 342
0150 WRITE(6p341) ZJI
0151 341 FOEMAT(' ADVANCE RATIO TOO HIGH = ' 9 F8.4)
0152 GO TO 600
0 153 342 CONTlNUE
C ITERAIION ON CT OR CP TO GET 50 PFRCENT STALL TIPSPEED
0154
0155
IF IN=O
IF (STALIT(ICJ.LE..501 GO TO 399
01 56
IUSV=IW
0157
0158
I u=3
CALL PERFY ~3~CP~ZJI~AFT~BLAOT~CLl~Cl~Z~S~77lO~
0159
I d=I HSV
0160
IF(Ik.EQ.2) GO TU 712
0161 711 BHPClITSI=2.0*TIPSDG~ITS)**3*DIA*~2~69hb.~CP/~l~.El~~R~R~~IC~l
0 162 IF(ABSIBHPlIC)-BHPG(lTS)).GE..OO5*BHP(IC~) GO TO 705
0163 THHUSTIIC~=CT*TIPSPD**2*~IA~~2/~1.515E06*RORO~IC~~~~64.76
0164
T4 1C=l.
0165
0 166
0 167
GO TO 720
705 IF(ITS.EQ.11 GO TO 7000
TIPSDG~ITS+1~~~ALOG~BHPlIC~~-ALOGlBHPGllTS-l~
))*(TIPSDG(lTS)-
1TIPSDGllTS-1l~/~ALOG~BHPG(ITS~~-ALOG~~HPG~lTS-l~~~+TlPSDG~ITS-l~
0 168 GO TO 709
0 169 7000 TIPSDG121=400.
0170 TIPSPO=TIPSDGIITS+l)
0171 GO TO 600
0172 712 THRSTG(ITS)=TIPSDG(ITS)**2*DlA**2*364.76*CT/l 1.515E06*RDRO(lC)~
0173 IFlABS(THRUST(ICI-THRSTG(iTS~~.~E..O~5*THRUST~IC)~ GO TO722
0174 TIPSPD=TIPSCG(ITS)
0175 BHPlICI=CP*2.0*TIPSP0**7*DIA**2/(10.E10*ROR~~IC~I*6966.
01 76
TRIG=l.
0177
GO TO 720
0178 722 lFlITS.EQ.1) GO TO 7000
0179 TIPSDG~ITS+L)=(ALOGITHRUST(IC))-ALOG(THRSTG(ITS-~~I~*(TIPS~G(ITS~oleo
lTIPSDG~ITS-1~I/~ALOG~THRSTG~ITSI~-ALOG~THRSTGll~S-lJ~~+TIPSDG
2( I is-I 1
709 TIPSPD=TIPSOG(ITS+l)
0101 1FlNTS.NE.ITS) GO TO 600
0182
0183
WRITE I 615981
598 FORMAT {//' FAILED STALL ITERATION '// I
0184 GO TO 700
C END OF T IPSPD ITERATION 50 PERCENT STALL
C CALCULATION OF REQUIRED CP OR CT
0185 399 IFllU-1)400~400~430
FIGURE 3A. FORTRAN IV L-ISTING (CONTINUED)
51