149 lines
9.8 KiB
Fortran
149 lines
9.8 KiB
Fortran
subroutine tiopf(t,pf)
|
|
c ======================
|
|
c
|
|
c TiO partition function (data from Kurucz web site)
|
|
c
|
|
INCLUDE 'IMPLIC.FOR'
|
|
dimension pf0(800)
|
|
data pf0/
|
|
* 29.107, 55.425, 82.417, 111.190, 142.564, 176.916,
|
|
* 214.340, 254.774, 298.065, 344.021, 392.431, 443.089,
|
|
* 495.795, 550.365, 606.632, 664.449, 723.686, 784.230,
|
|
* 845.981, 908.862, 972.800, 1037.739, 1103.636, 1170.451,
|
|
* 1238.155, 1306.723, 1376.144, 1446.403, 1517.492, 1589.409,
|
|
* 1662.152, 1735.724, 1810.122, 1885.352, 1961.428, 2038.351,
|
|
* 2116.119, 2194.758, 2274.260, 2354.633, 2435.907, 2518.063,
|
|
* 2601.125, 2685.096, 2769.992, 2855.809, 2942.560, 3030.257,
|
|
* 3118.897, 3208.496, 3299.067, 3390.598, 3483.106, 3576.598,
|
|
* 3671.095, 3766.569, 3863.048, 3960.522, 4059.035, 4158.545,
|
|
* 4259.074, 4360.642, 4463.259, 4566.905, 4671.582, 4777.321,
|
|
* 4884.105, 4991.937, 5100.852, 5210.813, 5321.838, 5433.972,
|
|
* 5547.154, 5661.417, 5776.789, 5893.211, 6010.774, 6129.422,
|
|
* 6249.173, 6370.026, 6491.973, 6615.042, 6739.240, 6864.542,
|
|
* 6990.959, 7118.533, 7247.214, 7377.053, 7508.012, 7640.121,
|
|
* 7773.370, 7907.764, 8043.309, 8180.032, 8317.835, 8456.861,
|
|
* 8597.055, 8738.396, 8880.926, 9024.672, 9169.570, 9315.610,
|
|
* 9462.927, 9611.339, 9760.963, 9911.798, 10063.900, 10217.148,
|
|
* 10371.572, 10527.253, 10684.109, 10842.173, 11001.469, 11161.970,
|
|
* 11323.751, 11486.758, 11650.978, 11816.415, 11983.159, 12151.134,
|
|
* 12320.243, 12490.668, 12662.333, 12835.234, 13009.470, 13184.926,
|
|
* 13361.601, 13539.660, 13718.891, 13899.456, 14081.252, 14264.326,
|
|
* 14448.643, 14634.341, 14821.225, 15009.476, 15199.021, 15389.829,
|
|
* 15581.955, 15775.377, 15970.188, 16166.239, 16363.513, 16562.006,
|
|
* 16761.930, 16963.301, 17165.906, 17369.881, 17575.236, 17781.814,
|
|
* 17989.816, 18198.996, 18409.707, 18621.680, 18835.068, 19049.715,
|
|
* 19265.768, 19483.375, 19702.006, 19922.209, 20143.668, 20366.555,
|
|
* 20590.742, 20816.402, 21043.338, 21271.672, 21501.369, 21732.563,
|
|
* 21965.119, 22199.068, 22434.432, 22671.266, 22909.307, 23148.898,
|
|
* 23389.893, 23632.322, 23875.969, 24121.160, 24367.707, 24615.848,
|
|
* 24865.471, 25116.320, 25368.604, 25622.342, 25877.512, 26134.055,
|
|
* 26392.404, 26651.764, 26912.826, 27175.250, 27439.197, 27704.539,
|
|
* 27971.287, 28239.572, 28509.373, 28780.707, 29053.516, 29327.602,
|
|
* 29603.338, 29880.539, 30159.105, 30439.322, 30721.055, 31004.254,
|
|
* 31288.818, 31575.061, 31862.693, 32151.781, 32442.586, 32734.619,
|
|
* 33027.777, 33323.023, 33619.535, 33917.707, 34217.711, 34518.996,
|
|
* 34821.676, 35126.195, 35432.141, 35739.602, 36048.926, 36359.488,
|
|
* 36672.023, 36985.633, 37300.863, 37617.965, 37936.469, 38256.309,
|
|
* 38578.074, 38901.668, 39226.461, 39552.969, 39880.852, 40210.785,
|
|
* 40541.852, 40874.691, 41209.359, 41545.535, 41883.602, 42222.715,
|
|
* 42563.895, 42906.508, 43250.656, 43596.902, 43944.355, 44293.695,
|
|
* 44644.504, 44997.621, 45351.590, 45707.242, 46065.008, 46424.367,
|
|
* 46785.605, 47148.023, 47512.496, 47878.418, 48246.426, 48615.895,
|
|
* 48987.336, 49360.082, 49734.758, 50111.004, 50489.383, 50868.996,
|
|
* 51250.250, 51633.691, 52018.945, 52405.715, 52794.090, 53184.340,
|
|
* 53576.375, 53970.605, 54366.176, 54763.148, 55162.430, 55563.215,
|
|
* 55966.391, 56371.000, 56777.176, 57185.570, 57596.074, 58007.617,
|
|
* 58421.418, 58837.172, 59254.539, 59673.418, 60094.066, 60517.410,
|
|
* 60941.844, 61368.660, 61797.395, 62227.590, 62659.789, 63094.238,
|
|
* 63529.695, 63967.488, 64407.887, 64849.496, 65292.867, 65735.922,
|
|
* 66182.000, 66631.266, 67082.055, 67534.391, 67988.992, 68446.117,
|
|
* 68904.789, 69365.180, 69827.914, 70292.781, 70759.352, 71228.500,
|
|
* 71699.375, 72171.672, 72647.086, 73123.984, 73603.023, 74083.516,
|
|
* 74566.359, 75050.555, 75537.758, 76027.258, 76518.125, 77012.008,
|
|
* 77507.063, 78003.813, 78503.977, 79006.125, 79509.320, 80015.375,
|
|
* 80522.461, 81031.938, 81544.164, 82058.313, 82574.352, 83093.914,
|
|
* 83614.367, 84136.820, 84662.211, 85188.867, 85719.375, 86249.977,
|
|
* 86783.781, 87319.219, 87857.180, 88396.797, 88939.805, 89484.266,
|
|
* 90032.023, 90580.930, 91132.563, 91686.148, 92242.742, 92799.406,
|
|
* 93360.016, 93923.453, 94488.313, 95055.211, 95625.297, 96197.477,
|
|
* 96771.531, 97348.156, 97926.922, 98507.453, 99091.563, 99677.938,
|
|
*100267.234,100856.438,101449.828,102045.750,102643.094,103244.117,
|
|
*103846.969,104450.313,105057.641,105667.188,106279.516,106894.937,
|
|
*107512.789,108133.117,108754.758,109377.687,110005.039,110634.602,
|
|
*111266.141,111902.133,112537.984,113178.891,113819.766,114464.312,
|
|
*115110.969,115760.687,116412.469,117068.055,117724.547,118384.383,
|
|
*119047.469,119712.469,120380.187,121051.336,121724.102,122399.250,
|
|
*123076.266,123756.977,124441.195,125126.406,125816.453,126506.766,
|
|
*127202.367,127899.086,128598.266,129299.969,130004.969,130712.016,
|
|
*131409.266,132117.719,132828.969,133544.016,134262.750,134986.344,
|
|
*135712.891,136439.937,137170.969,137905.562,138641.578,139380.266,
|
|
*140122.937,140868.641,141615.484,142366.703,143123.078,143880.000,
|
|
*144638.484,145401.594,146168.125,146935.359,147707.484,148482.641,
|
|
*149256.578,150037.281,150821.953,151606.750,152396.094,153188.766,
|
|
*153983.391,154782.141,155582.203,156387.234,157192.719,158003.156,
|
|
*158815.125,159632.437,160450.766,161274.750,162098.172,162926.000,
|
|
*163756.609,164593.141,165430.859,166270.937,167114.750,167960.797,
|
|
*168811.562,169663.906,170517.203,171376.531,172239.469,173105.891,
|
|
*173975.250,174847.203,175721.453,176597.250,177480.984,178366.094,
|
|
*179253.828,180145.734,181038.000,181936.031,182837.969,183739.922,
|
|
*184645.937,185558.281,186470.844,187387.422,188307.234,189232.281,
|
|
*190156.000,191088.234,192022.062,192957.250,193899.328,194842.984,
|
|
*195788.391,196736.156,197687.828,198645.719,199603.422,200569.234,
|
|
*201536.437,202508.641,203481.000,204459.016,205438.750,206424.312,
|
|
*207409.953,208398.734,209393.234,210391.047,211390.984,212395.516,
|
|
*213401.547,214420.141,215431.812,216453.453,217476.734,218501.266,
|
|
*219530.219,220560.719,221597.891,222637.875,223677.750,224725.500,
|
|
*225777.406,226829.297,227893.125,228954.547,230020.969,231086.453,
|
|
*232157.469,233233.047,234315.406,235395.625,236480.953,237572.125,
|
|
*238666.484,239765.125,240863.281,241969.750,243079.250,244191.719,
|
|
*245304.812,246427.937,247548.234,248673.562,249804.984,250942.781,
|
|
*252078.953,253222.812,254369.641,255519.359,256671.406,257827.906,
|
|
*258988.859,260154.734,261322.281,262458.781,263606.437,264770.625,
|
|
*265947.750,267125.156,268314.125,269507.687,270702.344,271905.156,
|
|
*273110.156,274318.937,275531.687,276751.344,277970.781,279198.531,
|
|
*280425.750,281663.250,282897.469,284138.906,285383.594,286637.031,
|
|
*287891.156,289147.625,290413.312,291678.719,292946.031,294225.875,
|
|
*295501.344,296782.656,298070.094,299363.875,300652.250,301953.750,
|
|
*303260.062,304563.781,305874.375,307191.437,308517.031,309835.750,
|
|
*311159.375,312490.937,313827.469,315166.781,316511.031,317860.406,
|
|
*319214.969,320565.875,321929.344,323296.906,324660.219,326035.687,
|
|
*327413.844,328794.406,330173.156,331566.156,332953.469,334356.187,
|
|
*335757.625,337165.562,338566.094,339984.750,341402.937,342828.125,
|
|
*344257.562,345686.750,347123.125,348564.250,350008.906,351453.219,
|
|
*352908.062,354361.469,355828.000,357292.500,358765.719,360233.687,
|
|
*361713.562,363200.187,364685.656,366174.500,367673.594,369174.906,
|
|
*370678.969,372191.125,373708.937,375225.281,376743.719,378270.406,
|
|
*379804.500,381334.250,382879.125,384420.812,385969.531,387519.812,
|
|
*389078.937,390639.781,392213.875,393782.437,395359.156,396943.625,
|
|
*398527.625,400110.937,401711.750,403310.344,404908.937,406513.875,
|
|
*408125.781,409741.906,411356.875,412979.500,414613.125,416245.500,
|
|
*417889.094,419530.000,421179.906,422831.531,424484.344,426153.187,
|
|
*427816.406,429489.094,431161.312,432840.656,434517.000,436215.281,
|
|
*437896.000,439602.594,441300.625,443016.156,444722.906,446445.437,
|
|
*448164.812,449885.937,451615.094,453351.594,455090.125,456833.281,
|
|
*458582.719,460335.344,462094.844,463857.094,465629.906,467402.781,
|
|
*469178.406,470963.750,472745.906,474539.594,476333.312,478131.125,
|
|
*479934.000,481740.750,483557.844,485376.625,487202.937,489033.562,
|
|
*490868.031,492709.281,494547.375,496401.094,498249.594,500110.250,
|
|
*501966.594,503836.062,505704.437,507580.687,509469.187,511349.781,
|
|
*513239.000,515137.187,517038.812,518942.906,520858.156,522767.094,
|
|
*524610.625,526433.812,528331.062,530253.437,532185.500,534127.875,
|
|
*536073.937,538028.312,539983.375,541954.687,543916.312,545902.500,
|
|
*547874.812,549857.125,551850.937,553839.937,555836.625,557838.500,
|
|
*559849.937,561859.375,563880.625,565889.875,567916.000,569953.625,
|
|
*571990.375,574034.937,576085.062,578127.375,580188.937,582251.000,
|
|
*584328.812,586385.562,588464.062,590551.875,592644.625,594722.250,
|
|
*596829.937,598931.375,601029.687,603142.812,605262.812,607384.625,
|
|
*609513.125,611644.000,613775.875,615930.375,618073.750,620218.437,
|
|
*622381.937,624524.312,626697.500,628869.000,631040.937,633223.562,
|
|
*635409.187,637597.562,639800.187,642002.125,644212.562,646416.250,
|
|
*648633.562,650864.187,653083.687,655315.312,657549.687,659795.500,
|
|
*662032.250,664292.875,666542.312,668806.250,671071.312,673340.937,
|
|
*675626.938,677898.750/
|
|
c
|
|
it=int(t/10.)
|
|
if(it.gt.800) it=800
|
|
pf=pf0(it)
|
|
return
|
|
end
|