->

(0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) ( 35 ) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (35)

goto 36i 361 irtNN.ne.dgoto 360

urite(6,3U3)

363 f0rmat(3x, ».-/3x.

» ""10.08-0.111 0.15-0.2510.35-0.4110.51-0.93!cb.l. /

------------------------------------------------------ ,/>

« ! 1.3 ! 1.25,, ! 1.2 1 1.1 ! 1.05 )

read(6.4)

l..l=a»xn»dopt

366 360 tt=xn»nOPT«»2/ll

urite<6.362) 342 porkat( ",12.5!

urite(6.15)

read(6,9)k1

if(n1.E0.1)g0t0 100 366 write(6r364)ll

264 format( pabha"cm"e12.5) write(6>15) reaii<6.9)k2 if<n2.E0.1)gdt0 100 call OnQ<doft,xn,f,d,xl>q) mrite(6,365) 365 f0rmat<3x.b :>

write<6,368)tt .368 format (. ,e12.5, "cm")

write(6,367)xnrd0pt.ll>q 7 format(3x,HEOFX0aHM0E .e12.5,/, « d -e12.5,/, « "",/12.5,/, / ,12.5) 3 format(, )

4 format(f10.5) 5 format (, 1-") £55 stop-end



SUBROUTINE EKRAN(UD,XL,XL1) URITE(6>200)

200 FDRHT<5X. i/5X> « D i/5Xr [l3KF/D=1.65...2.5./5Xi » D3KP/Iel.45...2.2>/5Xi » [l3KP/n=2.S 1/SXr « /) READ (6,202) DDEI

202 F0RMAT<F6.3) [F(L[l.GT.1.0)G0 205 ZN=0.i+0.9»LB

SO TO 203 205 ZN=1.4

203 XK=ZN»((1/D[ID)»«3) XH=XL/(1-XIV) IJRITE<6,222)XL1

222 rORMAT(3X,BAM L ./. « =fE12.3rmTH) RETURN END

SUBROUTINE IlOPTl<tlrF,XNiDOPT) MRITE<6>103)

103 FORMAT(3X. ». "*(.199))

READ(6,104) 102 ZZ=0.106«SaRT<F«1000000.> IO* FORMAT(F6.3)

PSI=(<XK»XN)/(2»ZZ«D))«»2

IF(PSI.GT.1.0)GD TO 100

Z0PT=3.2-1.2»PS1

GO TO 110

100 IF(PSI.GT.10.0)G0 TO 101 Z0FI=2.1-0.09«PSI

GO TO no

101 IF<PSI.GT.100.0)00 TO 105 ZOPT=1.2S-0.0043«PSI

GO TO no 105 Z0PT=0.9-0.00044«PSI 110 DOFT=ZOPT/ZZ

URrTEt6>lU)D0Pr 11; FORMAT(3Xf =,7.3)

RETURN



subroutine oqq(doptixn,f>d,xl.q)

real lp

write,<a.29)

29 f0r«at(3x, 1« ) read(4.30)roo

30 (.) lp=3.14»d»xn«0.01 r0=lp»roo WRnE(*f31)

31 fdrmat<3x, . "(.199)) read(6>32)xk

32 f0rhat(f6.3) z=0.106«sqrt<fi(1000ooo.) write(6.33)z

33 f0rmat(3x, . (2)" z=> f4.3,/, » *. "(2)" z (.200))

read(6,34)f0 read(634)

34 format(f6.3) mrite(6,35)

35 format(. ?>) read(6,36>it

36 format(il) ifdt.eq.dgoto 37 urite(6,38)

38 f0rmat(3x, too ) read(6,39)tgd

39 f0rmat(f6.3> write(6,40>

40 format(3x, xe » » )

read(6,41>xe

41 f0rhat(f6.3) urite(6>42)

n format(3x, "" . «. (. .197)) read(6,43)b

43 f0rmat(f6.3) write(6.44)

44 format («. xkl (xk1=f (/) .. 195)) read(6,4s)xkl

45 f0rmat(f6.3> urite(6.46)

46 f0rmat(3x, xk2 (xk2=f(ll/d) ..195)) rean(6.47)xk2

47 f0rmat(f6.3)

rd=(0.25«tgd»xk1»xk2«d«xl»xl«(f*»3))/((1.+1./(b»xe))«1000.) goto 49

37 wr1te(6.50)

50 format(3x. tgd ) reab(6.51)tgd

51 format(f6.3) write(6.52)

52 f0rmat(3x. ck :/ «. (.197);/

«. (=15.,.30) f «. ) read(6.53)ck 53f0rmat(f6.3)

rd=(0.25«tgd«0.9«ck«xl»xl«(f«»3))/1000, 49 rf=r0» (f0+ (xk4fxn»d0pt/ (2»d)) »»2»g) q=6.28«f»xl/(r0+rf+rd) urite(6.54)q

54 format(3x,aobf;qth0ctb 3kpftha=siFiOiS)

return £



(0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) ( 35 ) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68)