Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/665/719 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 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_ (nr__) ! Source sequence number (1)
character*18 CXOUJ (nr__) ! Source IAU name (HHMMSS.ss+DDMMSS.s)
real*8 RAdeg (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
real*4 ePos (nr__) ! (arcsec) Position error (2)
real*4 theta (nr__) ! (arcmin) Off-axis angle
real*4 CNet (nr__) ! (ct) Estimated (0.5-8keV) net counts
real*4 e_CNet (nr__) ! (ct) Averaged upper and lower 1{sigma} error in CNet
real*4 Back (nr__) ! (ct) Background (0.5-8keV) counts (3)
real*4 CHrd (nr__) ! (ct) Estimated hard (2-8keV) net counts
real*4 PSF (nr__) ! PSF fraction enclosed within extraction
* region (4)
real*4 Sign (nr__) ! Photometric significance
character*1 l_logP (nr__) ! Limit flag on logP
real*4 logP (nr__) ! ([-]) Log probability extracted (0.5-8keV) counts
* are solely from background
character*1 Anom (nr__) ! [geps] Source anomalies (5)
character*1 Var (nr__) ! [abc] (0.5-8keV) variability characterization
* from K-S statistic (6)
real*4 Eexp (nr__) ! (ks) Effective exposure time (7)
real*4 Emed (nr__) ! (keV) Background-corrected (0.5-8keV) median photon
* energy
real*4 FX (nr__) ! (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(nr__1) ! Source sequence number (1)
character*18 CXOUJ_1 (nr__1) ! Source IAU name (HHMMSS.ss+DDMMSS.s)
real*4 CNet_1 (nr__1) ! (ct) Estimated (0.5-8keV) net counts
real*4 Sign_1 (nr__1) ! Photometric significance
real*4 e_logNH (nr__1) ! ([cm-2]) ? Lower 90% confidence interval in logNH (2)
real*4 logNH (nr__1) ! ([cm-2]) ? Log of column density from spectral fit (3)
real*4 E_logNH_1 (nr__1) ! ([cm-2]) ? Upper 90% confidence interval in logNH (2)
real*4 e_kT (nr__1) ! (keV) ? Lower 90% confidence interval in kT (2)
real*4 kT (nr__1) ! (keV) ? Plasma temperature from spectral fit (3)
real*4 E_kT_1 (nr__1) ! (keV) ? Upper 90% confidence interval in kT (2)
real*4 e_logEM (nr__1) ! ([cm-3]) ? Lower 90% confidence interval in logEM (2)
real*4 logEM (nr__1) ! ([cm-3]) ? Log of emission measure from spectral
* fit (3)
real*4 E_logEM_1 (nr__1) ! ([cm-3]) ? Upper 90% confidence interval in logEM (2)
real*4 logLs (nr__1) ! ([10-7W]) ? Log of soft (0.5-2keV) luminosity (4)
real*4 logLh (nr__1) ! ([10-7W]) ? Log of hard (2-8keV) luminosity (4)
real*4 logLhc (nr__1) ! ([10-7W]) ? Log of absorption-corrected (2-8keV)
* luminosity (4)
real*4 logLt (nr__1) ! ([10-7W]) ? Log of (0.5-8keV) luminosity (4)
real*4 logLtc (nr__1) ! ([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(nr__2) ! Source sequence number (1)
character*12 NOMAD (nr__2) ! The NOMAD identification (Cat. I/297)
character*16 v2MASS (nr__2) ! The 2MASS identification (Cat. II/246)
character*17 IRAC (nr__2) ! The IRAC identification
character*1 NIR (nr__2) ! [TS] The Near-IR photometry flag (2)
real*4 Rmag (nr__2) ! (mag) ? R-band magnitude
real*4 Jmag (nr__2) ! (mag) ? J-band magnitude
real*4 Hmag (nr__2) ! (mag) ? H-band magnitude
real*4 Kmag (nr__2) ! (mag) ? Ks-band magnitude
real*4 v3_6mag (nr__2) ! (mag) ? IRAC 3.6 micron band magnitude
real*4 v4_5mag (nr__2) ! (mag) ? IRAC 4.5 micron band magnitude
real*4 v5_8mag (nr__2) ! (mag) ? IRAC 5.8 micron band magnitude
real*4 v8_0mag (nr__2) ! (mag) ? IRAC 8.0 micron band magnitude
character*10 Flags (nr__2) ! Flags (3)
character*12 OID (nr__2) ! 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_(i__),CXOUJ(i__),RAdeg(i__),DEdeg(i__),ePos(i__),
+ theta(i__),CNet(i__),e_CNet(i__),Back(i__),CHrd(i__),PSF(i__),
+ Sign(i__),l_logP(i__),logP(i__),Anom(i__),Var(i__),Eexp(i__),
+ Emed(i__),FX(i__)
c ..............Just test output...........
write(6,1)
+ v_TFT2007_(i__),CXOUJ(i__),RAdeg(i__),DEdeg(i__),ePos(i__),
+ theta(i__),CNet(i__),e_CNet(i__),Back(i__),CHrd(i__),PSF(i__),
+ Sign(i__),l_logP(i__),logP(i__),Anom(i__),Var(i__),Eexp(i__),
+ Emed(i__),FX(i__)
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(i__),CXOUJ_1(i__),CNet_1(i__),Sign_1(i__),
+ e_logNH(i__),logNH(i__),E_logNH_1(i__),e_kT(i__),kT(i__),
+ E_kT_1(i__),e_logEM(i__),logEM(i__),E_logEM_1(i__),logLs(i__),
+ logLh(i__),logLhc(i__),logLt(i__),logLtc(i__)
if(ar__1(35:38) .EQ. '') e_logNH(i__) = rNULL__
if(ar__1(40:43) .EQ. '') logNH(i__) = rNULL__
if(ar__1(45:48) .EQ. '') E_logNH_1(i__) = rNULL__
if(ar__1(50:53) .EQ. '') e_kT(i__) = rNULL__
if(ar__1(55:58) .EQ. '') kT(i__) = rNULL__
if(ar__1(60:63) .EQ. '') E_kT_1(i__) = rNULL__
if(ar__1(65:68) .EQ. '') e_logEM(i__) = rNULL__
if(ar__1(70:73) .EQ. '') logEM(i__) = rNULL__
if(ar__1(75:78) .EQ. '') E_logEM_1(i__) = rNULL__
if(ar__1(80:84) .EQ. '') logLs(i__) = rNULL__
if(ar__1(86:90) .EQ. '') logLh(i__) = rNULL__
if(ar__1(92:96) .EQ. '') logLhc(i__) = rNULL__
if(ar__1(98:102) .EQ. '') logLt(i__) = rNULL__
if(ar__1(104:108) .EQ. '') logLtc(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ v_TFT2007__1(i__),CXOUJ_1(i__),CNet_1(i__),Sign_1(i__),
+ e_logNH(i__),logNH(i__),E_logNH_1(i__),e_kT(i__),kT(i__),
+ E_kT_1(i__),e_logEM(i__),logEM(i__),E_logEM_1(i__),logLs(i__),
+ logLh(i__),logLhc(i__),logLt(i__),logLtc(i__)
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(i__),NOMAD(i__),v2MASS(i__),IRAC(i__),NIR(i__),
+ Rmag(i__),Jmag(i__),Hmag(i__),Kmag(i__),v3_6mag(i__),
+ v4_5mag(i__),v5_8mag(i__),v8_0mag(i__),Flags(i__),OID(i__)
if(ar__2(55:58) .EQ. '') Rmag(i__) = rNULL__
if(ar__2(60:63) .EQ. '') Jmag(i__) = rNULL__
if(ar__2(65:68) .EQ. '') Hmag(i__) = rNULL__
if(ar__2(70:73) .EQ. '') Kmag(i__) = rNULL__
if(ar__2(75:78) .EQ. '') v3_6mag(i__) = rNULL__
if(ar__2(80:83) .EQ. '') v4_5mag(i__) = rNULL__
if(ar__2(85:88) .EQ. '') v5_8mag(i__) = rNULL__
if(ar__2(90:93) .EQ. '') v8_0mag(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ v_TFT2007__2(i__),NOMAD(i__),v2MASS(i__),IRAC(i__),NIR(i__),
+ Rmag(i__),Jmag(i__),Hmag(i__),Kmag(i__),v3_6mag(i__),
+ v4_5mag(i__),v5_8mag(i__),v8_0mag(i__),Flags(i__),OID(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end