Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/721/1 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-Jun-10
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/721/1 Keck ESI observations of Ly{alpha} systems (Penprase+, 2010)
*================================================================================
*Keck Echellette Spectrograph and Imager observations of metal-poor damped
*Ly{alpha} systems.
* Penprase B.E., Prochaska J.X., Sargent W.L.W., Toro-Martinez I.,
* Beeler D.J.
* <Astrophys. J., 721, 1-25 (2010)>
* =2010ApJ...721....1P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! SDSS-DR5 metal-poor DLA candidates
integer*4 nr__
parameter (nr__=405) ! Number of records
character*31 ar__ ! Full-size record
character*19 Name ! QSO identification (JHMMSS.ss+DDMMSS or
* JHHMMSS.ss+DDMMSS) (1)
real*4 zabs ! Absorption redshift
real*4 logNHI ! ([cm-2]) Log of the HI column density
*Note (1): List of all candidate metal-poor DLAs flagged by Prochaska and
* collaborators from their analysis of the SDSS-DR5 (See
* Cat. J/ApJ/635/123 and 2009ApJ...696.1543P). We caution that
* many of these will not be bona-fide DLAs.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Targets observed, magnitudes, and exposure times
integer*4 nr__1
parameter (nr__1=44) ! 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*9 QSO ! Quasar designation (HHMM+DDMM; J2000) (G1)
real*4 zDLA ! DLA redshift
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 declination (J2000)
integer*4 DEd ! (deg) Degree of declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs ! (arcsec) Arcsecond of declination (J2000)
real*4 Rmag ! (mag) R-band magnitude
integer*4 Exp08 ! (s) ? Exposure time (2008)
integer*4 Exp07 ! (s) ? Exposure time (2007)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! Elemental abundances for DLA systems
integer*4 nr__2
parameter (nr__2=35) ! Number of records
character*85 ar__2 ! Full-size record
character*9 QSO_1 ! Quasar designation (HHMM+DDMM; J2000) (G1)
real*4 zDLA_1 ! DLA redshift
real*4 logNHI_1 ! ([cm-2]) Log of HI column density
real*4 e_logNHI ! ([cm-2]) logNHI uncertainty
character*1 l__C_H_ ! Limit flag on [C/H]
real*4 v_C_H_ ! ([-]) ? Abundance [C/H]
character*1 f__C_H_ ! [abc] Flag on [C/H]
character*1 l__O_H_ ! Limit flag on [O/H]
real*4 v_O_H_ ! ([-]) ? Abundance [O/H]
character*1 f__O_H_ ! [abc] Flag on [O/H]
character*1 l__Si_H_ ! Limit flag on [Si/H]
real*4 v_Si_H_ ! ([-]) ? Abundance [Si/H]
character*1 f__Si_H_ ! [abc] Flag on [Si/H]
character*1 l__Al_H_ ! Limit flag on [Al/H]
real*4 v_Al_H_ ! ([-]) ? Abundance [Al/H]
character*1 f__Al_H_ ! [abc] Flag on [Al/H]
character*1 l__Fe_H_ ! Limit flag on [Fe/H]
real*4 v_Fe_H_ ! ([-]) ? Abundance [Fe/H]
character*1 f__Fe_H_ ! [abc] Flag on [Fe/H]
character*1 l__Mg_H_ ! Limit flag on [Mg/H]
real*4 v_Mg_H_ ! ([-]) ? Abundance [Mg/H]
character*1 f__Mg_H_ ! [abc] Flag on [Mg/H]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table14.dat' ! Column densities for the species CIV and SiIV
integer*4 nr__3
parameter (nr__3=41) ! Number of records
character*52 ar__3 ! Full-size record
character*9 QSO_2 ! Quasar designation (HHMM+DDMM; J2000) (G1)
real*4 zDLA_2 ! DLA redshift
real*4 logNHI_2 ! ([cm-2]) Log of HI column density
character*1 l_logNCIV ! Limit flag on logNCIV
real*4 logNCIV ! ([cm-2]) Log of CIV column density
real*4 e_logNCIV ! ([cm-2]) ? logNCIV uncertainty
character*1 l_logNSiIV ! Limit flag on logNSiIV
real*4 logNSiIV ! ([cm-2]) ? Log of SiIV column density
real*4 e_logNSiIV ! ([cm-2]) ? logNSiIV uncertainty
C=============================================================================
C Loading file 'table1.dat' ! SDSS-DR5 metal-poor DLA candidates
C Format for file interpretation
1 format(A19,1X,F5.3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,405
read(1,'(A31)')ar__
read(ar__,1)Name,zabs,logNHI
c ..............Just test output...........
write(6,1)Name,zabs,logNHI
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Targets observed, magnitudes, and exposure times
C Format for file interpretation
2 format(
+ 4X,A9,1X,F5.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,I2,1X,
+ F5.2,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,44
read(1,'(A57)')ar__1
read(ar__1,2)
+ QSO,zDLA,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Rmag,Exp08,Exp07
if(ar__1(49:52) .EQ. '') Exp08 = iNULL__
if(ar__1(54:57) .EQ. '') Exp07 = 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)
+ QSO,zDLA,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Rmag,Exp08,Exp07
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! Elemental abundances for DLA systems
C Format for file interpretation
3 format(
+ 4X,A9,1X,F6.4,1X,F5.2,1X,F4.2,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,
+ A1,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,1X,A1,
+ F5.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,35
read(1,'(A85)')ar__2
read(ar__2,3)
+ QSO_1,zDLA_1,logNHI_1,e_logNHI,l__C_H_,v_C_H_,f__C_H_,l__O_H_,
+ v_O_H_,f__O_H_,l__Si_H_,v_Si_H_,f__Si_H_,l__Al_H_,v_Al_H_,
+ f__Al_H_,l__Fe_H_,v_Fe_H_,f__Fe_H_,l__Mg_H_,v_Mg_H_,f__Mg_H_
if(ar__2(34:38) .EQ. '') v_C_H_ = rNULL__
if(ar__2(43:47) .EQ. '') v_O_H_ = rNULL__
if(ar__2(52:56) .EQ. '') v_Si_H_ = rNULL__
if(ar__2(61:65) .EQ. '') v_Al_H_ = rNULL__
if(ar__2(70:74) .EQ. '') v_Fe_H_ = rNULL__
if(ar__2(79:83) .EQ. '') v_Mg_H_ = rNULL__
c ..............Just test output...........
write(6,3)
+ QSO_1,zDLA_1,logNHI_1,e_logNHI,l__C_H_,v_C_H_,f__C_H_,l__O_H_,
+ v_O_H_,f__O_H_,l__Si_H_,v_Si_H_,f__Si_H_,l__Al_H_,v_Al_H_,
+ f__Al_H_,l__Fe_H_,v_Fe_H_,f__Fe_H_,l__Mg_H_,v_Mg_H_,f__Mg_H_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table14.dat' ! Column densities for the species CIV and SiIV
C Format for file interpretation
4 format(
+ 4X,A9,1X,F6.4,1X,F5.2,1X,A1,F5.2,1X,F5.3,1X,A1,F5.2,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table14.dat')
write(6,*) '....Loading file: table14.dat'
do i__=1,41
read(1,'(A52)')ar__3
read(ar__3,4)
+ QSO_2,zDLA_2,logNHI_2,l_logNCIV,logNCIV,e_logNCIV,l_logNSiIV,
+ logNSiIV,e_logNSiIV
if(ar__3(35:39) .EQ. '') e_logNCIV = rNULL__
if(ar__3(42:46) .EQ. '') logNSiIV = rNULL__
if(ar__3(48:52) .EQ. '') e_logNSiIV = rNULL__
c ..............Just test output...........
write(6,4)
+ QSO_2,zDLA_2,logNHI_2,l_logNCIV,logNCIV,e_logNCIV,l_logNSiIV,
+ logNSiIV,e_logNSiIV
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end