10 'Program to calculate Steam Properties
20 'Ref : Program computes turbine steam rates and properties; Hydrocarbon Proces   sing,Nov 1988,105-108
30 DIM A(12),B(12),C(12),D(12),E(12),F(12),G(12),PROP(12)
40 FOR I=1 TO 12:READ A(I),B(I),C(I),D(I),E(I),F(I),G(I):NEXT
50 DATA -.17724,3.83986,11.48345,31.1311,8.762969e-5,-2.78794e-8,86.594,-5.280126e-7,2.99461e-5,1.521874e-4,6.62512e-5,8.408856e-10,1.86401e-14,.01596,-.48799,304.717614,9.8299035,-16.455274,9.474745e-4,-1.363366e-6,19.53953
60 DATA 2.662e-3,457.5802,-.176959,.826862,-4.601876e-7,6.3181e-11,-2.3928
70 DATA -.15115567,3.671404,11.622558,30.832667,8.74117e-5,-2.62306e-8,54.55
80 DATA .008676153,-1.3049844,-8.2137368,-16.37649,-4.3043e-5,9.763e-9,1045.81
90 DATA -.14129,2.258225,3.4014802,14.438078,4.222624e-5,-1.569916e-8,1100.5
100 DATA -1.67772e-4,4.272688e-3,.01048048,.05801509,9.101291e-8,-2.7592e-11,.11801
110 DATA 3.454439e-5,-2.75287e-3,-7.33044e-3,-.14263733,-3.49366e-8,7.433711e-12,1.85565
120 DATA -1.476933e-4,1.2617946e-3,3.44201e-3,-.08494128,6.89138e-8,-2.4941e-11,1.97364
130 DATA -.1549439,3.662121,11.632628,30.82137,8.76248e-5,-2.646533e-8,54.56
140 DATA -.0993951,1.93961,2.428354,10.9818864,2.737201e-5,-1.057475e-8,1040.03
150 KEY OFF:SCREEN 0,0,0:CLS
160 COLOR 10
170 LOCATE 2,25:PRINT "Program to calculate Steam Properties"
180 LOCATE 5,15:INPUT "Enter the Steam Pressure(Kg/cm2 G)  :";P
190 P1=(P+1)*14.7
200 SATTEMP=A(1)*P1+B(1)/P1+C(1)*(P1^.5)+D(1)*LOG(P1)+E(1)*P1*P1+F(1)*(P1^3)+G(1)
210 SATTEMP1=(SATTEMP-32)/1.8  'øC
220 LOCATE 7,20:COLOR 12:PRINT "Saturated temp. of this steam is :";SATTEMP1;"øC"
230 COLOR 10:LOCATE 9,15:INPUT "Whether your steam is saturated or not(y/n) ";ANS$
240 IF ANS$="y" OR ANS$="Y" THEN 270
250 IF ANS$="n" OR ANS$="N" THEN 550
260 GOTO 230
270 FOR I=2 TO 12
280 PROP(I)=A(I)*P1+B(I)/P1+C(I)*(P1^.5)+D(I)*LOG(P1)+E(I)*P1*P1+F(I)*(P1^3)+G(I)
290 NEXT
300 SPVOLLIQ=PROP(2):SPVOLVAPUPTO200=PROP(3):SPVOLVAP200TO500=PROP(4)  'ft3/lb
310 ENTHLIQ=PROP(5):VAPNENTH=PROP(6):ENTHVAP=PROP(7)  'Btu/lb
320 ENTRLIQ=PROP(8):VAPNENTR=PROP(9):ENTRVAP=PROP(10)  'Btu/lb/øR
330 INTENERLIQ=PROP(11):INTENERVAP=PROP(12)   'Btu/lb
340 'conversion of units
350 SPVOLLIQ1=SPVOLLIQ*6.243386E-02:SPVOLVAPUPTO2001=SPVOLVAPUPTO200*6.243386E-02::SPVOLVAP200TO5001=SPVOLVAP200TO500*6.243386E-02  'm3/Kg
360 ENTHLIQ1=ENTHLIQ*.5556:VAPNENTH1=VAPNENTH*.5556:ENTHVAP1=ENTHVAP*.5556  'Kcal/kg
370 ENTRLIQ1=ENTRLIQ:VAPNENTR1=VAPNENTR:ENTRVAP1=ENTRVAP  'Kcal/Kg/øC
380 INTENERLIQ1=INTENERLIQ*.5556:INTENERVAP1=INTENERVAP*.5556   'Kcal/Kg
390 CLS
400 LOCATE 1,25:PRINT "Program to calculate Steam Properties"
410 LOCATE 3,10:PRINT "Steam Pressure         (Kg/cm2 g) :";P
420 LOCATE 4,10:PRINT "Steam Saturation temperature (øC) :";SATTEMP1
430 LOCATE 5,10:PRINT "Liquid Specific Volume   (m3/kg)  :";SPVOLLIQ1
440 IF P1<200 THEN DUMMY=SPVOLVAPUPTO2001 ELSE DUMMY=SPVOLVAP200TO5001
450 LOCATE 6,10:PRINT "Vapour Specific Volume  (m3/kg)   :";DUMMY
460 LOCATE 7,10:PRINT "Liquid Enthalpy         (Kcal/Kg) :";ENTHLIQ1
470 LOCATE 8,10:PRINT "Vaporization Enthalpy   (Kcal/Kg) :";VAPNENTH1
480 LOCATE 9,10:PRINT "Vapour Enthalpy         (Kcal/Kg) :";ENTHVAP1
490 LOCATE 10,10:PRINT "Liquid Entropy        (Kcal/Kg/øC):";ENTRLIQ1
500 LOCATE 11,10:PRINT "Vaporization Entropy (Kcal/Kg/øC) :";VAPNENTR1
510 LOCATE 12,10:PRINT "Vapour Entropy       (Kcal/Kg/øC) :";ENTRVAP1
520 LOCATE 13,10:PRINT "Liquid Internal Energy (Kcal/Kg)  :";INTENERLIQ1
530 LOCATE 14,10:PRINT "Vapour Internal Energy (Kcal/Kg)  :";INTENERVAP1
540 GOTO 850
550 'routine for superheated steam
560 LOCATE 10,15:INPUT "Enter the temp. of your superheated steam(øC) :";T
570 P1=P+1 'atm
580 T1=T+273 'øK
590 C1=80870!/T1/T1
600 C2=(10^C1)-(2641.62/T1)
610 C3=1.89+C2
620 C4=C3*(P1^2)/(T1^2)
630 C5=372420!/T1/T1
640 C6=C5*C2
650 C7=1.89+C6
660 C8=.21878*T1-(126970!/T1)
670 C9=2*C8*C7-((C3/T1)*126970!)
680 C10=82.546-(162460!/T1)
690 C11=2*C10*C7-(C3*162460!/T1)
700 V=((((C8*C4*C3+C10)*C4/P1)+1)*C3+(4.55504*T1/P1))*.016018  'ft3/lb
710 H=775.596+.63296*T1+1.62467E-04*T1*T1+47.3635*LOG(T1)+.043557*(C7*P1+.5*C4*(C11+C3*(C10+C9*C4)))  'Btu/lb
720 S=((C8*C3-2*C9)*C3*C4/2-C11)*C4/2+(C3-C7)*P1/T1*(-.0241983)-.355579-11.4276/T1+1.8052E-04*T1-.253801*LOG(P1)+.809691*LOG(T1)  'Btu/lb/øF
730 V1=V*6.243386E-02 'm3/Kg
740 H1=H*.5556       'Kcal/Kg
750 S1=S            'Kcal/Kg/øC
760 CLS
770 LOCATE 2,25:PRINT "Program to calculate Steam Properties"
780 LOCATE 4,10:PRINT "Steam Pressure     (Kg/cm2 g)   :";P
790 LOCATE 5,10:PRINT "Steam Saturation Temp (øC)      :";SATTEMP1
800 LOCATE 6,10:PRINT "Steam Temperature  (øC)         :";T
810 LOCATE 7,10:PRINT "Degree of superheat(øC)         :";T-SATTEMP1
820 LOCATE 8,10:PRINT "Specific Volume (m3/kg)         :";V1
830 LOCATE 9,10:PRINT "Steam Enthalpy  (Kcal/Kg)       :";H1
840 LOCATE 10,10:PRINT "Steam Entropy   (Kcal/Kg/øC)    :";S1
850 'common
860 LOCATE 20,20:COLOR 25:PRINT "Press SHIFT+PrintScreen to take the print"
870 COLOR 12:LOCATE 17,15:INPUT "Want to rerun the program or quit (r/q) :";ANS$
880 IF ANS$="q" OR ANS$="Q" THEN SYSTEM
890 IF ANS$="r" OR ANS$="R" THEN 150
900 GOTO 870