Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/699/1610 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.5, on 2013-Jun-20
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/699/1610 Spitzer survey of submillimeter galaxies (Hainline+, 2009)
*================================================================================
*A mid-infrared imaging survey of submillimeter-selected galaxies with the
*Spitzer Space Telescope.
* Hainline L.J., Blain A.W., Smail I., Frayer D.T., Chapman S.C.,
* Ivison R.J., Alexander D.M.
* <Astrophys. J., 699, 1610-1632 (2009)>
* =2009ApJ...699.1610H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! IRAC Flux Densities for Chapman et al.
(2005, Cat. J/ApJ/622/772) SMG Sample
integer*4 nr__
parameter (nr__=76) ! Number of records
character*115 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension J2000
real*8 Dec (nr__) ! (deg) Declination J2000
character*19 SMM (nr__) ! Chapman et al. 2005, Cat. J/ApJ/622/772 ID
* (<[CBS2005] SMM JHHMMSS.ss+DDMMSS.s> in Simbad)
character*3 n_SMM (nr__) ! Individual notes (1)
real*4 zsp (nr__) ! Spectroscopic redshift
integer*4 RAh (nr__) ! (h) IRAC hour of Right Ascension (J2000) (2)
integer*4 RAm (nr__) ! (min) IRAC minute of Right Ascension (J2000)
real*4 RAs (nr__) ! (s) IRAC second of Right Ascension (J2000)
character*1 DE_ (nr__) ! IRAC sign of the Declination (J2000) (2)
integer*4 DEd (nr__) ! (deg) IRAC degree of Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) IRAC arcminute of Declination (J2000)
real*4 DEs (nr__) ! (arcsec) IRAC arcsecond of Declination (J2000)
character*1 s (nr__) ! [s] indicates position from SCUBA name (2)
character*1 l_S3_6 (nr__) ! Limit flag on S3.6
real*4 S3_6 (nr__) ! (uJy) ? Spitzer/IRAC 3.6 micron band flux density
real*4 e_S3_6 (nr__) ! (uJy) ? Uncertainty in S3.6
character*1 f_S3_6 (nr__) ! [d] Flag on S3.6 (3)
character*1 l_S4_5 (nr__) ! Limit flag on S4.5
real*4 S4_5 (nr__) ! (uJy) ? Spitzer/IRAC 4.5 micron band flux density
real*4 e_S4_5 (nr__) ! (uJy) ? Uncertainty in S4.5
character*1 f_S4_5 (nr__) ! [d] Flag on S4.5 (3)
character*1 l_S5_8 (nr__) ! Limit flag on S5.8
real*4 S5_8 (nr__) ! (uJy) ? Spitzer/IRAC 5.8 micron band flux density
real*4 e_S5_8 (nr__) ! (uJy) ? Uncertainty in S5.8
character*1 f_S5_8 (nr__) ! [d] Flag on S5.8 (3)
character*1 l_S8_0 (nr__) ! Limit flag on S8.0
real*4 S8_0 (nr__) ! (uJy) ? Spitzer/IRAC 8.0 micron band flux density
real*4 e_S8_0 (nr__) ! (uJy) ? Uncertainty in S8.0
*Note (1): Flag as follows:
* a = This galaxy was identified by Chapman et al. (2005, Cat. J/ApJ/622/772)
* as a low-z lens unassociated with the SMG. It has been excluded from the
* analysis presented here.
* b = The redshift of this SMG has been revised in Menendez-Delmestre et al.
* (2009ApJ...699..667M).
* c = This SMG lies on the edges of the IRAC images in all fields but is
* clearly detected.
* e = The redshift of this SMG has been revised in Pope et al.
* (2008ApJ...675.1171P).
* f = The optical counterpart for this SMG is offset from the radio
* counterpart, so we have allowed a larger IRAC position offset than the
* nominal tolerance of 1.1".
* g = This SMG is detected in the 3.6 micron channel of IRAC but falls within
* a column affected by column pull-down and with resulting large
* uncertainties cause it to be formally undetected.
* h = There is an ambiguity about object SMMJ105219.15+571858.4 which should
* be SMMJ105225.79+571906.4 as in Chapman et al. (2005,
* Cat. J/ApJ/622/772). Object not corrected in CDS tables.
* i = Misprint on object name in Chapman et al. (2005, Cat. J/ApJ/622/772)
* which gives SMMJ141800.40+512820.3 and not SMMJ141800.40+522820.3.
*Note (2): Unless flagged, position of the IRAC counterpart coming from the
* bluest wavelength of detection, since for the 3.6 and 4.5um IRAC bands
* absolute pointing refinement was successful, and we expect the
* uncertainties in the positions to be smaller.
* s = position of that particular galaxy is not covered in our IRAC
* image mosaics (in that case, given position is from the SMM name).
*Note (3):
* d = This SMG has a near neighbor contaminating the 4" photometry aperture.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! MIPS 24 micron and 70 micron Results for Chapman
et al. (2005, Cat. J/ApJ/622/772) SMG Sample
integer*4 nr__1
parameter (nr__1=65) ! Number of records
character*103 ar__1 ! Full-size record
real*8 RA__1 (nr__1) ! (deg) Right Ascension J2000
real*8 Dec_1 (nr__1) ! (deg) Declination J2000
character*19 SMM_1 (nr__1) ! Chapman et al. 2005, Cat. J/ApJ/622/772 ID
character*5 n_SMM_1 (nr__1) ! Individual notes (1)
integer*4 RA24h (nr__1) ! (h) ? 24um image hour of Right Ascension (J2000)
integer*4 RA24m (nr__1) ! (min) ? 24um image minute of Right Ascension (J2000)
real*4 RA24s (nr__1) ! (s) ? 24um image second of Right Ascension (J2000)
character*1 DE24_ (nr__1) ! 24um image sign of the Declination (J2000)
integer*4 DE24d (nr__1) ! (deg) ? 24um image degree of Declination (J2000)
integer*4 DE24m (nr__1) ! (arcmin) ? 24um image arcminute of Declination (J2000)
real*4 DE24s (nr__1) ! (arcsec) ? 24um image arcsecond of Declination (J2000)
integer*4 RA70h (nr__1) ! (h) ? 70um image hour of Right Ascension (J2000)
integer*4 RA70m (nr__1) ! (min) ? 70um image minute of Right Ascension (J2000)
real*4 RA70s (nr__1) ! (s) ? 70um image second of Right Ascension (J2000)
character*1 DE70_ (nr__1) ! 70um image sign of the Declination (J2000)
integer*4 DE70d (nr__1) ! (deg) ? 70um image degree of Declination (J2000)
integer*4 DE70m (nr__1) ! (arcmin) ? 70um image arcminute of Declination (J2000)
real*4 DE70s (nr__1) ! (arcsec) ? 70um image arcsecond of Declination (J2000)
character*1 l_S24 (nr__1) ! Limit flag on S24
real*4 S24 (nr__1) ! (uJy) ? Spitzer/MIPS 24 micron band flux density
real*4 e_S24 (nr__1) ! (uJy) ? Uncertainty in S24
character*1 l_S70 (nr__1) ! Limit flag on S70
real*4 S70 (nr__1) ! (mJy) ? Spitzer/MIPS 70 micron band flux density
real*4 e_S70 (nr__1) ! (mJy) ? Uncertainty in S70
*Note (1): Flag as follows:
* a = This galaxy was identified by Chapman et al. (2005, Cat. J/ApJ/622/772)
* as a low-z lens unassociated with the SMG. It has been excluded from the
* analysis presented here.
* b = This SMG has a mid-IR spectrum from IRS published in Menendez-Delmestre
* et al. (2007ApJ...655L..65M).
* c = This SMG has a mid-IR spectrum from IRS published in Menendez-Delmestre
* et al. (2009ApJ...699..667M).
* d = This SMG has a mid-IR spectrum from IRS published in Valiante et al.
* (2007ApJ...660.1060V).
* e = This SMG is detected at 24 micron but is blended with another nearby
* source.
* f = The 24 micron measurement for this SMG has been taken from Ivison et al.
* (2007, Cat. J/MNRAS/380/199).
* g = This SMG has a mid-IR spectrum from IRS published in Pope et al.
* (2008ApJ...675.1171P).
* h = The 24 micron measurement for this SMG has been taken from Pope et al.
* (2006MNRAS.370.1185P).
* i = The optical counterpart for this SMG is offset from the radio
* counterpart so we have allowed a larger offset from the radio position
* than the nominal 2" at 24 micron.
* j = This SMG has been allowed a larger offset from the radio position at
* 24 micron due to the difference in resolution of the radio images and
* the 24 micron images.
* k = There is an ambiguity about object SMMJ105219.15+571858.4 which should
* be SMMJ105225.79+571906.4 as in Chapman et al. (2005,
* Cat. J/ApJ/622/772). Object not corrected in CDS tables.
* l = Misprint on object name in Chapman et al. (2005, Cat. J/ApJ/622/772)
* which gives SMMJ141800.40+512820.3 and not SMMJ141800.40+522820.3.
C=============================================================================
C Loading file 'table4.dat' ! IRAC Flux Densities for Chapman et al.
* (2005, Cat. J/ApJ/622/772) SMG Sample
C Format for file interpretation
1 format(
+ 3X,A19,1X,A3,1X,F5.3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,
+ F4.1,A1,1X,A1,F5.1,1X,F4.1,1X,A1,1X,A1,F5.1,1X,F4.1,1X,A1,1X,
+ A1,F6.1,1X,F5.1,1X,A1,1X,A1,F6.1,1X,F5.1)
C Effective file loading
open(unit=1,file='table4.dat', status='old')
write(6,*) '....Loading file: table4.dat'
do i__=1,76
read(1,'(A115)')ar__
read(ar__,1)
+ SMM(i__),n_SMM(i__),zsp(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),s(i__),l_S3_6(i__),
+ S3_6(i__),e_S3_6(i__),f_S3_6(i__),l_S4_5(i__),S4_5(i__),
+ e_S4_5(i__),f_S4_5(i__),l_S5_8(i__),S5_8(i__),e_S5_8(i__),
+ f_S5_8(i__),l_S8_0(i__),S8_0(i__),e_S8_0(i__)
if(ar__(60:64) .EQ. '') S3_6(i__) = rNULL__
if(ar__(66:69) .EQ. '') e_S3_6(i__) = rNULL__
if(ar__(74:78) .EQ. '') S4_5(i__) = rNULL__
if(ar__(80:83) .EQ. '') e_S4_5(i__) = rNULL__
if(ar__(88:93) .EQ. '') S5_8(i__) = rNULL__
if(ar__(95:99) .EQ. '') e_S5_8(i__) = rNULL__
if(ar__(104:109) .EQ. '') S8_0(i__) = rNULL__
if(ar__(111:115) .EQ. '') e_S8_0(i__) = rNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_(i__) = RAh(i__)
if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
Dec(i__) = DEd(i__)
if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,1)
+ SMM(i__),n_SMM(i__),zsp(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),s(i__),l_S3_6(i__),
+ S3_6(i__),e_S3_6(i__),f_S3_6(i__),l_S4_5(i__),S4_5(i__),
+ e_S4_5(i__),f_S4_5(i__),l_S5_8(i__),S5_8(i__),e_S5_8(i__),
+ f_S5_8(i__),l_S8_0(i__),S8_0(i__),e_S8_0(i__)
write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! MIPS 24 micron and 70 micron Results for Chapman
* et al. (2005, Cat. J/ApJ/622/772) SMG Sample
C Format for file interpretation
2 format(
+ 3X,A19,1X,A5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,I2,
+ 1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A1,1X,F6.1,1X,F5.1,1X,
+ A1,1X,F5.1,1X,F3.1)
C Effective file loading
open(unit=1,file='table5.dat', status='old')
write(6,*) '....Loading file: table5.dat'
do i__=1,65
read(1,'(A103)')ar__1
read(ar__1,2)
+ SMM_1(i__),n_SMM_1(i__),RA24h(i__),RA24m(i__),RA24s(i__),
+ DE24_(i__),DE24d(i__),DE24m(i__),DE24s(i__),RA70h(i__),
+ RA70m(i__),RA70s(i__),DE70_(i__),DE70d(i__),DE70m(i__),
+ DE70s(i__),l_S24(i__),S24(i__),e_S24(i__),l_S70(i__),S70(i__),
+ e_S70(i__)
if(ar__1(30:31) .EQ. '') RA24h(i__) = iNULL__
if(ar__1(33:34) .EQ. '') RA24m(i__) = iNULL__
if(ar__1(36:40) .EQ. '') RA24s(i__) = rNULL__
if(ar__1(43:44) .EQ. '') DE24d(i__) = iNULL__
if(ar__1(46:47) .EQ. '') DE24m(i__) = iNULL__
if(ar__1(49:52) .EQ. '') DE24s(i__) = rNULL__
if(ar__1(54:55) .EQ. '') RA70h(i__) = iNULL__
if(ar__1(57:58) .EQ. '') RA70m(i__) = iNULL__
if(ar__1(60:64) .EQ. '') RA70s(i__) = rNULL__
if(ar__1(67:68) .EQ. '') DE70d(i__) = iNULL__
if(ar__1(70:71) .EQ. '') DE70m(i__) = iNULL__
if(ar__1(73:76) .EQ. '') DE70s(i__) = rNULL__
if(ar__1(80:85) .EQ. '') S24(i__) = rNULL__
if(ar__1(87:91) .EQ. '') e_S24(i__) = rNULL__
if(ar__1(95:99) .EQ. '') S70(i__) = rNULL__
if(ar__1(101:103) .EQ. '') e_S70(i__) = rNULL__
c Derive coordinates RA__1 and Dec_1 from input data
c (RA__1 and Dec_1 are set to rNULL_ when unknown)
RA__1(i__) = RA24h(i__)
if(RA24h(i__) .GE. 0) RA__1(i__)=RA24h(i__)*15.
if(RA24m(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RA24m(i__)/4.
if(RA24s(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RA24s(i__)/240.
Dec_1(i__) = DE24d(i__)
if(DE24m(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DE24m(i__)/60.
if(DE24s(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DE24s(i__)/3600.
if(DE24_(i__).EQ.'-'.AND.Dec_1(i__).GE.0) Dec_1(i__)=-Dec_1(i__)
c ..............Just test output...........
write(6,2)
+ SMM_1(i__),n_SMM_1(i__),RA24h(i__),RA24m(i__),RA24s(i__),
+ DE24_(i__),DE24d(i__),DE24m(i__),DE24s(i__),RA70h(i__),
+ RA70m(i__),RA70s(i__),DE70_(i__),DE70d(i__),DE70m(i__),
+ DE70s(i__),l_S24(i__),S24(i__),e_S24(i__),l_S70(i__),S70(i__),
+ e_S70(i__)
write(6,'(6H Pos: 2F8.4)') RA__1(i__),Dec_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end