Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/665/719 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 2024-Apr-19
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/665/719 X-ray imaging of RCW 49 (Tsujimoto+, 2007)
*================================================================================
*An X-ray imaging study of the stellar population in RCW 49.
* Tsujimoto M., Feigelson E.D., Townsley L.K., Broos P.S., Getman K.V.,
* Wang J., Garmire G.P., Baba D., Nagayama T., Tamura M., Churchwell E.B.
* <Astrophys. J., 665, 719-735 (2007)>
* =2007ApJ...665..719T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Chandra catalog: basic source properties
integer*4 nr__
parameter (nr__=468) ! Number of records
character*109 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 v_TFT2007_ ! Source sequence number (1)
character*18 CXOUJ ! Source IAU name (HHMMSS.ss+DDMMSS.s)
real*8 RAdeg ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000)
real*4 ePos ! (arcsec) Position error (2)
real*4 theta ! (arcmin) Off-axis angle
real*4 CNet ! (ct) Estimated (0.5-8keV) net counts
real*4 e_CNet ! (ct) Averaged upper and lower 1{sigma} error in CNet
real*4 Back ! (ct) Background (0.5-8keV) counts (3)
real*4 CHrd ! (ct) Estimated hard (2-8keV) net counts
real*4 PSF ! PSF fraction enclosed within extraction
* region (4)
real*4 Sign ! Photometric significance
character*1 l_logP ! Limit flag on logP
real*4 logP ! ([-]) Log probability extracted (0.5-8keV) counts
* are solely from background
character*1 Anom ! [geps] Source anomalies (5)
character*1 Var ! [abc] (0.5-8keV) variability characterization
* from K-S statistic (6)
real*4 Eexp ! (ks) Effective exposure time (7)
real*4 Emed ! (keV) Background-corrected (0.5-8keV) median photon
* energy
real*4 FX ! (mW/m2) Photometric (0.5-8keV) flux estimate (8)
*Note (1): Sorted by RA.
*Note (2): Estimated standard deviation of the random component of the
* position error, {sqrt}[{sigma}_x_^2^+{sigma}_y_^2^]. The single-axis
* position errors, {sigma}_x_ and {sigma}_y_, are estimated from the
* single-axis standard deviations of the PSF inside the extraction
* region and the number of counts extracted.
*Note (3): Scaled to source extraction area.
*Note (4): At 1.497keV. A reduced PSF fraction (significantly below 90%) may
* indicate that the source is in a crowded region.
*Note (5): Source anomalies flag as follows:
* g = fractional time that source was on a detector
* (FRACEXPO from mkarf) is <0.9
* e = source on field edge
* p = source piled up
* s = source on readout streak
*Note (6): Variability characterization flag as follows:
* a = no evidence for variability (0.05<P_KS_)
* b = possibly variable (0.005<P_KS_<0.05)
* c = definitely variable (P_KS_<0.005)
* No value is reported for sources in chip gaps or on field edges.
*Note (7): Approximate time the source would have to be observed on-axis to
* obtain the reported number of counts.
*Note (8): See discussion in Section 4.1.3.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! X-ray spectroscopic fits
integer*4 nr__1
parameter (nr__1=228) ! Number of records
character*108 ar__1 ! Full-size record
integer*4 v_TFT2007__1 ! Source sequence number (1)
character*18 CXOUJ_1 ! Source IAU name (HHMMSS.ss+DDMMSS.s)
real*4 CNet_1 ! (ct) Estimated (0.5-8keV) net counts
real*4 Sign_1 ! Photometric significance
real*4 e_logNH ! ([cm-2]) ? Lower 90% confidence interval in logNH (2)
real*4 logNH ! ([cm-2]) ? Log of column density from spectral fit (3)
real*4 E_logNH_1 ! ([cm-2]) ? Upper 90% confidence interval in logNH (2)
real*4 e_kT ! (keV) ? Lower 90% confidence interval in kT (2)
real*4 kT ! (keV) ? Plasma temperature from spectral fit (3)
real*4 E_kT_1 ! (keV) ? Upper 90% confidence interval in kT (2)
real*4 e_logEM ! ([cm-3]) ? Lower 90% confidence interval in logEM (2)
real*4 logEM ! ([cm-3]) ? Log of emission measure from spectral
* fit (3)
real*4 E_logEM_1 ! ([cm-3]) ? Upper 90% confidence interval in logEM (2)
real*4 logLs ! ([10-7W]) ? Log of soft (0.5-2keV) luminosity (4)
real*4 logLh ! ([10-7W]) ? Log of hard (2-8keV) luminosity (4)
real*4 logLhc ! ([10-7W]) ? Log of absorption-corrected (2-8keV)
* luminosity (4)
real*4 logLt ! ([10-7W]) ? Log of (0.5-8keV) luminosity (4)
real*4 logLtc ! ([10-7W]) ? Log of absorption-corrected (0.5-8keV)
* luminosity (4)
*Note (1): Sorted by RA.
*Note (2): More significant digits are used for uncertainties <=0.1 in order
* to avoid large rounding errors; for consistency, the same number of
* significant digits is used for both lower and upper uncertainties.
* Uncertainties are missing when XSPEC was unable to compute them or
* when their values were so large that the parameter is effectively
* unconstrained. Fits lacking uncertainties, fits with large
* uncertainties, and fits with frozen parameters should be viewed merely
* as splines to the data to obtain rough estimates of luminosities; the
* listed parameter values are unreliable.
*Note (3): Sources with photometric significance of larger than 2 were fit
* with an absorbed thin-thermal plasma model. The abundance is fixed to
* be 0.3 times solar value. The emission measure for the model spectrum
* assumes a distance of 4.2kpc.
*Note (4): Absorption-corrected luminosities are omitted when logN>22.5
* since the soft band emission is essentially unmeasurable. Luminosities
* are derived assuming a distance of 4.2kpc.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Optical, NIR and MIR counterparts of ACIS sources
integer*4 nr__2
parameter (nr__2=379) ! Number of records
character*117 ar__2 ! Full-size record
integer*4 v_TFT2007__2 ! Source sequence number (1)
character*12 NOMAD ! The NOMAD identification (Cat. I/297)
character*16 v2MASS ! The 2MASS identification (Cat. II/246)
character*17 IRAC ! The IRAC identification
character*1 NIR ! [TS] The Near-IR photometry flag (2)
real*4 Rmag ! (mag) ? R-band magnitude
real*4 Jmag ! (mag) ? J-band magnitude
real*4 Hmag ! (mag) ? H-band magnitude
real*4 Kmag ! (mag) ? Ks-band magnitude
real*4 v3_6mag ! (mag) ? IRAC 3.6 micron band magnitude
real*4 v4_5mag ! (mag) ? IRAC 4.5 micron band magnitude
real*4 v5_8mag ! (mag) ? IRAC 5.8 micron band magnitude
real*4 v8_0mag ! (mag) ? IRAC 8.0 micron band magnitude
character*10 Flags ! Flags (3)
character*12 OID ! Literature identifications (4)
*Note (1): Sorted by RA.
*Note (2): Near-IR photometry flag as follows:
* T = NIR photometry is from 2MASS
* S = NIR photometry is SIRIUS data
* NIR magnitudes are given only for good photometry; i.e., flags A, B,
* or C for 2MASS sources and the uncertainty of <0.1 mag for SIRIUS sources.
*Note (3): Three types of flags (NE, ME, and ET) with the following meaning:
* NE = NIR excess sources identified in the NIR color-color diagram (Fig. 9a)
* ME = MIR excess sources (Whitney et al., 2004ApJS..154..315W)
* ET = early-type star candidates that have a brighter K-band magnitude
* than a main-sequence B2V star in the NIR color-magnitude diagram
* (Fig 9b)
*Note (4): Also given are the names in the literature:
* MSP = Moffat et al. (1991AJ....102..642M)
* BM = Belloni & Mereghetti (1994A&A...286..935B, [BM94] in Simbad)
C=============================================================================
C Loading file 'table1.dat' ! Chandra catalog: basic source properties
C Format for file interpretation
1 format(
+ I3,1X,A18,1X,F9.5,1X,F9.5,1X,F3.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,
+ F3.1,1X,F5.1,1X,F4.2,1X,F4.1,1X,A1,F4.1,1X,A1,1X,A1,1X,F4.1,
+ 1X,F3.1,1X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,468
read(1,'(A109)')ar__
read(ar__,1)
+ v_TFT2007_,CXOUJ,RAdeg,DEdeg,ePos,theta,CNet,e_CNet,Back,CHrd,
+ PSF,Sign,l_logP,logP,Anom,Var,Eexp,Emed,FX
c ..............Just test output...........
write(6,1)
+ v_TFT2007_,CXOUJ,RAdeg,DEdeg,ePos,theta,CNet,e_CNet,Back,CHrd,
+ PSF,Sign,l_logP,logP,Anom,Var,Eexp,Emed,FX
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! X-ray spectroscopic fits
C Format for file interpretation
2 format(
+ I3,1X,A18,1X,F5.1,1X,F4.1,1X,F4.2,1X,F4.1,1X,F4.2,1X,F4.1,1X,
+ F4.1,1X,F4.1,1X,F4.2,1X,F4.1,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,228
read(1,'(A108)')ar__1
read(ar__1,2)
+ v_TFT2007__1,CXOUJ_1,CNet_1,Sign_1,e_logNH,logNH,E_logNH_1,
+ e_kT,kT,E_kT_1,e_logEM,logEM,E_logEM_1,logLs,logLh,logLhc,
+ logLt,logLtc
if(ar__1(35:38) .EQ. '') e_logNH = rNULL__
if(ar__1(40:43) .EQ. '') logNH = rNULL__
if(ar__1(45:48) .EQ. '') E_logNH_1 = rNULL__
if(ar__1(50:53) .EQ. '') e_kT = rNULL__
if(ar__1(55:58) .EQ. '') kT = rNULL__
if(ar__1(60:63) .EQ. '') E_kT_1 = rNULL__
if(ar__1(65:68) .EQ. '') e_logEM = rNULL__
if(ar__1(70:73) .EQ. '') logEM = rNULL__
if(ar__1(75:78) .EQ. '') E_logEM_1 = rNULL__
if(ar__1(80:84) .EQ. '') logLs = rNULL__
if(ar__1(86:90) .EQ. '') logLh = rNULL__
if(ar__1(92:96) .EQ. '') logLhc = rNULL__
if(ar__1(98:102) .EQ. '') logLt = rNULL__
if(ar__1(104:108) .EQ. '') logLtc = rNULL__
c ..............Just test output...........
write(6,2)
+ v_TFT2007__1,CXOUJ_1,CNet_1,Sign_1,e_logNH,logNH,E_logNH_1,
+ e_kT,kT,E_kT_1,e_logEM,logEM,E_logEM_1,logLs,logLh,logLhc,
+ logLt,logLtc
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Optical, NIR and MIR counterparts of ACIS sources
C Format for file interpretation
3 format(
+ I3,1X,A12,1X,A16,1X,A17,1X,A1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,
+ 1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,A10,1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,379
read(1,'(A117)')ar__2
read(ar__2,3)
+ v_TFT2007__2,NOMAD,v2MASS,IRAC,NIR,Rmag,Jmag,Hmag,Kmag,
+ v3_6mag,v4_5mag,v5_8mag,v8_0mag,Flags,OID
if(ar__2(55:58) .EQ. '') Rmag = rNULL__
if(ar__2(60:63) .EQ. '') Jmag = rNULL__
if(ar__2(65:68) .EQ. '') Hmag = rNULL__
if(ar__2(70:73) .EQ. '') Kmag = rNULL__
if(ar__2(75:78) .EQ. '') v3_6mag = rNULL__
if(ar__2(80:83) .EQ. '') v4_5mag = rNULL__
if(ar__2(85:88) .EQ. '') v5_8mag = rNULL__
if(ar__2(90:93) .EQ. '') v8_0mag = rNULL__
c ..............Just test output...........
write(6,3)
+ v_TFT2007__2,NOMAD,v2MASS,IRAC,NIR,Rmag,Jmag,Hmag,Kmag,
+ v3_6mag,v4_5mag,v5_8mag,v8_0mag,Flags,OID
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end