Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/831/7 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 2026-Jul-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/831/7 SDSS-RM project: peak velocities of QSOs (Shen+, 2016)
*================================================================================
*The Sloan Digital Sky Survey Reverberation Mapping project: velocity shifts
*of quasar emission lines.
* Shen Y., Brandt W.N., Richards G.T., Denney K.D., Greene J.E., Grier C.J.,
* Ho L.C., Peterson B.M., Petitjean P., Schneider D.P., Tao C., Trump J.R.
* <Astrophys. J., 831, 7-7 (2016)>
* =2016ApJ...831....7S (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! Peak velocity measurements for each of the 849
quasars in the Sloan Digital Sky Survey
Reverberation Mapping (SDSS-RM) sample
integer*4 nr__
parameter (nr__=849) ! Number of records
character*267 ar__ ! Full-size record
integer*4 RMID ! [0/848]? Unique RMID number (1)
real*4 z ! [0.1/4.4] Pipeline redshift (Z)
real*8 logL1350 ! ([10-7W]) [44.3/47]?=0 Continuum luminosity
* at rest 1350{AA}, in erg/s (LOGL1350)
real*4 e_logL1350 ! ([10-7W]) [0.0003/0.03]? logL1350 uncertainty
real*8 logL1700 ! ([10-7W]) [44.1/47]?=0 Continuum luminosity
* at rest 1350{AA}, in erg/s (LOGL1700)
real*4 e_logL1700 ! ([10-7W]) [0.0002/0.02]? logL1700 uncertainty
real*8 logL3000 ! ([10-7W]) [42.2/46.8]?=0 Continuum luminosity
* at rest 3000{AA}, in erg/s (LOGL3000)
real*4 e_logL3000 ! ([10-7W]) [0.0001/0.04]? logL3000 uncertainty
real*8 logL5100 ! ([10-7W]) [42.6/45.7]?=0 Continuum luminosity
* at rest 5100{AA}, in erg/s (LOGL5100)
real*4 e_logL5100 ! ([10-7W]) [0.0002/0.03]? logL5100 uncertainty
real*8 VHb ! (km/s) [-3025/3088]?=0 Peak velocity (wrt pipeline z)
* for broad H{beta} (HBETA_BR)
real*8 e_VHb ! (km/s) [2/2546]?=-1 VHb uncertainty (2)
real*8 VOIIIc ! (km/s) [-319/341]?=0 Peak velocity (wrt pipeline z)
* for the core [OIII] (OIII5008C)
real*4 e_VOIIIc ! (km/s) [0.9/83]?=-1 VOIIIc uncertainty (2)
real*8 VOIIIf ! (km/s) [-906/760]?=0 Peak velocity (wrt pipeline z)
* for the full [OIII] (OIII5008A)
real*8 e_VOIIIf ! (km/s) [0.8/1227]?=-1 VOIIIf uncertainty (2)
real*8 VCaII ! (km/s) [-1351/802]?=0 Peak velocity (wrt pipeline z)
* for CaII (CAII3934)
real*4 e_VCaII ! (km/s) [6/992]?=-1 VCaII uncertainty (2)
real*8 VOII ! (km/s) [-1094/2085] Peak velocity (wrt pipeline z)
* for [OII] (OII3728)
real*4 e_VOII ! (km/s) [1/858]?=-1 VOII uncertainty (2)
real*8 VNeV ! (km/s) [-2983/3013]?=0 Peak velocity (wrt pipeline z)
* for [NeV]3426 (NEV3426)
real*8 e_VNeV ! (km/s) [0.9/1035]?=-1 VNeV uncertainty (2)
real*8 VMgII ! (km/s) [-2699/3292]?=0 Peak velocity (wrt pipeline z)
* for MgII (MGII)
real*8 e_VMgII ! (km/s) [8/2003]?=-1 VMgII uncertainty (2)
real*8 VCIII ! (km/s) [-4513/4573]?=0 Peak velocity (wrt pipeline z)
* for CIII] (CIII)
real*8 e_VCIII ! (km/s) [8/4426]?=-1 VCIII uncertainty (2)
real*8 VComp ! (km/s) [-8937/4573]?=0 Peak velocity (wrt pipeline z)
* for the CIII]/SiIII]/AlIII complex (CIIIA)
real*8 e_VComp ! (km/s) [8/4426]?=-1 VComp uncertainty (2)
real*8 VHeII ! (km/s) [-2390/1963]?=0 Peak velocity (wrt pipeline z)
* for HeII1640 (HEII1640)
real*8 e_VHeII ! (km/s) [0.7/1634]?=-1 VHeII uncertainty (2)
real*8 VCIV ! (km/s) [-4316/2226]?=0 Peak velocity (wrt pipeline z)
* for CIV (CIV)
real*8 e_VCIV ! (km/s) [6/2140]?=-1 VCIV uncertainty (2)
real*8 VSi_OIV ! (km/s) [-2882/2744]?=0 Peak velocity (wrt pipeline z)
* for SiIV/OIV] (SIIV)
real*4 e_VSi_OIV ! (km/s) [11/843]?=-1 VSi/OIV uncertainty (2)
*Note (1): More information about the SDSS-RM sample can be found in the
* technical overview paper (Shen+, 2015, J/ApJS/216/4).
*Note (2): The uncertainty is -1 if the measurement is unavailable
* (e.g., line not covered).
* In the paper, we only used measurements with velocity uncertainties
* <500km/s.
C=============================================================================
C Loading file 'catalog.dat' ! Peak velocity measurements for each of the 849
* quasars in the Sloan Digital Sky Survey
* Reverberation Mapping (SDSS-RM) sample
C Format for file interpretation
1 format(
+ I3,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,
+ F7.4,1X,F6.4,1X,F8.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F7.2,1X,F7.2,
+ 1X,F8.2,1X,F6.2,1X,F8.2,1X,F6.2,1X,F8.2,1X,F7.2,1X,F8.2,1X,
+ F7.2,1X,F8.2,1X,F7.2,1X,F8.2,1X,F7.2,1X,F8.2,1X,F7.2,1X,F8.2,
+ 1X,F7.2,1X,F8.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,849
read(1,'(A267)')ar__
read(ar__,1)
+ RMID,z,logL1350,e_logL1350,logL1700,e_logL1700,logL3000,
+ e_logL3000,logL5100,e_logL5100,VHb,e_VHb,VOIIIc,e_VOIIIc,
+ VOIIIf,e_VOIIIf,VCaII,e_VCaII,VOII,e_VOII,VNeV,e_VNeV,VMgII,
+ e_VMgII,VCIII,e_VCIII,VComp,e_VComp,VHeII,e_VHeII,VCIV,e_VCIV,
+ VSi_OIV,e_VSi_OIV
if(ar__(1:3) .EQ. '') RMID = iNULL__
if(ar__(20:25) .EQ. '') e_logL1350 = rNULL__
if(ar__(35:40) .EQ. '') e_logL1700 = rNULL__
if(ar__(50:55) .EQ. '') e_logL3000 = rNULL__
if(ar__(65:70) .EQ. '') e_logL5100 = rNULL__
c ..............Just test output...........
write(6,1)
+ RMID,z,logL1350,e_logL1350,logL1700,e_logL1700,logL3000,
+ e_logL3000,logL5100,e_logL5100,VHb,e_VHb,VOIIIc,e_VOIIIc,
+ VOIIIf,e_VOIIIf,VCaII,e_VCaII,VOII,e_VOII,VNeV,e_VNeV,VMgII,
+ e_VMgII,VCIII,e_VCIII,VComp,e_VComp,VHeII,e_VHeII,VCIV,e_VCIV,
+ VSi_OIV,e_VSi_OIV
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end