Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/199/23 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-Aug-15
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/ApJS/199/23 IR and UV star formation in ACCEPT BCGs (Hoffer+, 2012)
*================================================================================
*Infrared and ultraviolet star formation in brightest cluster galaxies in the
*ACCEPT sample.
* Hoffer A.S., Donahue M., Hicks A., Barthelemy R.S.
* <Astrophys. J. Suppl. Ser., 199, 23 (2012)>
* =2012ApJS..199...23H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Brightest Cluster Galaxy identifications
integer*4 nr__
parameter (nr__=258) ! Number of records
character*164 ar__ ! Full-size record
character*18 Name ! Cluster name
character*1 f_Name ! [g] Flag on name (1)
character*19 GALEX ! First GALEX IAU name(s) (mostly HHMMSS.s+DDMMSS)
character*1 f_GALEX ! Flag on GALEX (1)
character*16 GALEX2 ! Second GALEX IAU name(s) (HHMMSS.s+DDMMSS)
character*1 f_GALEX2 ! Flag on GALEX2 (1)
character*18 v2MASX ! 2MASX identifier (mostly HHMMSSss+DDMMSSs)
character*37 IRAC ! Spitzer/IRAC AOR number(s) (1)
character*48 MIPS ! Spitzer/MIPS AOR number(s) (1)
*Note (1): Comments for the notes following an identifier are:
* a = MIPS observation includes a 24 micron observation.
* b = MIPS observation includes a 70 micron observation.
* c = GALEX observation only includes FUV measurement.
* d = GALEX observation only includes NUV measurement.
* e = GALEX observation includes both NUV and FUV measurement.
* f = MIPS observation includes a 160 micron observation.
* g = Spitzer observations were made by Donahue as part of the
* DDT program 488.
* h = Due to an anomalous point source in the pbcd reduction, the images
* needed to be reduced from the bcd frames.
* i = GALEX observation was taken as a Guest Investigator for
* Hicks et al. (2010ApJ...719.1844H).
* j = IRAC observation only has bands 1 and 3 on target.
* k = IRAC observation only has bands 2 and 4 on target.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Physical properties
integer*4 nr__1
parameter (nr__1=243) ! Number of records
character*72 ar__1 ! Full-size record
character*18 Name_1 ! Cluster name
character*1 n_Name ! [l] Note on Abell 644 (1)
real*4 z ! [0.0031/0.89] ACCEPT database redshift
real*4 Entr ! (keV.cm2) ACCEPT database central entropy (2)
real*8 K100 ! (keV.cm2) ACCEPT database entropy profile (2)
real*4 alpha ! ACCEPT database best-fit power law index (2)
real*4 Scale ! (kpc/arcsec) Size scale (3)
real*4 IRrad ! (arcsec) Spitzer radius (4)
real*4 UVrad ! (arcsec) ? GALEX radius (5)
real*8 Offset ! (kpc) ? Centroid offset
*Note (1): The UV data for Abell 644 is taken from Hicks et al.
* (2010ApJ...719.1844H).
*Note (2): Quantities defined in Cavagnolo et al. (2009, Cat. J/ApJS/182/12)
* are from radially fit entropy profiles with a functional form
* K(r)=K_0_+K_x_(r/r_x_)^{alpha}^, where K_0_ is the central
* entropy in excess above the power law fit, K100 is the entropy
* profile normalization at 100h_70_^-1^kpc from the X-ray centroid,
* and alpha is the best-fit power law index.
*Note (3): The size scale is calculated as the angular distance size assuming
* the standard cosmology in the paper.
*Note (4): The IR radii are set at 14.3h_70_^-1^kpc in size and are used for
* 2MASS and IRAC aperture measurements except in the case where the
* aperture is below 5". In this case the aperture is set to be 5"
* to minimize large aperture corrections.
*Note (5): The UV radius is set by the aperture photometry in the GALEX database
* which most closely matches the GALEXView total flux measurement.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Fluxes matched to UV aperture
integer*4 nr__2
parameter (nr__2=243) ! Number of records
character*86 ar__2 ! Full-size record
character*18 Name_2 ! Cluster name
character*2 R ! [ab ] in Hicks et al. sample (6)
real*4 NUV ! (mag) ? GALEX/NUV band magnitude
real*4 e_NUV ! (mag) ? Error in NUV (7)
real*4 FUV ! (mag) ? GALEX/FUV band magnitude
real*4 e_FUV ! (mag) ? Error in FUV (7)
real*4 Jmag ! (mag) ? 2MASS J band magnitude
real*4 e_Jmag ! (mag) ? Error in Jmag (7)
real*4 Hmag ! (mag) ? 2MASS H band magnitude
real*4 e_Hmag ! (mag) ? Error in Hmag (7)
real*4 Ksmag ! (mag) ? 2MASS Ks band magnitude
real*4 e_Ksmag ! (mag) ? Error in Ksmag (7)
*Note (6): Flag as follows:
* a = BCGs are also in the Hicks et al. (2010ApJ...719.1844H) sample.
* b = BCG GALEX fluxes are taken from Hicks et al. (2010ApJ...719.1844H).
*Note (7): Fluxes reported with errors equal to 0 are 3{sigma} upper limits.
* For NUV upper limits, the 2MASS flux is matched with a 7" aperture
* such that it is similar in size to the GALEX PSF.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Spitzer aperture flux
integer*4 nr__3
parameter (nr__3=243) ! Number of records
character*132 ar__3 ! Full-size record
character*18 Name_3 ! Cluster name
character*1 S ! [e] IRAC flux remosaicked (8)
real*8 S3_6 ! (mJy) ? Spizter/IRAC 3.6 micron band flux density
real*8 e_S3_6 ! (mJy) ? Error in S3.6
real*8 S4_5 ! (mJy) ? Spizter/IRAC 4.5 micron band flux density
real*4 e_S4_5 ! (mJy) ? Error in S4.5
real*4 S5_8 ! (mJy) ? Spizter/IRAC 5.8 micron band flux density
real*4 e_S5_8 ! (mJy) ? Error in S5.8
real*4 S8_0 ! (mJy) ? Spizter/IRAC 8.0 micron band flux density
real*4 e_S8_0 ! (mJy) ? Error in S8.0
character*1 f_S8_0 ! [f] Flag on S8.0 (9)
real*8 S24 ! (mJy) ? Spitzer/MIPS 24 micron band flux density
character*1 f_S24 ! [bcf] Flag on S24 (9)
real*4 e_S24 ! (mJy) ? Error in S24
real*8 S70 ! (mJy) ? Spitzer/MIPS 70 micron band flux density
character*1 f_S70 ! [ad] Flag on S70 (9)
real*8 e_S70 ! (mJy) ? Error in S70
real*8 S160 ! (mJy) ? Spitzer/MIPS 160 micron band flux density
character*1 f_S160 ! [ad] Flag on S160 (9)
real*4 e_S160 ! (mJy) ? Error in S160
*Note (8):
* e = IRAC flux measurements derived from a Basic Calibrated Data (BCD) image
* that was remosaicked.
*Note (9): Flag as follows:
* a = Originally, aperture flux measurement indicated a detection. However,
* 24 micron flux measurement and visual inspection indicated
* contamination in the aperture where flux is likely from an unrelated
* source. The reported measurement is now an upper limit computed using
* the point source estimate at 16".
* b = Source not extended but has significant contamination. Point source
* measurement at 35" radius is greater than 10% error margin.
* c = Source not extended but has mild source contamination. Point source
* measurement at 35" radius is within 10% error.
* d = Flux measurement is a filtered detection.
* f = Source is extended in the MIPS 24 micron image. 24 micron flux measured
* within the aperture (see Table 3). Note an aperture radius of 35" was
* used for NGC4636 because of significant point source contamination
* outside of 35".
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! 2MASS aperture flux
integer*4 nr__4
parameter (nr__4=243) ! Number of records
character*76 ar__4 ! Full-size record
character*18 Name_4 ! Cluster Name
real*8 SJ ! (mJy) ? 2MASS J band flux density
real*4 e_SJ ! (mJy) ? Error in SJ
real*8 SH ! (mJy) ? 2MASS H band flux density
real*4 e_SH ! (mJy) ? Error in SH
real*8 SKs ! (mJy) ? 2MASS Ks band flux density
real*4 e_SKs ! (mJy) ? Error in SKs
real*4 SK24 ! (mJy) ? The K_24_ band flux density (24um aperture
* for K band)
real*4 e_SK24 ! (mJy) ? Error in SK24
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Star Formation Rates
integer*4 nr__5
parameter (nr__5=243) ! Number of records
character*68 ar__5 ! Full-size record
character*18 Name_5 ! Cluster name
real*4 SFRUV ! (Msun/yr) ? Star formation rate from UV data (10)
real*4 e_SFRUV ! (Msun/yr) ? Error in SFRUV (11)
real*4 SFRIR ! (Msun/yr) ? Star formation rate from IR data (12)
real*4 SFR70 ! (Msun/yr) ? Star formation rate from 70 micron data (13)
real*4 e_SFR70 ! (Msun/yr) ? Error in SFR70 (11)
real*4 SFR24 ! (Msun/yr) ? Star formation rate from 24 micron data (13)
real*4 e_SFR24 ! (Msun/yr) ? Error in SFR24 (11)
real*4 Mass ! (10+10Msun) ? Stellar mass
*Note (10): UV SFR are calculated for all objects, including those in high
* K_0_ systems. Upper limits are then calculated using the uncertainties
* on the NUV, K, and inert BCG color. If a star formation rate is
* consistent with a SFR of zero within 3{sigma}, a 3{sigma} upper limit
* is reported.
*Note (11): A star formation rate uncertainty of 0 identifies the quoted
* rate as a 3{sigma} upper limit.
*Note (12): The IR SFR is estimated by a fit to the Groves et al.
* (2008ApJS..176..438G) models.
*Note (13): The 24 and 70 micron SFRs are estimated from empirical
* relationships reported in Calzetti et al. (2010ApJ...714.1256C) and
* the uncertainties are calculated using the 24 and 70 micron flux
* uncertainties, respectively.
C=============================================================================
C Loading file 'table2.dat' ! Brightest Cluster Galaxy identifications
C Format for file interpretation
1 format(A18,A1,1X,A19,A1,1X,A16,A1,1X,A18,1X,A37,1X,A48)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,258
read(1,'(A164)')ar__
read(ar__,1)
+ Name,f_Name,GALEX,f_GALEX,GALEX2,f_GALEX2,v2MASX,IRAC,MIPS
c ..............Just test output...........
write(6,1)
+ Name,f_Name,GALEX,f_GALEX,GALEX2,f_GALEX2,v2MASX,IRAC,MIPS
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Physical properties
C Format for file interpretation
2 format(
+ A18,A1,1X,F6.4,1X,F6.2,1X,F8.2,1X,F4.2,1X,F4.2,1X,F6.2,1X,
+ F4.1,1X,F7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,243
read(1,'(A72)')ar__1
read(ar__1,2)
+ Name_1,n_Name,z,Entr,K100,alpha,Scale,IRrad,UVrad,Offset
if(ar__1(61:64) .EQ. '') UVrad = rNULL__
if(ar__1(66:72) .EQ. '') Offset = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,n_Name,z,Entr,K100,alpha,Scale,IRrad,UVrad,Offset
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Fluxes matched to UV aperture
C Format for file interpretation
3 format(
+ A18,1X,A2,1X,F5.2,1X,F6.4,1X,F5.2,1X,F6.4,1X,F5.2,1X,F6.4,1X,
+ F5.2,1X,F6.4,1X,F5.2,1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,243
read(1,'(A86)')ar__2
read(ar__2,3)
+ Name_2,R,NUV,e_NUV,FUV,e_FUV,Jmag,e_Jmag,Hmag,e_Hmag,Ksmag,
+ e_Ksmag
if(ar__2(23:27) .EQ. '') NUV = rNULL__
if(ar__2(29:34) .EQ. '') e_NUV = rNULL__
if(ar__2(36:40) .EQ. '') FUV = rNULL__
if(ar__2(42:47) .EQ. '') e_FUV = rNULL__
if(ar__2(49:53) .EQ. '') Jmag = rNULL__
if(ar__2(55:60) .EQ. '') e_Jmag = rNULL__
if(ar__2(62:66) .EQ. '') Hmag = rNULL__
if(ar__2(68:73) .EQ. '') e_Hmag = rNULL__
if(ar__2(75:79) .EQ. '') Ksmag = rNULL__
if(ar__2(81:86) .EQ. '') e_Ksmag = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,R,NUV,e_NUV,FUV,e_FUV,Jmag,e_Jmag,Hmag,e_Hmag,Ksmag,
+ e_Ksmag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Spitzer aperture flux
C Format for file interpretation
4 format(
+ A18,1X,A1,1X,F7.2,1X,F9.5,1X,F7.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,
+ F6.2,1X,F5.2,1X,A1,1X,F7.2,1X,A1,1X,F6.2,1X,F7.2,1X,A1,1X,
+ F7.2,1X,F7.2,1X,A1,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,243
read(1,'(A132)')ar__3
read(ar__3,4)
+ Name_3,S,S3_6,e_S3_6,S4_5,e_S4_5,S5_8,e_S5_8,S8_0,e_S8_0,
+ f_S8_0,S24,f_S24,e_S24,S70,f_S70,e_S70,S160,f_S160,e_S160
if(ar__3(22:28) .EQ. '') S3_6 = rNULL__
if(ar__3(30:38) .EQ. '') e_S3_6 = rNULL__
if(ar__3(40:46) .EQ. '') S4_5 = rNULL__
if(ar__3(48:52) .EQ. '') e_S4_5 = rNULL__
if(ar__3(54:59) .EQ. '') S5_8 = rNULL__
if(ar__3(61:65) .EQ. '') e_S5_8 = rNULL__
if(ar__3(67:72) .EQ. '') S8_0 = rNULL__
if(ar__3(74:78) .EQ. '') e_S8_0 = rNULL__
if(ar__3(82:88) .EQ. '') S24 = rNULL__
if(ar__3(92:97) .EQ. '') e_S24 = rNULL__
if(ar__3(99:105) .EQ. '') S70 = rNULL__
if(ar__3(109:115) .EQ. '') e_S70 = rNULL__
if(ar__3(117:123) .EQ. '') S160 = rNULL__
if(ar__3(127:132) .EQ. '') e_S160 = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3,S,S3_6,e_S3_6,S4_5,e_S4_5,S5_8,e_S5_8,S8_0,e_S8_0,
+ f_S8_0,S24,f_S24,e_S24,S70,f_S70,e_S70,S160,f_S160,e_S160
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! 2MASS aperture flux
C Format for file interpretation
5 format(
+ A18,1X,F7.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F6.2,
+ 1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,243
read(1,'(A76)')ar__4
read(ar__4,5)Name_4,SJ,e_SJ,SH,e_SH,SKs,e_SKs,SK24,e_SK24
if(ar__4(20:26) .EQ. '') SJ = rNULL__
if(ar__4(28:33) .EQ. '') e_SJ = rNULL__
if(ar__4(35:41) .EQ. '') SH = rNULL__
if(ar__4(43:48) .EQ. '') e_SH = rNULL__
if(ar__4(50:56) .EQ. '') SKs = rNULL__
if(ar__4(58:63) .EQ. '') e_SKs = rNULL__
if(ar__4(65:70) .EQ. '') SK24 = rNULL__
if(ar__4(72:76) .EQ. '') e_SK24 = rNULL__
c ..............Just test output...........
write(6,5)Name_4,SJ,e_SJ,SH,e_SH,SKs,e_SKs,SK24,e_SK24
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Star Formation Rates
C Format for file interpretation
6 format(
+ A18,1X,F5.2,1X,F4.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,
+ 1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,243
read(1,'(A68)')ar__5
read(ar__5,6)
+ Name_5,SFRUV,e_SFRUV,SFRIR,SFR70,e_SFR70,SFR24,e_SFR24,Mass
if(ar__5(20:24) .EQ. '') SFRUV = rNULL__
if(ar__5(26:29) .EQ. '') e_SFRUV = rNULL__
if(ar__5(31:36) .EQ. '') SFRIR = rNULL__
if(ar__5(38:43) .EQ. '') SFR70 = rNULL__
if(ar__5(45:49) .EQ. '') e_SFR70 = rNULL__
if(ar__5(51:56) .EQ. '') SFR24 = rNULL__
if(ar__5(58:62) .EQ. '') e_SFR24 = rNULL__
if(ar__5(64:68) .EQ. '') Mass = rNULL__
c ..............Just test output...........
write(6,6)
+ Name_5,SFRUV,e_SFRUV,SFRIR,SFR70,e_SFR70,SFR24,e_SFR24,Mass
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end