17 December 2023 2:04:41.662 PM steam_test(): FORTRAN90 version Test steam(). TEST01 THERM computes CP, the specific heat capacity at constant pressure. P T Rho CP CP (bar) (C) g/cm3 tabulated computed 1.000000 0.000000 0.999828 4.22800000 4.22822199 1.000000 100.000000 0.000590 2.04200000 2.04184907 1.000000 200.000000 0.000460 1.97500000 1.97520198 1.000000 300.000000 0.000379 2.01300000 2.01293308 1.000000 350.000000 0.000348 2.04000000 2.04026979 1.000000 400.000000 0.000322 2.07000000 2.07013339 1.000000 500.000000 0.000280 2.13500000 2.13455122 1.000000 600.000000 0.000248 2.20300000 2.20288857 1.000000 850.000000 0.000193 2.37800000 2.37809639 1.000000 1100.000000 0.000158 2.54100000 2.54070979 1.000000 1600.000000 0.000116 2.79200000 2.79230214 1.000000 2000.000000 0.000095 2.93100000 2.93149885 5.000000 0.000000 1.000032 4.22600000 4.22601245 10.000000 0.000000 1.000286 4.22300000 4.22326535 50.000000 0.000000 1.002311 4.20200000 4.20186900 100.000000 0.000000 1.004813 4.17700000 4.17651887 200.000000 0.000000 1.009728 4.13000000 4.13013704 300.000000 0.000000 1.014530 4.08900000 4.08901036 400.000000 0.000000 1.019226 4.05300000 4.05261620 500.000000 0.000000 1.023819 4.02100000 4.02048977 1000.000000 0.000000 1.045306 3.90900000 3.90916764 1500.000000 0.000000 1.064366 3.84400000 3.84428059 2000.000000 0.000000 1.081131 3.78600000 3.78637972 5000.000000 0.000000 1.152606 2.89000000 2.89459734 Maximum absolute error was 0.459734E-02 TEST02 DENSE computes fluid density. Pressure is fixed at 0.100000 MegaPascals. The saturation temperature is 372.782 degrees Kelvin T (C) T (K) RHO (KG/M3) dPdR 0.00000 273.150 999.82809 1.9611900 5.00000 278.150 999.99261 2.0338820 10.0000 283.150 999.74023 2.0942830 15.0000 288.150 999.13712 2.1442288 20.0000 293.150 998.23194 2.1851042 25.0000 298.150 997.06136 2.2179888 30.0000 303.150 995.65354 2.2437491 35.0000 308.150 994.03051 2.2630976 40.0000 313.150 992.20979 2.2766327 45.0000 318.150 990.20557 2.2848670 50.0000 323.150 988.02950 2.2882461 55.0000 328.150 985.69128 2.2871634 60.0000 333.150 983.19910 2.2819701 65.0000 338.150 980.55998 2.2729827 70.0000 343.150 977.77995 2.2604895 75.0000 348.150 974.86426 2.2447543 80.0000 353.150 971.81749 2.2260203 85.0000 358.150 968.64361 2.2045129 90.0000 363.150 965.34613 2.1804414 95.0000 368.150 961.92805 2.1540009 100.000 373.150 0.58958161 0.16695330 110.000 383.150 0.57313852 0.17208280 120.000 393.150 0.55769897 0.17713767 130.000 403.150 0.54315537 0.18213200 140.000 413.150 0.52941819 0.18707668 150.000 423.150 0.51641181 0.19198028 160.000 433.150 0.50407152 0.19684956 170.000 443.150 0.49234136 0.20168996 180.000 453.150 0.48117245 0.20650584 190.000 463.150 0.47052175 0.21130078 200.000 473.150 0.46035108 0.21607769 300.000 573.150 0.37896022 0.26323817 400.000 673.150 0.32230020 0.30986592 500.000 773.150 0.28046167 0.35628420 600.000 873.150 0.24826823 0.40260240 700.000 973.150 0.22271800 0.44886646 800.000 1073.15 0.20194253 0.49509848 900.000 1173.15 0.18471575 0.54131019 1000.00 1273.15 0.17019901 0.58750831 TEST03 DIELEC computes the dielectric constant. P T Rho EPS EPS (bar) (C) g/cm3 tabulated computed 100.000000 0.000000 1.004813 88.29000000 88.28851185 500.000000 0.000000 1.023819 90.07000000 90.06662931 1000.000000 0.000000 1.045306 92.02000000 92.02142668 2000.000000 0.000000 1.081131 95.14000000 95.13976079 5000.000000 0.000000 1.152606 100.77000000 100.77176538 100.000000 25.000000 1.001475 78.85000000 78.84656060 100.000000 50.000000 0.992305 70.27000000 70.26716075 100.000000 75.000000 0.979217 62.60000000 62.59810640 100.000000 100.000000 0.962980 55.78000000 55.77714444 100.000000 150.000000 0.922375 44.31000000 44.31246493 100.000000 200.000000 0.871032 35.11000000 35.10808595 100.000000 300.000000 0.715582 20.40000000 20.39699056 100.000000 400.000000 0.037867 1.17000000 1.17122010 100.000000 500.000000 0.030503 1.11000000 1.11356056 100.000000 600.000000 0.026068 1.08000000 1.08870726 Maximum absolute error = 0.870726E-02 TEST04 PRANDTL computes the Prandtl number given the pressure and temperature. P T PR PR (bar) (C) tabulated computed 1.000000 0.000000 13.50000000 13.50352690 5.000000 0.000000 13.48000000 13.48395658 10.000000 0.000000 13.46000000 13.45961514 25.000000 0.000000 13.39000000 13.38739260 50.000000 0.000000 13.27000000 13.26964555 75.000000 0.000000 13.15000000 13.15509364 100.000000 0.000000 13.04000000 13.04364396 125.000000 0.000000 12.93000000 12.93520835 150.000000 0.000000 12.83000000 12.82970308 175.000000 0.000000 12.73000000 12.72704837 200.000000 0.000000 12.63000000 12.62716802 225.000000 0.000000 12.53000000 12.52998945 250.000000 0.000000 12.43000000 12.43544281 275.000000 0.000000 12.34000000 12.34346116 300.000000 0.000000 12.25000000 12.25398012 350.000000 0.000000 12.08000000 12.08227380 400.000000 0.000000 11.92000000 11.91985236 450.000000 0.000000 11.77000000 11.76627356 500.000000 0.000000 11.62000000 11.62112011 550.000000 0.000000 11.48000000 11.48399608 600.000000 0.000000 11.36000000 11.35452397 650.000000 0.000000 11.23000000 11.23234239 700.000000 0.000000 11.12000000 11.11710419 800.000000 0.000000 10.91000000 10.90613199 900.000000 0.000000 10.72000000 10.71906685 1000.000000 0.000000 10.55000000 10.55353223 1.000000 25.000000 6.13700000 6.13710111 1.000000 50.000000 3.55500000 3.55462527 1.000000 75.000000 2.37800000 2.37794267 1.000000 100.000000 1.00000000 0.99960846 1.000000 150.000000 0.97400000 0.97373693 1.000000 200.000000 0.96000000 0.96028283 1.000000 400.000000 0.92400000 0.92434641 VISCOSITY - Warning (once only)! The input temperature T should be no more than 800.000 The input value was T = 873.150 1.000000 600.000000 0.89900000 0.89932214 1.000000 800.000000 0.88200000 0.88158114 Maximum absolute error = 0.547603E-02 TEST05 PSAT computes the saturation pressure given the temperature. T P P (C) (bar) (bar) tabulated computed 0.010000 0.00611730 0.00611732 1.000000 0.00657160 0.00657156 5.000000 0.00872600 0.00872604 50.000000 0.12344000 0.12344470 100.000000 1.01320000 1.01321994 125.000000 2.32010000 2.32014417 150.000000 4.75720000 4.75716897 200.000000 15.53700000 15.53650004 250.000000 39.73700000 39.73649419 300.000000 85.83800000 85.83784638 350.000000 165.21000000 165.21127636 373.976000 220.55000000 220.59772620 Maximum absolute error = 0.477262E-01 TEST06 SECVIR computes the second virial coefficient as a function of temperature. T VIR VIR (C) tabulated computed 0.000000 -98.96000000 -98.96364064 5.000000 -90.08000000 -90.08483145 10.000000 -82.29000000 -82.29403175 20.000000 -69.36000000 -69.36483095 30.000000 -59.19000000 -59.18835975 40.000000 -51.07000000 -51.06674386 60.000000 -39.13000000 -39.12547788 90.000000 -27.81000000 -27.80930215 120.000000 -20.83000000 -20.82778985 150.000000 -16.21000000 -16.20744958 180.000000 -12.98000000 -12.97900648 210.000000 -10.63000000 -10.62587300 240.000000 -8.85000000 -8.85270570 300.000000 -6.39000000 -6.39454283 400.000000 -4.03000000 -4.03226704 500.000000 -2.71000000 -2.70564339 700.000000 -1.32000000 -1.31848524 1000.000000 -0.39000000 -0.39280130 2000.000000 0.53000000 0.53481651 Maximum absolute error = 0.483145E-02 TEST07 SOUND computes the speed of sound. P T C C (bar) (C) tabulated computed 1.000000 0.000000 1401.00000000 1400.99622226 1.000000 100.000000 472.80000000 472.83836607 1.000000 200.000000 533.70000000 533.71269122 1.000000 300.000000 585.70000000 585.65694023 1.000000 350.000000 609.50000000 609.47882785 1.000000 400.000000 632.20000000 632.15734219 1.000000 500.000000 674.60000000 674.63137322 1.000000 600.000000 713.90000000 713.93144211 1.000000 850.000000 802.00000000 801.99270196 1.000000 1100.000000 880.10000000 880.04633430 1.000000 1600.000000 1017.80000000 1017.75218294 1.000000 2000.000000 1115.90000000 1115.93283773 5.000000 0.000000 1401.70000000 1401.70421558 10.000000 0.000000 1402.60000000 1402.58704023 50.000000 0.000000 1409.60000000 1409.56938512 100.000000 0.000000 1418.10000000 1418.12638234 250.000000 0.000000 1443.10000000 1443.10393822 500.000000 0.000000 1484.60000000 1484.63536363 1000.000000 0.000000 1577.10000000 1577.14652411 2500.000000 0.000000 1913.40000000 1913.37379474 Maximum absolute error = 0.536657E-01 TEST08 SURTEN computes the surface tension given the temperature. T SIGMA SIGMA (C) tabulated computed 10.000000 74.22000000 74.22448389 20.000000 72.74000000 72.73984283 30.000000 71.20000000 71.19812265 40.000000 69.60000000 69.60055920 50.000000 67.95000000 67.94843377 100.000000 58.92000000 58.91786932 150.000000 48.75000000 48.74887099 200.000000 37.68000000 37.68352563 250.000000 26.05000000 26.05327715 300.000000 14.37000000 14.37056572 325.000000 8.78000000 8.78523918 350.000000 3.67000000 3.67520935 370.000000 0.40000000 0.39484620 373.976000 0.00000000 0.00064188 Maximum absolute error = 0.523918E-02 TEST09 THERCON computes thermal conductivity. P T LAMBDA LAMBDA (bar) (C) tabulated computed 1.000000 0.000000 561.00000000 561.02950035 5.000000 0.000000 561.30000000 561.25078490 10.000000 0.000000 561.50000000 561.52742144 25.000000 0.000000 562.40000000 562.35749299 50.000000 0.000000 563.70000000 563.74121511 75.000000 0.000000 565.10000000 565.12479451 100.000000 0.000000 566.50000000 566.50768943 125.000000 0.000000 567.90000000 567.88937061 150.000000 0.000000 569.30000000 569.26932082 175.000000 0.000000 570.60000000 570.64703447 200.000000 0.000000 572.00000000 572.02201730 225.000000 0.000000 573.40000000 573.39378616 250.000000 0.000000 574.80000000 574.76186885 275.000000 0.000000 576.10000000 576.12580397 300.000000 0.000000 577.50000000 577.48514085 350.000000 0.000000 580.20000000 580.18827092 400.000000 0.000000 582.90000000 582.86786856 450.000000 0.000000 585.50000000 585.52071675 500.000000 0.000000 588.10000000 588.14377501 550.000000 0.000000 590.70000000 590.73418338 600.000000 0.000000 593.30000000 593.28926678 650.000000 0.000000 595.80000000 595.80653572 700.000000 0.000000 598.30000000 598.28370663 800.000000 0.000000 603.10000000 603.10956144 900.000000 0.000000 607.80000000 607.75247767 1000.000000 0.000000 612.20000000 612.20133542 1.000000 25.000000 607.20000000 607.17353903 1.000000 50.000000 643.60000000 643.55651739 1.000000 75.000000 666.80000000 666.78588466 1.000000 100.000000 25.08000000 25.07665477 1.000000 150.000000 28.85000000 28.85257097 1.000000 200.000000 33.28000000 33.28003963 1.000000 400.000000 54.76000000 54.75549185 1.000000 600.000000 79.89000000 79.88981795 1.000000 800.000000 107.30000000 107.32449185 Maximum absolute error = 0.492151E-01 TEST10 TSAT computes the saturation temperature given the pressure. P T T (bar) (C) (C) tabulated computed 0.006117 0.01000000 0.00996289 0.012000 9.65500000 9.65456141 0.025000 21.08000000 21.08049199 0.055000 34.58900000 34.58876805 0.080000 41.51800000 41.51824122 0.110000 47.69500000 47.69455640 0.160000 55.32700000 55.32696030 0.250000 64.98000000 64.97989850 0.500000 81.33900000 81.33882463 0.750000 91.78300000 91.78263369 1.000000 99.63200000 99.63238208 1.500000 111.37800000 111.37825185 2.000000 120.44300000 120.24059079 5.000000 151.86600000 151.86640634 10.000000 179.91600000 179.91631089 20.000000 212.41700000 212.41660361 50.000000 263.97700000 263.97732634 100.000000 311.03100000 311.03063189 200.000000 365.80000000 365.79975211 220.550000 373.97600000 373.95813643 Maximum absolute error = 0.202409 TEST11 VISCOSITY computes the viscosity. P T RHO ETA ETA (bar) (C) g/cm3 tabulated computed 1.000000 0.000000 0.999828 1792.00000000 1791.74058525 5.000000 0.000000 1.000032 1791.00000000 1790.78535511 10.000000 0.000000 1.000286 1790.00000000 1789.59699714 25.000000 0.000000 1.001048 1786.00000000 1786.06951039 50.000000 0.000000 1.002311 1780.00000000 1780.31397681 75.000000 0.000000 1.003566 1775.00000000 1774.71003032 100.000000 0.000000 1.004813 1769.00000000 1769.25445157 125.000000 0.000000 1.006053 1764.00000000 1763.94418185 150.000000 0.000000 1.007285 1759.00000000 1758.77631018 175.000000 0.000000 1.008510 1754.00000000 1753.74806154 200.000000 0.000000 1.009728 1749.00000000 1748.85678111 225.000000 0.000000 1.010938 1744.00000000 1744.09994446 250.000000 0.000000 1.012142 1739.00000000 1739.47511801 275.000000 0.000000 1.013339 1735.00000000 1734.97997119 300.000000 0.000000 1.014530 1731.00000000 1730.61226354 350.000000 0.000000 1.016891 1722.00000000 1722.25061235 400.000000 0.000000 1.019226 1714.00000000 1714.37378676 450.000000 0.000000 1.021535 1707.00000000 1706.96645028 500.000000 0.000000 1.023819 1700.00000000 1700.01414904 550.000000 0.000000 1.026078 1694.00000000 1693.50317669 600.000000 0.000000 1.028313 1687.00000000 1687.42045988 650.000000 0.000000 1.030522 1682.00000000 1681.75346156 700.000000 0.000000 1.032708 1676.00000000 1676.49009977 800.000000 0.000000 1.037005 1667.00000000 1667.12784174 900.000000 0.000000 1.041204 1659.00000000 1659.24386525 1000.000000 0.000000 1.045306 1653.00000000 1652.75246271 1.000000 25.000000 0.997061 890.80000000 890.79821418 1.000000 50.000000 0.988030 547.10000000 547.08093048 1.000000 75.000000 0.974864 378.40000000 378.41428910 1.000000 100.000000 0.958392 12.28000000 282.31580357 1.000000 200.000000 0.863624 16.18000000 133.24483037 1.000000 400.000000 0.000322 24.45000000 24.44916952 1.000000 600.000000 0.000248 32.61000000 32.61475996 1.000000 800.000000 0.000202 40.38000000 40.37452760 Maximum absolute error = 270.036 TEST12 VOLUME computes the specific volume. T(C) = 0.0000000 T(K) = 273.15000 RHO V dVdT dVdRHO 0.999828 1.0001719 -0.80236982E-04 -1.0003439 1.00003 0.99996806 -0.78550925E-04 -0.99993612 1.00029 0.99971362 -0.76451762E-04 -0.99942733 1.00105 0.99895306 -0.70209390E-04 -0.99790722 1.00231 0.99769447 -0.59983636E-04 -0.99539425 1.00357 0.99644689 -0.49971396E-04 -0.99290641 1.00481 0.99521009 -0.40163045E-04 -0.99044313 1.00605 0.99398383 -0.30550008E-04 -0.98800385 1.00728 0.99276788 -0.21124661E-04 -0.98558806 1.00851 0.99156204 -0.11880249E-04 -0.98319528 1.00973 0.99036612 -0.28107971E-05 -0.98082505 1.01094 0.98917995 0.60889186E-05 -0.97847696 1.01214 0.98800335 0.14823495E-04 -0.97615062 1.01334 0.98683618 0.23396922E-04 -0.97384564 1.01453 0.98567829 0.31812646E-04 -0.97156169 1.01689 0.98338984 0.48182338E-04 -0.96705557 1.01923 0.98113702 0.63951061E-04 -0.96262986 1.02153 0.97891901 0.79131437E-04 -0.95828244 1.02382 0.97673508 0.93731375E-04 -0.95401142 1.02608 0.97458458 0.10775508E-03 -0.94981510 1.02831 0.97246695 0.12120396E-03 -0.94569197 1.03052 0.97038169 0.13407738E-03 -0.94164062 1.03271 0.96832834 0.14637340E-03 -0.93765978 1.03700 0.96431581 0.16922240E-03 -0.92990498 1.04120 0.96042642 0.18973003E-03 -0.92241891 1.04531 0.95665751 0.20788516E-03 -0.91519359 T(C) = 25.000000 T(K) = 298.15000 RHO V dVdT dVdRHO 0.997061 1.0029473 0.26019115E-03 -1.0059033 0.997242 1.0027660 0.26057424E-03 -1.0055397 0.997467 1.0025397 0.26105220E-03 -1.0050859 0.998140 1.0018630 0.26248025E-03 -1.0037295 0.999258 1.0007423 0.26484118E-03 -1.0014851 1.00037 0.99963023 0.26717883E-03 -0.99926059 1.00148 0.99852677 0.26949382E-03 -0.99705570 1.00257 0.99743171 0.27178666E-03 -0.99487001 1.00367 0.99634489 0.27405774E-03 -0.99270315 1.00476 0.99526619 0.27630738E-03 -0.99055479 1.00584 0.99419545 0.27853579E-03 -0.98842459 1.00691 0.99313254 0.28074310E-03 -0.98631225 1.00799 0.99207735 0.28292937E-03 -0.98421746 1.00905 0.99102973 0.28509458E-03 -0.98213993 1.01011 0.98998959 0.28723866E-03 -0.98007939 1.01222 0.98793127 0.29146284E-03 -0.97600820 1.01430 0.98590155 0.29560025E-03 -0.97200187 1.01636 0.98389966 0.29964853E-03 -0.96805855 1.01841 0.98192489 0.30360481E-03 -0.96417648 1.02043 0.97997657 0.30746577E-03 -0.96035407 1.02244 0.97805409 0.31122781E-03 -0.95658980 1.02443 0.97615688 0.31488708E-03 -0.95288226 1.02639 0.97428442 0.31843964E-03 -0.94923012 1.03028 0.97061173 0.32520875E-03 -0.94208714 1.03409 0.96703239 0.33150409E-03 -0.93515165 1.03784 0.96354312 0.33729714E-03 -0.92841534 steam_test(): Normal end of execution. 17 December 2023 2:04:41.681 PM