PROGRAM READLOW C SAMPLE PROGRAM FOR READING LOWLINES.DAT C C These data were collected or computed by Kurucz and described at the C IAU General Assembly in Baltimore in 1988. Up to now they have been C distributed on 30 magnetic tapes. The data are resident in the mass C storage of the San Diego Supercomputer Center and have been used by C Kurucz to compute opacities. From those files we have extracted the C absolute minimum data for each line that is needed for computing the C line opacity and we have packed it into 16 bytes. The atomic data C can be fitted onto one CD-ROM. For detailed spectrum synthesis it is C necessary to have much more information about each line, the classi- C fication, exact energy levels, etc. All these data will be made C available on future CD-ROMs. Much of these data can be improved using C newer laboratory data and the bigger, faster computers now available. C As significant progress is made, newer versions of this CD will be C published. C Robert Kurucz 8 February 1993 C C1988 Semiempirical calculation of gf values for the iron group. Presented C at the meeting of Commission 14 on Atomic and Molecular Spectroscopy C at the IAU General Assembly in Baltimore, August 1988. in Transactions C of the International Astronomical Union, Volume XXB. M.McNally, ed., C Dordrecht: Kluwer, 168-172. C C1991 New opacity calculations. Presented at the NATO Advanced Research Work- C shop, Stellar Atmospheres: Beyond Classical Models, Marina d'Aurisina C (Trieste) Italy, 3-7 September 1990. pp.440-448 in Stellar Atmospheres: C Beyond Classical Models (ed. by L. Crivellari, I. Hubeny, and D.G. Hummer), C NATO ASI Series, Kluwer, Dordrecht. C C1992 Atomic and molecular data for opacity calculations. Presented at the C Workshop on Astrophysical Opacities, Caracas, 15-19 July 1991. Revista C Mexicana de Astronomia y Astrofisica, vol. 23, pp. 45-48. C C1992 "Finding" the "missing" solar ultraviolet opacity. Presented at the C Workshop on Astrophysical Opacities, Caracas, 15-19 July 1991. Revista C Mexicana de Astronomia y Astrofisica, vol. 23, pp. 181-186. C C1992 Remaining line opacity problems for the solar spectrum. Presented at the C Workshop on Astrophysical Opacities, Caracas, 15-19 July 1991. Revista C Mexicana de Astronomia y Astrofisica, vol. 23, pp. 187-194. C C1992 A new opacity-sampling model atmosphere program for arbitrary abundances. C Presented at IAU Colloquium 138, Peculiar versus Normal Phenomena C in A-type and Related Stars, Trieste, Italy, 6-10 July 1992. To be C published in IAU Colloquium 138, Peculiar versus Normal Phenomena in C A-type and Related Stars (eds. M.M. Dworetsky, F. Castelli, and R. C Faraggiana), A.S.P. Conference Series, Astronomical Society of the C Pacific, San Francisco. C REAL*8 RESOLU,RATIO,RATIOLG REAL*8 WLVAC,START,STOP REAL*4 TABLOG(32768) INTEGER*2 IELION,IELO,IGFLOG,IGR,IGS,IGW COMMON /IIIIIII/IWL,IELION,IELO,IGFLOG,IGR,IGS,IGW INTEGER*4 IIIIIII(4) EQUIVALENCE (IIIIIII(1),IWL) BYTE IIBYTE(16),ONEBYTE EQUIVALENCE (IIIIIII(1),IIBYTE(1)) DIMENSION ELEM(840) DIMENSION ELEMA(119),ELEMB(156),ELEMC(159),ELEMD(61) DIMENSION ELEME(95),ELEMF(95),ELEMG(95),ELEMH(60) EQUIVALENCE (ELEM( 1),ELEMA(1)),(ELEM(120),ELEMB(1)) EQUIVALENCE (ELEM(276),ELEMC(1)),(ELEM(435),ELEMD(1)) EQUIVALENCE (ELEM(496),ELEME(1)),(ELEM(591),ELEMF(1)) EQUIVALENCE (ELEM(686),ELEMG(1)),(ELEM(781),ELEMH(1)) DIMENSION I100000(317) DIMENSION I100000A(54),I100000B(54),I100000C(54) DIMENSION I100000D(54),I100000E(54),I100000F(47) EQUIVALENCE (I100000( 1),I100000A(1)) EQUIVALENCE (I100000( 55),I100000B(1)) EQUIVALENCE (I100000(109),I100000C(1)) EQUIVALENCE (I100000(163),I100000D(1)) EQUIVALENCE (I100000(217),I100000E(1)) EQUIVALENCE (I100000(271),I100000F(1)) DATA I100000A/ 1 1824404, 5562510, 5753867, 5892965, 6059758, 6240560, 2 6357512, 6447033, 6523911, 6602547, 6680540, 6763683, 3 6843602, 6918299, 6984962, 7048369, 7107200, 7164745, 4 7220908, 7275195, 7328276, 7380007, 7432694, 7482725, 5 7529234, 7575114, 7620053, 7662936, 7704228, 7743754, 6 7782868, 7822468, 7859858, 7895448, 7931629, 7969049, 7 8007695, 8045971, 8083602, 8119895, 8155315, 8190043, 8 8224088, 8257935, 8291286, 8324839, 8358054, 8389984, 9 8420984, 8451429, 8480402, 8507602, 8532831, 8559171/ DATA I100000B/ 1 8585057, 8609681, 8634371, 8658295, 8683163, 8709114, 2 8734649, 8760093, 8785346, 8811071, 8836026, 8860308, 3 8884321, 8909551, 8935323, 8960959, 8986179, 9010581, 4 9034534, 9057973, 9081067, 9104264, 9127266, 9150082, 5 9172529, 9195183, 9217143, 9238406, 9260301, 9282902, 6 9306487, 9330156, 9353782, 9376444, 9398163, 9418327, 7 9438824, 9459836, 9480500, 9500751, 9520426, 9539296, 8 9557608, 9576754, 9595781, 9615176, 9634293, 9653520, 9 9672132, 9690468, 9708534, 9726772, 9743851, 9761371/ DATA I100000C/ 1 9779750, 9799094, 9818982, 9838780, 9858055, 9877387, 2 9896925, 9915888, 9934516, 9952460, 9969808, 9985833, 3 10002634, 10020826, 10039751, 10058794, 10077860, 10097083, 4 10116197, 10134981, 10153827, 10172623, 10191590, 10210761, 5 10230321, 10250698, 10271057, 10291616, 10311886, 10332309, 6 10352883, 10373217, 10393005, 10412224, 10431911, 10451520, 7 10472072, 10493332, 10514832, 10536144, 10556601, 10576172, 8 10594857, 10613731, 10633085, 10652023, 10671757, 10691809, 9 10712097, 10731969, 10752112, 10773186, 10794268, 10814818/ DATA I100000D/ 1 10834916, 10854205, 10872871, 10889859, 10907076, 10925650, 2 10943904, 10962575, 10981709, 11002566, 11023742, 11045250, 3 11066864, 11088538, 11109937, 11131129, 11152360, 11173491, 4 11194440, 11215201, 11235648, 11255643, 11274791, 11292941, 5 11310681, 11327595, 11344622, 11362941, 11382251, 11402653, 6 11424377, 11446729, 11469396, 11494588, 11520571, 11546998, 7 11573428, 11600186, 11626097, 11651951, 11677482, 11702707, 8 11727717, 11752834, 11776550, 11799396, 11822470, 11846309, 9 11871752, 11897863, 11924876, 11952587, 11979394, 12006915/ DATA I100000E/ 1 12033968, 12059671, 12085645, 12111339, 12136347, 12160165, 2 12183906, 12208073, 12233170, 12258275, 12282888, 12306500, 3 12331097, 12358173, 12385877, 12415029, 12445544, 12475290, 4 12503337, 12532860, 12563887, 12594838, 12628481, 12663446, 5 12698565, 12733735, 12768464, 12802812, 12836253, 12869774, 6 12904079, 12938737, 12973517, 13009866, 13046094, 13081625, 7 13115750, 13149532, 13184840, 13222071, 13258728, 13295121, 8 13331700, 13368916, 13405474, 13441611, 13480841, 13522517, 9 13564681, 13606259, 13647523, 13690632, 13734617, 13780762/ DATA I100000F/ 1 13827522, 13874879, 13924448, 13976425, 14027411, 14075924, 2 14124969, 14173772, 14223660, 14277151, 14335259, 14396198, 3 14457484, 14522736, 14588410, 14656501, 14727084, 14800333, 4 14874624, 14947410, 15022147, 15101157, 15185806, 15278350, 5 15375114, 15479650, 15588185, 15703838, 15818869, 15947717, 6 16077228, 16224174, 16381290, 16552267, 16727667, 16917649, 7 17118363, 17346959, 17590897, 17875438, 18183158, 18502692, 8 18894577, 19374281, 20010771, 20935370, 22342326/ C 9 23026377/=last line= 31627892th line DATA ELEMA/ 1 1.00, 1.01, 2 2.00, 2.01, 2.02, 3 3.00, 3.01, 3.02, 3.03, 4 4.00, 4.01, 4.02, 4.03, 4.04, 5 5.00, 5.01, 5.02, 5.03, 5.04, 5.05, 6 6.00, 6.01, 6.02, 6.03, 6.04, 6.05, 6.06, 7 7.00, 7.01, 7.02, 7.03, 7.04, 7.05, 7.06, 7.07, 8 8.00, 8.01, 8.02, 8.03, 8.04, 8.05, 8.06, 8.07, 8.08, 9 9.00, 9.01, 9.02, 9.03, 9.04, 9.05, 9.06, 9.07, 9.08, 9.09, A 10.00,10.01,10.02,10.03,10.04,10.05,10.06,10.07,10.08,10.09, 1 10.10, 2 11.00,11.01,11.02,11.03,11.04,11.05,11.06,11.07,11.08,11.09, 3 11.10,11.11, 4 12.00,12.01,12.02,12.03,12.04,12.05,12.06,12.07,12.08,12.09, 5 12.10,12.11,12.12, 6 13.00,13.01,13.02,13.03,13.04,13.05,13.06,13.07,13.08,13.09, 7 13.10,13.11,13.12,13.13, 8 14.00,14.01,14.02,14.03,14.04,14.05,14.06,14.07,14.08,14.09, 9 14.10,14.11,14.12,14.13,14.14/ DATA ELEMB/ 1 15.00,15.01,15.02,15.03,15.04,15.05,15.06,15.07,15.08,15.09, 2 15.10,15.11,15.12,15.13,15.14,15.15, 3 16.00,16.01,16.02,16.03,16.04,16.05,16.06,16.07,16.08,16.09, 4 16.10,16.11,16.12,16.13,16.14,16.15,16.16, 5 17.00,17.01,17.02,17.03,17.04,17.05,17.06,17.07,17.08,17.09, 6 17.10,17.11,17.12,17.13,17.14,17.15,17.16,17.17, 7 18.00,18.01,18.02,18.03,18.04,18.05,18.06,18.07,18.08,18.09, 8 18.10,18.11,18.12,18.13,18.14,18.15,18.16,18.17,18.18, 9 19.00,19.01,19.02,19.03,19.04,19.05,19.06,19.07,19.08,19.09, A 19.10,19.11,19.12,19.13,19.14,19.15,19.16,19.17,19.18,19.19, 1 20.00,20.01,20.02,20.03,20.04,20.05,20.06,20.07,20.08,20.09, 2 20.10,20.11,20.12,20.13,20.14,20.15,20.16,20.17,20.18,20.19, 3 20.20, 4 21.00,21.01,21.02,21.03,21.04,21.05,21.06,21.07,21.08,21.09, 5 21.10,21.11,21.12,21.13,21.14,21.15,21.16,21.17,21.18,21.19, 6 21.20,21.21, 7 22.00,22.01,22.02,22.03,22.04,22.05,22.06,22.07,22.08,22.09, 8 22.10,22.11,22.12,22.13,22.14,22.15,22.16,22.17,22.18,22.19, 9 22.20,22.21,22.22/ DATA ELEMC/ 1 23.00,23.01,23.02,23.03,23.04,23.05,23.06,23.07,23.08,23.09, 2 23.10,23.11,23.12,23.13,23.14,23.15,23.16,23.17,23.18,23.19, 3 23.20,23.21,23.22,23.23, 4 24.00,24.01,24.02,24.03,24.04,24.05,24.06,24.07,24.08,24.09, 5 24.10,24.11,24.12,24.13,24.14,24.15,24.16,24.17,24.18,24.19, 6 24.20,24.21,24.22,24.23,24.24, 7 25.00,25.01,25.02,25.03,25.04,25.05,25.06,25.07,25.08,25.09, 8 25.10,25.11,25.12,25.13,25.14,25.15,25.16,25.17,25.18,25.19, 9 25.20,25.21,25.22,25.23,25.24,25.25, A 26.00,26.01,26.02,26.03,26.04,26.05,26.06,26.07,26.08,26.09, 1 26.10,26.11,26.12,26.13,26.14,26.15,26.16,26.17,26.18,26.19, 2 26.20,26.21,26.22,26.23,26.24,26.25,26.26, 3 27.00,27.01,27.02,27.03,27.04,27.05,27.06,27.07,27.08,27.09, 4 27.10,27.11,27.12,27.13,27.14,27.15,27.16,27.17,27.18,27.19, 5 27.20,27.21,27.22,27.23,27.24,27.25,27.26,27.27, 6 28.00,28.01,28.02,28.03,28.04,28.05,28.06,28.07,28.08,28.09, 7 28.10,28.11,28.12,28.13,28.14,28.15,28.16,28.17,28.18,28.19, 8 28.20,28.21,28.22,28.23,28.24,28.25,28.26,28.27,28.28/ DATA ELEMD/ 1 29.00,29.01,29.02,29.03,29.04,29.05,29.06,29.07,29.08,29.09, 2 29.10,29.11,29.12,29.13,29.14,29.15,29.16,29.17,29.18,29.19, 3 29.20,29.21,29.22,29.23,29.24,29.25,29.26,29.27,29.28,29.29, 4 30.00,30.01,30.02,30.03,30.04,30.05,30.06,30.07,30.08,30.09, 5 30.10,30.11,30.12,30.13,30.14,30.15,30.16,30.17,30.18,30.19, 6 30.20,30.21,30.22,30.23,30.24,30.25,30.26,30.27,30.28,30.29, 7 30.30/ DATA ELEME/ 1 31.00,31.01,31.02,31.03,31.04, 2 32.00,32.01,32.02,32.03,32.04, 3 33.00,33.01,33.02,33.03,33.04, 4 34.00,34.01,34.02,34.03,34.04, 5 35.00,35.01,35.02,35.03,35.04, 6 36.00,36.01,36.02,36.03,36.04, 7 37.00,37.01,37.02,37.03,37.04, 8 38.00,38.01,38.02,38.03,38.04, 9 39.00,39.01,39.02,39.03,39.04, A 40.00,40.01,40.02,40.03,40.04, 1 41.00,41.01,41.02,41.03,41.04, 2 42.00,42.01,42.02,42.03,42.04, 3 43.00,43.01,43.02,43.03,43.04, 4 44.00,44.01,44.02,44.03,44.04, 5 45.00,45.01,45.02,45.03,45.04, 6 46.00,46.01,46.02,46.03,46.04, 7 47.00,47.01,47.02,47.03,47.04, 8 48.00,48.01,48.02,48.03,48.04, 9 49.00,49.01,49.02,49.03,49.04/ DATA ELEMF/ 1 50.00,50.01,50.02,50.03,50.04, 2 51.00,51.01,51.02,51.03,51.04, 3 52.00,52.01,52.02,52.03,52.04, 4 53.00,53.01,53.02,53.03,53.04, 5 54.00,54.01,54.02,54.03,54.04, 6 55.00,55.01,55.02,55.03,55.04, 7 56.00,56.01,56.02,56.03,56.04, 8 57.00,57.01,57.02,57.03,57.04, 9 58.00,58.01,58.02,58.03,58.04, A 59.00,59.01,59.02,59.03,59.04, 1 60.00,60.01,60.02,60.03,60.04, 2 61.00,61.01,61.02,61.03,61.04, 3 62.00,62.01,62.02,62.03,62.04, 4 63.00,63.01,63.02,63.03,63.04, 5 64.00,64.01,64.02,64.03,64.04, 6 65.00,65.01,65.02,65.03,65.04, 7 66.00,66.01,66.02,66.03,66.04, 8 67.00,67.01,67.02,67.03,67.04, 9 68.00,68.01,68.02,68.03,68.04/ DATA ELEMG/ 1 69.00,69.01,69.02,69.03,69.04, 2 70.00,70.01,70.02,70.03,70.04, 3 71.00,71.01,71.02,71.03,71.04, 4 72.00,72.01,72.02,72.03,72.04, 5 73.00,73.01,73.02,73.03,73.04, 6 74.00,74.01,74.02,74.03,74.04, 7 75.00,75.01,75.02,75.03,75.04, 8 76.00,76.01,76.02,76.03,76.04, 9 77.00,77.01,77.02,77.03,77.04, A 78.00,78.01,78.02,78.03,78.04, 1 79.00,79.01,79.02,79.03,79.04, 2 80.00,80.01,80.02,80.03,80.04, 3 81.00,81.01,81.02,81.03,81.04, 4 82.00,82.01,82.02,82.03,82.04, 5 83.00,83.01,83.02,83.03,83.04, 6 84.00,84.01,84.02,84.03,84.04, 7 85.00,85.01,85.02,85.03,85.04, 8 86.00,86.01,86.02,86.03,86.04, 9 87.00,87.01,87.02,87.03,87.04/ DATA ELEMH/ 1 88.00,88.01,88.02,88.03,88.04, 2 89.00,89.01,89.02,89.03,89.04, 3 90.00,90.01,90.02,90.03,90.04, 4 91.00,91.01,91.02,91.03,91.04, 5 92.00,92.01,92.02,92.03,92.04, 6 93.00,93.01,93.02,93.03,93.04, 7 94.00,94.01,94.02,94.03,94.04, 8 95.00,95.01,95.02,95.03,95.04, 9 96.00,96.01,96.02,96.03,96.04, A 97.00,97.01,97.02,97.03,97.04, 1 98.00,98.01,98.02,98.03,98.04, 2 99.00,99.01,99.02,99.03,99.04/ DO 1 I=1,32768 1 TABLOG(I)=10.**((I-16384)*.001) RATIOLOG=LOG(1.D0+1.D0/2000000.D0) OPEN(UNIT=11,STATUS='OLD',READONLY,SHARED,FORM='UNFORMATTED', 1RECORDTYPE='FIXED',BLOCKSIZE=8000,RECORDSIZE=4,access='direct') C LENTAP=31627892 NSKIP=0 START=500. STOP=500.1 ISTART=DLOG(START)/RATIOLOG+.5 ISTOP=DLOG(STOP)/RATIOLOG+.5 DO 2 I100=2,317 IF(ISTART.LT.I100000(I100))GO TO 5 2 NSKIP=NSKIP+100000 5 NSKIP1=NSKIP+1 DO 20 N=NSKIP1,LENTAP READ(11,rec=n)IIIIIII C 12429219 -3550 13005 22110 25757 12087 8909 C IF COMPUTER REQUIRES BYTE ROTATION C DO 18 I=1,16,2 C ONEBYTE=IIBYTE(I) C IIBYTE(I)=IIBYTE(I+1) C 18 IIBYTE(I+1)=ONEBYTE C ONEBYTE=IIBYTE(1) C IIBYTE(1)=IIBYTE(3) C IIBYTE(3)=ONEBYTE C ONEBYTE=IIBYTE(2) C IIBYTE(2)=IIBYTE(4) C IIBYTE(4)=ONEBYTE IF(IWL.LT.ISTART)GO TO 20 IF(IWL.GT.ISTOP)GO TO 21 C IELION NEGATIVE = LINE HAS PREDICTED ENERGY LEVELS AND WAVELENGTH C IELION POSITIVE = LINE HAS MEASURED ENERGY LEVELS AND REAL WAVELENGTH C UNITS DIGIT IS FOR ISOTOPE, NOT USED AT PRESENT NELION=ABS(IELION)/10 CODE=ELEM(NELION) WLVAC=EXP(IWL*RATIOLOG) IF(IELION.LT.0)WLVAC=-WLVAC GF=TABLOG(IGFLOG) ELO=TABLOG(IELO) GAMMAR=TABLOG(IGR) GAMMAS=TABLOG(IGS) GAMMAW=TABLOG(IGW) PRINT 19,N,WLVAC,CODE,GF,ELO,GAMMAR,GAMMAS,GAMMAW, 1 IWL,IELION,IGFLOG,IELO,IGR,IGS,IGW 19 FORMAT(I9,F10.4,F7.2,1PE11.3,0PF12.3,1P3E10.2,I12,6I6) C23146629 -500.0001 26.04 4.178E-04 532108.625 2.36E+09 5.05E-05 3.35E-08 C 12429219 -3550 13005 22110 25757 12087 8909 20 CONTINUE 21 CALL EXIT END