SUBROUTINE PFSPEC(IAT,IZI,T,ANE,U) C ================================== C Non-standard evaluation of the partition function C user supplied procedure C C Input: C IAT - atomic number C IZI - ionic charge (=1 for neutrals, =1 for once ionized, etc) C T - temperature C ANE - electron density C XMAX - principal quantum number of the last bound level C C Output: C U - partition function C * * Modified from the ATMOS related programme 5-April-1990 * as an addition to TLUSTY to allow high ionisation states * of C, N and O * * M.A.Barstow - University of Leicester, Dept of Physics & Astronomy * INCLUDE 'PARAMS.FOR' real nvii PARAMETER (MH=100,MHEI=100,MHEII=100,MCI=135, + MCII=157,MCIII=156,MCIV=55,MCV=15,MCVI=100,MNI=228,MNII=122, + MNIII=133,MNIV=73,MNV=51,MNVI=8,MNVII=100,MOI=174,MOII=191, + MOIII=168,MOIV=166,MOV=115,MOVI=52,MOVII=16,MOVIII=100) DIMENSION GHYD(MH),SHYD(MH),ENHYD(MH), + GHEL(MH),ENHEL(MH),SHEL(MH), + GCI(MCI),ENCI(MCI),SCI(MCI), + GCII(MCII),ENCII(MCII),SCII(MCII), + GCIII(MCIII),ENCIII(MCIII),SCIII(MCIII), + GCIV(MCIV),ENCIV(MCIV),SCIV(MCIV), + GCV(MCV),ENCV(MCV),SCV(MCV), + GNI(MNI),ENNI(MNI),SNI(MNI), + GNII(MNII),ENNII(MNII),SNII(MNII), + GNIII(MNIII),ENNIII(MNIII),SNIII(MNIII), + GNIV(MNIV),ENNIV(MNIV),SNIV(MNIV), + GNV(MNV),ENNV(MNV),SNV(MNV), + GNVI(MNVI),ENNVI(MNVI),SNVI(MNVI), + GOI(MOI),ENOI(MOI),SOI(MOI), + GOII(MOII),ENOII(MOII),SOII(MOII), + GOIII(MOIII),ENOIII(MOIII),SOIII(MOIII), + GOIV(MOIV),ENOIV(MOIV),SOIV(MOIV), + GOV(MOV),ENOV(MOV),SOV(MOV), + GOVI(MOVI),ENOVI(MOVI),SOVI(MOVI), + GOVII(MOVII),ENOVII(MOVII),SOVII(MOVII) INTEGER NHYD(MH),NHEL(MHEI),NCI(MCI),NCII(MCII), + NCIII(MCIII),NCIV(MCIV),NCV(MCV),NNI(MNI),NNII(MNII), + NNIII(MNIII),NNIV(MNIV),NNV(MNV),NNVI(MNVI),NOI(MOI), + NOII(MOII),NOIII(MOIII),NOIV(MOIV),NOV(MOV),NOVI(MOVI), + NOVII(MOVII) PARAMETER (HI=13.5878,HEI=24.587,HEII=54.416,CVI=489.84, + NVII=666.83,OVIII=871.12) PARAMETER (ZH=1.0,ZHE=2.0,ZC=6.0,ZN=7.0,ZO=8.0) C N***=QUANTUM NO. OF LEVEL C DATA FOR IONS G***=STATISTICAL WEIGHT OF LEVEL C EN***=ENERGY OF LEVEL C S*=SCREENING NO. OF LEVEL DATA NHYD/ 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,69,70,71,72, + 73,74,75,76,77,78, + 79,80,81,82,83,84, + 85,86,87,88,89,90, + 91,92,93,94,95,96, + 97,98,99, 100 / DATA GHYD/ 2.000000, 8.000000, 18.00000, + 32.00000, 50.00000, 72.00000, + 98.00000, 128.0000, 162.0000, + 200.0000, 242.0000, 288.0000, + 338.0000, 392.0000, 450.0000, + 512.0000, 578.0000, 648.0000, + 722.0000, 800.0000, 882.0000, + 968.0000, 1058.000, 1152.000, + 1250.000, 1352.000, 1458.000, + 1568.000, 1682.000, 1800.000, + 1922.000, 2048.000, 2178.000, + 2312.000, 2450.000, 2592.000, + 2738.000, 2888.000, 3042.000, + 3200.000, 3362.000, 3528.000, + 3698.000, 3872.000, 4050.000, + 4232.000, 4418.000, 4608.000, + 4802.000, 5000.000, 5202.000, + 5408.000, 5618.000, 5832.000, + 6050.000, 6272.000, 6498.000, + 6728.000, 6962.000, 7200.000, + 7442.000, 7688.000, 7938.000, + 8192.000, 8450.000, 8712.000, + 8978.000, 9248.000, 9522.000, + 9800.000, 10082.00, 10368.00, + 10658.00, 10952.00, 11250.00, + 11552.00, 11858.00, 12168.00, + 12482.00, 12800.00, 13122.00, + 13448.00, 13778.00, 14112.00, + 14450.00, 14792.00, 15138.00, + 15488.00, 15842.00, 16200.00, + 16562.00, 16928.00, 17298.00, + 17672.00, 18050.00, 18432.00, + 18818.00, 19208.00, 19602.00, + 20000.00/ DATA ENHYD /0.0000000E+00,10.19085000000000,12.07804444444444, + 12.73856250000000,13.04428800000000,13.21036111111111, + 13.31049795918367,13.37549062500000,13.42004938271605, + 13.45192200000000,13.47550413223140,13.49344027777778, + 13.50739881656805,13.51847448979592,13.52740977777778, + 13.53472265625000,13.54078339100346,13.54586234567901, + 13.55016066481994,13.55383050000000,13.55698866213152, + 13.55972603305785,13.56211417769376,13.56421006944444, + 13.56605952000000,13.56769970414201,13.56916104252401, + 13.57046862244898,13.57164328180737,13.57270244444444, + 13.57366077003122,13.57453066406250,13.57532268135905, + 13.57604584775087,13.57670791836735,13.57731558641975, + 13.57787465303141,13.57839016620499,13.57886653517423, + 13.57930762500000,13.57971683521713,13.58009716553288, + 13.58045127095727,13.58078150826446,13.58108997530864, + 13.58137854442344,13.58164889090086,13.58190251736111, + 13.58214077467722,13.58236488000000,13.58257593233372, + 13.58277492603550,13.58296276254895,13.58314026063100, + 13.58330816528926,13.58346715561225,13.58361785164666, + 13.58376082045184,13.58389658144211,13.58402561111111, + 13.58414834721849,13.58426519250780,13.58437651801461, + 13.58448266601563,13.58458395266272,13.58468067033976, + 13.58477308977501,13.58486146193772,13.58494601974375, + 13.58502697959184,13.58510454274945,13.58517889660494, + 13.58525021580034,13.58531866325785,13.58538439111111, + 13.58544754155125,13.58550824759656,13.58556663379356, + 13.58562281685627,13.58567690625000,13.58572900472489, + 13.58577920880428,13.58582760923211,13.58587429138322, + 13.58591933564014,13.58596281773932,13.58600480908971, + 13.58604537706612,13.58608458527964,13.58612249382716, + 13.58615915952180,13.58619463610586,13.58622897444791, + 13.58626222272522,13.58629442659280,13.58632562934028, + 13.58635587203741,13.58638519366930,13.58641363126212, + 13.58644122000000/ DATA SHYD/100*0.0D0/ DATA NHEL/1,2,2,2,2,2,2,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,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,69, + 70,71,72,73,74,75, + 76,77,78,79,80,81/ DATA GHEL/1.0D0,3.0D0,1.0D0,5.0D0,3.0D0,1.0D0,3.0D0, + 3.0D0,1.0D0,5.0D0,3.0D0, + 1.0D0,15.0D0,5.0D0,3.0D0,3.0D0,1.0D0,9.0D0, + 15.0D0,5.0D0,21.0D0,7.0D0, + 3.0D0,100.0D0,144.0D0,196.0D0,256.0D0,324.0D0, + 400.0D0,484.0D0, + 576.0D0,676.0D0,784.0D0,900.0D0,1024.0D0,1156.0D0, + 1296.0D0,1444.0D0,1600.0D0,1764.0D0,1936.0D0, + 2116.0D0,2304.0D0,2500.0D0,2704.0D0,3136.0D0, + 3136.000000000000,3364.000000000000,3600.000000000000, + 3844.000000000000,4096.000000000000,4356.000000000000, + 4624.000000000000,4900.000000000000,5184.000000000000, + 5476.000000000000,5776.000000000000,6084.000000000000, + 6400.000000000000,6724.000000000000,7056.000000000000, + 7396.000000000000,7744.000000000000,8100.000000000000, + 8464.000000000000,8836.000000000000,9216.000000000000, + 9604.000000000000,10000.00000000000,10404.00000000000, + 10816.00000000000,11236.00000000000,11664.00000000000, + 12100.00000000000,12544.00000000000,12996.00000000000, + 13456.00000000000,13924.00000000000,14400.00000000000, + 14884.00000000000,15376.00000000000,15876.00000000000, + 16384.00000000000,16900.00000000000,17424.00000000000, + 17956.00000000000,18496.00000000000,19044.00000000000, + 19600.00000000000,20164.00000000000,20736.00000000000, + 21316.00000000000,21904.00000000000,22500.00000000000, + 23104.00000000000,23716.00000000000,24336.00000000000, + 24964.00000000000,25600.00000000000,26244.00000000000/ DATA ENHEL/0.0D0,19.819D0,20.615D0,20.964D0, + 20.964D0,20.964D0,21.218D0, + 22.718D0,22.920D0,23.007D0,23.007D0, + 23.007D0,23.073D0,23.074D0, + 23.087D0,23.593D0,23.673D0,23.707D0, + 23.736D0,23.736D0,23.737D0, + 23.737D0,23.742D0,24.028D0,24.201D0, + 24.304D0,24.371D0,24.417D0, + 24.449D0,24.473D0,24.491D0,24.506D0, + 24.517D0,24.526D0,24.534D0, + 24.540D0,24.545D0,24.549D0,24.553D0, + 24.556D0,24.559D0,24.562D0, + 24.564D0,24.566D0,24.568D0,24.570D0, + 24.57131951530612,24.57238228299643,24.57334055555556, + 24.57420759625390,24.57499462890625,24.57571120293848, + 24.57636548442907,24.57696448979592,24.57751427469136, + 24.57802008765522,24.57848649584488,24.57891748849441, + 24.57931656250000,24.57968679357525,24.58003089569161, + 24.58035127095727,24.58065005165289,24.58092913580247, + 24.58119021739130,24.58143481213219,24.58166427951389, + 24.58187984173261,24.58208260000000,24.58227354863514, + 24.58245358727811,24.58262353150587,24.58278412208505, + 24.58293603305785,24.58307987882653,24.58321622037550, + 24.58334557074911,24.58346839988509,24.58358513888889, + 24.58369618382155,24.58380189906348,24.58390262030738, + 24.58399865722656,24.58409029585799,24.58417780073462, + 24.58426141679661,24.58434137110727,24.58441787439614, + 24.58449112244898,24.58456129736163,24.58462856867284, + 24.58469309438919,24.58475502191381,24.58481448888889, + 24.58487162396122,24.58492654747850,24.58497937212360, + 24.58503020349303,24.58507914062500,24.58512627648224/ DATA SHEL/0.375D0,0.622D0,0.622D0,0.842D0, + 0.842D0,0.842D0,0.842D0,0.747D0, + 0.747D0,0.912D0,0.912D0,0.912D0, + 0.993D0,0.993D0,0.912D0,0.810D0, + 0.810D0,0.937D0,0.995D0,0.995D0, + 1.000D0,1.000D0,0.937D0,0.949D0, + 0.958D0,75*1.000D0/ DATA NCI/2,2,2,2,2,2,3,3,3,3,2,2,2,3,3,3,3,3, + 3,3,3,3,3,2,3,4,4,4,3,3,3,3,3,3,4,3, + 3,3,3,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4, + 4,4,4,5,4,4,4,4,4,5,5,5,5,5,5,5,5,5, + 5,5,5,5,6,5,5,5,5,5,6,6,6,6,6,6,6,7, + 6,6,6,6,6,7,7,7,7,7,7,7,7,7,7,7,8,8, + 8,8,8,8,8,8,8,8,9,9,9,9,9,9,9,10,10, + 10,11,11,11,2,3,3,3,2,2/ DATA GCI/1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 5.0D0,1.0D0,3.0D0,5.0D0,3.0D0, + 7.0D0,5.0D0,3.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,3.0D0,1.0D0,3.0D0, + 5.0D0,5.0D0,1.0D0,9.0D0,5.0D0, + 1.0D0,3.0D0,5.0D0,5.0D0,7.0D0, + 9.0D0,3.0D0,5.0D0,7.0D0,3.0D0, + 3.0D0,3.0D0,5.0D0,3.0D0,1.0D0, + 3.0D0,5.0D0,7.0D0,3.0D0,3.0D0, + 1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 5.0D0,5.0D0,7.0D0,9.0D0,3.0D0, + 5.0D0,7.0D0,3.0D0,7.0D0,3.0D0, + 5.0D0,3.0D0,1.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,5.0D0,1.0D0,5.0D0, + 5.0D0,7.0D0,9.0D0,3.0D0,5.0D0, + 7.0D0,3.0D0,7.0D0,3.0D0,5.0D0, + 3.0D0,1.0D0,5.0D0,5.0D0,7.0D0, + 9.0D0,3.0D0,5.0D0,7.0D0,3.0D0, + 7.0D0,5.0D0,3.0D0,1.0D0,3.0D0, + 5.0D0,7.0D0,9.0D0,3.0D0,5.0D0, + 7.0D0,7.0D0,3.0D0,5.0D0,3.0D0, + 1.0D0,9.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,7.0D0,5.0D0,3.0D0, + 1.0D0,9.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,7.0D0,3.0D0,5.0D0, + 7.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 3.0D0,5.0D0,7.0D0,3.0D0,3.0D0/ DATA ENCI/0.0D0,2.0333605D-03,5.3933649D-03,1.263870,2.684086, + 4.182672,7.480511,7.482891,7.487915,7.684888, + 7.946046,7.946620,7.946474,8.537387,8.640516, + 8.643146,8.647287,8.771255,8.846707,8.848247, + 8.850785,9.002712,9.171972,9.330682,9.631248, + 9.683908,9.685375,9.689256,9.695577,9.697620, + 9.701885,9.708156,9.708925,9.710041,9.712769, + 9.714380,9.761111,9.833419,9.834406,9.834934, + 9.940317,9.942698,9.946449,9.988707,10.05592, + 10.08144,10.08328,10.08553,10.13833,10.19809, + 10.35278,10.38514,10.38514,10.38514,10.39370, + 10.39456,10.39580,10.40021,10.40845,10.41874, + 10.42750,10.42990,10.42990,10.52043,10.52041, + 10.52041,10.53705,10.58840,10.61635,10.67973, + 10.70230,10.70328,10.70328,10.70878,10.70878, + 10.71184,10.71407,10.71854,10.72362,10.72523, + 10.72684,10.72684,10.86509,10.87426,10.87513, + 10.87513,10.87997,10.87997,10.88257,10.88533, + 10.88679,10.88964,10.89075,10.89075,10.88980, + 10.97789,10.97854,10.97854,10.98597,10.98597, + 10.98597,10.98808,10.98913,10.98994,10.98994, + 10.98994,11.04474,11.04474,11.04487,11.05280, + 11.05280,11.05280,11.05392,11.05429,11.05429, + 11.05429,11.09049,11.09049,11.09049,11.09843, + 11.09843,11.09843,11.09880,11.13129,11.13129, + 11.13129,11.15477,11.15477,11.15477,12.13544, + 12.83767,12.84024,12.84331,13.11772,14.86312/ DATA NCII/2,2,2,2,2,2,2,2,2,2,3,3,3,2,3,3,2,2,4,4,4,3,3,3, + 4,4,2,2,4,4,5,5,5,3,3,5,5,5,5,6,3,3,3,3,3,3,6,6, + 6,6,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, + 3,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4, + 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4, + 4,4,4,4,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5, + 5,5,5,5,5,5,6,6,6,6,6,6,6/ DATA GCII/2.0D0,4.0D0,2.0D0,4.0D0,6.0D0, + 6.0D0,4.0D0,2.0D0,2.0D0,4.0D0, + 2.0D0,2.0D0,4.0D0,4.0D0,4.0D0, + 6.0D0,6.0D0,4.0D0,2.0D0,2.0D0, + 4.0D0,2.0D0,4.0D0,6.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 2.0D0,2.0D0,4.0D0,2.0D0,4.0D0, + 4.0D0,6.0D0,6.0D0,8.0D0,2.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,6.0D0,8.0D0, + 4.0D0,2.0D0,4.0D0,6.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 10.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 4.0D0,6.0D0,6.0D0,4.0D0,2.0D0, + 6.0D0,8.0D0,4.0D0,2.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,4.0D0,6.0D0,4.0D0, + 6.0D0,8.0D0,10.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,4.0D0,6.0D0,6.0D0, + 4.0D0,2.0D0,6.0D0,8.0D0,4.0D0, + 6.0D0,8.0D0,10.0D0,6.0D0,8.0D0, + 6.0D0,8.0D0,10.0D0,12.0D0,8.0D0, + 10.0D0,8.0D0,6.0D0,4.0D0,2.0D0, + 6.0D0,4.0D0,4.0D0,2.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,6.0D0,4.0D0, + 2.0D0,6.0D0,8.0D0,4.0D0,6.0D0, + 8.0D0,10.0D0,6.0D0,8.0D0,10.0D0, + 12.0D0,8.0D0,6.0D0,4.0D0,2.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,6.0D0, + 4.0D0,2.0D0/ DATA ENCII/0.0D0,7.9350658D-03,5.331397,5.334075,5.337658, + 9.290338,9.290624,11.96386,13.71590,13.72101, + 14.44900,16.33194,16.33332,17.60895,18.04607, + 18.04625,18.65519,18.65582,19.49478,20.14995, + 20.15068,20.70119,20.70413,20.70971,20.84491, + 20.84496,20.92025,20.92256,20.95094,20.95094, + 21.49265,21.73314,21.73405,22.09347,22.13075, + 22.13075,22.13075,22.18799,22.18799,22.47211, + 22.52747,22.52929,22.53239,22.53689,22.56844, + 22.57086,22.82136,22.82136,22.85996,22.85996, + 22.89870,23.11398,23.11600,23.11878,23.38108, + 23.38522,24.12408,24.27024,24.27201,24.27444, + 24.27787,24.37010,24.37079,24.37187,24.37315, + 24.60198,24.60332,24.65351,24.65617,24.65793, + 24.78982,24.79512,25.06741,25.07039,25.98117, + 25.98415,25.98986,26.58329,26.58615,26.62689, + 26.62867,26.63139,26.63554,26.75178,26.82771, + 26.82771,26.83016,26.89454,26.89578,27.22147, + 27.22329,27.22585,27.22930,27.29263,27.29263, + 27.29378,27.29509,27.35131,27.35294,27.37703, + 27.37957,27.38104,27.41188,27.41302,27.41395, + 27.41395,27.41395,27.41409,27.46301,27.46301, + 27.46810,27.46936,27.47200,27.47561,27.47330, + 27.47864,27.48713,27.49096,27.49330,27.49330, + 27.48854,27.49412,27.55688,27.56022,27.99752, + 27.99752,27.99752,28.25640,28.25640,28.61124, + 28.61124,28.61124,28.61124,28.64683,28.64683, + 28.64683,28.66803,26.43629,28.66875,28.66875, + 28.66875,28.66875,28.70253,28.70253,28.70253, + 28.70253,28.70515,28.70515,28.70515,28.70515, + 29.31561,29.31561,29.31561,29.31561,29.33557, + 29.33557,29.33557/ DATA NCIII/2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,4, + 3,4,4,4,4,3,4,4,4,4,4,4,4,4,3,3,3,4,3,3,3,3,3,3, + 3,3,3,3,3,3,5,3,3,3,3,5,5,5,5,3,5,5,5,5,5,5,5,5, + 5,5,5,5,5,6,6,6,6,6,6,6,6,6,6,6,6,6,6,7,7,7,7,7, + 7,8,8,8,8,9,9,9,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4, + 4,4,4,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,6,6,6,6,6,6, + 6,6,6,6,6,6,7,7,7,7,7,7/ DATA GCIII/1.0D0,1.0D0,3.0D0,5.0D0,3.0D0, + 1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 3.0D0,1.0D0,3.0D0,1.0D0,3.0D0, + 5.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 1.0D0,3.0D0,5.0D0,3.0D0,3.0D0, + 1.0D0,1.0D0,3.0D0,5.0D0,3.0D0, + 3.0D0,5.0D0,7.0D0,5.0D0,7.0D0, + 9.0D0,3.0D0,7.0D0,3.0D0,5.0D0, + 7.0D0,5.0D0,3.0D0,1.0D0,3.0D0, + 5.0D0,5.0D0,5.0D0,5.0D0,7.0D0, + 9.0D0,3.0D0,5.0D0,7.0D0,3.0D0, + 5.0D0,3.0D0,1.0D0,7.0D0,3.0D0, + 1.0D0,3.0D0,5.0D0,1.0D0,3.0D0, + 5.0D0,7.0D0,7.0D0,9.0D0,11.0D0, + 9.0D0,5.0D0,3.0D0,5.0D0,7.0D0, + 9.0D0,7.0D0,3.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,7.0D0,9.0D0,11.0D0, + 9.0D0,5.0D0,5.0D0,7.0D0,9.0D0, + 7.0D0,3.0D0,3.0D0,3.0D0,5.0D0, + 7.0D0,5.0D0,3.0D0,3.0D0,5.0D0, + 7.0D0,3.0D0,5.0D0,7.0D0,1.0D0, + 3.0D0,5.0D0,3.0D0,3.0D0,5.0D0, + 7.0D0,1.0D0,3.0D0,5.0D0,5.0D0, + 5.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,7.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,1.0D0,3.0D0,5.0D0, + 5.0D0,5.0D0,3.0D0,5.0D0,7.0D0, + 5.0D0,3.0D0,1.0D0,3.0D0,5.0D0, + 7.0D0,1.0D0,3.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,5.0D0,3.0D0,1.0D0, + 3.0D0,5.0D0,7.0D0,1.0D0,3.0D0,5.0D0/ DATA ENCIII/0.0D0,6.486296,6.489148,6.496191,12.69008, + 17.03237,17.03602,17.04185,18.08638,22.62984, + 29.52845,30.64541,32.10371,32.19328,32.19396, + 32.19555,33.47080,33.45866,33.47146,34.27982, + 38.20770,38.21183,38.22034,38.36164,38.43612, + 38.64882,39.39549,39.39549,39.39611,39.64054, + 39.84380,39.84582,39.84874,39.91699,39.91782, + 39.91892,39.97328,40.01022,40.05026,40.05341, + 40.05822,40.19756,40.57121,40.86969,40.87231, + 40.87686,41.24874,41.30157,41.32848,41.33158, + 41.33611,41.85783,41.80309,41.86202,42.14028, + 42.16117,42.16444,42.16623,42.32471,42.55869, + 42.67342,42.67342,42.67342,42.78661,42.83001, + 42.83001,42.83001,42.96405,42.96405,42.96416, + 42.96405,42.98029,42.98736,43.03527,43.03550, + 43.03579,43.25349,43.98952,44.27370,44.39248, + 44.39248,44.39248,44.46592,44.46592,44.46600, + 44.47219,44.47673,44.48596,44.48596,44.48596, + 44.52591,45.07626,45.24178,45.32720,45.32720, + 45.32720,45.38200,45.86543,45.92891,45.92891, + 45.92891,46.33929,46.33929,46.33929,46.69749, + 46.69749,46.69749,47.25143,47.35238,47.35238, + 47.35722,47.64920,47.64920,47.65379,47.81342, + 47.83558,48.06245,48.06245,48.06245,48.16114, + 48.16114,48.16114,48.20208,50.51542,50.55803, + 50.55803,50.55803,50.69428,50.69428,50.69428, + 50.77264,50.79460,50.90022,50.90022,50.90022, + 50.93829,50.93829,50.93829,52.24497,52.24497, + 52.24497,52.31775,52.31775,52.31775,52.43107, + 52.43107,52.43107,52.45302,52.45302,52.45302, + 53.23251,53.23251,53.23251,53.27802,53.27802, + 53.27802/ DATA NCIV/2,2,2,3,3,3,3,3,4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,5, + 6,6,6,6,6,6,6,6,6,6,6,7,7,7,7,7,7,7,7,7,7,7,8,8, + 8,8,8,8,8,8,8/ DATA GCIV/2.0D0,2.0D0,4.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,6.0D0,8.0D0, + 2.0D0,2.0D0,4.0D0,4.0D0,6.0D0, + 6.0D0,8.0D0,8.0D0,10.0D0,2.0D0, + 2.0D0,4.0D0,4.0D0,6.0D0,6.0D0, + 8.0D0,8.0D0,10.0D0,10.0D0,12.0D0, + 2.0D0,2.0D0,4.0D0,4.0D0,6.0D0, + 6.0D0,8.0D0,8.0D0,10.0D0,10.0D0, + 12.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 8.0D0,10.0D0,12.0D0,14.0D0,16.0D0/ DATA ENCIV/0.0D0,7.995100,8.008378,37.54872,39.68134, + 39.68525,40.28040,40.28173,49.76113,50.62434, + 50.62599,50.87540,50.87595,50.88784,50.88784, + 55.21889,55.65134,55.65221,55.77947,55.77947, + 55.78577,55.78578,55.78703,55.78703,58.12002, + 58.36774,58.36774,58.44275,58.44275,58.44709, + 58.44709,58.44764,58.44764,58.44770,58.44770, + 59.84267,60.00038,60.00038,60.04725,60.04725, + 60.05156,60.05156,60.05191,60.05191,60.05194, + 60.05194,61.05946,61.05946,61.09294,61.09294, + 61.09319,61.09319,61.09319,61.09319,61.09319/ DATA NCV/1,2,2,2,2,2,3,3,3,3,4,5,6,7,8/ DATA GCV/1.0D0,3.0D0,1.0D0,3.0D0,5.0D0, + 3.0D0,3.0D0,5.0D0,7.0D0,3.0D0, + 3.0D0,3.0D0,3.0D0,3.0D0,3.0D0/ DATA ENCV/0.0D0,298.9618,304.4046,304.4030,304.4199, + 307.8855,354.2645,354.2645,354.2645,354.5177, + 370.9247,378.5349,382.6710,385.1917,386.6807/ DATA NNI/2,2,2,2,2,3,3,3,3,3,2,2,2,3,3,3,3,3,3,3,3,3,3,3, + 3,3,3,3,4,4,4,4,4,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, + 3,3,4,4,4,4,4,4,4,4,4,5,5,5,5,5,4,4,4,4,4,4,4,4, + 4,4,4,4,4,4,4,4,4,3,3,3,3,6,6,6,6,6,5,5,5,5,5,5, + 5,5,5,5,5,5,5,5,5,5,5,7,7,7,7,7,6,6,6,6,6,6,6,6, + 6,6,6,6,6,6,6,6,6,8,8,8,8,8,7,7,7,7,7,7,7,7,7,7, + 7,7,7,9,9,9,9,9,8,8,8,8,8,8,8,8,8,8,8,8,8,10,10,10, + 10,10,9,9,9,9,9,9,9,9,9,9,9,9,9,11,11,11,11,11,10, + 10,10,10,10,10,10,10,10,10,10,10,10,12,12,12,12,12, + 11,11,11,11,11,11,11,11,11,11,11,11,11,13,13,12,12, + 12,12,12,12,12/ DATA GNI/4.0D0,6.0D0,4.0D0,4.0D0,2.0D0, + 2.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 6.0D0,4.0D0,2.0D0,2.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,2.0D0,4.0D0, + 6.0D0,4.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,4.0D0,2.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,10.0D0,6.0D0, + 8.0D0,2.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,4.0D0,6.0D0, + 2.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 2.0D0,4.0D0,6.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,4.0D0, + 6.0D0,8.0D0,10.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,4.0D0,2.0D0,6.0D0, + 8.0D0,2.0D0,4.0D0,6.0D0,4.0D0, + 6.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 4.0D0,6.0D0,8.0D0,10.0D0,4.0D0, + 2.0D0,6.0D0,8.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,2.0D0,4.0D0,6.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,4.0D0,6.0D0,8.0D0, + 10.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 4.0D0,2.0D0,6.0D0,8.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,6.0D0,8.0D0, + 4.0D0,2.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,6.0D0, + 8.0D0,4.0D0,2.0D0,6.0D0,8.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,2.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,4.0D0, + 2.0D0,6.0D0,8.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,4.0D0,2.0D0, + 6.0D0,8.0D0,2.0D0,4.0D0,6.0D0, + 8.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,2.0D0,4.0D0, + 6.0D0,4.0D0,2.0D0,6.0D0,8.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,4.0D0,2.0D0,2.0D0,4.0D0, + 6.0D0,4.0D0,6.0D0/ DATA ENNI/0.0D0,2.383371,2.384363,3.575739,3.575739, + 10.32619,10.33038,10.33617,10.67904,10.69042, + 10.92429,10.92973,10.93217,11.60284,11.75037, + 11.75317,11.75780,11.76412,11.83769,11.83997, + 11.84472,11.99580,12.00032,12.00975,12.12207, + 12.12649,12.35701,12.35614,12.84713,12.85333, + 12.86185,12.91211,12.92268,12.97078,12.97568, + 12.97693,12.97929,12.98350,12.98958,12.99502, + 13.00392,13.00161,13.00483,13.00074,13.01686, + 13.01822,13.01983,13.02095,13.03344,13.03636, + 13.20179,13.23674,13.23917,13.24364,13.25041, + 13.26429,13.26623,13.27127,13.32189,13.61527, + 13.62076,13.62945,13.64202,13.65185,13.66270, + 13.66493,13.66914,13.67609,13.66580,13.67249, + 13.67410,13.68043,13.66588,13.66872,13.67695, + 13.68464,13.67869,13.68191,13.68836,13.69398, + 13.69673,13.70310,13.70607,13.92292,13.92614, + 13.95653,13.96207,13.97100,13.97749,13.98841, + 13.97948,13.98097,13.98543,13.99324,13.98568, + 13.98754,13.98803,13.99674,13.98865,13.98865, + 13.98865,13.99696,13.99237,13.99473,13.99944, + 14.00155,14.00384,14.13620,14.14326,14.15244, + 14.15045,14.15455,14.15417,14.15417,14.15417, + 14.15417,14.15690,14.15690,14.15690,14.16508, + 14.15827,14.16025,14.15864,14.16843,14.16313, + 14.17035,14.16645,14.16645,14.16831,14.23464, + 14.24468,14.25113,14.25212,14.25212,14.25683, + 14.25683,14.25683,14.25683,14.25882,14.25882, + 14.26043,14.26043,14.26545,14.27073,14.27109, + 14.27109,14.27109,14.36247,14.36247,14.31821, + 14.31821,14.31821,14.32329,14.32329,14.32329, + 14.32329,14.32403,14.32403,14.32465,14.32465, + 14.33234,14.33544,14.33494,14.33494,14.33494, + 14.36272,14.36272,14.36433,14.36433,14.36433, + 14.36830,14.36830,14.36830,14.36830,14.36854, + 14.36854,14.37016,14.37016,14.37896,14.38119, + 14.38107,14.38107,14.38107,14.39557,14.39557, + 14.39768,14.39768,14.39768,14.40152,14.40152, + 14.40202,14.40202,14.40264,14.40264,14.40264, + 14.40264,14.41206,14.41206,14.41442,14.41442, + 14.41442,14.42012,14.42012,14.42099,14.42099, + 14.42099,14.42583,14.42583,14.42682,14.42682, + 14.42781,14.42781,14.42781,14.42781,14.43636, + 14.43636,14.43698,14.43698,14.43698,14.46253, + 14.44021,14.44455,14.44455,14.45434,14.45434, + 14.45434,14.45980,14.45980/ DATA NNII/2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,2,3,3,3,3,2,3, + 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,4,4,4,3,4,4,4, + 4,4,4,4,4,4,3,3,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4, + 4,4,4,4,4,4,4,4,4,4,3,3,3,5,5,5,5,5,5,5,5,5,5,5, + 5,5,5,5,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3/ DATA GNII/1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 5.0D0,7.0D0,5.0D0,3.0D0,5.0D0, + 3.0D0,1.0D0,5.0D0,1.0D0,3.0D0, + 5.0D0,3.0D0,3.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,3.0D0,3.0D0,1.0D0, + 3.0D0,5.0D0,5.0D0,1.0D0,5.0D0, + 7.0D0,9.0D0,5.0D0,3.0D0,5.0D0, + 7.0D0,5.0D0,3.0D0,1.0D0,7.0D0, + 3.0D0,1.0D0,3.0D0,5.0D0,3.0D0, + 3.0D0,3.0D0,5.0D0,7.0D0,1.0D0, + 3.0D0,5.0D0,3.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,1.0D0,5.0D0,7.0D0, + 9.0D0,5.0D0,3.0D0,5.0D0,7.0D0, + 5.0D0,3.0D0,1.0D0,7.0D0,5.0D0, + 7.0D0,9.0D0,7.0D0,7.0D0,9.0D0, + 11.0D0,3.0D0,9.0D0,7.0D0,5.0D0, + 3.0D0,5.0D0,1.0D0,3.0D0,5.0D0, + 1.0D0,3.0D0,5.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,5.0D0,7.0D0,9.0D0, + 7.0D0,7.0D0,9.0D0,11.0D0,9.0D0, + 1.0D0,3.0D0,5.0D0,7.0D0,9.0D0, + 3.0D0,5.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,9.0D0,11.0D0,7.0D0, + 5.0D0,3.0D0,1.0D0,3.0D0,5.0D0, + 7.0D0,9.0D0/ DATA ENNII/0.0D0,6.0876831D-03,1.6279284D-02,1.898923,4.052723, + 5.848106,11.43604,11.43781,11.43801,13.54146, + 13.54146,13.54228,17.87734,18.46259,18.46651, + 18.48341,18.49722,19.23384,20.40944,20.64636, + 20.65389,20.66582,20.67651,20.94027,21.14861, + 21.15298,21.16022,21.59986,22.10340,23.12481, + 23.13218,23.14229,23.19670,23.23962,23.24260, + 23.24636,23.41565,23.42207,23.42555,23.47490, + 23.57225,24.36823,24.37465,24.38944,24.53166, + 25.06612,25.13369,25.14001,25.15193,25.18946, + 25.19245,25.20124,25.23510,25.46049,25.53877, + 25.54572,25.55447,25.58160,25.99668,26.00464, + 26.01527,26.02787,26.06667,26.06994,26.07548, + 26.12440,26.13011,26.13327,26.16475,26.16510, + 26.16800,26.16849,26.17391,26.19663,26.19758, + 26.20937,26.20252,26.21087,26.21191,26.21252, + 26.22134,26.22182,26.25393,26.25770,26.26368, + 26.55921,26.56489,26.58065,26.63554,27.36569, + 27.36569,27.36569,27.40948,27.40948,27.40999, + 27.41783,27.42901,27.42963,27.43824,27.43947, + 27.77609,27.77805,27.78169,27.78704,27.79372, + 28.01910,28.02209,28.02755,28.54429,30.17253, + 30.17448,30.17763,30.18179,30.18682,30.34387, + 30.34864,30.35188,30.41607,30.41652,30.41750, + 30.41894,30.42068/ DATA NNIII/2,2,2,2,2,2,2,2,2,2,2,2,2,3,2,2,3,3,3,3,3,3,3,3, + 3,4,3,3,3,3,3,3,4,4,3,3,3,3,4,4,4,4,3,3,3,3,3,3, + 3,3,3,3,3,5,3,3,3,3,3,3,3,5,5,3,3,5,5,5,5,6,6,6, + 6,6,6,4,4,4,3,3,4,4,4,4,4,4,3,3,4,4,4,4,4,4,4,4, + 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4, + 4,4,4,4,4,4,4,3,3,5,5,5,5/ DATA GNIII/2.0D0,4.0D0,2.0D0,4.0D0,6.0D0, + 6.0D0,4.0D0,2.0D0,2.0D0,4.0D0, + 4.0D0,6.0D0,4.0D0,2.0D0,2.0D0, + 4.0D0,2.0D0,4.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 2.0D0,2.0D0,4.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,2.0D0,4.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,4.0D0,6.0D0, + 6.0D0,8.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,10.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,2.0D0,4.0D0, + 6.0D0,6.0D0,4.0D0,2.0D0,6.0D0, + 8.0D0,4.0D0,6.0D0,4.0D0,2.0D0, + 6.0D0,8.0D0,8.0D0,10.0D0,4.0D0, + 6.0D0,6.0D0,8.0D0,8.0D0,10.0D0, + 2.0D0,4.0D0,6.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,2.0D0,4.0D0,6.0D0, + 8.0D0,2.0D0,4.0D0,4.0D0,6.0D0, + 4.0D0,2.0D0,4.0D0,6.0D0,4.0D0, + 6.0D0,8.0D0,10.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,6.0D0, + 4.0D0,2.0D0,6.0D0,8.0D0,4.0D0, + 6.0D0,8.0D0,10.0D0,6.0D0,8.0D0, + 6.0D0,8.0D0,10.0D0,12.0D0,8.0D0, + 10.0D0,8.0D0,6.0D0,4.0D0,2.0D0, + 6.0D0,4.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0/ DATA ENNIII/0.0D0,2.1635452D-02,7.180255,7.098413,7.108480, + 12.52548,12.52643,16.24252,18.08651,18.10019, + 23.16076,25.17799,25.18006,27.43827,28.56680, + 28.56730,30.45896,30.46342,33.13367,33.13441, + 35.65022,35.65797,35.67233,36.84229,36.85629, + 38.44641,38.32793,38.33453,38.39367,38.39807, + 38.40689,38.41771,38.64517,38.64825,38.95919, + 39.34056,39.34595,39.35325,39.39646,39.40031, + 39.71098,39.71098,39.79651,39.80747,40.55027, + 40.94474,40.94909,40.95552,40.96437,41.26192, + 41.26358,41.26631,41.26982,41.37555,41.47835, + 41.48166,41.68555,41.69232,41.69667,42.12335, + 42.13715,42.39634,42.39655,42.48893,42.49769, + 42.49625,42.49625,42.54757,42.54757,43.95493, + 43.95493,44.00932,44.00932,44.04135,44.04135, + 45.69180,45.69957,45.71402,46.28896,46.29317, + 46.46321,46.47039,46.71232,46.71811,46.72555, + 46.73671,46.81577,46.81788,46.85206,46.86286, + 46.92110,47.02857,47.03412,47.04068,47.61238, + 47.61238,47.61845,47.62763,47.75000,47.75000, + 47.77108,47.77108,49.01428,47.77802,47.88887, + 47.88887,47.88887,47.97657,47.97913,47.98245, + 47.98245,47.98363,47.98760,48.07270,48.08297, + 48.11119,48.11662,48.12305,48.13089,48.12993, + 48.14229,48.14024,48.14488,48.15087,48.15427, + 48.15307,48.16119,49.16950,49.17073,50.71214, + 50.71214,50.71214,50.71214/ DATA NNIV/2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3, + 3,3,3,3,3,3,4,3,3,3,3,3,4,4,4,3,3,3,3,4,4,4,4,3, + 3,3,4,4,4,4,3,4,5,5,5,5,5,5,5,6,6,6,4,4,4,4,5,5,4/ DATA GNIV/1.0D0,1.0D0,3.0D0,7.0D0,3.0D0, + 1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 3.0D0,1.0D0,1.0D0,3.0D0,5.0D0, + 3.0D0,5.0D0,7.0D0,5.0D0,1.0D0, + 3.0D0,5.0D0,3.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,3.0D0,1.0D0,3.0D0, + 5.0D0,5.0D0,5.0D0,5.0D0,7.0D0, + 9.0D0,1.0D0,3.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,7.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,5.0D0,3.0D0,1.0D0, + 5.0D0,5.0D0,7.0D0,9.0D0,3.0D0, + 7.0D0,3.0D0,3.0D0,5.0D0,7.0D0, + 7.0D0,9.0D0,11.0D0,3.0D0,5.0D0, + 7.0D0,5.0D0,3.0D0,5.0D0,7.0D0, + 3.0D0,5.0D0,7.0D0/ DATA ENNIV/0.0D0,8.323934,8.331770,8.349648,16.20427, + 21.75491,21.76399,21.77946,23.41898,29.18244, + 46.76804,50.15470,50.32483,50.32679,50.33118, + 52.06988,52.07031,52.07132,53.20933,57.68086, + 57.69048,57.71067,58.64906,59.62210,60.05779, + 60.05779,60.07403,60.44809,61.27855,61.27855, + 61.29070,61.78379,61.95650,61.97423,61.97423, + 61.97423,62.44215,62.44215,62.44215,62.67301, + 62.67685,62.68218,62.77282,62.86333,63.40415, + 63.40415,63.40415,63.41109,63.41767,63.41767, + 63.80760,64.05482,64.05569,64.05706,64.39976, + 64.70402,68.21900,68.53058,68.53058,68.53058, + 68.73986,68.73986,68.73986,71.28416,71.28416, + 71.28416,73.28070,73.60580,73.60580,73.61063, + 78.63129,78.63129,78.63129/ DATA NNV/2,2,2,3,3,3,3,3,4,4,4,4,4,5,5,5,5,5,6,6,6,6,6,6, + 6,6,6,6,7,7,7,7,7,7,7,7,7,7,7,8,8,8,8,8,8,8,8,8,8,8,8/ DATA GNV/2.0D0,2.0D0,4.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,6.0D0,8.0D0, + 8.0D0,10.0D0,12.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,6.0D0,8.0D0, + 8.0D0,10.0D0,12.0D0,14.0D0,2.0D0, + 2.0D0,4.0D0,4.0D0,6.0D0,6.0D0, + 8.0D0,8.0D0,10.0D0,12.0D0,14.0D0,16.0D0/ DATA ENNV/0.0D0,9.976473,10.00851,56.55396,59.23740, + 59.24660,60.05890,60.06188,75.17694,76.26962, + 76.26962,76.61120,76.61120,83.55153,84.09893, + 84.09893,84.27598,84.27598,88.02306,88.33514, + 88.33514,88.43854,88.43742,88.44214,88.44214, + 88.44313,88.44313,88.44313,90.68689,90.88043, + 90.88043,90.94527,90.94527,90.94912,90.94912, + 90.94974,90.94974,90.94974,90.94974,92.40136, + 92.53167,92.53167,92.57358,92.57358,92.57618, + 92.57618,92.57668,92.57668,92.57668,92.57668, + 92.57668/ DATA NNVI/1,2,2,2,2,2,3,4/ DATA GNVI/1.0D0,3.0D0,1.0D0,3.0D0,5.0D0,3.0D0,3.0D0,3.0D0/ DATA ENNVI/0.0D0,419.8009,426.2953,426.2965,426.3325, + 425.7398,497.9737,521.5830/ DATA NOI/2,2,2,2,2,3,3,3,3,3,3,3,3,4,4,3,3,3,3,3,3,3,3,3, + 4,4,4,4,4,4,3,3,3,5,5,3,4,4,4,4,4,4,4,4,5,5,5,6, + 6,5,5,5,5,5,5,5,5,6,6,6,7,7,6,6,6,6,6,6,6,6,8,8, + 7,7,7,7,7,7,7,7,9,9,8,8,8,8,8,8,8,8,10,10,9,9,9,9, + 9,9,9,9,11,11,10,10,10,10,10,10,10,10,3,3,3,3,3,3, + 3,3,3,3,3,3,4,3,3,3,3,3,3,3,3,3,3,3,4,4,4,2,2,2,3, + 3,3,3,3,5,4,4,4,4,4,4,4,4,4,4,4,3,6,5,5,5,5,5,5,5, + 5,5,5,7,6,6,6,2/ DATA GOI/5.0D0,3.0D0,1.0D0,5.0D0,1.0D0, + 5.0D0,3.0D0,3.0D0,5.0D0,7.0D0, + 5.0D0,3.0D0,1.0D0,5.0D0,3.0D0, + 9.0D0,7.0D0,5.0D0,5.0D0,3.0D0, + 1.0D0,7.0D0,5.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,5.0D0,3.0D0,1.0D0, + 7.0D0,5.0D0,3.0D0,5.0D0,3.0D0, + 5.0D0,9.0D0,7.0D0,5.0D0,3.0D0, + 1.0D0,7.0D0,5.0D0,3.0D0,5.0D0, + 3.0D0,1.0D0,5.0D0,3.0D0,9.0D0, + 7.0D0,5.0D0,3.0D0,1.0D0,7.0D0, + 5.0D0,3.0D0,5.0D0,3.0D0,1.0D0, + 5.0D0,3.0D0,9.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,3.0D0,9.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,3.0D0,9.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,3.0D0,9.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,3.0D0,9.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,7.0D0,5.0D0,3.0D0, + 7.0D0,5.0D0,3.0D0,9.0D0,7.0D0, + 5.0D0,5.0D0,3.0D0,1.0D0,7.0D0, + 3.0D0,5.0D0,5.0D0,5.0D0,3.0D0, + 1.0D0,9.0D0,7.0D0,5.0D0,9.0D0, + 11.0D0,9.0D0,7.0D0,7.0D0,7.0D0, + 5.0D0,3.0D0,5.0D0,3.0D0,1.0D0, + 7.0D0,5.0D0,3.0D0,3.0D0,5.0D0, + 5.0D0,9.0D0,7.0D0,5.0D0,9.0D0, + 11.0D0,9.0D0,7.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,1.0D0,5.0D0,9.0D0, + 7.0D0,5.0D0,9.0D0,11.0D0,9.0D0, + 7.0D0,5.0D0,3.0D0,1.0D0,5.0D0, + 5.0D0,3.0D0,1.0D0,3.0D0/ DATA ENOI/0.0D0,01.9651687D-02,2.8082693D-02,1.967363,0.4206081, + 9.146132,9.521420,10.74028,10.74053,10.74098, + 10.98893,10.98886,10.98895,11.83768,11.93056, + 12.07869,12.07870,12.07870,12.07872,12.07872, + 12.07872,12.08711,12.08711,12.08711,12.28604, + 12.28612,12.28627,12.35891,12.35891,12.35891, + 12.53927,12.54078,12.54176,12.66092,12.69755, + 12.72854,12.75377,12.75377,12.75377,12.75377, + 12.75377,12.75911,12.75911,12.75911,12.87829, + 12.87829,12.87829,13.02082,13.03891,13.06624, + 13.06624,13.06624,13.06624,13.06624,13.06913, + 13.06913,13.06913,13.13145,13.13145,13.13145, + 13.21004,13.22030,13.23559,13.23559,13.23559, + 13.23559,13.23559,13.23740,13.23740,13.23740, + 13.32166,13.32807,13.33749,13.33749,13.33749, + 13.33749,13.33749,13.33869,13.33869,13.33869, + 13.39308,13.39756,13.40353,385.3597,13.40353, + 13.40353,13.40353,13.40488,13.40488,13.40488, + 13.44262,13.44449,13.44872,13.44872,13.44872, + 13.44872,13.44872,13.44966,13.44966,13.44966, + 13.47577,13.47812,13.48112,13.48112,13.48112, + 13.48112,13.48112,13.48148,13.48148,13.48148, + 14.04685,14.04687,14.04730,14.09888,14.09975, + 14.10046,14.12320,14.12450,14.12526,14.13382, + 14.37218,14.46048,15.22525,15.28698,15.29424, + 15.29817,15.40062,15.40062,15.40062,15.40372, + 15.40390,15.40622,15.40550,15.41465,15.59420, + 15.59514,15.59577,15.65520,15.66431,15.66970, + 15.78109,15.78181,15.78222,15.82895,15.94391, + 16.01073,16.07676,16.07676,16.07676,16.07836, + 16.07844,16.08080,16.08005,16.08545,16.11433, + 16.11550,16.11614,16.23505,16.35702,16.35702, + 16.35702,16.35702,16.39057,16.39063,16.39308, + 16.39308,16.40451,16.40451,16.40451,16.54127, + 16.56668,16.56668,16.56668,23.53702/ DATA NOII/2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,2,3,3,3,3,3,3,3,3, + 3,3,3,3,3,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, + 3,3,3,3,3,3,3,3,3,4,4,4,4,4,3,4,4,4,4,4,4,4,4,3, + 3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,3,4,4,4,4,4,4, + 4,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,5,5, + 5,5,5,4,4,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5, + 5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,3,3,3,4,4,4,4, + 4,4,4,4,4,4,3,3,4,4,4,4,4,4,4,5,5,3,3,3,3,3,4/ DATA GOII/4.0D0,6.0D0,4.0D0,4.0D0,2.0D0, + 6.0D0,4.0D0,2.0D0,6.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 2.0D0,2.0D0,2.0D0,4.0D0,6.0D0, + 8.0D0,6.0D0,4.0D0,2.0D0,4.0D0, + 6.0D0,4.0D0,6.0D0,4.0D0,4.0D0, + 2.0D0,2.0D0,4.0D0,2.0D0,6.0D0, + 8.0D0,6.0D0,4.0D0,4.0D0,6.0D0, + 8.0D0,10.0D0,6.0D0,4.0D0,2.0D0, + 2.0D0,4.0D0,2.0D0,4.0D0,6.0D0, + 8.0D0,6.0D0,8.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,2.0D0,4.0D0,8.0D0,6.0D0, + 10.0D0,8.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,8.0D0,10.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,4.0D0, + 6.0D0,4.0D0,2.0D0,4.0D0,2.0D0, + 6.0D0,8.0D0,2.0D0,6.0D0,4.0D0, + 8.0D0,5.80D0,4.0D0,2.0D0,6.0D0, + 8.0D0,10.0D0,12.0D0,8.0D0,10.0D0, + 4.0D0,6.0D0,4.0D0,6.0D0,8.0D0, + 10.0D0,6.0D0,8.0D0,2.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,6.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,2.0D0, + 4.0D0,6.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,8.0D0,6.0D0,4.0D0, + 2.0D0,6.0D0,8.0D0,8.0D0,6.0D0, + 4.0D0,2.0D0,6.0D0,8.0D0,10.0D0, + 12.0D0,8.0D0,10.0D0,4.0D0,6.0D0, + 4.0D0,6.0D0,8.0D0,10.0D0,6.0D0, + 8.0D0,4.0D0,6.0D0,8.0D0,6.0D0, + 8.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 8.0D0,10.0D0,6.0D0,8.0D0,6.0D0, + 4.0D0,2.0D0,4.0D0,6.0D0,10.0D0, + 12.0D0,2.0D0,4.0D0,4.0D0,6.20D0, + 10.0D0,8.0D0,6.0D0,4.0D0,2.0D0, + 6.0D0/ DATA ENOII/0.0D0,3.323850,3.326454,5.017305,5.017491, + 14.85813,14.87838,14.88860,20.58005,20.57736, + 22.96648,22.97954,23.001876,23.41940,23.44172, + 24.26523,25.28586,25.63160,25.63849,25.64984, + 25.66529,25.66142,25.66154,25.83188,25.83760, + 25.84900,26.22564,26.24928,26.30498,26.35845, + 26.37943,26.55392,26.56133,28.12621,28.35835, + 28.36128,28.51330,28.51270,28.67733,28.68403, + 28.69369,28.70637,28.82200,28.83108,28.83932, + 28.82414,28.82992,28.85285,28.85711,28.85729, + 28.85808,28.86334,28.88355,28.94193,28.95606, + 29.06249,29.06893,29.58618,29.59923,29.61924, + 29.79726,29.82051,30.42546,30.47162,30.47763, + 30.48836,30.50400,30.74951,30.77135,30.80112, + 30.81214,31.02747,31.02747,31.14773,31.14812, + 31.31967,31.31982,31.37404,31.37430,31.46620, + 31.46649,31.55199,31.55199,31.55199,31.56553, + 31.61407,31.61407,31.61407,31.61407,31.61407, + 31.62925,31.63375,31.63644,31.63766,31.65117, + 31.65364,31.67396,31.69345,31.70178,31.71699, + 31.70200,31.71709,31.72948,31.72935,31.70999, + 31.71043,31.71889,31.73747,31.71911,31.73823, + 31.72081,31.72752,31.75062,31.75112,31.75553, + 31.75715,31.75586,31.75803,31.95026,31.96318, + 31.98375,32.03889,32.06284,32.14771,32.14780, + 32.35511,32.35511,32.36540,32.38251,32.39264, + 32.39264,32.40412,32.44667,32.46798,32.88345, + 32.88345,32.88345,32.88345,32.90963,32.91418, + 32.91418,32.92780,32.92780,32.93536,32.94354, + 32.95061,32.96264,32.93858,32.94181,32.95049, + 32.97082,32.95073,32.97146,32.96227,32.96227, + 32.97119,32.97528,32.97826,32.97999,32.97863, + 32.97999,33.19875,33.19968,33.20123,34.06365, + 34.06901,34.08607,34.08607,34.17174,34.17174, + 34.20029,34.20029,34.20504,34.20504,34.21390, + 34.21390,34.21960,34.22819,34.22819,34.23350, + 34.23350,34.25269,34.25269,34.48530,34.48530, + 36.19083,36.18759,36.19109,36.19123,36.19131,37.05294/ DATA NOIII/2,2,2,2,2,2,2,3,3,2,2,2,2,2,2,3,3,3,3,2,2,2,3,3, + 3,3,3,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, + 2,3,3,3,4,4,4,4,3,3,3,3,3,3,4,4,4,4,4,3,3,3,4,4, + 4,4,4,3,3,3,3,4,4,4,4,3,3,3,4,4,4,4,4,4,4,4,5,5, + 5,5,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3, + 5,5,5,5,5,5,5,5,5,3,3,3,6,6,6,6,7,3,3,4,4,4,3,4, + 4,4,4,4,4,4,4,4,4,4,4,4,4,4,3,3,3,3,3,3,3,3,3,5/ DATA GOIII/1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 5.0D0,7.0D0,5.0D0,3.0D0,5.0D0, + 3.0D0,1.0D0,5.0D0,3.0D0,3.0D0, + 1.0D0,3.0D0,5.0D0,3.0D0,5.0D0, + 3.0D0,1.0D0,3.0D0,3.0D0,5.0D0, + 7.0D0,3.0D0,5.0D0,1.0D0,3.0D0, + 5.0D0,5.0D0,1.0D0,5.0D0,7.0D0, + 9.0D0,5.0D0,3.0D0,5.0D0,7.0D0, + 5.0D0,3.0D0,1.0D0,7.0D0,3.0D0, + 3.0D0,5.0D0,7.0D0,1.0D0,1.0D0, + 3.0D0,5.0D0,1.0D0,3.0D0,5.0D0, + 3.0D0,3.0D0,1.0D0,3.0D0,5.0D0, + 7.0D0,9.0D0,3.0D0,3.0D0,5.0D0, + 7.0D0,3.0D0,3.0D0,5.0D0,7.0D0, + 1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 3.0D0,5.0D0,7.0D0,5.0D0,5.0D0, + 7.0D0,9.0D0,5.0D0,5.0D0,3.0D0, + 1.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 3.0D0,1.0D0,7.0D0,3.0D0,1.0D0, + 3.0D0,5.0D0,3.0D0,3.0D0,5.0D0, + 7.0D0,3.0D0,5.0D0,7.0D0,9.0D0, + 11.0D0,1.0D0,3.0D0,5.0D0,7.0D0, + 9.0D0,7.0D0,5.0D0,3.0D0,5.0D0, + 3.0D0,1.0D0,5.0D0,7.0D0,9.0D0, + 5.0D0,7.0D0,9.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,7.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,5.0D0,3.0D0,5.0D0, + 7.0D0,7.0D0,7.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,3.0D0,3.0D0,1.0D0, + 3.0D0,5.0D0,7.0D0,9.0D0,3.0D0, + 5.0D0,7.0D0,3.0D0,5.0D0,7.0D0, + 7.0D0,5.0D0,3.0D0,5.0D0,7.0D0, + 9.0D0,3.0D0,5.0D0,7.0D0,1.0D0, + 3.0D0,5.0D0,3.0D0/ DATA ENOIII/0.0D0,1.4059945D-02,3.8038719D-02,2.513308,5.354124, + 7.477820,14.88140,14.88477,14.88550,17.65325, + 17.65339,17.65514,23.19140,24.43587,26.09378, + 33.13600,33.15068,33.18253,33.85794,35.18196, + 35.20895,35.22094,36.07438,36.43500,36.45190, + 36.47919,36.89279,36.98353,37.22392,37.23410, + 37.25028,38.01204,38.90675,40.22861,40.25288, + 40.27497,40.26230,40.57149,40.57759,40.58673, + 40.84922,40.86335,40.87098,41.14086,41.25951, + 41.97723,41.99266,42.14902,42.56451,43.39812, + 43.41013,43.43237,44.22956,44.24270,44.27655, + 44.46952,45.03978,45.31862,45.32294,45.33144, + 45.34384,45.35962,45.34443,45.43903,45.45230, + 45.47797,45.62070,45.69189,45.69899,45.71153, + 45.91510,45.92614,45.93959,45.98626,46.25228, + 46.44183,45.21283,46.46955,46.62690,46.78899, + 46.78899,46.78899,46.82767,46.91713,46.91867, + 46.92080,47.01923,47.02679,47.03461,47.20199, + 47.20199,47.20199,47.21141,47.24910,48.62968, + 48.62968,48.62968,48.69874,48.86141,48.86587, + 48.87442,48.91428,48.91908,48.92621,48.93560, + 48.94701,49.36293,49.36248,49.36198,49.36323, + 49.37332,49.40500,49.41368,49.41845,49.63815, + 49.65178,49.65844,49.76514,49.77709,49.79367, + 49.78386,49.78386,49.78386,49.81572,49.78386, + 49.78386,49.78386,50.01249,50.03133,50.31391, + 50.31750,50.32357,51.41365,51.47638,51.47638, + 51.47638,52.44297,52.69355,52.85969,53.12613, + 53.14089,53.16110,53.31682,54.18348,54.33549, + 54.33549,54.34320,54.35460,54.36977,54.46407, + 54.47044,54.48261,54.88958,54.88958,54.88958, + 55.81414,55.82281,55.82951,56.14741,56.14741, + 56.14741,56.31095,56.31095,56.31095,56.73994, + 56.73994,56.73994,58.73808/ DATA NOIV/2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3, + 3,3,3,3,3,3,3,3,3,3,3,4,3,3,3,3,3,3,3,3,3,3,3,3, + 3,3,4,4,3,3,3,3,3,3,5,3,3,3,3,5,5,5,5,3,4,4,4,3, + 3,4,4,6,6,4,4,3,3,3,3,3,3,3,4,4,7,7,4,4,4,4,4,4, + 4,4,4,4,4,4,4,4,4,4,3,8,8,4,4,3,3,3,3,3,3,3,3,3, + 3,3,3,3,3,5,5,3,3,5,5,3,3,5,5,5,5,5,5,5,5,5,5,5, + 3,3,3,3,3,3,3,3,3,6,6,6,6,4,4,3,4,4,7,7,7,7/ DATA GOIV/2.0D0,4.0D0,2.0D0,4.0D0,6.0D0, + 6.0D0,4.0D0,2.0D0,2.0D0,6.0D0, + 4.0D0,6.0D0,4.0D0,2.0D0,4.0D0, + 2.0D0,2.0D0,4.0D0,2.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,2.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,4.0D0,6.0D0, + 2.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 10.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 4.0D0,6.0D0,6.0D0,4.0D0,2.0D0, + 4.0D0,6.0D0,6.0D0,8.0D0,4.0D0, + 2.0D0,2.0D0,4.0D0,2.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,4.0D0,6.0D0, + 6.0D0,8.0D0,2.0D0,2.0D0,4.0D0, + 6.0D0,6.0D0,8.0D0,2.0D0,4.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,4.0D0, + 6.0D0,2.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,6.0D0,8.0D0, + 2.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 6.0D0,4.0D0,2.0D0,4.0D0,6.0D0, + 6.0D0,8.0D0,4.0D0,6.0D0,6.0D0, + 8.0D0,2.0D0,6.0D0,8.0D0,4.0D0, + 2.0D0,4.0D0,6.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,2.0D0,4.0D0,6.0D0, + 6.0D0,4.0D0,4.0D0,6.0D0,8.0D0, + 2.0D0,4.0D0,6.0D0,8.0D0,4.0D0, + 6.0D0,6.0D0,4.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,6.0D0,4.0D0,2.0D0, + 6.0D0,8.0D0,4.0D0,2.0D0,6.0D0, + 4.0D0,2.0D0,4.0D0,6.0D0,6.0D0, + 8.0D0,4.0D0,2.0D0,2.0D0,4.0D0, + 6.0D0,8.0D0,4.0D0,6.0D0,2.0D0, + 4.0D0,6.0D0,2.0D0,4.0D0,6.0D0,8.0D0/ DATA ENOIV/0.0D0,4.7920357D-02,8.824909,8.841201,8.864076, + 15.73825,15.73998,20.37910,22.37705,22.40721, + 28.67474,31.63571,31.63934,35.83378,35.83476, + 44.33902,48.37428,48.38508,54.37857,54.39532, + 54.42593,56.14158,56.17444,57.92984,57.94415, + 58.03452,58.04428,58.06108,58.08709,58.79609, + 59.33789,59.34961,59.36561,59.84372,59.87542, + 60.23497,61.10992,61.36131,61.37108,61.38501, + 61.40412,61.93150,61.93509,61.94088,61.94888, + 62.18008,62.18691,62.46812,62.48219,62.49133, + 63.30199,63.30286,63.32506,63.35387,63.75540, + 63.77412,64.30924,64.30999,66.87376,67.85857, + 67.86167,68.16618,68.17400,68.44416,68.44416, + 68.50069,68.50069,68.74507,70.50282,70.51955, + 70.55017,70.76975,70.76975,71.12993,71.15609, + 71.21387,71.21387,71.31690,71.33785,71.39315, + 71.39737,71.48887,71.50672,71.53300,72.12492, + 72.12764,72.47591,72.50269,72.88482,72.88482, + 73.16019,73.37047,73.37047,73.37047,73.37047, + 73.52322,73.52322,73.52322,73.60108,73.61112, + 73.64819,73.65725,73.68911,73.71453,73.93237, + 73.95444,74.05078,74.06293,74.06293,74.10930, + 74.12628,74.40265,74.40438,74.76035,74.76035, + 74.76035,74.76035,75.18896,75.18896,75.18896, + 76.30446,76.30806,76.44791,77.47625,77.47625, + 77.92433,77.92433,78.12258,78.12258,78.19797, + 78.21979,78.41159,78.43242,78.59385,78.59385, + 78.59385,78.59385,78.63718,78.63718,78.63718, + 78.85769,78.88398,78.91572,78.91572,78.96023, + 78.97250,78.98019,80.20107,80.20107,80.72665, + 80.72900,81.00314,81.01343,81.37509,81.37509, + 81.37509,81.37509,81.42716,81.42716,81.83012, + 82.88895,82.88895,83.03365,83.03365,83.03365,83.03365/ DATA NOV/2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3, + 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,4,4,4,4,4, + 4,4,4,4,4,4,5,5,5,5,5,5,4,4,4,4,4,4,4,4,4,4,4,6, + 6,6,6,6,4,4,4,6,4,4,4,4,4,7,7,7,7,7,8,8,8,8,5,5, + 5,5,5,5,5,5,5,5,5,5,5,6,6,6,6,6,6,6,6/ DATA GOV/1.0D0,1.0D0,3.0D0,5.0D0,3.0D0, + 1.0D0,3.0D0,5.0D0,5.0D0,1.0D0, + 3.0D0,1.0D0,3.0D0,1.0D0,3.0D0, + 5.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 1.0D0,3.0D0,5.0D0,3.0D0,3.0D0, + 3.0D0,5.0D0,7.0D0,3.0D0,1.0D0, + 3.0D0,5.0D0,5.0D0,5.0D0,3.0D0, + 5.0D0,7.0D0,1.0D0,5.0D0,3.0D0, + 1.0D0,7.0D0,3.0D0,3.0D0,1.0D0, + 1.0D0,3.0D0,5.0D0,3.0D0,3.0D0, + 5.0D0,7.0D0,5.0D0,7.0D0,3.0D0, + 3.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 3.0D0,3.0D0,3.0D0,5.0D0,7.0D0, + 3.0D0,1.0D0,3.0D0,5.0D0,5.0D0, + 5.0D0,3.0D0,7.0D0,3.0D0,5.0D0, + 7.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 5.0D0,3.0D0,1.0D0,7.0D0,3.0D0, + 3.0D0,3.0D0,5.0D0,7.0D0,5.0D0, + 3.0D0,3.0D0,5.0D0,7.0D0,3.0D0, + 3.0D0,5.0D0,7.0D0,1.0D0,3.0D0, + 5.0D0,5.0D0,5.0D0,3.0D0,5.0D0, + 7.0D0,7.0D0,3.0D0,3.0D0,5.0D0, + 7.0D0,1.0D0,3.0D0,5.0D0,5.0D0/ DATA ENOV/0.0D0,10.18183,10.19878,10.23674,19.68863, + 26.48845,26.50776,26.54108,28.73015,35.69651, + 67.83862,69.59028,72.01395,72.28146,72.28596, + 72.29554,74.50599,74.50733,74.50979,75.95557, + 80.97483,80.99497,81.03748,82.38657,83.40436, + 83.97941,84.00407,84.04314,84.82139,85.49855, + 85.51269,85.53633,86.12596,86.43890,87.33036, + 87.33829,87.35107,87.73579,87.80076,87.81837, + 87.82866,88.39750,89.17985,89.60004,90.71603, + 91.26665,91.26665,91.26888,91.48672,92.04689, + 92.04763,92.04937,92.59937,92.97132,98.72523, + 99.49233,106.7049,106.7049,106.7049,100.2237, + 102.1987,102.8568,103.0377,103.0583,103.0944, + 103.1870,103.5465,103.5465,103.5676,103.8792, + 103.8829,104.1001,104.2509,104.2990,104.2990, + 104.2990,104.3064,104.3181,104.3333,104.4087, + 104.5556,104.5689,104.5754,105.0316,105.0733, + 106.7358,106.9963,106.9963,106.9963,106.9274, + 108.4187,108.5325,108.5325,108.5325,111.4108, + 111.5461,111.5461,111.5461,111.7535,111.7535, + 111.7535,111.8898,111.9082,112.1444,112.1444, + 112.1444,112.3809,115.9379,116.0435,116.0435, + 116.0435,116.1501,116.1501,116.1501,116.2166/ DATA NOVI/2,2,2,3,3,3,3,3,4,4,4,4,4,4,4,5,5,5,5,5,6,6,6,6, + 6,6,6,6,6,7,7,7,7,7,7,7,7,7,7,7,8,8,8,8,8,8,8,8,8,8,8,8/ DATA GOVI/2.0D0,2.0D0,4.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,2.0D0,2.0D0, + 4.0D0,4.0D0,6.0D0,6.0D0,8.0D0, + 2.0D0,2.0D0,4.0D0,4.0D0,6.0D0, + 2.0D0,4.0D0,4.0D0,6.0D0,6.0D0, + 8.0D0,8.0D0,10.0D0,12.0D0,2.0D0, + 2.0D0,4.0D0,4.0D0,6.0D0,6.0D0, + 8.0D0,8.0D0,10.0D0,12.0D0,14.0D0, + 2.0D0,2.0D0,4.0D0,6.0D0,8.0D0, + 8.0D0,10.0D0,12.0D0,14.0D0,16.0D0,4.0D0,6.0D0/ DATA ENOVI/0.0D0,11.94909,12.01505,79.35559,82.58831, + 82.60773,83.64374,83.65008,105.7219,107.0408, + 107.0487,107.4805,107.4831,107.5050,107.5062, + 117.6237,118.2920,118.2920,118.5122,118.5122, + 124.3735,124.3735,124.5034,124.5034,124.5142, + 124.5142,124.5156,124.5156,124.5156,127.8017, + 128.0311,128.0311,128.1171,128.1171,128.1243, + 128.1243,128.1252,128.1252,128.1252,128.1252, + 130.2520,130.3984,130.3984,130.4674,130.4674, + 130.4680,130.4680,130.4680,130.4680,130.4680, + 130.4693,130.4693/ DATA NOVII/1,2,2,2,2,2,3,3,3,3,3,3,3,4,5,6/ DATA GOVII/1.0D0,3.0D0,1.0D0,3.0D0,5.0D0, + 3.0D0,1.0D0,3.0D0,5.0D0,7.0D0, + 5.0D0,3.0D0,3.0D0,3.0D0,3.0D0,3.0D0/ DATA ENOVII/0.0D0,561.0761,568.6182,568.6255,568.6938, + 573.9532,664.1129,664.1129,664.1129,665.1804, + 665.1804,665.1804,665.6218,697.8022,712.7239,720.8449/ DATA SCI/4.179704,4.179868,4.180140,4.284864,4.411317, + 4.556712,4.417538,4.418036,4.419087,4.460873, + 5.012059,5.012145,5.012123,4.656621,4.682271, + 4.682931,4.683973,4.715525,4.735114,4.735517, + 4.736181,4.776610,4.823287,5.245868,4.960446, + 4.636463,4.637096,4.638772,4.981131,4.981795, + 4.983181,4.985224,4.985476,4.985839,4.648973, + 4.987257,5.002644,5.026932,5.027268,5.027448, + 4.751991,4.753114,4.754886,4.775015,4.807735, + 4.820394,4.821310,4.822433,4.849118,4.880081, + 4.964531,4.983084,4.983084,4.983084,4.988046, + 4.988550,4.989272,4.739796,4.996660,5.002712, + 5.007890,5.009317,5.009317,4.830776,4.830763, + 4.830763,4.843919,4.885500,4.908801,4.963564, + 4.983777,4.984663,4.984663,4.989659,4.989659, + 4.992449,4.793379,4.998577,5.003253,5.004741, + 5.006231,5.006231,4.972326,4.984214,4.985345, + 4.985345,4.991673,4.991673,4.995098,4.831870, + 5.000666,5.004451,5.005935,5.005935,5.004665, + 4.984616,4.985763,4.985763,4.999064,4.999064, + 4.999064,5.002865,5.004758,5.006233,5.006233, + 5.006233,4.984145,4.984145,4.984432,5.002990, + 5.002990,5.002990,5.005628,5.006508,5.006508, + 5.006508,4.983365,4.983365,4.983365,5.006886, + 5.006886,5.006886,5.008002,5.012074,5.012074, + 5.012074,5.014099,5.014099,5.014099,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000/ DATA SCII/3.322208,3.322644,3.633090,3.633257,3.633480, + 3.893420,3.893440,4.089184,4.229172,4.229597, + 3.436720,3.692591,3.692789,4.589103,3.953149, + 3.953178,4.702748,4.702820,3.603431,3.770060, + 3.770254,4.440431,4.441056,4.442240,3.961645, + 3.961659,4.991754,4.992091,3.992479,3.992479, + 3.697579,3.795690,3.796068,4.770876,4.780956, + 3.968260,3.968260,3.994325,3.994325,3.754887, + 4.893884,4.894430,4.895358,4.896706,4.906213, + 4.906944,3.971236,3.971236,3.996578,3.996578, + 5.011171,5.086056,5.086787,5.087796,5.188515, + 5.190205,5.591661,5.735427,5.737656,5.740739, + 5.745144,5.937513,5.941305,5.947732,5.956571, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000/ DATA SCIII/2.247678,2.511178,2.511298,2.511595,2.783334, + 2.988424,2.988601,2.988887,3.040348,3.275480, + 2.516355,2.624130,2.770250,2.779440,2.779510, + 2.779673,2.913505,2.912204,2.913576,3.001503, + 3.471912,3.472452,3.473566,2.656192,3.501991, + 2.707107,2.843331,2.843331,2.843447,3.667003, + 2.928022,2.928408,2.928967,2.942070,2.942230, + 2.942442,2.952919,2.960061,3.725865,3.726323, + 3.727023,2.996535,3.802981,3.848412,3.848814, + 3.849514,3.907525,3.915897,3.920174,3.920667, + 3.921389,4.006182,3.997117,4.006877,2.756046, + 4.057183,4.057739,4.058045,4.085242,2.876864, + 2.910818,2.910818,2.910818,4.166811,2.957773, + 2.957773,2.957773,2.998549,2.998549,2.998583, + 2.998549,3.003525,3.005696,3.020441,3.020510, + 3.020601,3.088544,2.797248,2.916939,2.968366, + 2.968366,2.968366,3.000603,3.000603,3.000640, + 3.003372,3.005378,3.009464,3.009464,3.009464, + 3.027199,2.830505,2.926040,2.976524,2.976524, + 2.976524,3.009360,2.932986,2.982088,2.982088, + 2.982088,2.986294,2.986294,2.986294,4.828427, + 4.828427,4.828427,5.151012,5.224114,5.224114, + 5.227788,5.497263,5.497263,5.502663,5.756044, + 5.817122,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000,6.000000,6.000000,6.000000,6.000000, + 6.000000/ DATA SCIV/1.644934,1.923884,1.924364,1.778341,1.948965, + 1.949284,1.998202,1.998312,1.838941,1.962835, + 1.963075,1.999589,1.999669,2.001418,2.001418, + 1.874532,1.972046,1.972243,2.001394,2.001394, + 2.002842,2.002845,2.003133,2.003133,1.897882, + 1.978616,1.978616,2.003384,2.003384,2.004821, + 2.004821,2.005002,2.005002,2.005023,2.005023, + 1.913888,1.984033,1.984033,2.005113,2.005113, + 2.007061,2.007061,2.007217,2.007217,2.007234, + 2.007234,1.989961,1.989961,2.009654,2.009654, + 2.009800,2.009800,2.009800,2.009800,2.009800/ DATASCV/0.6309066,0.7688928,0.9242349,0.9241881,0.9246764, + 1.026124,1.003322,1.003322,1.003322,1.020118, + 1.021842,1.027042,1.033988,1.051912,1.003001/ DATA SNI/4.931870,5.108953,5.109031,5.204087,5.204087, + 5.329969,5.330800,5.331948,5.401419,5.403777, + 5.968683,5.969459,5.969806,5.605717,5.641185, + 5.641868,5.642995,5.644538,5.662621,5.663186, + 5.664362,5.702335,5.703489,5.705897,5.734946, + 5.736104,5.797976,5.797737,5.588642,5.591228, + 5.594790,5.615995,5.620492,5.980873,5.982464, + 5.982872,5.983638,5.985012,5.986994,5.988774, + 5.991692,5.990931,5.991989,5.990646,5.995945, + 5.996395,5.996926,5.997295,6.001428,6.002394, + 5.745163,5.761658,5.762813,5.764937,5.768167, + 5.774817,5.775745,5.778173,5.802792,5.696104, + 5.699982,5.706142,5.715098,5.722151,5.983984, + 5.985277,5.987724,5.991767,5.985780,5.989671, + 5.990610,5.994303,5.985830,5.987479,5.992273, + 5.996771,5.993288,5.995173,5.998955,6.002261, + 6.003886,6.255743,6.257061,6.360914,6.362586, + 5.757127,5.763044,5.772635,5.779661,5.791555, + 5.984846,5.986194,5.990250,5.997386,5.990476, + 5.992170,5.992623,6.000597,5.993189,5.993189, + 5.993189,6.000802,5.996590,5.998751,6.003087, + 6.005032,6.007153,5.793718,5.804320,5.818227, + 5.815205,5.821445,5.989322,5.989322,5.989322, + 5.989322,5.992901,5.992901,5.992901,6.003714, + 5.994695,5.997311,5.995185,6.008173,6.001116, + 6.010741,6.005529,6.005529,6.008008,5.801159, + 5.821038,5.833978,5.835981,5.835981,5.989852, + 5.989852,5.989852,5.989852,5.993398,5.993398, + 5.996287,5.996287,6.005342,6.014956,6.015613, + 6.015613,6.015613,5.971639,5.971639,5.850568, + 5.850568,5.850568,5.990061,5.990061,5.990061, + 5.990061,5.991796,5.991796,5.993244,5.993244, + 6.011374,6.018783,6.017594,6.017594,6.017594, + 5.858175,5.858175,5.863377,5.863377,5.863377, + 5.988659,5.988659,5.988659,5.988659,5.989390, + 5.989390,5.994151,5.994151,6.020563,6.027375, + 6.026996,6.026996,6.026996,5.866343,5.866343, + 5.874646,5.874646,5.874646,5.990859,5.990859, + 5.992667,5.992667,5.994933,5.994933,5.994933, + 5.994933,6.030020,6.030020,6.038991,6.038991, + 6.038991,5.873279,5.873279,5.877365,5.877365, + 5.877365,5.992039,5.992039,5.996431,5.996431, + 6.000838,6.000838,6.000838,6.000838,6.039686, + 6.039686,6.042562,6.042562,6.042562,6.018730, + 5.886326,5.994597,5.994597,6.047575,6.047575, + 6.047575,6.078406,6.078406/ DATA SNII/4.048939,4.049242,4.049750,4.145151,4.258360, + 4.356432,4.688145,4.688257,4.688270,4.826217, + 4.826217,4.826272,5.142618,4.284333,4.284811, + 4.286872,4.288557,5.253337,4.532960,4.564950, + 4.565974,4.567596,5.379367,4.605227,4.634192, + 4.634804,4.635817,4.698180,4.771750,4.928997, + 4.930175,4.931792,4.940517,4.947424,4.947905, + 4.948512,4.976005,4.977055,4.977624,4.985716, + 5.001774,4.517682,4.519204,4.522714,5.167543, + 4.688999,4.706267,4.707887,4.710950,4.720615, + 4.721386,4.723658,4.732426,4.791679,5.359476, + 5.360877,5.362645,4.824183,4.939473,4.941747, + 4.944790,4.948400,4.959554,4.960499,4.962098, + 4.976268,4.977930,4.978850,4.988034,4.988136, + 4.988983,4.989127,4.990716,4.997378,4.997656, + 5.001124,4.999108,5.001566,5.001872,5.002052, + 5.004650,5.004791,5.510713,5.511550,5.512880, + 4.633615,4.635821,4.641956,4.663452,4.970949, + 4.970949,4.970949,4.990886,4.990886,4.991119, + 4.994713,4.999842,5.000126,5.004091,5.004655, + 5.899771,5.900361,5.901458,5.903069,5.905087, + 5.975470,5.976437,5.978201,6.162114,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000/ DATA SNIII/3.264886,3.265738,3.559230,3.555733,3.556163, + 3.795859,3.795903,3.971288,4.062202,4.062887, + 4.328298,4.441761,4.441879,3.362787,4.644640, + 4.644671,3.648880,3.649321,3.924340,3.924418, + 4.208216,4.209135,4.210838,4.353292,4.355043, + 3.749472,4.546074,4.546964,4.554955,4.555551, + 4.556746,4.558211,3.785649,3.786211,4.632736, + 4.686664,4.687436,4.688480,3.926233,3.926969, + 3.987034,3.987034,4.752837,4.754452,4.866728, + 4.928826,4.929522,4.930549,4.931964,4.980142, + 4.980413,4.980860,4.981436,3.664742,5.015918, + 5.016470,5.050785,5.051935,5.052674,5.126587, + 5.129026,3.959079,3.959144,5.192320,5.193925, + 3.989434,3.989434,4.005147,4.005147,3.968565, + 3.968565,3.992409,3.992409,4.006539,4.006539, + 5.571515,5.574721,5.580696,6.132490,6.134099, + 5.935632,5.939608,6.083616,6.087341,6.092149, + 6.099411,6.364476,6.365573,6.178216,6.185985, + 6.229222,6.316157,6.320952,6.326658,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000,7.000000,7.000000, + 7.000000,7.000000,7.000000/ DATA SNIV/2.226836,2.490623,2.490878,2.491461,2.755431, + 2.952339,2.952669,2.953231,3.013266,3.231892, + 2.493613,2.749589,2.762857,2.763010,2.763353, + 2.901417,2.901452,2.901533,2.994477,3.382731, + 3.383611,3.385459,3.472422,3.564919,3.607152, + 3.607152,3.608737,3.645438,3.728391,3.728391, + 2.639493,3.779903,3.797702,3.799535,3.799535, + 3.799535,2.797721,2.797721,2.797721,3.872625, + 3.873032,3.873596,3.883205,2.857107,2.934635, + 2.934635,2.934635,3.951730,3.952443,3.952443, + 2.993446,3.029914,3.030044,3.030246,4.061023, + 3.127314,2.880353,2.950477,2.950477,2.950477, + 2.998267,2.998267,2.998267,2.959708,2.959708, + 2.959708,4.785085,4.873190,4.873190,4.874527, + 7.000000,7.000000,7.000000/ DATA SNV/1.634565,1.915400,1.916327,1.771111,1.943796, + 1.944398,1.997854,1.998051,1.833396,1.959357, + 1.959357,1.999384,1.999384,1.870467,1.969524, + 1.969524,2.001982,2.001982,1.895972,1.977561, + 1.977561,2.004889,2.004593,2.005843,2.005843, + 2.006106,2.006106,2.006106,1.914798,1.983841, + 1.983841,2.007186,2.007186,2.008574,2.008574, + 2.008797,2.008797,2.008797,2.008797,1.929893, + 1.990742,1.990742,2.010469,2.010469,2.011696, + 2.011696,2.011930,2.011930,2.011930,2.011930, + 2.011930/ DATA SNVI/0.6290283,0.7657156,0.9208641,0.9208946,0.9217644, + 0.9074407,1.024314,1.024867/ DATA SOI/5.998809,6.000254,6.000874,6.149045,6.029965, + 6.280362,6.354168,6.620857,6.620917,6.621027, + 6.681873,6.681856,6.681878,6.554275,6.592577, + 6.991942,6.991948,6.991948,6.991953,6.991953, + 6.991953,6.994710,6.994710,6.994710,6.749977, + 6.750016,6.750087,6.784759,6.784759,6.784759, + 7.156593,7.157186,7.157570,6.676266,6.701954, + 7.234456,6.993919,6.993919,6.993919,6.993919, + 6.993919,6.997045,6.997045,6.997045,6.836974, + 6.836974,6.836974,6.746833,6.766088,6.996467, + 6.996467,6.996467,6.996467,6.996467,6.999115, + 6.999115,6.999115,6.869720,6.869720,6.869720, + 6.793480,6.808908,6.999084,6.999084,6.999084, + 6.999084,6.999084,7.001479,7.001479,7.001479, + 6.826997,6.839929,7.001803,7.001803,7.001803, + 7.001803,7.001803,7.003955,7.003955,7.003955, + 6.852827,6.864539,7.004705,8.000000,7.004705, + 7.004705,7.004705,7.007905,7.007905,7.007905, + 6.877349,6.883498,7.007761,7.007761,7.007761, + 7.007761,7.007761,7.010591,7.010591,7.010591, + 6.890946,6.900387,7.011453,7.011453,7.011453, + 7.011453,7.011453,7.012788,7.012788,7.012788, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000/ DATA SOII/4.784610,4.940430,4.940555,5.022952,5.022962, + 5.557054,5.558274,5.558889,5.930025,5.929834, + 5.160761,5.162283,5.164577,5.214052,5.216705, + 6.210938,5.445367,5.490556,5.491464,5.492962, + 5.495002,5.494491,5.494507,5.517109,5.517870, + 5.519391,5.570158,5.573380,5.580988,6.392210, + 6.394130,5.615287,5.616316,5.844496,5.880436, + 5.880893,5.904250,5.904674,5.930839,5.931911, + 5.933457,5.935489,5.954107,5.955576,5.956912, + 5.954453,5.955388,5.959104,5.959794,5.959825, + 5.959952,5.960805,5.964087,5.973599,5.975908, + 5.993385,5.994448,5.442262,5.445265,5.449878, + 5.491284,5.496742,6.232405,5.654758,5.656267, + 5.658962,5.662895,5.725537,5.731195,5.738926, + 5.741796,6.348958,6.348958,6.373241,6.373322, + 6.408604,6.408635,6.419951,6.420006,6.439373, + 6.439435,5.943564,5.943564,5.943564,5.947441, + 5.961402,5.961402,5.961402,5.961402,6.471051, + 5.965786,5.967088,5.967867,5.968222,5.972136, + 5.972852,5.978758,6.488330,5.986873,5.991322, + 5.986938,5.991353,5.994984,5.994948,5.989273, + 5.989404,5.991879,5.997332,5.991945,5.997554, + 5.992443,5.994410,6.001196,6.001346,6.002642, + 6.003121,6.002740,6.003380,5.576062,5.580966, + 5.588796,5.609913,5.619139,6.121709,6.121740, + 5.734796,5.734796,5.738977,5.745944,5.750079, + 5.750079,5.754774,5.772264,5.781077,5.960446, + 5.960446,5.960446,5.960446,5.972282,5.974347, + 5.974347,5.980535,5.980535,5.983981,5.987715, + 5.990945,5.996456,5.985451,5.986923,5.990890, + 6.000214,5.991003,6.000511,5.996286,5.996286, + 6.000386,6.002267,6.003636,6.004436,6.003808, + 6.004436,6.864734,6.865004,6.865458,6.871479, + 6.874276,6.883227,6.883227,6.929313,6.929313, + 6.945119,6.945119,6.947771,6.947771,7.214550, + 7.214550,6.955942,6.960794,6.960794,6.963802, + 6.963802,6.974759,6.974759,6.897856,6.897856, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000/ DATA SOIII/3.980091,3.980605,3.981483,4.073126,4.181011, + 4.263698,4.567496,2.851461,2.851508,4.688399, + 4.688406,4.688483,4.944256,5.004756,5.087306, + 4.201648,4.202927,4.205704,4.265077,5.589531, + 5.591178,5.591910,4.466920,4.500863,4.502461, + 4.505044,4.544430,5.702087,4.576288,4.577272, + 4.578837,4.653335,4.743010,4.880211,4.882788, + 4.885133,4.883788,4.916797,4.917453,4.918434, + 4.946753,4.948285,4.949112,4.978528,4.991553, + 5.071567,5.073311,5.091046,6.092469,5.236801, + 5.238239,5.240906,4.450987,4.453166,4.458786, + 4.490992,5.440957,5.477274,5.477840,5.478956, + 5.480584,5.482658,4.640882,4.657492,4.659830, + 4.664354,4.689624,5.526725,5.527675,5.529354, + 4.742366,4.744359,4.746791,4.755242,4.803841, + 5.629193,5.463434,5.633066,5.655169,4.904211, + 4.904211,4.904211,4.911571,5.696494,5.696715, + 5.697021,4.948280,4.949739,4.951245,4.983719, + 4.983719,4.983719,4.985557,4.992922,4.595488, + 4.595488,4.595488,4.614187,5.995187,5.995924, + 5.997336,6.003933,6.004729,6.005913,6.007472, + 6.009368,6.079757,6.079680,6.079593,6.079808, + 6.081549,6.087021,6.088523,6.089349,6.127790, + 6.130199,6.131379,6.150372,6.152512,6.155484, + 4.922875,4.922875,4.922875,4.932409,4.922875, + 4.922875,4.922875,4.991952,4.997716,6.251313, + 6.251994,6.253141,4.947119,4.974443,4.974443, + 4.974443,5.003925,6.782259,6.828280,6.541492, + 6.547457,6.555665,6.965416,7.060263,7.160808, + 7.160808,7.166230,7.174317,7.185196,7.256399, + 7.261456,7.271212,7.771374,7.771374,7.771374, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000/ DATA SOIV/3.228562,3.230039,3.508826,3.509360,3.510109, + 3.741247,3.741307,3.904661,3.977057,3.978160, + 4.214302,4.331145,4.331291,4.503491,4.503533, + 3.322590,3.617383,3.618198,4.097019,4.098440, + 4.101037,4.249484,4.252384,4.410741,4.412061, + 4.420406,4.421309,4.422863,4.425270,4.491520, + 4.543003,4.544125,4.545658,4.591770,4.594849, + 3.506632,4.717019,4.742457,4.743450,4.744866, + 4.746809,4.800908,4.801279,4.801878,4.802707, + 4.826727,4.827440,4.856910,4.858391,4.859354, + 3.927959,3.928085,4.948471,4.951597,4.995503, + 4.997566,5.057139,5.057223,3.602068,5.487784, + 5.488193,5.528638,5.529685,3.943577,3.943577, + 3.956409,3.956409,5.607411,5.152442,5.155901, + 5.162243,5.906104,5.906104,5.285100,5.290775, + 3.955026,3.955026,5.325924,5.330538,6.007065, + 6.007765,6.023023,6.026014,6.030425,6.132527, + 6.133009,5.594399,5.600957,3.969004,3.969004, + 5.768017,5.824149,5.824149,5.824149,5.824149, + 5.865850,5.865850,5.865850,5.887424,5.890223, + 5.900586,5.903125,5.912084,5.919259,5.981793, + 5.988238,6.512458,4.040435,4.040435,6.034046, + 6.039135,6.592915,6.593322,6.679720,6.679720, + 6.679720,6.679720,6.791924,6.791924,6.791924, + 7.150804,7.152208,7.208681,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000, + 8.000000/ DATA SOV/2.212299,2.477108,2.477559,2.478570,2.736373, + 2.929941,2.930501,2.931468,2.995396,3.204502, + 2.480139,2.586177,2.736415,2.753261,2.753545, + 2.754149,2.895501,2.895588,2.895747,2.990361, + 3.333698,3.335127,3.338143,3.434917,3.509306, + 3.551885,3.553720,3.556629,3.614975,3.666382, + 3.667461,3.669268,3.714561,3.738797,3.808601, + 3.809227,3.810236,3.840735,3.845908,3.847311, + 3.848131,3.893723,3.957266,2.655746,2.780047, + 2.842480,2.842480,2.842734,2.867645,2.932266, + 2.932352,2.932553,2.996815,3.040747,2.722722, + 2.858081,4.369750,4.369750,4.369750,2.990545, + 4.293693,4.399677,4.429359,4.432752,4.438707, + 4.454041,4.514207,4.514207,4.517767,4.570812, + 4.571450,2.913395,2.952782,2.965416,2.965416, + 2.965416,4.644915,4.646958,4.649636,2.994349, + 4.688903,4.691261,4.692408,4.774586,4.782194, + 2.928606,3.022013,3.022013,3.022013,2.997125, + 2.933282,2.986425,2.986425,2.986425,5.872365, + 5.931636,5.931636,5.931636,6.025977,6.025977, + 6.025977,6.090481,6.099398,6.217293,6.217293, + 6.217293,6.343698,8.000000,8.000000,8.000000, + 8.000000,8.000000,8.000000,8.000000,8.000000/ DATA SOVI/1.626749,1.908750,1.910343,1.765577,1.939605, + 1.940666,1.997515,1.997865,1.829541,1.956604, + 1.957376,1.999562,1.999822,2.001964,2.002083, + 1.867336,1.968342,1.968342,2.001995,2.001995, + 1.976059,1.976059,2.004681,2.004681,2.007063, + 2.007063,2.007365,2.007365,2.007365,1.914098, + 1.982390,1.982390,2.008208,2.008208,2.010369, + 2.010369,2.010631,2.010631,2.010631,2.010631, + 1.930105,1.987143,1.987143,2.014184,2.014184, + 2.014431,2.014431,2.014431,2.014431,2.014431, + 2.014965,2.014965/ DATA SOVII/0.6273875,0.7631111,0.9180546,0.9182081,0.9196253, + 1.029738,0.9543552,0.9543552,0.9543552,1.004676, + 1.004676,1.004676,1.025589,1.027917,1.034432, + 1.045346/ * * Find index for atom and ion, 10*IAT+IZI * c IF(IAT.EQ.26.AND.IZI.GE.6.AND.IZI.LE.9) GO TO 260 IF(IAT.GT.2.AND.IAT.LT.6)GO TO 9999 IF(IAT.LT.1.OR.IAT.GT.8)GO TO 9999 IND=10*IAT+IZI IF(IND.EQ.11) GO TO 11 IF(IND.EQ.21) GO TO 21 IF(IND.EQ.22) GO TO 22 IF(IND.EQ.61) GO TO 61 IF(IND.EQ.61) GO TO 62 IF(IND.EQ.63) GO TO 63 IF(IND.EQ.64) GO TO 64 IF(IND.EQ.65) GO TO 65 IF(IND.EQ.66) GO TO 66 IF(IND.EQ.71) GO TO 71 IF(IND.EQ.72) GO TO 72 IF(IND.EQ.73) GO TO 73 IF(IND.EQ.74) GO TO 74 IF(IND.EQ.75) GO TO 75 IF(IND.EQ.76) GO TO 76 IF(IND.EQ.77) GO TO 77 IF(IND.EQ.81) GO TO 81 IF(IND.EQ.82) GO TO 82 IF(IND.EQ.83) GO TO 83 IF(IND.EQ.84) GO TO 84 IF(IND.EQ.85) GO TO 85 IF(IND.EQ.86) GO TO 86 IF(IND.EQ.87) GO TO 87 IF(IND.EQ.88) GO TO 88 * * CALCULATING PARTITION FUNCTIONS FOR HYDROGEN * 11 CALL PARTDV(T,ANE,ZH,MH,NHYD,GHYD,ENHYD,SHYD,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR HEI * 21 CALL PARTDV(T,ANE,ZHE,MHEI,NHEL,GHEL,ENHEL,SHEL,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR HEII * 22 CALL PARTDV(T,ANE,ZHE,MHEII,NHYD,GHYD,ENHYD,SHYD,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR CI * 61 CALL PARTDV(T,ANE,ZC,MCI,NCI,GCI,ENCI,SCI,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR CII * 62 CALL PARTDV(T,ANE,ZC,MCII,NCII,GCII,ENCII,SCII,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR CIII * 63 CALL PARTDV(T,ANE,ZC,MCIII,NCIII,GCIII,ENCIII,SCIII,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR CIV * 64 CALL PARTDV(T,ANE,ZC,MCIV,NCIV,GCIV,ENCIV,SCIV,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR CV * 65 CALL PARTDV(T,ANE,ZC,MCV,NCV,GCV,ENCV,SCV,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR CVI * 66 CALL PARTDV(T,ANE,ZC,MH,NHYD,GHYD,ENHYD,SHYD,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR NI * 71 CALL PARTDV(T,ANE,ZN,MNI,NNI,GNI,ENNI,SNI,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR NII * 72 CALL PARTDV(T,ANE,ZN,MNII,NNII,GNII,ENNII,SNII,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR NIII * 73 CALL PARTDV(T,ANE,ZN,MNIII,NNIII,GNIII,ENNIII,SNIII,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR NIV * 74 CALL PARTDV(T,ANE,ZN,MNIV,NNIV,GNIV,ENNIV,SNIV,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR NV * 75 CALL PARTDV(T,ANE,ZN,MNV,NNV,GNV,ENNV,SNV,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR NVI * 76 CALL PARTDV(T,ANE,ZN,MNVI,NNVI,GNVI,ENNVI,SNVI,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR NVII * 77 CALL PARTDV(T,ANE,ZN,MH,NHYD,GHYD,ENHYD,SHYD,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OI * 81 CALL PARTDV(T,ANE,ZO,MOI,NOI,GOI,ENOI,SOI,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OII * 82 CALL PARTDV(T,ANE,ZO,MOII,NOII,GOII,ENOII,SOII,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OIII * 83 CALL PARTDV(T,ANE,ZO,MOIII,NOIII,GOIII,ENOIII,SOIII,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OIV * 84 CALL PARTDV(T,ANE,ZO,MOIV,NOIV,GOIV,ENOIV,SOIV,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OV * 85 CALL PARTDV(T,ANE,ZO,MOV,NOV,GOV,ENOV,SOV,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OVI * 86 CALL PARTDV(T,ANE,ZO,MOVI,NOVI,GOVI,ENOVI,SOVI,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OVII * 87 CALL PARTDV(T,ANE,ZO,MOVII,NOVII,GOVII,ENOVII,SOVII,U) GO TO 8888 * * CALCULATING PARTITION FUNCTIONS FOR OVIII * 88 CALL PARTDV(T,ANE,ZO,MH,NHYD,GHYD,ENHYD,SHYD,U) GO TO 8888 C C C CALCULATING PARTITION FUNCTIONS FOR FE VI - FE IX C C260 CALL PFFE(IZI,T,ANE,U) 8888 CONTINUE RETURN 9999 U=0 WRITE(*,*)!! INVALID ATOM IN USER SUPPLIED ROUTINE PARTFUN !! STOP END