Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/613/682 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-19
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/613/682 AGN central masses and broad-line region sizes (Peterson+, 2004)
*================================================================================
*Central masses and broad-line region sizes of active galactic nuclei.
*II. A homogeneous analysis of a large reverberation-mapping database.
* Peterson B.M., Ferrarese L., Gilbert K.M., Kaspi S., Malkan M.A., Maoz D.,
* Merritt D., Netzer H., Onken C.A., Pogge R.W., Vestergaard M., Wandel A.
* <Astrophys. J., 613, 682-699 (2004)>
* =2004ApJ...613..682P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Object List
integer*4 nr__
parameter (nr__=37) ! Number of records
character*129 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*11 Name ! Name
character*18 r_Name ! Reference on Name
integer*4 RAh ! (h) Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*8 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! Sign of the declination (J2000)
integer*4 DEd ! (deg) Degree of declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of declination (J2000)
real*8 z ! (km/s) Redshift
real*4 Vmag ! (mag) V magnitude
real*4 AB ! (mag) B-band extinction
character*16 OName1 ! Veron-Cetty & Veron Catalog (Cat. <VII/224>) Name
character*30 OName2 ! Alternative names
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Wavelength Windows for Line Width Measurements
integer*4 nr__1
parameter (nr__1=117) ! Number of records
character*88 ar__1 ! Full-size record
character*11 Name_1 ! Common object name
character*5 r_Name_1 ! Reference on Name
character*18 Line ! Individual emission line
real*4 Res ! (0.1nm) Spectral resolution in Angstroms
character*7 Date1 ! (d) Julian date of start of observation
character*7 Date2 ! (d) Julian date of end of observation
integer*4 RWin1 ! (0.1nm) Red continuum window in Angstroms (low)
integer*4 RWin2 ! (0.1nm) Red continuum window in Angstroms (up)
integer*4 BWin1 ! (0.1nm) Blue continuum window in Angstroms (low)
integer*4 BWin2 ! (0.1nm) Blue continuum window in Angstroms (up)
integer*4 Lim1 ! (0.1nm) Line limit (low)
integer*4 Lim2 ! (0.1nm) Line limit (up)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__2
parameter (nr__2=27) ! Number of records
character*76 ar__2 ! Full-size record
integer*4 Ref ! Reference number
character*19 Bibcode ! Bibcode
character*24 Aut ! Author's name
character*25 Com ! Comments
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Time-Series Analysis
integer*4 nr__3
parameter (nr__3=117) ! Number of records
character*115 ar__3 ! Full-size record
character*11 Name_2 ! Common object name
character*18 Line_1 ! Individual emission line
character*7 Date1_1 ! (d) Julian date of start of observation
character*7 Date2_1 ! (d) Julian date of end of observation
real*4 Rmax ! Cross-correlation function peak amplitude
real*4 e_Rmax ! The RMS uncertainty in Rmax
real*4 FvarC ! Continuum noise-corrected fractional variation
* (1)
real*4 FvarL ! Line noise-corrected fractional variation (1)
character*1 u_tcent ! [:] Uncertainty flag in tcent
real*4 tcent ! (d) Cross-correlation function centroid peak value
real*4 E_tcent ! (d) Upper error limit on tcent
real*4 e_tcent_1 ! (d) Lower error limit on tcent
character*1 u_tpeak ! [:] Uncertainty flag in tpeak
real*4 tpeak ! (d) Cross-correlation function centroid peak
* location
real*4 E_tpeak ! (d) Upper error limit on tpeak
real*4 e_tpeak_1 ! (d) Lower error limit on tpeak
*Note (1): cf. Rodriguez-Pascual et al. (1997ApJS..110....9R).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Rest-Frame Lags, Line Widths, and Virial
Products
integer*4 nr__4
parameter (nr__4=121) ! Number of records
character*116 ar__4 ! Full-size record
character*11 Name_3 ! Common object name
character*18 Line_2 ! Individual emission line
character*1 u_tcent_1 ! [:] Uncertainty flag in tcent
real*4 tcent_1 ! (d) Cross-correlation function centroid peak
* value (1)
real*4 E_tcent_2 ! (d) Upper error limit on tcent
real*4 e_tcent_3 ! (d) Lower error limit on tcent
character*1 u_tpeak_1 ! [:] Uncertainty flag in tpeak
real*4 tpeak_1 ! (d) Cross-correlation function centroid peak
* location (1)
real*4 E_tpeak_2 ! (d) Upper error limit on tpeak
real*4 e_tpeak_3 ! (d) Lower error limit on tpeak
character*1 u_LDisp ! [:] Uncertainty flag in LDisp
integer*4 LDisp ! (km/s) Line dispersion
integer*4 e_LDisp ! (km/s) Uncertainty in LDisp
character*1 f_Line ! Flag on Line (2)
character*1 u_FWHM ! [:] Uncertainty flag in FWHM
integer*4 FWHM ! (km/s) Line Full-Width at Half-Maximum
integer*4 e_FWHM ! (km/s) Uncertainty in FWHM
character*1 f_FWHM ! Flag on FWHM (2)
character*1 u_Mass ! [:]Uncertainty in Mass
real*4 Mass ! (10+6solMass) Virial mass
real*4 E_Mass ! (10+6solMass) Upper error limit on Mass
real*4 e_Mass_1 ! (10+6solMass) Lower error limit on Mass
*Note (1): Corrected for time dilation by dividing by 1+z.
*Note (2): Flag on Line and on FWHM:
* 1 = Measurements for blended line;
* 2 = Measurements based on unblended half of line.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Adopted Virial Products and Derived Black Hole
Masses
integer*4 nr__5
parameter (nr__5=35) ! Number of records
character*78 ar__5 ! Full-size record
character*11 Name_4 ! Name
real*8 VP ! (10+6solMass) Mean virial product (1)
character*1 f_VP ! Uncertainty flag on VP
real*8 E_VP ! (10+6solMass) Upper limit error on VP
real*8 e_VP_1 ! (10+6solMass) RMS uncertainty on VP or lower limit error
* if E_VP exists
real*8 MBH ! (10+6solMass) Black hole mass
character*1 u_MBH ! Uncertainty flag on MBH
real*8 E_MBH ! (10+6solMass) Upper limit on MBH
real*8 e_MBH_1 ! (10+6solMass) RMS uncertainty on MBH or lower limit error
* if E_MBH exists
real*4 logLum ! ([10-7W]) Optical luminosity log({lambda}L_{lambda}_)
* in units of ergs/s
real*4 e_logLum ! ([10-7W]) Error in logLum
*Note (1): The mean virial product is <c{tau}_cent_{sigma}_line_^2^/G>, where
* * c is the speed of light,
* * {tau}_cent_ is the emission-line lag for the central object relative
* to continuum variations and
* * G is the gravitational constant.
C=============================================================================
C Loading file 'table1.dat' ! Object List
C Format for file interpretation
1 format(
+ A11,5X,A18,1X,I2,1X,I2,1X,F7.4,1X,A1,I2,1X,I2,1X,F6.3,1X,F7.5,
+ 1X,F4.1,1X,F5.3,1X,A16,1X,A30)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,37
read(1,'(A129)')ar__
read(ar__,1)
+ Name,r_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,Vmag,AB,OName1,
+ OName2
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,r_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,z,Vmag,AB,OName1,
+ OName2
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Wavelength Windows for Line Width Measurements
C Format for file interpretation
2 format(
+ A11,1X,A5,1X,A18,1X,F4.1,1X,A7,1X,A7,1X,I4,1X,I4,1X,I4,1X,I4,
+ 1X,I4,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,117
read(1,'(A88)')ar__1
read(ar__1,2)
+ Name_1,r_Name_1,Line,Res,Date1,Date2,RWin1,RWin2,BWin1,BWin2,
+ Lim1,Lim2
c ..............Just test output...........
write(6,2)
+ Name_1,r_Name_1,Line,Res,Date1,Date2,RWin1,RWin2,BWin1,BWin2,
+ Lim1,Lim2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
3 format(I2,2X,A19,2X,A24,2X,A25)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,27
read(1,'(A76)')ar__2
read(ar__2,3)Ref,Bibcode,Aut,Com
c ..............Just test output...........
write(6,3)Ref,Bibcode,Aut,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Time-Series Analysis
C Format for file interpretation
4 format(
+ A11,1X,A18,1X,A7,1X,A7,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,A1,
+ F5.1,1X,F4.1,1X,F5.1,1X,A1,F5.1,1X,F5.1,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,117
read(1,'(A115)')ar__3
read(ar__3,4)
+ Name_2,Line_1,Date1_1,Date2_1,Rmax,e_Rmax,FvarC,FvarL,u_tcent,
+ tcent,E_tcent,e_tcent_1,u_tpeak,tpeak,E_tpeak,e_tpeak_1
c ..............Just test output...........
write(6,4)
+ Name_2,Line_1,Date1_1,Date2_1,Rmax,e_Rmax,FvarC,FvarL,u_tcent,
+ tcent,E_tcent,e_tcent_1,u_tpeak,tpeak,E_tpeak,e_tpeak_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Rest-Frame Lags, Line Widths, and Virial
* Products
C Format for file interpretation
5 format(
+ A11,1X,A18,1X,A1,F5.1,1X,F4.1,1X,F5.1,1X,A1,F5.1,1X,F5.1,1X,
+ F5.1,1X,A1,I4,1X,I3,1X,A1,1X,A1,I5,1X,I4,1X,A1,1X,A1,F6.2,1X,
+ F6.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,121
read(1,'(A116)')ar__4
read(ar__4,5)
+ Name_3,Line_2,u_tcent_1,tcent_1,E_tcent_2,e_tcent_3,u_tpeak_1,
+ tpeak_1,E_tpeak_2,e_tpeak_3,u_LDisp,LDisp,e_LDisp,f_Line,
+ u_FWHM,FWHM,e_FWHM,f_FWHM,u_Mass,Mass,E_Mass,e_Mass_1
c ..............Just test output...........
write(6,5)
+ Name_3,Line_2,u_tcent_1,tcent_1,E_tcent_2,e_tcent_3,u_tpeak_1,
+ tpeak_1,E_tpeak_2,e_tpeak_3,u_LDisp,LDisp,e_LDisp,f_Line,
+ u_FWHM,FWHM,e_FWHM,f_FWHM,u_Mass,Mass,E_Mass,e_Mass_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Adopted Virial Products and Derived Black Hole
* Masses
C Format for file interpretation
6 format(
+ A11,1X,F8.3,A1,F8.3,1X,F8.3,1X,F8.3,A1,F8.3,1X,F8.3,2X,F5.2,
+ 2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,35
read(1,'(A78)')ar__5
read(ar__5,6)
+ Name_4,VP,f_VP,E_VP,e_VP_1,MBH,u_MBH,E_MBH,e_MBH_1,logLum,
+ e_logLum
c ..............Just test output...........
write(6,6)
+ Name_4,VP,f_VP,E_VP,e_VP_1,MBH,u_MBH,E_MBH,e_MBH_1,logLum,
+ e_logLum
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end