Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/788/39 into FORTRAN code for reading data files line by line.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe C============================================================================= C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-28 C============================================================================= * This code was generated from the ReadMe file documenting a catalogue * according to the "Standard for Documentation of Astronomical Catalogues" * currently in use by the Astronomical Data Centers (CDS, ADC, A&A) * (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx) * Please report problems or questions to C============================================================================= implicit none * Unspecified or NULL values, generally corresponding to blank columns, * are assigned one of the following special values: * rNULL__ for unknown or NULL floating-point values * iNULL__ for unknown or NULL integer values real*4 rNULL__ integer*4 iNULL__ parameter (rNULL__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. J/ApJ/788/39 Hot Jupiter exoplanets host stars EW and abundances (Teske+, 2014) *================================================================================ *C/O ratios of stars with transiting hot Jupiter exoplanets. * Teske J.K., Cunha K., Smith V.V., Schuler S.C., Griffith C.A. * <Astrophys. J., 788, 39 (2014)> * =2014ApJ...788...39T (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'stars.dat' ! Stars integer*4 nr__ parameter (nr__=17) ! Number of records character*48 ar__ ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg ! (deg) Right Ascension J2000 real*8 DEdeg ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*9 Name ! Host star's name integer*4 RAh ! (h) Simbad Hour of Right Ascension (J2000) integer*4 RAm ! (min) Simbad Minute of Right Ascension (J2000) real*4 RAs ! (s) Simbad Second of Right Ascension (J2000) character*1 DE_ ! Simbad Sign of the Declination (J2000) integer*4 DEd ! (deg) Simbad Degree of Declination (J2000) integer*4 DEm ! (arcmin) Simbad Arcminute of Declination (J2000) real*4 DEs ! (arcsec) Simbad Arcsecond of Declination (J2000) real*4 Vmag ! (mag) V band magnitude from Table 1 character*8 Name1 ! Host star's name as written in table4.dat file c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Lines Measured, Equivalent Widths, and Abundances integer*4 nr__1 parameter (nr__1=106) ! Number of records character*246 ar__1 ! Full-size record character*5 Ion ! Ion designation (1) real*8 Wave ! (0.1nm) Wavelength, in Angstroms real*4 ExPot ! (eV) Excitation Potential real*4 log_gf ! ([-]) Log of oscillator strength real*4 EWSolHDSI ! (10-13m) ? Solar equivalent width (2) real*4 EWSolHDSSP ! (10-13m) ? Solar equivalent width (3) real*4 EWSolHIRSP ! (10-13m) ? Solar equivalent width (4) real*4 logNSolHDSI ! ([-]) ? Log of solar number abundance (2) real*4 logNSolHDSS ! ([-]) ? Log of solar number abundance (3) real*4 logNSolHIRS ! ([-]) ? Log of solar number abundance (4) real*4 EW_COROT2 ! (10-13m) ? CoRoT-2 equivalent width (4) real*4 logN_COROT2 ! ([-]) ? Log of CoRoT-2 number abundance (4) real*4 EW_TRES4 ! (10-13m) ? TrES-4 equivalent width (4) real*4 logN_TRES4 ! ([-]) ? Log of TrES-4 number abundance (4) real*4 EW_TRES2 ! (10-13m) ? TrES-2 equivalent width (4) real*4 logN_TRES2 ! ([-]) ? Log of TrES-2 number abundance (4) real*4 EW_WASP2 ! (10-13m) ? WASP-2 equivalent width (4) real*4 logN_WASP2 ! ([-]) ? Log of WASP-2 number abundance (4) real*4 EW_WASP12 ! (10-13m) ? WASP-12 equivalent width (3) real*4 logN_WASP12 ! ([-]) ? Log of WASP-12 number abundance (3) real*4 EW_HD149026 ! (10-13m) ? HD 149026 equivalent width (3) real*4 logN_HD149026 ! ([-]) ? Log of HD 149026 number abundance (3) character*1 f_EW_HD149026 ! [*] Flag on EW-HD149026 (5) real*4 EW_HATP1 ! (10-13m) ? HAT-P-1 equivalent width (4) real*4 logN_HATP1 ! ([-]) ? Log of HAT-P-1 number abundance (4) real*4 EW_XO2S ! (10-13m) ? XO-2S equivalent width (2) real*4 logN_XO2S ! ([-]) ? Log of XO-2S number abundance (2) real*4 EW_XO2N ! (10-13m) ? XO-2N equivalent width (2) real*4 logN_XO2N ! ([-]) ? Log of XO-2N number abundance (2) real*4 EW_XO1 ! (10-13m) ? XO-1 equivalent width (2) real*4 logN_XO1 ! ([-]) ? Log of XO-1 number abundance (2) real*4 EW_TRES3 ! (10-13m) ? TRES-3 equivalent width (3) real*4 logN_TRES3 ! ([-]) ? Log of TRES-3 number abundance (3) character*1 f_EW_TRES3 ! [*] Flag on EW-TRES3 (5) real*4 EW_HD189733 ! (10-13m) ? HD 189733 equivalent width (2) real*4 logN_HD189733 ! ([-]) ? Log of HD 189733 number abundance (2) character*1 f_EW_HD189733 ! [*] Flag on EW-HD189733 (5) real*4 EW_HD80606 ! (10-13m) ? HD 80606 equivalent width (3) real*4 logN_HD80606 ! ([-]) ? Log of HD 80606 number abundance (3) real*4 EW_HATP7 ! (10-13m) ? HAT-P-7 equivalent width (3) real*4 logN_HATP7 ! ([-]) ? Log of HAT-P-7 number abundance (3) character*1 f_EW_HATP7 ! [*] Flag on EW-HATP7 (5) real*4 EW_HATP13 ! (10-13m) ? HAT-P-13 equivalent width (2) real*4 logN_HATP13 ! ([-]) ? Log of HAT-P-13 number abundance (2) real*4 EW_HATP16 ! (10-13m) ? HAT-P-16 equivalent width (4) real*4 logN_HATP16 ! ([-]) ? Log of HAT-P-16 number abundance (4) real*4 EW_WASP32 ! (10-13m) ? WASP-32 equivalent width (4) real*4 logN_WASP32 ! ([-]) ? Log of WASP-32 number abundance (4) *Note (1): For [O I] (6300.30), the log N values represent those derived from * synthesis fitting, as these are the values we use in calculating the * final [O/H] for each object. The reported EWs refer to the total EW * of the 6300.3 Angstrom blend. *Note (2): Measured with IRAF in Subaru/HDS data. EWs in milliAngstroms. *Note (3): Measured with SPECTRE in Subaru/HDS data. EWs in milliAngstroms. *Note (4): Measured with SPECTRE in Keck/HIRES data. EWs in milliAngstroms. *Note (5): Flag as follows: * * = Measured with SPECTRE in Keck/HIRES data. C============================================================================= C Loading file 'stars.dat' ! Stars C Format for file interpretation 1 format( + A9,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,A8) C Effective file loading open(unit=1,status='old',file= +'stars.dat') write(6,*) '....Loading file: stars.dat' do i__=1,17 read(1,'(A48)')ar__ read(ar__,1)Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,Name1 RAdeg = rNULL__ DEdeg = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh .GT. -180) RAdeg=RAh*15. if(RAm .GT. -180) RAdeg=RAdeg+RAm/4. if(RAs .GT. -180) RAdeg=RAdeg+RAs/240. if(DEd .GE. 0) DEdeg=DEd if(DEm .GE. 0) DEdeg=DEdeg+DEm/60. if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600. if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg c ..............Just test output........... write(6,1)Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,Name1 write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Lines Measured, Equivalent Widths, * and Abundances C Format for file interpretation 2 format( + A5,1X,F7.2,1X,F4.2,1X,F6.3,1X,F4.1,1X,F5.2,1X,F4.1,1X,F4.2,1X, + F4.2,1X,F4.2,1X,F4.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F4.1,1X,F4.2, + 1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,A1,1X,F4.1, + 1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F4.1,1X,F4.2,1X, + F4.1,1X,F4.2,1X,A1,1X,F5.1,1X,F4.2,1X,A1,1X,F5.1,1X,F4.2,1X, + F5.1,1X,F4.2,1X,A1,1X,F5.1,1X,F4.2,1X,F5.1,1X,F4.2,1X,F5.1,1X, + F4.2) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,106 read(1,'(A246)')ar__1 read(ar__1,2) + Ion,Wave,ExPot,log_gf,EWSolHDSI,EWSolHDSSP,EWSolHIRSP, + logNSolHDSI,logNSolHDSS,logNSolHIRS,EW_COROT2,logN_COROT2, + EW_TRES4,logN_TRES4,EW_TRES2,logN_TRES2,EW_WASP2,logN_WASP2, + EW_WASP12,logN_WASP12,EW_HD149026,logN_HD149026,f_EW_HD149026, + EW_HATP1,logN_HATP1,EW_XO2S,logN_XO2S,EW_XO2N,logN_XO2N, + EW_XO1,logN_XO1,EW_TRES3,logN_TRES3,f_EW_TRES3,EW_HD189733, + logN_HD189733,f_EW_HD189733,EW_HD80606,logN_HD80606,EW_HATP7, + logN_HATP7,f_EW_HATP7,EW_HATP13,logN_HATP13,EW_HATP16, + logN_HATP16,EW_WASP32,logN_WASP32 if(ar__1(27:30) .EQ. '') EWSolHDSI = rNULL__ if(ar__1(32:36) .EQ. '') EWSolHDSSP = rNULL__ if(ar__1(38:41) .EQ. '') EWSolHIRSP = rNULL__ if(ar__1(43:46) .EQ. '') logNSolHDSI = rNULL__ if(ar__1(48:51) .EQ. '') logNSolHDSS = rNULL__ if(ar__1(53:56) .EQ. '') logNSolHIRS = rNULL__ if(ar__1(58:61) .EQ. '') EW_COROT2 = rNULL__ if(ar__1(63:66) .EQ. '') logN_COROT2 = rNULL__ if(ar__1(68:72) .EQ. '') EW_TRES4 = rNULL__ if(ar__1(74:77) .EQ. '') logN_TRES4 = rNULL__ if(ar__1(79:82) .EQ. '') EW_TRES2 = rNULL__ if(ar__1(84:87) .EQ. '') logN_TRES2 = rNULL__ if(ar__1(89:93) .EQ. '') EW_WASP2 = rNULL__ if(ar__1(95:98) .EQ. '') logN_WASP2 = rNULL__ if(ar__1(100:104) .EQ. '') EW_WASP12 = rNULL__ if(ar__1(106:109) .EQ. '') logN_WASP12 = rNULL__ if(ar__1(111:115) .EQ. '') EW_HD149026 = rNULL__ if(ar__1(117:120) .EQ. '') logN_HD149026 = rNULL__ if(ar__1(124:127) .EQ. '') EW_HATP1 = rNULL__ if(ar__1(129:132) .EQ. '') logN_HATP1 = rNULL__ if(ar__1(134:138) .EQ. '') EW_XO2S = rNULL__ if(ar__1(140:143) .EQ. '') logN_XO2S = rNULL__ if(ar__1(145:149) .EQ. '') EW_XO2N = rNULL__ if(ar__1(151:154) .EQ. '') logN_XO2N = rNULL__ if(ar__1(156:159) .EQ. '') EW_XO1 = rNULL__ if(ar__1(161:164) .EQ. '') logN_XO1 = rNULL__ if(ar__1(166:169) .EQ. '') EW_TRES3 = rNULL__ if(ar__1(171:174) .EQ. '') logN_TRES3 = rNULL__ if(ar__1(178:182) .EQ. '') EW_HD189733 = rNULL__ if(ar__1(184:187) .EQ. '') logN_HD189733 = rNULL__ if(ar__1(191:195) .EQ. '') EW_HD80606 = rNULL__ if(ar__1(197:200) .EQ. '') logN_HD80606 = rNULL__ if(ar__1(202:206) .EQ. '') EW_HATP7 = rNULL__ if(ar__1(208:211) .EQ. '') logN_HATP7 = rNULL__ if(ar__1(215:219) .EQ. '') EW_HATP13 = rNULL__ if(ar__1(221:224) .EQ. '') logN_HATP13 = rNULL__ if(ar__1(226:230) .EQ. '') EW_HATP16 = rNULL__ if(ar__1(232:235) .EQ. '') logN_HATP16 = rNULL__ if(ar__1(237:241) .EQ. '') EW_WASP32 = rNULL__ if(ar__1(243:246) .EQ. '') logN_WASP32 = rNULL__ c ..............Just test output........... write(6,2) + Ion,Wave,ExPot,log_gf,EWSolHDSI,EWSolHDSSP,EWSolHIRSP, + logNSolHDSI,logNSolHDSS,logNSolHIRS,EW_COROT2,logN_COROT2, + EW_TRES4,logN_TRES4,EW_TRES2,logN_TRES2,EW_WASP2,logN_WASP2, + EW_WASP12,logN_WASP12,EW_HD149026,logN_HD149026,f_EW_HD149026, + EW_HATP1,logN_HATP1,EW_XO2S,logN_XO2S,EW_XO2N,logN_XO2N, + EW_XO1,logN_XO1,EW_TRES3,logN_TRES3,f_EW_TRES3,EW_HD189733, + logN_HD189733,f_EW_HD189733,EW_HD80606,logN_HD80606,EW_HATP7, + logN_HATP7,f_EW_HATP7,EW_HATP13,logN_HATP13,EW_HATP16, + logN_HATP16,EW_WASP32,logN_WASP32 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end