Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/679/239 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-Jun-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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJ/679/239 Associated MgII absorbers (Vanden+, 2008)
*================================================================================
*Average properties of a large sample of z_abs_~ z_em_Associated Mg II
*absorption line systems.
* Vanden Berk D., Khare P., York D.G., Richards G.T., Lundgren B.,
* Alsayyad Y., Kulkarni V.P., Subbarao M., Schneider D.P., Heckman T.,
* Anderson S., Crotts A.P.S., Frieman J., Stoughton C., Lauroesch J.T.,
* Hall P.B., Meiksin A., Steffing M., Vanlandingham J.
* <Astrophys. J., 679, 239-259 (2008)>
* =2008ApJ...679..239V
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! The full absorber and matching non-absorber samples
integer*4 nr__
parameter (nr__=415) ! Number of records
character*87 ar__ ! Full-size record
integer*4 Plate ! Absorber plate number
integer*4 Fiber ! Absorber fiber number
integer*4 MJD ! (d) Absorber Modified Julian Date
real*4 zem_A ! Absorber emission line redshift
real*8 Dg_i_A ! (mag) Absorber {Delta}(g-i) (1)
real*4 zab_A ! Absorber absorption line redshift
real*4 imag_A ! (mag) Absorber SDSS i band magnitude
real*4 beta_A ! Absorber {beta} value (2)
integer*4 Plate2 ! Non-absorber plate number
integer*4 Fiber2 ! Non-absorber fiber number
integer*4 MJD2 ! (d) Non-absorber Modified Julian Date
real*4 zem_N ! Non-absorber emission line redshift
real*8 Dg_i_N ! (mag) Non-absorber {Delta}(g-i) (1)
real*4 imag_N ! (mag) Non-absorber SDSS i band magnitude
*Note (1): Defined as the difference between the (g-i) color of a QSO and
* the median value of (g-i) of all other verified SDSS QSOs with nearly
* the same redshift (Richards et al., 2003AJ....126.1131R).
*Note (2): {beta} =
* ((1+z_em_)^2^-(1+z_abs_)^2^)/((1+z_em_)^2^+(1+z_abs_)^2^).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Equivalent widths of chosen lines
integer*4 nr__1
parameter (nr__1=415) ! Number of records
character*122 ar__1 ! Full-size record
integer*4 Plate_1 ! Absorber plate number
integer*4 Fiber_1 ! Absorber fiber number
integer*4 MJD_1 ! (d) Absorber Modified Julian Date
real*4 AlII ! (0.1nm) ?=-1.0 Equivalent width of AlII (1670{AA}) (1)
real*4 e_AlII ! (0.1nm) ?=-1.0 The 1{sigma} error in AlII (1)
real*4 CIVa ! (0.1nm) ?=-1.0 Equivalent width of CIV (1548{AA}) (1)
real*4 e_CIVa ! (0.1nm) ?=-1.0 The 1{sigma} error in CIVa (1)
real*4 CIVb ! (0.1nm) ?=-1.0 Equivalent width of CIV (1550{AA}) (1)
real*4 e_CIVb ! (0.1nm) ?=-1.0 The 1{sigma} error in CIVb (1)
real*4 MgI ! (0.1nm) ?=-1.0 Equivalent width of MgI (2852{AA}) (1)
real*4 e_MgI ! (0.1nm) ?=-1.0 The 1{sigma} error in MgI (1)
real*4 MgIIa ! (0.1nm) ?=-1.0 Equivalent width of MgII (2796{AA}) (1)
real*4 e_MgIIa ! (0.1nm) ?=-1.0 The 1{sigma} error in MgIIa (1)
real*4 MgIIb ! (0.1nm) ?=-1.0 Equivalent width of MgII (2803{AA}) (1)
real*4 e_MgIIb ! (0.1nm) ?=-1.0 The 1{sigma} error in MgIIb (1)
real*4 SiIVa ! (0.1nm) ?=-1.0 Equivalent width of SiIV (1396{AA}) (1)
real*4 e_SiIVa ! (0.1nm) ?=-1.0 The 1{sigma} error in SiIVa (1)
real*4 SiIVb ! (0.1nm) ?=-1.0 Equivalent width of SiIV (1403{AA}) (1)
real*4 e_SiIVb ! (0.1nm) ?=-1.0 The 1{sigma} error in SiIVb (1)
real*4 FeII ! (0.1nm) ?=-1.0 Equivalent width of FeII (2382{AA}) (1)
real*4 e_FeII ! (0.1nm) ?=-1.0 The 1{sigma} error in FeII (1)
*Note (1): In the absorber rest frame and in units of Angstroms.
* Equivalent widths and error values of "-1.0" indicate that the line
* was not detected and the line was not covered by the SDSS spectrum,
* respectively.
C=============================================================================
C Loading file 'table5.dat' ! The full absorber and matching non-absorber samples
C Format for file interpretation
1 format(
+ I4,1X,I3,1X,I5,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.3,1X,F6.3,1X,I4,
+ 1X,I3,1X,I5,1X,F6.4,1X,F7.4,1X,F6.3)
C Effective file loading
open(unit=1,file='table5.dat', status='old')
write(6,*) '....Loading file: table5.dat'
do i__=1,415
read(1,'(A87)')ar__
read(ar__,1)
+ Plate,Fiber,MJD,zem_A,Dg_i_A,zab_A,imag_A,beta_A,Plate2,
+ Fiber2,MJD2,zem_N,Dg_i_N,imag_N
c ..............Just test output...........
write(6,1)
+ Plate,Fiber,MJD,zem_A,Dg_i_A,zab_A,imag_A,beta_A,Plate2,
+ Fiber2,MJD2,zem_N,Dg_i_N,imag_N
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Equivalent widths of chosen lines
C Format for file interpretation
2 format(
+ I4,1X,I3,1X,I5,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='table6.dat', status='old')
write(6,*) '....Loading file: table6.dat'
do i__=1,415
read(1,'(A122)')ar__1
read(ar__1,2)
+ Plate_1,Fiber_1,MJD_1,AlII,e_AlII,CIVa,e_CIVa,CIVb,e_CIVb,MgI,
+ e_MgI,MgIIa,e_MgIIa,MgIIb,e_MgIIb,SiIVa,e_SiIVa,SiIVb,e_SiIVb,
+ FeII,e_FeII
c ..............Just test output...........
write(6,2)
+ Plate_1,Fiber_1,MJD_1,AlII,e_AlII,CIVa,e_CIVa,CIVb,e_CIVb,MgI,
+ e_MgI,MgIIa,e_MgIIa,MgIIb,e_MgIIb,SiIVa,e_SiIVa,SiIVb,e_SiIVb,
+ FeII,e_FeII
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end