Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/685/773 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-Apr-20 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/685/773 X-ray-UV relations in SDSS DR5 QSOs (Gibson+, 2008) *================================================================================ *Are optically selected quasars universally X-ray luminous? *X-ray-UV relations in Sloan Digital Sky Survey quasars. * Gibson R.R., Brandt W.N., Schneider D.P. * <Astrophys. J., 685, 773-786 (2008)> * =2008ApJ...685..773G C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Source information integer*4 nr__ parameter (nr__=536) ! Number of records character*111 ar__ ! Full-size record character*19 SDSS ! SDSS source name JHHMMSS.ss+DDMMSS.s real*4 z ! Spectroscopic redshift integer*8 ObsId ! Chandra or XMM observation identification * number (1) real*4 Expo ! (ks) X-ray exposure time (2) integer*4 Scts ! (ct) Observed-frame soft (0.5-2keV) counts * in source region integer*4 Hcts ! (ct) Observed-frame hard (2-8keV) counts * in source region integer*4 X ! [0/1] X-ray detected? (1=yes) integer*4 B ! [0/1] Source in sample B? (1=yes) character*1 l_logR_ ! Limit flag on logR* real*4 logR_ ! ([-]) Radio-loudness parameter (5) integer*4 CIV_0 ! (km/s) Rest-frame balnicity index * for CIV absorption (6) integer*4 Mg ! Number of narrow Mg II absorption systems (7) real*4 L1400 ! ([10-7W/Hz]) Log of monochromatic luminosity at * 1400{AA} (8) real*4 L2500 ! ([10-7W/Hz]) Log of monochromatic luminosity at * 2500{AA} (8) character*1 l_L2keV ! Limit flag on log2keV real*4 L2keV ! ([10-7W/Hz]) Log of monochromatic luminosity at 2keV (8) character*1 l_a_OX ! Limit flag on a.OX real*4 a_OX ! X-ray/optical spectral index (9) real*4 e_a_OX ! ? Uncertainty in a.OX (10) character*1 l_Da_OX ! Limit flag on Da.OX real*4 Da_OX ! Relative X-ray brightness (11) *Note (1): The Chandra or XMM-Newton observation identification number. * Ten-digit numbers correspond to XMM-Newton observations, while shorter * numbers correspond to Chandra observations. *Note (2): Reported by the CIAO or SAS toolchains for the extraction of the * source region. *Note (5): Where log(R*)=log[L_{nu}_(5GHz)/L_{nu}_(2500{AA})]. We classify * sources with log(R*)>=1 as "radio-loud" and all other sources as * "radio-quiet". *Note (6): Similar to the traditional "balnicity index", BI, of Weymann et * al. (1991ApJ...373...23W), except that it is integrated from velocity * offsets -25000km/s to 0km/s from the QSO rest frame. *Note (7): Found using the method described in Section 2.2. *Note (8): Where the luminosity is in units of erg/s/Hz. *Note (9): Where a.OX = {alpha}OX = 0.3838*log(L(2keV)/L(2500{AA})). *Note (10): Errors assumed to be dominated by the error on X-ray luminosity. *Note (11): Relative to "ordinary" QSOs. Where * Da.OX = {Delta}{alpha}OX = {alpha}OX - * (-0.217+/-0.036)*log(L_2500{AA}_)+(5.075+/-1.118). * (for instance, a value of -0.5 corresponds to X-ray weakness * by a factor of ~20 compared to QSOs with the same UV luminosity) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! UV emission line measurements integer*4 nr__1 parameter (nr__1=536) ! Number of records character*55 ar__1 ! Full-size record character*19 SDSS_1 ! SDSS source identification, JHHMMSS.ss+DDMMSS.s real*4 C4EW ! (0.1nm) Equivalent width of CIV (1549{AA}) real*4 C4Wave ! (0.1nm) Central wavelength of CIV (1549{AA}) real*4 C4FWHM ! (0.1nm) Full-Width at Half-Maximum of CIV (1549{AA}) real*4 C3EW ! (0.1nm) ? Equivalent width of CIII] (1909{AA}) (1) real*4 C3Wave ! (0.1nm) ? Central wavelength of CIII] (1909{AA}) (1) real*4 C3FWHM ! (0.1nm) ? Full-Width at Half-Maximum of * CIII] (1909{AA}) (1) *Note (1): Blank entries indicate that the C III] emission line could not be * fit because the spectrum was bad in that region. C============================================================================= C Loading file 'table1.dat' ! Source information C Format for file interpretation 1 format( + A19,1X,F4.2,1X,I10,1X,F5.1,1X,I4,1X,I3,1X,I1,1X,I1,1X,A1,1X, + F5.2,1X,I5,1X,I1,1X,F5.2,1X,F5.2,1X,A1,1X,F5.2,1X,A1,1X,F5.2, + 1X,F4.2,1X,A1,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,536 read(1,'(A111)')ar__ read(ar__,1) + SDSS,z,ObsId,Expo,Scts,Hcts,X,B,l_logR_,logR_,CIV_0,Mg,L1400, + L2500,l_L2keV,L2keV,l_a_OX,a_OX,e_a_OX,l_Da_OX,Da_OX if(ar__(100:103) .EQ. '') e_a_OX = rNULL__ c ..............Just test output........... write(6,1) + SDSS,z,ObsId,Expo,Scts,Hcts,X,B,l_logR_,logR_,CIV_0,Mg,L1400, + L2500,l_L2keV,L2keV,l_a_OX,a_OX,e_a_OX,l_Da_OX,Da_OX c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! UV emission line measurements C Format for file interpretation 2 format(A19,1X,F5.1,1X,F6.1,1X,F4.1,1X,F4.1,1X,F6.1,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,536 read(1,'(A55)')ar__1 read(ar__1,2)SDSS_1,C4EW,C4Wave,C4FWHM,C3EW,C3Wave,C3FWHM if(ar__1(39:42) .EQ. '') C3EW = rNULL__ if(ar__1(44:49) .EQ. '') C3Wave = rNULL__ if(ar__1(51:55) .EQ. '') C3FWHM = rNULL__ c ..............Just test output........... write(6,2)SDSS_1,C4EW,C4Wave,C4FWHM,C3EW,C3Wave,C3FWHM c .......End.of.Just test output........... end do close(1) C============================================================================= stop end