Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/713/330 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.5, on 2013-May-21
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/713/330 Spitzer observations of major-merger galaxies (Xu+, 2010)
*================================================================================
*Local benchmarks for the evolution of major-merger galaxies--
*Spitzer observations of a K-band selected sample.
* Xu C.K., Domingue D., Cheng Y.-W., Lu N., Huang J., Gao Y.,
* Mazzarella J.M., Cutri R., Sun W.-H., Surace J.
* <Astrophys. J., 713, 330-355 (2010)>
* =2010ApJ...713..330X
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! KPAIR galaxy sample
integer*4 nr__
parameter (nr__=54) ! Number of records
character*83 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension J2000
real*8 Dec ! (deg) Declination J2000
character*10 KPAIR ! Pair identification (JHHMM+DDMM)
character*17 v2MASX ! Galaxy identification (JHHMMSSss+DDMMSSs)
integer*4 RAh ! (h) Hour of right ascension (J2000)
integer*4 RAm ! (min) Minute of right ascension (J2000)
real*4 RAs ! (s) Second of right ascension (J2000)
character*1 DE_ ! Sign of declination (J2000)
integer*4 DEd ! (deg) Degree of declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs ! (arcsec) Arcsecond of declination (J2000)
integer*4 Vz ! (km/s) SDSS-DR3 recession velocity
real*4 Ksmag ! (mag) 2MASX Ks-band (K_20_) magnitude
real*4 logM ! ([Msun]) Logarithm of the stellar mass of the galaxy
character*1 Type ! [SE] Spiral (S) or Elliptical (E) (1)
character*1 n_Type ! [12R]? AGN-type activity (2)
real*4 F60 ! (Jy) ? IRAS 60um flux density of the total pair
real*4 logL60 ! ([Lsun]) ? Log of the 60um luminosity ({nu}L_{nu}_) of
* the total pair
*Note (1): Galaxies are classified as either "S" or "E" using the following
* scheme: first, classical "eyeball" classifications were carried out.
* Independently, two of us (C.K.X. and D.D.) inspected the SDSS optical
* image of every galaxy in the sample and assigned a type to it
* according to its morphology. In addition, we also run an automatic
* classification script. It classifies a galaxy as "E" if u-r>2.22 and
* R_50_/R_90_<0.35, where R_50_ and R_90_ are the radii containing 50%
* and 90% of the Petrosian flux (Park & Choi 2005ApJ...635L..29P).
* Otherwise, the galaxy is classified as "S" type. The final
* classification is the median of the two eyeball results and the result
* of the automatic method.
*Note (2): Galaxies containing known active galactic nuclei (AGNs; via NED)
* are flagged with 1, 2, or R to indicate type 1, type 2, or radio AGNs.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! The Spitzer observations of the 27 pairs
integer*4 nr__1
parameter (nr__1=29) ! Number of records
character*91 ar__1 ! Full-size record
real*8 RA__1 ! (deg) Right Ascension J2000
real*8 Dec_1 ! (deg) Declination J2000
character*10 KPAIR_1 ! Pair identification
character*1 m_KPAIR ! [AB]? Multiple observation flag
integer*4 RAh_1 ! (h) Hour of right ascension (J2000)
integer*4 RAm_1 ! (min) Minute of right ascension (J2000)
real*4 RAs_1 ! (s) Second of right ascension (J2000)
character*1 DE__1 ! Sign of declination (J2000)
integer*4 DEd_1 ! (deg) Degree of declination (J2000)
integer*4 DEm_1 ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs_1 ! (arcsec) Arcsecond of declination (J2000)
character*8 AOR ! IRAC AOR key
character*8 Date ! ("YY/MM/DD") IRAC observation date
character*8 AORM ! MIPS AOR key
character*8 DateM ! ("YY/MM/DD") MIPS observation date
real*4 IRAC ! (s) IRAC on-source time
real*4 v24um ! (s) MIPS 24um on-source time
real*4 v70um ! (s) MIPS 70um on-source time
real*4 v160um ! (s) MIPS 160um on-source time
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! IR emission of paired galaxies
integer*4 nr__2
parameter (nr__2=162) ! Number of records
character*137 ar__2 ! Full-size record
character*17 v2MASX_1 ! Galaxy identification (JHHMMSSss+DDMMSSs)
character*5 Ap ! (kpc) Diameter of IRAC Aperture (4kpc, 10kpc or
* Total)
real*4 F3_6 ! (mJy) Spitzer/IRAC 3.6um flux density
real*4 e_F3_6 ! (mJy) F3.6 uncertainty
real*4 F4_5 ! (mJy) Spitzer/IRAC 4.5um flux density
real*4 e_F4_5 ! (mJy) F4.5 uncertainty
real*4 F5_8 ! (mJy) Spitzer/IRAC 5.8um flux density
real*4 e_F5_8 ! (mJy) F5.8 uncertainty
real*4 F8_0 ! (mJy) Spitzer/IRAC 8.0um flux density
real*4 e_F8_0 ! (mJy) F8.0 uncertainty
character*1 l_F24 ! Limit flag on F24
real*4 F24 ! (mJy) ? Spitzer/MIPS 24um flux density
real*4 e_F24 ! (mJy) ? F24 uncertainty
character*1 l_F70 ! Limit flag on F70
real*8 F70 ! (mJy) ? Spitzer/MIPS 70um flux density
real*4 e_F70 ! (mJy) ? F70 uncertainty
character*1 l_F160 ! Limit flag on F160
real*8 F160 ! (mJy) ? Spitzer/MIPS 160um flux density
real*4 e_F160 ! (mJy) ? F160 uncertainty
character*3 Flag ! Photometry flag (1)
character*1 l_LTIR ! Limit flag on LTIR
real*4 LTIR ! ([Lsun]) ? Total IR luminosity log(L_TIR_) (2)
character*1 l_LIR ! Limit flag on LIR
real*4 LIR ! ([Lsun]) ? IR luminosity log(L_IR_) (3)
*Note (1): Photometry methods for measurements of MIPS flux densities. "A"
* stands for aperture photometry, "W" for weak source, and "P" for PRF
* (Point Response Function) fitting. The three letters correspond to the
* 24, 70, and 160um bands, respectively. For example, "AAP" means that
* flux densities in both the 24 and 70um bands were measured by aperture
* photometry while the 160um flux density was measured by PRF fitting.
*Note (2): 3-1000um IR luminosity, calculated from L_24um_, L_70um_, and
* L_160um_ using the formula of Dale & Helou (2002ApJ...576..159D).
*Note (3): Second estimate of 3-1000um IR luminosity, calculated from L_24um_
* and L_8um_. L_IR_ is an unbiased estimator of L_TIR_
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Galaxies in the control sample and their KPAIR
matches (non-AGN spirals)
integer*4 nr__3
parameter (nr__3=39) ! Number of records
character*144 ar__3 ! Full-size record
real*8 RA__2 ! (deg) Right Ascension J2000
real*8 Dec_2 ! (deg) Declination J2000
character*10 Gal ! Galaxy identification in the control sample
integer*4 RAh_2 ! (h) Hour of right ascension (J2000)
integer*4 RAm_2 ! (min) Minute of right ascension (J2000)
real*4 RAs_2 ! (s) Second of right ascension (J2000)
character*1 DE__2 ! Sign of declination (J2000)
integer*4 DEd_2 ! (deg) Degree of declination (J2000)
integer*4 DEm_2 ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs_2 ! (arcsec) Arcsecond of declination (J2000)
real*4 z ! SDSS redshift
real*4 Ksmag_1 ! (mag) 2MASX Ks-band (K_20_) magnitude
character*1 l_F8_0 ! Limit flag on F8.0
real*8 F8_0_1 ! (mJy) Spitzer/IRAC 8.0um flux density from SWIRE-DR2
real*4 e_F8_0_1 ! (mJy) ? F8.0 uncertainty
character*1 l_F24_1 ! Limit flag on F24
real*8 F24_1 ! (mJy) Spitzer/MIPS 24um flux density from SWIRE-DR2
real*4 e_F24_1 ! (mJy) ? F24 uncertainty
real*4 logM_1 ! ([Msun]) Logarithm of the mass of the KPAIR galaxy
character*1 l_logSFR ! Limit flag on logSFR
real*4 logSFR ! ([yr-1]) Specific star formation rate (log(SFR/M))
character*17 v2MASX_2 ! Matched KPAIR galaxy identification
* (JHHMMSSss+DDMMSSs)
real*4 z2 ! Redshift of the KPAIR galaxy
real*4 logM2 ! ([Msun]) Log of the mass of the KPAIR galaxy
character*1 l_logSFR2 ! Limit flag on logSFRK
real*4 logSFR2 ! ([yr-1]) Specific SFR of the KPAIR galaxy (log scale)
character*1 l_e ! Limit flag on e
real*4 e ! SFR enhancement of the KPAIR galaxy (1)
character*3 CAT ! Category of the KPAIR galaxy (2)
real*4 SEP ! Scale free separation parameter (3)
*Note (1): SFR enhancement: {epsilon}=log((SFR/M)_KPAIR-S_)-log((SFR/M)_control_)
*Note (2): Category of the KPAIR galaxy. For example, "SE1" indicates that
* the galaxy is the primary of an S+E pair, and "SS2" denotes a
* secondary in an S+S pair.
*Note (3): Scale free separation parameter: SEP=s/(r_1_+r_2_), where s is
* the pair separation, r_1_ the K-band Kron radius of the primary, and
* r_2_ the K band Kron radius of the secondary.
C=============================================================================
C Loading file 'table1.dat' ! KPAIR galaxy sample
C Format for file interpretation
1 format(
+ A10,1X,A17,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,I5,1X,
+ F5.2,1X,F5.2,1X,A1,1X,A1,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,54
read(1,'(A83)')ar__
read(ar__,1)
+ KPAIR,v2MASX,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vz,Ksmag,logM,Type,
+ n_Type,F60,logL60
if(ar__(73:77) .EQ. '') F60 = rNULL__
if(ar__(79:83) .EQ. '') logL60 = rNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAs .GE. 0) RA_=RA_+RAs/240.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DEs .GE. 0) Dec=Dec+DEs/3600.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)
+ KPAIR,v2MASX,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vz,Ksmag,logM,Type,
+ n_Type,F60,logL60
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! The Spitzer observations of the 27 pairs
C Format for file interpretation
2 format(
+ A10,A1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A8,1X,A8,
+ 1X,A8,1X,A8,1X,F5.1,1X,F5.1,1X,F5.1,1X,F4.1)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,29
read(1,'(A91)')ar__1
read(ar__1,2)
+ KPAIR_1,m_KPAIR,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,AOR,
+ Date,AORM,DateM,IRAC,v24um,v70um,v160um
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 = RAh_1
if(RAh_1 .GE. 0) RA__1=RAh_1*15.
if(RAm_1 .GE. 0) RA__1=RA__1+RAm_1/4.
if(RAs_1 .GE. 0) RA__1=RA__1+RAs_1/240.
Dec_1 = DEd_1
if(DEm_1 .GE. 0) Dec_1=Dec_1+DEm_1/60.
if(DEs_1 .GE. 0) Dec_1=Dec_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.Dec_1.GE.0) Dec_1=-Dec_1
c ..............Just test output...........
write(6,2)
+ KPAIR_1,m_KPAIR,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,AOR,
+ Date,AORM,DateM,IRAC,v24um,v70um,v160um
write(6,'(6H Pos: 2F8.4)') RA__1,Dec_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! IR emission of paired galaxies
C Format for file interpretation
3 format(
+ A17,1X,A5,1X,F6.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,
+ F6.2,1X,F5.2,1X,A1,F6.2,1X,F5.2,2X,A1,F7.2,1X,F6.2,1X,A1,F7.2,
+ 1X,F6.2,1X,A3,1X,A1,F5.2,1X,A1,F5.2)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,162
read(1,'(A137)')ar__2
read(ar__2,3)
+ v2MASX_1,Ap,F3_6,e_F3_6,F4_5,e_F4_5,F5_8,e_F5_8,F8_0,e_F8_0,
+ l_F24,F24,e_F24,l_F70,F70,e_F70,l_F160,F160,e_F160,Flag,
+ l_LTIR,LTIR,l_LIR,LIR
if(ar__2(75:80) .EQ. '') F24 = rNULL__
if(ar__2(82:86) .EQ. '') e_F24 = rNULL__
if(ar__2(90:96) .EQ. '') F70 = rNULL__
if(ar__2(98:103) .EQ. '') e_F70 = rNULL__
if(ar__2(106:112) .EQ. '') F160 = rNULL__
if(ar__2(114:119) .EQ. '') e_F160 = rNULL__
if(ar__2(126:130) .EQ. '') LTIR = rNULL__
if(ar__2(133:137) .EQ. '') LIR = rNULL__
c ..............Just test output...........
write(6,3)
+ v2MASX_1,Ap,F3_6,e_F3_6,F4_5,e_F4_5,F5_8,e_F5_8,F8_0,e_F8_0,
+ l_F24,F24,e_F24,l_F70,F70,e_F70,l_F160,F160,e_F160,Flag,
+ l_LTIR,LTIR,l_LIR,LIR
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Galaxies in the control sample and their KPAIR
* matches (non-AGN spirals)
C Format for file interpretation
4 format(
+ A10,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F6.4,1X,F5.2,
+ 1X,A1,F7.2,1X,F6.2,1X,A1,F7.2,1X,F5.2,1X,F5.2,1X,A1,F6.2,1X,
+ A17,1X,F6.4,1X,F5.2,1X,A1,F6.2,1X,A1,F5.2,1X,A3,1X,F4.2)
C Effective file loading
open(unit=1,file='table4.dat', status='old')
write(6,*) '....Loading file: table4.dat'
do i__=1,39
read(1,'(A144)')ar__3
read(ar__3,4)
+ Gal,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,z,Ksmag_1,
+ l_F8_0,F8_0_1,e_F8_0_1,l_F24_1,F24_1,e_F24_1,logM_1,l_logSFR,
+ logSFR,v2MASX_2,z2,logM2,l_logSFR2,logSFR2,l_e,e,CAT,SEP
if(ar__3(55:60) .EQ. '') e_F8_0_1 = rNULL__
if(ar__3(71:75) .EQ. '') e_F24_1 = rNULL__
c Derive coordinates RA__2 and Dec_2 from input data
c (RA__2 and Dec_2 are set to rNULL_ when unknown)
RA__2 = RAh_2
if(RAh_2 .GE. 0) RA__2=RAh_2*15.
if(RAm_2 .GE. 0) RA__2=RA__2+RAm_2/4.
if(RAs_2 .GE. 0) RA__2=RA__2+RAs_2/240.
Dec_2 = DEd_2
if(DEm_2 .GE. 0) Dec_2=Dec_2+DEm_2/60.
if(DEs_2 .GE. 0) Dec_2=Dec_2+DEs_2/3600.
if(DE__2.EQ.'-'.AND.Dec_2.GE.0) Dec_2=-Dec_2
c ..............Just test output...........
write(6,4)
+ Gal,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,z,Ksmag_1,
+ l_F8_0,F8_0_1,e_F8_0_1,l_F24_1,F24_1,e_F24_1,logM_1,l_logSFR,
+ logSFR,v2MASX_2,z2,logM2,l_logSFR2,logSFR2,l_e,e,CAT,SEP
write(6,'(6H Pos: 2F8.4)') RA__2,Dec_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end