Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/732/110 into FORTRAN code for loading all data files into arrays.
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/732/110 SDSS 2175{AA} extinction bump candidates (Jiang+, 2011)
*================================================================================
*Toward detecting the 2175 {AA} dust feature associated with strong high-redshift
*Mg II absorption lines.
* Jiang P., Ge J., Zhou H., Wang J., Wang T.
* <Astrophys. J., 732, 110 (2011)>
* =2011ApJ...732..110J
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Median composite DR7 quasar spectrum
integer*4 nr__
parameter (nr__=7751) ! Number of records
character*20 ar__ ! Full-size record
real*4 lambda (nr__) ! (0.1nm) [800/8550] Wavelength; in Angstroms
real*8 Flux (nr__) ! [1.909/122.729] Monochromatic flux; arbitrary
* units
real*4 e_Flux (nr__) ! Uncertainty in Flux
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! The list of rejected candidates
integer*4 nr__1
parameter (nr__1=220) ! Number of records
character*73 ar__1 ! Full-size record
character*19 SDSS (nr__1) ! SDSS identification (JHHMMSS.ss+DDMMSS.s)
real*4 zem (nr__1) ! [0.651/2.099] Schneider et al. 2010,
* Cat. VII/260, emission line redshift
real*4 zab (nr__1) ! [1.00/1.84] Prochter et al. 2006,
* Cat. J/ApJ/639/766, absorption line redshift
real*4 EW (nr__1) ! (0.1nm) MgII (2796) equivalent width; in Angstroms (1)
real*4 D_g_i (nr__1) ! (mag) Relative (g-i) color of background quasar
real*4 e_D_g_i (nr__1) ! (mag) Uncertainty in D(g-i)
character*1 f_gamma (nr__1) ! [mu] measurement or upper limit (2)
real*4 gamma (nr__1) ! (um-1) FWHM of Drude profile {gamma}
real*4 e_gamma (nr__1) ! (um-1) ? Uncertainty in gamma
character*1 f_Abump (nr__1) ! [mu] measurement or upper limit (2)
real*4 Abump (nr__1) ! [0.04/6.13] Bump strength = {pi} c_3_/(2{gamma})
real*4 e_Abump (nr__1) ! ? Uncertainty in Abump
character*1 Note (nr__1) ! [a-c] Additional note (3)
*Note (1): The expected typical error is of 0.3 Angstroms.
*Note (2): Flag as follows:
* m = indicates gamma/Abump is a measurement;
* u = indicates gamma/Abump is an upper limit.
*Note (3): Note as follows:
* a = at a statistical confidence level of <3{sigma};
* b = broad absorption line quasar;
* c = the quasar spectrum was trimmed a lot.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Twelve high-confidence 2175{AA} absorber candidates
integer*4 nr__2
parameter (nr__2=12) ! Number of records
character*108 ar__2 ! Full-size record
character*19 SDSS_1 (nr__2) ! SDSS identification (JHHMMSS.ss+DDMMSS.s)
real*4 zem_1 (nr__2) ! [1.18/2.04] Schneider et al. 2010,
* Cat. VII/260, emission line redshift
real*4 zab_1 (nr__2) ! [1.00/1.85] Prochter et al. 2006,
* Cat. J/ApJ/639/766, absorption line redshift
real*4 EW_1 (nr__2) ! (0.1nm) MgII (2796) equivalent width; in Angstroms (1)
real*4 D_g_i_1 (nr__2) ! (mag) Relative (g-i) color of background quasar
real*4 e_D_g_i_1 (nr__2) ! (mag) Uncertainty in D(g-i)
real*4 c1 (nr__2) ! (mag) c_1_ parameter (2)
real*4 e_c1 (nr__2) ! (mag) c1 uncertainty
real*4 c2 (nr__2) ! (mag) c_2_ parameter (2)
real*4 e_c2 (nr__2) ! (mag) c2 uncertainty
real*4 c3 (nr__2) ! (mag) c_3_ parameter (2)
real*4 e_c3 (nr__2) ! (mag) c3 uncertainty
real*4 x0 (nr__2) ! (um-1) Peak position of Drude profile (x_0_) (2)
real*4 e_x0 (nr__2) ! (um-1) x0 uncertainty
real*4 gamma_1 (nr__2) ! (um-1) FWHM of Drude profile ({gamma}) (2)
real*4 e_gamma_1 (nr__2) ! (um-1) ? Uncertainty in gamma
real*4 Chi2 (nr__2) ! {chi}^2^_{nu}_
real*4 Sig (nr__2) ! Significance ({sigma}) (3)
*Note (1): The expected typical error is of 0.3 Angstroms.
*Note (2): We use a parameterized extinction curve (Fitzpatrick & Massa
* 1990ApJS...72..163F) constituted by a linear component and a Drude
* component to describe the optical/UV extinction curve in the rest
* frame of an absorber on SDSS spectra for QAL (quasar absorption line)
* systems. The linear component is used to model the underlying
* extinction, while the Drude component is used to model the possible
* 2175{AA} extinction bump. The parameterized extinction curve is
* written as: A({lambda})=c_1_+c_2_x+c_3_D(x,x_0_,{gamma}).
* See section 2 for further explanations.
*Note (3): The simulation technique developed by Jiang et al.
* 2010ApJ...720..328J, is used to derive the statistical significance
* of the candidate 2175{AA} bumps.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Ten median-confidence 2175{AA} absorber candidates
integer*4 nr__3
parameter (nr__3=10) ! Number of records
character*108 ar__3 ! Full-size record
character*19 SDSS_2 (nr__3) ! SDSS identification (JHHMMSS.ss+DDMMSS.s)
real*4 zem_2 (nr__3) ! [1.18/2.04] Schneider et al. 2010,
* Cat. VII/260, emission line redshift
real*4 zab_2 (nr__3) ! [1.00/1.85] Prochter et al. 2006,
* Cat. J/ApJ/639/766, absorption line redshift
real*4 EW_2 (nr__3) ! (0.1nm) MgII (2796) equivalent width; in Angstroms (1)
real*4 D_g_i_1 (nr__3) ! (mag) Relative (g-i) color of background quasar
real*4 e_D_g_i_1 (nr__3) ! (mag) Uncertainty in D(g-i)
real*4 c1_1 (nr__3) ! (mag) c_1_ parameter (2)
real*4 e_c1_1 (nr__3) ! (mag) c1 uncertainty
real*4 c2_1 (nr__3) ! (mag) c_2_ parameter (2)
real*4 e_c2_1 (nr__3) ! (mag) c2 uncertainty
real*4 c3_1 (nr__3) ! (mag) c_3_ parameter (2)
real*4 e_c3_1 (nr__3) ! (mag) c3 uncertainty
real*4 x0_1 (nr__3) ! (um-1) Peak position of Drude profile (x_0_) (2)
real*4 e_x0_1 (nr__3) ! (um-1) x0 uncertainty
real*4 gamma_2 (nr__3) ! (um-1) FWHM of Drude profile ({gamma}) (2)
real*4 e_gamma_2 (nr__3) ! (um-1) ? Uncertainty in gamma
real*4 Chi2_1 (nr__3) ! {chi}^2^_{nu}_
real*4 Sig_1 (nr__3) ! Significance ({sigma}) (3)
*Note (1): The expected typical error is of 0.3 Angstroms.
*Note (2): We use a parameterized extinction curve (Fitzpatrick & Massa
* 1990ApJS...72..163F) constituted by a linear component and a Drude
* component to describe the optical/UV extinction curve in the rest
* frame of an absorber on SDSS spectra for QAL (quasar absorption line)
* systems. The linear component is used to model the underlying
* extinction, while the Drude component is used to model the possible
* 2175{AA} extinction bump. The parameterized extinction curve is
* written as: A({lambda})=c_1_+c_2_x+c_3_D(x,x_0_,{gamma}).
* See section 2 for further explanations.
*Note (3): The simulation technique developed by Jiang et al.
* 2010ApJ...720..328J, is used to derive the statistical significance
* of the candidate 2175{AA} bumps.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Seventeen low-confidence 2175{AA} absorber candidates
integer*4 nr__4
parameter (nr__4=17) ! Number of records
character*108 ar__4 ! Full-size record
character*19 SDSS_3 (nr__4) ! SDSS identification (JHHMMSS.ss+DDMMSS.s)
real*4 zem_3 (nr__4) ! [1.18/2.04] Schneider et al. 2010,
* Cat. VII/260, emission line redshift
real*4 zab_3 (nr__4) ! [1.00/1.85] Prochter et al. 2006,
* Cat. J/ApJ/639/766, absorption line redshift
real*4 EW_3 (nr__4) ! (0.1nm) MgII (2796) equivalent width; in Angstroms (1)
real*4 D_g_i_1 (nr__4) ! (mag) Relative (g-i) color of background quasar
real*4 e_D_g_i_1 (nr__4) ! (mag) Uncertainty in D(g-i)
real*4 c1_2 (nr__4) ! (mag) c_1_ parameter (2)
real*4 e_c1_2 (nr__4) ! (mag) c1 uncertainty
real*4 c2_2 (nr__4) ! (mag) c_2_ parameter (2)
real*4 e_c2_2 (nr__4) ! (mag) c2 uncertainty
real*4 c3_2 (nr__4) ! (mag) c_3_ parameter (2)
real*4 e_c3_2 (nr__4) ! (mag) c3 uncertainty
real*4 x0_2 (nr__4) ! (um-1) Peak position of Drude profile (x_0_) (2)
real*4 e_x0_2 (nr__4) ! (um-1) x0 uncertainty
real*4 gamma_3 (nr__4) ! (um-1) FWHM of Drude profile ({gamma}) (2)
real*4 e_gamma_3 (nr__4) ! (um-1) ? Uncertainty in gamma
real*4 Chi2_2 (nr__4) ! {chi}^2^_{nu}_
real*4 Sig_2 (nr__4) ! Significance ({sigma}) (3)
*Note (1): The expected typical error is of 0.3 Angstroms.
*Note (2): We use a parameterized extinction curve (Fitzpatrick & Massa
* 1990ApJS...72..163F) constituted by a linear component and a Drude
* component to describe the optical/UV extinction curve in the rest
* frame of an absorber on SDSS spectra for QAL (quasar absorption line)
* systems. The linear component is used to model the underlying
* extinction, while the Drude component is used to model the possible
* 2175{AA} extinction bump. The parameterized extinction curve is
* written as: A({lambda})=c_1_+c_2_x+c_3_D(x,x_0_,{gamma}).
* See section 2 for further explanations.
*Note (3): The simulation technique developed by Jiang et al.
* 2010ApJ...720..328J, is used to derive the statistical significance
* of the candidate 2175{AA} bumps.
C=============================================================================
C Loading file 'table1.dat' ! Median composite DR7 quasar spectrum
C Format for file interpretation
1 format(F6.1,1X,F7.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,7751
read(1,'(A20)')ar__
read(ar__,1)lambda(i__),Flux(i__),e_Flux(i__)
c ..............Just test output...........
write(6,1)lambda(i__),Flux(i__),e_Flux(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! The list of rejected candidates
C Format for file interpretation
2 format(
+ A19,1X,F5.3,1X,F5.3,1X,F4.2,1X,F5.2,1X,F4.2,1X,A1,1X,F4.2,1X,
+ F4.2,1X,A1,1X,F4.2,1X,F4.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,220
read(1,'(A73)')ar__1
read(ar__1,2)
+ SDSS(i__),zem(i__),zab(i__),EW(i__),D_g_i(i__),e_D_g_i(i__),
+ f_gamma(i__),gamma(i__),e_gamma(i__),f_Abump(i__),Abump(i__),
+ e_Abump(i__),Note(i__)
if(ar__1(56:59) .EQ. '') e_gamma(i__) = rNULL__
if(ar__1(68:71) .EQ. '') e_Abump(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ SDSS(i__),zem(i__),zab(i__),EW(i__),D_g_i(i__),e_D_g_i(i__),
+ f_gamma(i__),gamma(i__),e_gamma(i__),f_Abump(i__),Abump(i__),
+ e_Abump(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Twelve high-confidence 2175{AA} absorber candidates
C Format for file interpretation
3 format(
+ A19,1X,F5.3,1X,F5.3,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,12
read(1,'(A108)')ar__2
read(ar__2,3)
+ SDSS_1(i__),zem_1(i__),zab_1(i__),EW_1(i__),D_g_i_1(i__),
+ e_D_g_i_1(i__),c1(i__),e_c1(i__),c2(i__),e_c2(i__),c3(i__),
+ e_c3(i__),x0(i__),e_x0(i__),gamma_1(i__),e_gamma_1(i__),
+ Chi2(i__),Sig(i__)
if(ar__2(96:99) .EQ. '') e_gamma_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ SDSS_1(i__),zem_1(i__),zab_1(i__),EW_1(i__),D_g_i_1(i__),
+ e_D_g_i_1(i__),c1(i__),e_c1(i__),c2(i__),e_c2(i__),c3(i__),
+ e_c3(i__),x0(i__),e_x0(i__),gamma_1(i__),e_gamma_1(i__),
+ Chi2(i__),Sig(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Ten median-confidence 2175{AA} absorber candidates
C Format for file interpretation
4 format(
+ A19,1X,F5.3,1X,F5.3,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,10
read(1,'(A108)')ar__3
read(ar__3,4)
+ SDSS_2(i__),zem_2(i__),zab_2(i__),EW_2(i__),D_g_i_1(i__),
+ e_D_g_i_1(i__),c1_1(i__),e_c1_1(i__),c2_1(i__),e_c2_1(i__),
+ c3_1(i__),e_c3_1(i__),x0_1(i__),e_x0_1(i__),gamma_2(i__),
+ e_gamma_2(i__),Chi2_1(i__),Sig_1(i__)
if(ar__3(96:99) .EQ. '') e_gamma_2(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ SDSS_2(i__),zem_2(i__),zab_2(i__),EW_2(i__),D_g_i_1(i__),
+ e_D_g_i_1(i__),c1_1(i__),e_c1_1(i__),c2_1(i__),e_c2_1(i__),
+ c3_1(i__),e_c3_1(i__),x0_1(i__),e_x0_1(i__),gamma_2(i__),
+ e_gamma_2(i__),Chi2_1(i__),Sig_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Seventeen low-confidence 2175{AA} absorber candidates
C Format for file interpretation
5 format(
+ A19,1X,F5.3,1X,F5.3,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,17
read(1,'(A108)')ar__4
read(ar__4,5)
+ SDSS_3(i__),zem_3(i__),zab_3(i__),EW_3(i__),D_g_i_1(i__),
+ e_D_g_i_1(i__),c1_2(i__),e_c1_2(i__),c2_2(i__),e_c2_2(i__),
+ c3_2(i__),e_c3_2(i__),x0_2(i__),e_x0_2(i__),gamma_3(i__),
+ e_gamma_3(i__),Chi2_2(i__),Sig_2(i__)
if(ar__4(96:99) .EQ. '') e_gamma_3(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ SDSS_3(i__),zem_3(i__),zab_3(i__),EW_3(i__),D_g_i_1(i__),
+ e_D_g_i_1(i__),c1_2(i__),e_c1_2(i__),c2_2(i__),e_c2_2(i__),
+ c3_2(i__),e_c3_2(i__),x0_2(i__),e_x0_2(i__),gamma_3(i__),
+ e_gamma_3(i__),Chi2_2(i__),Sig_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end