Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/675/49 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 2026-Aug-12
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/675/49 High-redshift QSOs in the SWIRE survey (Siana+, 2008)
*================================================================================
*High-redshift QSOs in the SWIRE survey and the z~3 QSO luminosity function.
* Siana B., Polletta M.D.C., Smith H.E., Lonsdale C.J., Gonzalez-Solares E.,
* Farrah D., Babbedge T.S.R., Rowan-Robinson M., Surace J., Shupe D.,
* Fang F., Franceschini A., Oliver S.
* <Astrophys. J., 675, 49-70 (2008)>
* =2008ApJ...675...49S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! The broadband optical/IR template combined with
the Telfer et al. (2002ApJ...565..773T)/
Vandenberk et al. (2001AJ....122..549V)
composite spectra
integer*4 nr__
parameter (nr__=7762) ! Number of records
character*13 ar__ ! Full-size record
real*8 lambda (nr__) ! (um) Wavelength in microns
real*4 Flux (nr__) ! Normalized flux at lambda
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Optical/IR photometry of the 100 z~3 QSO candidates
integer*4 nr__1
parameter (nr__1=102) ! Number of records
character*114 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__1) ! Sequential number (1)
character*3 n_Seq (nr__1) ! Note on the object (2)
character*19 SWIRE (nr__1) ! SWIRE name (JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg (nr__1) ! (deg) Rigth ascension (J2000)
real*8 DEdeg (nr__1) ! (deg) Declination (J2000)
real*4 zsp (nr__1) ! ?=- Spectroscopic redshift
character*1 n_zsp (nr__1) ! [s?] 's' = redshift from SDSS
real*4 Umag (nr__1) ! (mag) ?=99.00 U magnitude (Vega)
character*1 u_Umag (nr__1) ! [)] ')' when U magnitude is 5{sigma} limit
real*4 g_mag (nr__1) ! (mag) g' magnitude (Vega) (4)
real*4 r_mag (nr__1) ! (mag) r' magnitude (Vega) (4)
real*4 i_mag (nr__1) ! (mag) ?=99.00 i' magnitude (Vega) (4)
real*4 Zmag (nr__1) ! (mag) ?=99.00 Z magnitude (Vega) (4)
integer*4 F3_6 (nr__1) ! (uJy) IRAC1 (3.6um) flux (5)
integer*4 F4_5 (nr__1) ! (uJy) IRAC2 (4.5um) flux (5)
character*1 l_F5_8 (nr__1) ! Limit flag on F5.8
integer*4 F5_8 (nr__1) ! (uJy) IRAC3 (5.8um) flux (5)
character*1 l_F8_0 (nr__1) ! Limit flag on F8.0
integer*4 F8_0 (nr__1) ! (uJy) IRAC4 (8.0um) magnitude (5)
character*1 l_F24 (nr__1) ! Limit flag on F24
integer*4 F24 (nr__1) ! (uJy) MIPS24 (24um) flux
*Note (1): In table6, initial sequential numbers (1-12) supplemented by 200.
* Objects named <[SPS2008] NNN> in Simbad.
*Note (2): Note as follows:
* a = Revised photometry moved object outside of color criteria.
* Two objects with spectroscopic confirmation were moved just
* outside of our color criteria (due to revised photometry) and are
* no longer part of our primary sample.
* b = Selected in the SWIRE Lockman field
* c = Low-redshift interlopers
* d = Could not see emission lines for redshift confirmation
*Note (4): Typical uncertainties are ~0.04mag.
*Note (5): Typical uncertainties are ~10%.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Redshifts and optical/IR photometry of the seven
spectroscopically confirmed z~4 (g'-dropout)
QSOs and two SDSS QSOs in our sample
integer*4 nr__2
parameter (nr__2=12) ! Number of records
character*114 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq_1 (nr__2) ! Sequential number (1)
character*3 n_Seq_1 (nr__2) ! Note on the object (2)
character*19 SWIRE_1 (nr__2) ! SWIRE name (JHHMMSS.ss+DDMMSS.s)
real*8 RAdeg_1 (nr__2) ! (deg) Rigth ascension (J2000)
real*8 DEdeg_1 (nr__2) ! (deg) Declination (J2000)
real*4 zsp_1 (nr__2) ! ?=- Spectroscopic redshift
character*1 n_zsp_1 (nr__2) ! [s?] 's' = redshift from SDSS
real*4 Umag_1 (nr__2) ! (mag) ?=99.00 U magnitude (Vega)
character*1 u_Umag_1 (nr__2) ! [)] ')' when U magnitude is 5{sigma} limit
real*4 g_mag_1 (nr__2) ! (mag) g' magnitude (Vega) (4)
real*4 r_mag_1 (nr__2) ! (mag) r' magnitude (Vega) (4)
real*4 i_mag_1 (nr__2) ! (mag) ?=99.00 i' magnitude (Vega) (4)
real*4 Zmag_1 (nr__2) ! (mag) ?=99.00 Z magnitude (Vega) (4)
integer*4 F3_6_1 (nr__2) ! (uJy) IRAC1 (3.6um) flux (5)
integer*4 F4_5_1 (nr__2) ! (uJy) IRAC2 (4.5um) flux (5)
character*1 l_F5_8_1 (nr__2) ! Limit flag on F5.8
integer*4 F5_8_1 (nr__2) ! (uJy) IRAC3 (5.8um) flux (5)
character*1 l_F8_0_1 (nr__2) ! Limit flag on F8.0
integer*4 F8_0_1 (nr__2) ! (uJy) IRAC4 (8.0um) magnitude (5)
character*1 l_F24_1 (nr__2) ! Limit flag on F24
integer*4 F24_1 (nr__2) ! (uJy) MIPS24 (24um) flux
*Note (1): In table6, initial sequential numbers (1-12) supplemented by 200.
* Objects named <[SPS2008] NNN> in Simbad.
*Note (2): Note as follows:
* a = Revised photometry moved object outside of color criteria.
* Two objects with spectroscopic confirmation were moved just
* outside of our color criteria (due to revised photometry) and are
* no longer part of our primary sample.
* b = Selected in the SWIRE Lockman field
* c = Low-redshift interlopers
* d = Could not see emission lines for redshift confirmation
*Note (4): Typical uncertainties are ~0.04mag.
*Note (5): Typical uncertainties are ~10%.
C=============================================================================
C Loading file 'table3.dat' ! The broadband optical/IR template combined with
* the Telfer et al. (2002ApJ...565..773T)/
* Vandenberk et al. (2001AJ....122..549V)
* composite spectra
C Format for file interpretation
1 format(F7.4,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,7762
read(1,'(A13)')ar__
read(ar__,1)lambda(i__),Flux(i__)
c ..............Just test output...........
write(6,1)lambda(i__),Flux(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Optical/IR photometry of the 100 z~3 QSO candidates
C Format for file interpretation
2 format(
+ I3,A3,7X,A19,1X,F9.5,1X,F8.5,1X,F5.3,A1,1X,F5.2,A1,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,I3,1X,A1,I3,1X,A1,I4,1X,A1,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,102
read(1,'(A114)')ar__1
read(ar__1,2)
+ Seq(i__),n_Seq(i__),SWIRE(i__),RAdeg(i__),DEdeg(i__),zsp(i__),
+ n_zsp(i__),Umag(i__),u_Umag(i__),g_mag(i__),r_mag(i__),
+ i_mag(i__),Zmag(i__),F3_6(i__),F4_5(i__),l_F5_8(i__),
+ F5_8(i__),l_F8_0(i__),F8_0(i__),l_F24(i__),F24(i__)
if (idig(ar__1(53:57)).EQ.0) zsp(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq(i__),n_Seq(i__),SWIRE(i__),RAdeg(i__),DEdeg(i__),zsp(i__),
+ n_zsp(i__),Umag(i__),u_Umag(i__),g_mag(i__),r_mag(i__),
+ i_mag(i__),Zmag(i__),F3_6(i__),F4_5(i__),l_F5_8(i__),
+ F5_8(i__),l_F8_0(i__),F8_0(i__),l_F24(i__),F24(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Redshifts and optical/IR photometry of the seven
* spectroscopically confirmed z~4 (g'-dropout)
* QSOs and two SDSS QSOs in our sample
C Format for file interpretation
3 format(
+ I3,A3,7X,A19,1X,F9.5,1X,F8.5,1X,F5.3,A1,1X,F5.2,A1,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,I3,1X,A1,I3,1X,A1,I4,1X,A1,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,12
read(1,'(A114)')ar__2
read(ar__2,3)
+ Seq_1(i__),n_Seq_1(i__),SWIRE_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),zsp_1(i__),n_zsp_1(i__),Umag_1(i__),
+ u_Umag_1(i__),g_mag_1(i__),r_mag_1(i__),i_mag_1(i__),
+ Zmag_1(i__),F3_6_1(i__),F4_5_1(i__),l_F5_8_1(i__),F5_8_1(i__),
+ l_F8_0_1(i__),F8_0_1(i__),l_F24_1(i__),F24_1(i__)
if (idig(ar__2(53:57)).EQ.0) zsp_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Seq_1(i__),n_Seq_1(i__),SWIRE_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),zsp_1(i__),n_zsp_1(i__),Umag_1(i__),
+ u_Umag_1(i__),g_mag_1(i__),r_mag_1(i__),i_mag_1(i__),
+ Zmag_1(i__),F3_6_1(i__),F4_5_1(i__),l_F5_8_1(i__),F5_8_1(i__),
+ l_F8_0_1(i__),F8_0_1(i__),l_F24_1(i__),F24_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end