Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/718/392 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.5, on 2013-May-18
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJ/718/392 SDSS Lyman limit systems at z~3.5 (Prochaska+, 2010)
*================================================================================
*A definitive survey for Lyman limit systems at z ~ 3.5 with the Sloan Digital
*Sky Survey.
* Prochaska J.X., O'Meara J.M., Worseck G.
* <Astrophys. J., 718, 392-416 (2010)>
* =2010ApJ...718..392P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! SDSS-DR7 quasar sample
integer*4 nr__
parameter (nr__=3759) ! Number of records
character*69 ar__ ! Full-size record
integer*4 Plate ! Plate number
integer*4 MJD ! (d) Modified Julian Date of observation
integer*4 Fiber ! Fiber identification
character*19 SDSS ! Object name (JHHMMSS.ss+DDMMSS.s)
real*4 z ! Emission line derived redshift (1)
integer*4 f_z ! Flag on z (0=normal) (2)
real*4 A ! ? Absorbed continuum fitting parameter A (3)
real*8 B ! ? Absorbed continuum fitting parameter B (3)
real*4 Scale ! ? Additional scaling factor (4)
real*4 S_N ! ? Estimated signal-to-noise (5)
character*1 f_S_N ! [I] I: indicates an infinity value for S/N
*Note (1): List of all objects tagged as QSO or HIZ_QSO in the SDSS Data
* Release 7 with estimated redshift 5{>=}z{>=}3.4.
*Note (2): Redshift flag as follows:
* 0 = Normal;
* 1 = Not at SDSS reported redshift and/or not a quasar (excluded);
* 2 = Too low S/N for evaluation (excluded);
* 3 = Strong BAL (excluded).
*Note (3): Absorbed continuum fitting parameters of the form:
* C=A+Blog_10_({lambda}/{AA}). We caution the reader that these
* models are not meant to describe the intrinsic spectral energy
* distributions of the quasars (see the text).
*Note (4): Imposed by the authors on the best-fit absorbed quasar continuum.
*Note (5): For the absorbed quasar continuum at
* {lambda}={lambda}_912_(1+z_em_).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! SDSS-DR7 intervening {tau}_912_>=2 LLS survey
integer*4 nr__1
parameter (nr__1=607) ! Number of records
character*54 ar__1 ! Full-size record
character*5 Range ! The z_em_ range (< 3.6 or >=3.6) (1)
character*19 SDSS_1 ! Quasar name (JHHMMSS.ss+DDMMSS.s)
real*4 z_1 ! Emission line derived redshift
real*4 zS_N2 ! Starting redshift at S/N=2 (2)
real*4 zS_N3 ! ? Starting redshift at S/N=3 (2)
real*4 zLLS ! ? Lyman Limit System redshift
*Note (1): The quasars with z_em_<3.6 are not included in the final analysis.
*Note (2): The starting redshifts correspond to the wavelength at which the
* absorbed continuum model, starting from z_em_, no longer exceeds the
* smoothed noise array by the specified S/N^A^_912_ threshold. This
* value is limited to a maximum offset from z_em_ of {delta}z=0.4.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! SDSS-DR7 non-statistical LLS systems and candidates
integer*4 nr__2
parameter (nr__2=2593) ! Number of records
character*41 ar__2 ! Full-size record
character*19 SDSS_2 ! Quasar name (JHHMMSS.ss+DDMMSS.s)
real*4 z_2 ! Emission line derived redshift
real*4 zLLS_1 ! Lyman Limit System redshift
real*4 logNHI ! ([cm-2]) Log of HI column density (1)
character*4 S_N_1 ! Signal-to-noise (2)
*Note (1): The logNHI values listed serve as a rough estimate. Typical
* uncertainties for systems with N_HI_<10^17.5^cm^-2^ are at
* least 0.2dex.
*Note (2): Estimate of the S/N of the absorbed continuum (ignoring the effects
* of pLLS) at {lambda}={lambda}_912_(1+z_LLS_).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! SDSS-DR7 proximate {tau}_912_>=2 LLS survey
integer*4 nr__3
parameter (nr__3=468) ! Number of records
character*31 ar__3 ! Full-size record
character*19 SDSS_3 ! Quasar name
real*4 z_3 ! Emission line derived redshift
real*4 zLLS_2 ! ? Lyman Limit System redshift
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! SDSS-DR7 quasars with KECK+LRIS comparison spectra
integer*4 nr__4
parameter (nr__4=32) ! Number of records
character*55 ar__4 ! Full-size record
character*19 SDSS_4 ! SDSS name (JHHMMSS.ss+DDMMSS.s)
real*4 z_4 ! Emission line redshift
real*4 S_N_2 ! Signal to noise of the absorbed continuum at
* the Lyman limit (S/N^A^_912_)
integer*4 f_SDSS ! [1/3] Flag indication (1=no LLS,
* 2={tau}_912_>=2LLS, 3={tau}_912_<2LLS) (1)
real*4 SLLS ! ? SDSS LLS modeling result
real*4 SpLLS ! ? SDSS proximate LLS modeling result
real*4 LLLS ! ? LRIS LLS modeling result
real*4 LpLLS ! ? LRIS proximate LLS modeling result
*Note (1): {tau}_912_ being the optical depth at the Lyman Limit (912{AA}),
* we define a "standard" Lyman Limit System (LLS) to be 1 (LLS=1) when
* {tau}_912_>=2, i.e. N(HI)>=10^17.5^cm^-2^. A flag=2 therefore indicates
* a "strong" LLS, a flag=3 a partial LLS.
C=============================================================================
C Loading file 'table1.dat' ! SDSS-DR7 quasar sample
C Format for file interpretation
1 format(
+ I4,1X,I5,1X,I3,1X,A19,1X,F5.3,1X,I1,1X,F6.2,1X,F7.2,1X,F4.2,
+ 1X,F4.1,1X,A1)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,3759
read(1,'(A69)')ar__
read(ar__,1)Plate,MJD,Fiber,SDSS,z,f_z,A,B,Scale,S_N,f_S_N
if(ar__(44:49) .EQ. '') A = rNULL__
if(ar__(51:57) .EQ. '') B = rNULL__
if(ar__(59:62) .EQ. '') Scale = rNULL__
if(ar__(64:67) .EQ. '') S_N = rNULL__
c ..............Just test output...........
write(6,1)Plate,MJD,Fiber,SDSS,z,f_z,A,B,Scale,S_N,f_S_N
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! SDSS-DR7 intervening {tau}_912_>=2 LLS survey
C Format for file interpretation
2 format(5X,A5,1X,A19,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,607
read(1,'(A54)')ar__1
read(ar__1,2)Range,SDSS_1,z_1,zS_N2,zS_N3,zLLS
if(ar__1(44:48) .EQ. '') zS_N3 = rNULL__
if(ar__1(50:54) .EQ. '') zLLS = rNULL__
c ..............Just test output...........
write(6,2)Range,SDSS_1,z_1,zS_N2,zS_N3,zLLS
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! SDSS-DR7 non-statistical LLS systems and candidates
C Format for file interpretation
3 format(A19,1X,F5.3,1X,F5.3,1X,F4.1,1X,A4)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,2593
read(1,'(A41)')ar__2
read(ar__2,3)SDSS_2,z_2,zLLS_1,logNHI,S_N_1
c ..............Just test output...........
write(6,3)SDSS_2,z_2,zLLS_1,logNHI,S_N_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! SDSS-DR7 proximate {tau}_912_>=2 LLS survey
C Format for file interpretation
4 format(A19,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,file='table6.dat', status='old')
write(6,*) '....Loading file: table6.dat'
do i__=1,468
read(1,'(A31)')ar__3
read(ar__3,4)SDSS_3,z_3,zLLS_2
if(ar__3(27:31) .EQ. '') zLLS_2 = rNULL__
c ..............Just test output...........
write(6,4)SDSS_3,z_3,zLLS_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! SDSS-DR7 quasars with KECK+LRIS comparison spectra
C Format for file interpretation
5 format(
+ A19,1X,F4.2,1X,F4.1,1X,I1,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,file='table8.dat', status='old')
write(6,*) '....Loading file: table8.dat'
do i__=1,32
read(1,'(A55)')ar__4
read(ar__4,5)SDSS_4,z_4,S_N_2,f_SDSS,SLLS,SpLLS,LLLS,LpLLS
if(ar__4(33:37) .EQ. '') SLLS = rNULL__
if(ar__4(39:43) .EQ. '') SpLLS = rNULL__
if(ar__4(45:49) .EQ. '') LLLS = rNULL__
if(ar__4(51:55) .EQ. '') LpLLS = rNULL__
c ..............Just test output...........
write(6,5)SDSS_4,z_4,S_N_2,f_SDSS,SLLS,SpLLS,LLLS,LpLLS
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end