Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/199/23 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-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 (nr__) ! Cluster name
character*1 f_Name (nr__) ! [g] Flag on name (1)
character*19 GALEX (nr__) ! First GALEX IAU name(s) (mostly HHMMSS.s+DDMMSS)
character*1 f_GALEX (nr__) ! Flag on GALEX (1)
character*16 GALEX2 (nr__) ! Second GALEX IAU name(s) (HHMMSS.s+DDMMSS)
character*1 f_GALEX2 (nr__) ! Flag on GALEX2 (1)
character*18 v2MASX (nr__) ! 2MASX identifier (mostly HHMMSSss+DDMMSSs)
character*37 IRAC (nr__) ! Spitzer/IRAC AOR number(s) (1)
character*48 MIPS (nr__) ! 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 (nr__1) ! Cluster name
character*1 n_Name (nr__1) ! [l] Note on Abell 644 (1)
real*4 z (nr__1) ! [0.0031/0.89] ACCEPT database redshift
real*4 Entr (nr__1) ! (keV.cm2) ACCEPT database central entropy (2)
real*8 K100 (nr__1) ! (keV.cm2) ACCEPT database entropy profile (2)
real*4 alpha (nr__1) ! ACCEPT database best-fit power law index (2)
real*4 Scale (nr__1) ! (kpc/arcsec) Size scale (3)
real*4 IRrad (nr__1) ! (arcsec) Spitzer radius (4)
real*4 UVrad (nr__1) ! (arcsec) ? GALEX radius (5)
real*8 Offset (nr__1) ! (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 (nr__2) ! Cluster name
character*2 R (nr__2) ! [ab ] in Hicks et al. sample (6)
real*4 NUV (nr__2) ! (mag) ? GALEX/NUV band magnitude
real*4 e_NUV (nr__2) ! (mag) ? Error in NUV (7)
real*4 FUV (nr__2) ! (mag) ? GALEX/FUV band magnitude
real*4 e_FUV (nr__2) ! (mag) ? Error in FUV (7)
real*4 Jmag (nr__2) ! (mag) ? 2MASS J band magnitude
real*4 e_Jmag (nr__2) ! (mag) ? Error in Jmag (7)
real*4 Hmag (nr__2) ! (mag) ? 2MASS H band magnitude
real*4 e_Hmag (nr__2) ! (mag) ? Error in Hmag (7)
real*4 Ksmag (nr__2) ! (mag) ? 2MASS Ks band magnitude
real*4 e_Ksmag (nr__2) ! (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 (nr__3) ! Cluster name
character*1 S (nr__3) ! [e] IRAC flux remosaicked (8)
real*8 S3_6 (nr__3) ! (mJy) ? Spizter/IRAC 3.6 micron band flux density
real*8 e_S3_6 (nr__3) ! (mJy) ? Error in S3.6
real*8 S4_5 (nr__3) ! (mJy) ? Spizter/IRAC 4.5 micron band flux density
real*4 e_S4_5 (nr__3) ! (mJy) ? Error in S4.5
real*4 S5_8 (nr__3) ! (mJy) ? Spizter/IRAC 5.8 micron band flux density
real*4 e_S5_8 (nr__3) ! (mJy) ? Error in S5.8
real*4 S8_0 (nr__3) ! (mJy) ? Spizter/IRAC 8.0 micron band flux density
real*4 e_S8_0 (nr__3) ! (mJy) ? Error in S8.0
character*1 f_S8_0 (nr__3) ! [f] Flag on S8.0 (9)
real*8 S24 (nr__3) ! (mJy) ? Spitzer/MIPS 24 micron band flux density
character*1 f_S24 (nr__3) ! [bcf] Flag on S24 (9)
real*4 e_S24 (nr__3) ! (mJy) ? Error in S24
real*8 S70 (nr__3) ! (mJy) ? Spitzer/MIPS 70 micron band flux density
character*1 f_S70 (nr__3) ! [ad] Flag on S70 (9)
real*8 e_S70 (nr__3) ! (mJy) ? Error in S70
real*8 S160 (nr__3) ! (mJy) ? Spitzer/MIPS 160 micron band flux density
character*1 f_S160 (nr__3) ! [ad] Flag on S160 (9)
real*4 e_S160 (nr__3) ! (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 (nr__4) ! Cluster Name
real*8 SJ (nr__4) ! (mJy) ? 2MASS J band flux density
real*4 e_SJ (nr__4) ! (mJy) ? Error in SJ
real*8 SH (nr__4) ! (mJy) ? 2MASS H band flux density
real*4 e_SH (nr__4) ! (mJy) ? Error in SH
real*8 SKs (nr__4) ! (mJy) ? 2MASS Ks band flux density
real*4 e_SKs (nr__4) ! (mJy) ? Error in SKs
real*4 SK24 (nr__4) ! (mJy) ? The K_24_ band flux density (24um aperture
* for K band)
real*4 e_SK24 (nr__4) ! (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 (nr__5) ! Cluster name
real*4 SFRUV (nr__5) ! (Msun/yr) ? Star formation rate from UV data (10)
real*4 e_SFRUV (nr__5) ! (Msun/yr) ? Error in SFRUV (11)
real*4 SFRIR (nr__5) ! (Msun/yr) ? Star formation rate from IR data (12)
real*4 SFR70 (nr__5) ! (Msun/yr) ? Star formation rate from 70 micron data (13)
real*4 e_SFR70 (nr__5) ! (Msun/yr) ? Error in SFR70 (11)
real*4 SFR24 (nr__5) ! (Msun/yr) ? Star formation rate from 24 micron data (13)
real*4 e_SFR24 (nr__5) ! (Msun/yr) ? Error in SFR24 (11)
real*4 Mass (nr__5) ! (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(i__),f_Name(i__),GALEX(i__),f_GALEX(i__),GALEX2(i__),
+ f_GALEX2(i__),v2MASX(i__),IRAC(i__),MIPS(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),f_Name(i__),GALEX(i__),f_GALEX(i__),GALEX2(i__),
+ f_GALEX2(i__),v2MASX(i__),IRAC(i__),MIPS(i__)
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(i__),n_Name(i__),z(i__),Entr(i__),K100(i__),alpha(i__),
+ Scale(i__),IRrad(i__),UVrad(i__),Offset(i__)
if(ar__1(61:64) .EQ. '') UVrad(i__) = rNULL__
if(ar__1(66:72) .EQ. '') Offset(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),n_Name(i__),z(i__),Entr(i__),K100(i__),alpha(i__),
+ Scale(i__),IRrad(i__),UVrad(i__),Offset(i__)
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(i__),R(i__),NUV(i__),e_NUV(i__),FUV(i__),e_FUV(i__),
+ Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),Ksmag(i__),
+ e_Ksmag(i__)
if(ar__2(23:27) .EQ. '') NUV(i__) = rNULL__
if(ar__2(29:34) .EQ. '') e_NUV(i__) = rNULL__
if(ar__2(36:40) .EQ. '') FUV(i__) = rNULL__
if(ar__2(42:47) .EQ. '') e_FUV(i__) = rNULL__
if(ar__2(49:53) .EQ. '') Jmag(i__) = rNULL__
if(ar__2(55:60) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__2(62:66) .EQ. '') Hmag(i__) = rNULL__
if(ar__2(68:73) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__2(75:79) .EQ. '') Ksmag(i__) = rNULL__
if(ar__2(81:86) .EQ. '') e_Ksmag(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),R(i__),NUV(i__),e_NUV(i__),FUV(i__),e_FUV(i__),
+ Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),Ksmag(i__),
+ e_Ksmag(i__)
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(i__),S(i__),S3_6(i__),e_S3_6(i__),S4_5(i__),
+ e_S4_5(i__),S5_8(i__),e_S5_8(i__),S8_0(i__),e_S8_0(i__),
+ f_S8_0(i__),S24(i__),f_S24(i__),e_S24(i__),S70(i__),
+ f_S70(i__),e_S70(i__),S160(i__),f_S160(i__),e_S160(i__)
if(ar__3(22:28) .EQ. '') S3_6(i__) = rNULL__
if(ar__3(30:38) .EQ. '') e_S3_6(i__) = rNULL__
if(ar__3(40:46) .EQ. '') S4_5(i__) = rNULL__
if(ar__3(48:52) .EQ. '') e_S4_5(i__) = rNULL__
if(ar__3(54:59) .EQ. '') S5_8(i__) = rNULL__
if(ar__3(61:65) .EQ. '') e_S5_8(i__) = rNULL__
if(ar__3(67:72) .EQ. '') S8_0(i__) = rNULL__
if(ar__3(74:78) .EQ. '') e_S8_0(i__) = rNULL__
if(ar__3(82:88) .EQ. '') S24(i__) = rNULL__
if(ar__3(92:97) .EQ. '') e_S24(i__) = rNULL__
if(ar__3(99:105) .EQ. '') S70(i__) = rNULL__
if(ar__3(109:115) .EQ. '') e_S70(i__) = rNULL__
if(ar__3(117:123) .EQ. '') S160(i__) = rNULL__
if(ar__3(127:132) .EQ. '') e_S160(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),S(i__),S3_6(i__),e_S3_6(i__),S4_5(i__),
+ e_S4_5(i__),S5_8(i__),e_S5_8(i__),S8_0(i__),e_S8_0(i__),
+ f_S8_0(i__),S24(i__),f_S24(i__),e_S24(i__),S70(i__),
+ f_S70(i__),e_S70(i__),S160(i__),f_S160(i__),e_S160(i__)
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(i__),SJ(i__),e_SJ(i__),SH(i__),e_SH(i__),SKs(i__),
+ e_SKs(i__),SK24(i__),e_SK24(i__)
if(ar__4(20:26) .EQ. '') SJ(i__) = rNULL__
if(ar__4(28:33) .EQ. '') e_SJ(i__) = rNULL__
if(ar__4(35:41) .EQ. '') SH(i__) = rNULL__
if(ar__4(43:48) .EQ. '') e_SH(i__) = rNULL__
if(ar__4(50:56) .EQ. '') SKs(i__) = rNULL__
if(ar__4(58:63) .EQ. '') e_SKs(i__) = rNULL__
if(ar__4(65:70) .EQ. '') SK24(i__) = rNULL__
if(ar__4(72:76) .EQ. '') e_SK24(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4(i__),SJ(i__),e_SJ(i__),SH(i__),e_SH(i__),SKs(i__),
+ e_SKs(i__),SK24(i__),e_SK24(i__)
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(i__),SFRUV(i__),e_SFRUV(i__),SFRIR(i__),SFR70(i__),
+ e_SFR70(i__),SFR24(i__),e_SFR24(i__),Mass(i__)
if(ar__5(20:24) .EQ. '') SFRUV(i__) = rNULL__
if(ar__5(26:29) .EQ. '') e_SFRUV(i__) = rNULL__
if(ar__5(31:36) .EQ. '') SFRIR(i__) = rNULL__
if(ar__5(38:43) .EQ. '') SFR70(i__) = rNULL__
if(ar__5(45:49) .EQ. '') e_SFR70(i__) = rNULL__
if(ar__5(51:56) .EQ. '') SFR24(i__) = rNULL__
if(ar__5(58:62) .EQ. '') e_SFR24(i__) = rNULL__
if(ar__5(64:68) .EQ. '') Mass(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ Name_5(i__),SFRUV(i__),e_SFRUV(i__),SFRIR(i__),SFR70(i__),
+ e_SFR70(i__),SFR24(i__),e_SFR24(i__),Mass(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end