CLS
PRINT "THERMODYNAMIC PROPERTIES OF AIR"
[ONONZE]
PRINT
PRINT "BAR PRESSURE (INCHES HG) = ";
INPUT " ";P
[ONFIZE]
PRINT "DRYBULB TEMP (DEG F) =     ";
INPUT " ";T1
PRINT "WET BULB TEMP(DEG F) =     ";
INPUT " ";T2
PRINT
PRINT
REM ==CHECK INPUT DATA==
IF T1>=T2 GOTO [TWSEZE]
PRINT "WET BULB MUST NOT EXCEED DRY BULB"
GOTO [ONFIZE]
REM ==CONVERT TO RANKIN==
[TWSEZE]
T4 = (T2 + 459.688)/1.8
REM ==CALCULATE SATURATION PRESSURE/WET BULB==
GOSUB [ONZETWZE]
P8 = P1
REM ==CONVERT TO RANKIN==
T4 = (T1 + 459.688)/1.8
REM ==CALC SATURATION PRESSURE/DRY BULB==
GOSUB [ONZETWZE]
P9 = P1
W1 = 0.622*P8/(P-P8)
IF T2 < 32 THEN [FOFIZE]
IF T2 = 32 THEN [FOFIZE]
H1 = 1093.049+0.441*T1-T2
Z1 = 0.24+0.441*W1
Z2 = W1-Z1*(T1-T2) /H1
P2 = P * Z2 / (0.622 + Z2)
GOTO [FOSEZE]
[FOFIZE]
P2 = P8 -(5.704/10^4) * P * (T1-T2) / 1.8
REM ==CALC HUMIDITY RATIO W==
[FOSEZE]
W = 0.622 * P2 /(P-P2)
REM ==CALC ENTHALPY H==
H = 0.24 * T1 + (1061 + 0.444 * T1) * W
REM ==CALC DEW POINT T3==
T3 = T2
IF T1 = T2 THEN [FINIZE]
L = LOG(P2)
IF P2 < .18036 THEN [FISEZE]
T3 = 79.047 + 30.579 * L + 1.8893 * L * L
GOTO [FINIZE]
[FISEZE]
T3 = 71.98 + 24.873 * L + 0.8927 * L * L
REM ==CALC VOLUME MOIST AIR V==
[FINIZE]
V = 0.754 * (T1 + 459.7) * (1 + 7000 * W / 4360) / P
REM ==CALC RELATIVE HUMIDITY==
R = (P2 / P9) * 100
PRINT "DEW POINT TEMP (DEG F) =     ";T3
PRINT "VAPOR PRESSURE (INCH HG) =   ";P2
PRINT "HUMIDTY RATIO  (#WAT/#AIR) = ";W
PRINT "VOL MOIST AIR  (FT^3/#) =    ";V
PRINT "ENTHALPY       (BTU/#)=      ";H
PRINT "REL HUMIDITY   (PERCENT)=    ";R
PRINT
PRINT "RERUN PROGRAM WITH NEW DATA (Y/N) ?  ";
INPUT " ";Q$
IF Q$ = "y" OR Q$="Y" THEN [ONONZE]
GOTO [ONTWSIZE]
REM ==SUBROUTINE VAPOR PRESSURE OF SATURATED AIR==
[ONZETWZE]
A1 = -7.90298
A2 = 5.02808
A3 = -1.3816/10^7
A4 = 11.344
A5 = 8.1328/10^3
A6 = -3.49149
B1 = -9.09718
B2 = -3.56654
B3 = .876793
B4 = 0.0060273
IF T4 < 273.16 THEN [ONONNIZE]
Z = 373.16 / T4
C1 = A1 * (Z-1)
C2 = A2*LOG(Z) /LOG(10)
C3 = A3 *(10^(A4*(1-1/Z))-1)
C4 = A5 *(10^(A6*(Z-1))-1)
GOTO [ONTWFOZE]
[ONONNIZE]
Z = 273.16/T4
C1 = B1*(Z-1)
C2 = B2*LOG(Z)/LOG(10)
C3 = B3*(1-1/Z)
C4 =LOG(84)/LOG(10)
[ONTWFOZE]
P1 = 29.921*(10^(C1 + C2 + C3 + C4))
RETURN
[ONTWSIZE]
END