Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/847/81 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/847/81 HeII photoionization rates of quasars (Schmidt+, 2017)
*================================================================================
*Statistical detection of the HeII transverse proximity effect: evidence for
*sustained quasar activity for >25 million years.
* Schmidt T.M., Worseck G., Hennawi J.F., Prochaska J.X., Crighton N.H.M.
* <Astrophys. J., 847, 81 (2017)>
* =2017ApJ...847...81S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Overview of the FUV spectra used for this work
integer*4 nr__
parameter (nr__=22) ! Number of records
character*55 ar__ ! Full-size record
character*16 HE2QS ! He II-transparent quasar ID
character*10 Inst ! HST instrument +grating
integer*4 R ! [1000/18000] Dispersion
integer*4 HST ! [9350/13013] HST Program Id
character*8 PI ! HST Program PI
character*6 Ref ! Reference(s) (1)
*Note (1): References as follows:
* 1 = Reimers et al. (2005A&A...442...63R);
* 2 = Shull et al. (2010ApJ...722.1312S);
* 3 = Syphers et al. (2011ApJ...742...99S);
* 4 = Syphers et al. (2012AJ....143..100S);
* 5 = Syphers et al. (2013ApJ...765..119S);
* 6 = Syphers et al. (2014ApJ...784...42S);
* 7 = Worseck et al., in prep
* 8 = Worseck et al. (2011ApJ...733L..24W);
* 9 = Worseck et al. (2016ApJ...825..144W);
* 10 = Zheng et al. (2015ApJ...806..142Z);
* 11 = this paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Overview of the HeII sightlines and the number of
foreground quasars
integer*4 nr__1
parameter (nr__1=22) ! Number of records
character*57 ar__1 ! 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*16 HE2QS_1 ! He II-transparent quasar ID
integer*4 RAh ! (h) Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*4 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*4 zBG ! [2.7/4] Spectroscopic redshift of HeII-transparent
* QSO
integer*4 N_VLT ! [0/2]? Number of quasars detected with
* VLT/VIMOS+LBT/LBC
integer*4 N_NTT ! [0/2]? Number of quasars detected with NTT/EFOSC2
integer*4 N_CAHA ! [0/5]? Number of quasars detected with CAHA/TWIN
integer*4 N_SDSS ! [0/7]? Number of quasars used from SDSS (1)
integer*4 N_Lit ! [1/2]? Number of quasars used from the
* Literature (2)
*Note (1): From SDSS DR12, Alam+ (2015, V/147); Paris+ (2017, VII/279).
*Note (2): Literature objects from Jakobsen+ (2003A&A...397..891J);
* Steidel+ (2003, J/ApJ/592/728); Hennawi+ (2006, J/ApJ/651/61);
* Worseck & Wisotzki (2006A&A...450..495W);
* Worseck+ (2007A&A...473..805W); Syphers & Shull (2013ApJ...765..119S).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Quasars discovered within our survey
integer*4 nr__2
parameter (nr__2=131) ! Number of records
character*86 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 HE2QS_2 ! He II-transparent quasar ID (1)
real*8 RAdeg_1 ! (deg) Right Ascension, foreground QSOs (J2000)
real*8 DEdeg_1 ! (deg) Declination, foreground QSOs (J2000)
real*4 zspec ! [0.6/4.2] Spectroscopic Redshift
real*4 rmag ! (mag) [17/24.5]? SDSS r band AB magnitude
real*4 M1450 ! (mag) [-27.8/-20.5]? Absolute M1450 magnitude
real*4 Sep ! (arcsec) [3.1/80.4] Angular separation
real*4 Dprop ! (Mpc) [1.5/39.8] Transverse separation in situation
* proper distances; in proper Mpc units
real*4 GamHeQSO ! (s-1) [6.8e-17/1.1e-06] HeII photoionization rate
character*13 Inst_1 ! Instrument used (CAHA 3.5/TWIN, VLT/VIMOS,
* NTT/EFOSC2 or Keck/LRIS; see the "Description"
* section above)
*Note (1): For the sightlines of SDSS J0139-0847, LBQS 1216+1656,
* SDSS J2200+0008 and J2251-0857 no HeII spectra are available.
* Despite that, we still list discovered forground quasars.
C=============================================================================
C Loading file 'table1.dat' ! Overview of the FUV spectra used for this work
C Format for file interpretation
1 format(A16,1X,A10,1X,I5,1X,I5,1X,A8,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,22
read(1,'(A55)')ar__
read(ar__,1)HE2QS,Inst,R,HST,PI,Ref
c ..............Just test output...........
write(6,1)HE2QS,Inst,R,HST,PI,Ref
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Overview of the HeII sightlines and the number of
* foreground quasars
C Format for file interpretation
2 format(
+ A16,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.3,1X,I1,
+ 1X,I1,1X,I1,1X,I1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,22
read(1,'(A57)')ar__1
read(ar__1,2)
+ HE2QS_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,zBG,N_VLT,N_NTT,N_CAHA,
+ N_SDSS,N_Lit
if(ar__1(48:48) .EQ. '') N_VLT = iNULL__
if(ar__1(50:50) .EQ. '') N_NTT = iNULL__
if(ar__1(52:52) .EQ. '') N_CAHA = iNULL__
if(ar__1(54:54) .EQ. '') N_SDSS = iNULL__
if(ar__1(56:57) .EQ. '') N_Lit = iNULL__
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,2)
+ HE2QS_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,zBG,N_VLT,N_NTT,N_CAHA,
+ N_SDSS,N_Lit
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Quasars discovered within our survey
C Format for file interpretation
3 format(
+ A16,1X,F9.5,1X,F9.5,1X,F5.3,1X,F4.1,1X,F5.1,1X,F4.1,1X,F4.1,
+ 1X,E8.2,1X,A13)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,131
read(1,'(A86)')ar__2
read(ar__2,3)
+ HE2QS_2,RAdeg_1,DEdeg_1,zspec,rmag,M1450,Sep,Dprop,GamHeQSO,
+ Inst_1
if(ar__2(44:47) .EQ. '') rmag = rNULL__
if(ar__2(49:53) .EQ. '') M1450 = rNULL__
c ..............Just test output...........
write(6,3)
+ HE2QS_2,RAdeg_1,DEdeg_1,zspec,rmag,M1450,Sep,Dprop,GamHeQSO,
+ Inst_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end