Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/673/331 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-18
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/673/331 X-ray observations of CG 12 (Getman+, 2008)
*================================================================================
*The stellar population and origin of the mysterious high-latitude star-forming
*cloud CG 12.
* Getman K.V., Feigelson E.D., Lawson W.A., Broos P.S., Garmire G.P.
* <Astrophys. J., 673, 331-353 (2008)>
* =2008ApJ...673..331G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Basic source properties
integer*4 nr__
parameter (nr__=128) ! Number of records
character*111 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*3 Fld (nr__) ! Chandra ACIS-I fields surveyed in the CG 12
* region (G1)
integer*4 Src (nr__) ! X-ray source number within the field,
* [GFL2008] "Field" "Source" in Simbad
character*18 CXOU (nr__) ! CXOUJ designation (HHMMSS.dd+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 e_Pos (nr__) ! (arcsec) Positional uncertainty
real*4 Theta (nr__) ! (arcmin) Off-axis angle
real*4 Fnet (nr__) ! (ct) Estimated net counts from extracted area in
* total energy band 0.5-8.0keV
real*4 e_Fnet (nr__) ! (ct) rms uncertainty on ctNetF
real*4 Fbg (nr__) ! (ct) Estimated background counts from extracted
* area in total energy band
real*4 Hnet (nr__) ! (ct) Estimated net counts from extracted area in
* hard energy band 2-8keV
real*4 Frac (nr__) ! Fraction of the PSF at the fiducial energy of
* 1.497keV enclosed within the extracted
* area (2)
real*4 Signi (nr__) ! Source significance in {sigma}
character*2 Anom (nr__) ! [ge ] Source anomalies (3)
character*1 Var (nr__) ! [abc] Variability characterization based on
* K-S statistics (4)
real*4 Texp (nr__) ! (ks) Source effective exposure
real*4 Emed (nr__) ! (keV) Background-corrected median photon energy in
* the total 0.5-8keV energy band
real*4 logPB (nr__) ! ([-]) Logarithm of K-S probability that background
* and source spectra are drawn from the same
* distribution (5)
*Note (2): Note that a reduced PSF fraction (significantly below 90%) may
* indicate that the source is in a crowded region or on the edge of the
* field.
*Note (3): Source anomalies flags 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 (4): Source variability flags as follows:
* a = no evidence for variability
* b = possibly variable
* c = definitely variable.
* blank = No test is performed for sources with fewer than 4 total full-band
* counts. No value is reported for sources in chip gaps
* or on field edges.
*Note (5): The value of logPB has been truncated at -4.0 for sources with
* significantly different background and source spectra.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! X-ray spectroscopy for photometrically selected
sources: thermal plasma fits
integer*4 nr__1
parameter (nr__1=45) ! Number of records
character*107 ar__1 ! Full-size record
character*3 Fld_1 (nr__1) ! Chandra ACIS-I fields surveyed in the CG 12
* region (G1)
integer*4 Src_1 (nr__1) ! X-ray source number within the field
character*18 CXOU_1 (nr__1) ! CXOUJ designation (HHMMSS.dd+DDMMSS.s)
real*4 Fnet_1 (nr__1) ! (ct) Estimated net counts from extracted area in
* total energy band 0.5-8.0keV
real*4 Emed_1 (nr__1) ! (keV) Background-corrected median photon energy in
* the total 0.5-8keV energy band
real*4 logNH (nr__1) ! ([cm-2]) Best-fit value or the column density (1)
real*4 E_logNH (nr__1) ! ([cm-2]) ? Error on logNH (upper limit)
real*4 e_logNH_1 (nr__1) ! ([cm-2]) ? Error on logNH (lower limit)
real*4 kT (nr__1) ! (keV) Best-fit value for plasma temperature (1)
real*4 E_kT (nr__1) ! (keV) ? Error on kT (upper limit)
real*4 e_kT_1 (nr__1) ! (keV) ? Error on kT (lower limit)
real*4 logLh (nr__1) ! ([10-7W]) Hard band (2-8keV) X-ray luminosity (2)
real*4 e_logLh (nr__1) ! ([10-7W]) rms uncertainty on logLh
real*4 logLhc (nr__1) ! ([10-7W]) Hard band (2-8keV) absorption-corrected
* X-ray luminosity (2)
real*4 e_logLhc (nr__1) ! ([10-7W]) rms uncertainty on logLhc
real*4 logLt (nr__1) ! ([10-7W]) Total band (-0.5-8keV) X-ray luminosity (2)
real*4 e_logLt (nr__1) ! ([10-7W]) rms uncertainty on logLt
real*4 logLtc (nr__1) ! ([10-7W]) Total band (-0.5-8keV) absorption-corrected
* X-ray luminosity (2)
real*4 e_logLtc (nr__1) ! ([10-7W]) rms uncertainty on logLtc
*Note (1): All fits used the "wabs(apec)" model in XSPEC and assumed
* 0.3Z_{sun}_ abundances. Uncertainties represent 90% confidence
* intervals. 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. Fitted values with logNH<20.0cm^-2^ are truncated at
* 20.0 because ACIS-I spectra are insensitive to differences in very
* low column densities. Fitted values with kT>15keV are truncated at
* 15keV because the data often cannot discriminate between very high
* temperature values.
*Note (2): X-ray luminosities are calculated assuming a distance of 550pc
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Optical and NIR counterparts
integer*4 nr__2
parameter (nr__2=128) ! Number of records
character*153 ar__2 ! Full-size record
character*3 Fld_2 (nr__2) ! Chandra ACIS-I fields surveyed in the CG 12
* region (G1)
integer*4 Src_2 (nr__2) ! X-ray source number within the field
character*18 CXOU_2 (nr__2) ! CXOUJ designation (HHMMSS.dd+DDMMSS.s)
character*12 NOMAD (nr__2) ! NOMAD identifier (NNNN-NNNNNNN)
real*4 Sep1 (nr__2) ! (arcsec) ? Chandra-NOMAD positional offset
character*5 v_MMB2004b_(nr__2) ! Source identifier from Maheswar et al.
* (2004, Cat. <J/MNRAS/355/1272>)
real*4 Vmag (nr__2) ! (mag) ? V magnitude (1)
real*4 e_Vmag (nr__2) ! (mag) ? rms uncertainty on Vmag
real*4 Rmag (nr__2) ! (mag) ? R magnitude (1)
real*4 e_Rmag (nr__2) ! (mag) ? rms uncertainty on Rmag
real*4 pmRA (nr__2) ! (mas/yr) ? Best proper motion along RA
* (pmRA*cosDE) (2)
real*4 e_pmRA (nr__2) ! (mas/yr) ? rms uncertainty on pmRA
real*4 pmDE (nr__2) ! (mas/yr) ? Best proper motion along DE (2)
real*4 e_pmDE (nr__2) ! (mas/yr) ? rms uncertainty on pmDE
character*16 v2MASS (nr__2) ! 2MASS designation (HHMMSSss+DDMMSSs)
real*4 Sep2 (nr__2) ! (arcsec) ? Chandra-2MASS positional offset
real*4 Jmag (nr__2) ! (mag) ? 2MASS J magnitude
real*4 e_Jmag (nr__2) ! (mag) ? rms uncertainty on Jmag
real*4 Hmag (nr__2) ! (mag) ? 2MASS H magnitude
real*4 e_Hmag (nr__2) ! (mag) ? rms uncertainty on Hmag
real*4 Ksmag (nr__2) ! (mag) ? 2MASS Ks magnitude
real*4 e_Ksmag (nr__2) ! (mag) ? rms uncertainty on Ksmag
character*6 Flags (nr__2) ! 2MASS photometry quality and
* confusion-contamination flag
*Note (1): When formal statistical errors are given, V and R magnitudes are
* from Maheswar et al. (2004, Cat. <J/MNRAS/355/1272>). V magnitudes for
* three B-type stars, I 57, I 58, and I 74, are from Williams et al.
* (1977MNRAS.180..709W). Otherwise, V and R magnitudes are from the
* NOMAD catalog, in which case all V magnitudes except for the Chandra
* source I 81 are from the YB6 (=USNO unpublished) catalog, and all R
* magnitudes except for the Chandra source I 57 are from the USNO-B1.0
* catalog (R magnitudes from USNO-B1.0 are second red magnitudes R_2_,
* except for the Chandra sources I 70 and II 9, for which it is a first
* red magnitude R_1_). For the Chandra source I 81 V magnitude is from
* the Tycho-2 catalog, and for the Chandra source I 57 R magnitude is
* from the UCAC2 catalog.
*Note (2): "Best" proper motions and their mean errors from catalogs
* compiled in NOMAD (Cat. I/297)
C=============================================================================
C Loading file 'table2.dat' ! Basic source properties
C Format for file interpretation
1 format(
+ A3,1X,I2,1X,A18,1X,F10.6,1X,F10.6,1X,F4.2,1X,F5.2,1X,F5.1,1X,
+ F4.1,1X,F4.1,1X,F5.1,1X,F4.2,1X,F4.1,1X,A2,1X,A1,1X,F4.1,1X,
+ F4.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,128
read(1,'(A111)')ar__
read(ar__,1)
+ Fld(i__),Src(i__),CXOU(i__),RAdeg(i__),DEdeg(i__),e_Pos(i__),
+ Theta(i__),Fnet(i__),e_Fnet(i__),Fbg(i__),Hnet(i__),Frac(i__),
+ Signi(i__),Anom(i__),Var(i__),Texp(i__),Emed(i__),logPB(i__)
c ..............Just test output...........
write(6,1)
+ Fld(i__),Src(i__),CXOU(i__),RAdeg(i__),DEdeg(i__),e_Pos(i__),
+ Theta(i__),Fnet(i__),e_Fnet(i__),Fbg(i__),Hnet(i__),Frac(i__),
+ Signi(i__),Anom(i__),Var(i__),Texp(i__),Emed(i__),logPB(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! X-ray spectroscopy for photometrically selected
* sources: thermal plasma fits
C Format for file interpretation
2 format(
+ A3,1X,I2,1X,A18,1X,F5.1,1X,F4.2,1X,F4.1,1X,F3.1,1X,F3.1,1X,
+ F5.2,F3.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,45
read(1,'(A107)')ar__1
read(ar__1,2)
+ Fld_1(i__),Src_1(i__),CXOU_1(i__),Fnet_1(i__),Emed_1(i__),
+ logNH(i__),E_logNH(i__),e_logNH_1(i__),kT(i__),E_kT(i__),
+ e_kT_1(i__),logLh(i__),e_logLh(i__),logLhc(i__),e_logLhc(i__),
+ logLt(i__),e_logLt(i__),logLtc(i__),e_logLtc(i__)
if(ar__1(43:45) .EQ. '') E_logNH(i__) = rNULL__
if(ar__1(47:49) .EQ. '') e_logNH_1(i__) = rNULL__
if(ar__1(56:58) .EQ. '') E_kT(i__) = rNULL__
if(ar__1(60:63) .EQ. '') e_kT_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Fld_1(i__),Src_1(i__),CXOU_1(i__),Fnet_1(i__),Emed_1(i__),
+ logNH(i__),E_logNH(i__),e_logNH_1(i__),kT(i__),E_kT(i__),
+ e_kT_1(i__),logLh(i__),e_logLh(i__),logLhc(i__),e_logLhc(i__),
+ logLt(i__),e_logLt(i__),logLtc(i__),e_logLtc(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Optical and NIR counterparts
C Format for file interpretation
3 format(
+ A3,1X,I2,1X,A18,1X,A12,1X,F3.1,1X,A5,1X,F5.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,A16,1X,F3.1,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,128
read(1,'(A153)')ar__2
read(ar__2,3)
+ Fld_2(i__),Src_2(i__),CXOU_2(i__),NOMAD(i__),Sep1(i__),
+ v_MMB2004b_(i__),Vmag(i__),e_Vmag(i__),Rmag(i__),e_Rmag(i__),
+ pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),v2MASS(i__),
+ Sep2(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),
+ Ksmag(i__),e_Ksmag(i__),Flags(i__)
if(ar__2(40:42) .EQ. '') Sep1(i__) = rNULL__
if(ar__2(50:54) .EQ. '') Vmag(i__) = rNULL__
if(ar__2(56:59) .EQ. '') e_Vmag(i__) = rNULL__
if(ar__2(61:65) .EQ. '') Rmag(i__) = rNULL__
if(ar__2(67:70) .EQ. '') e_Rmag(i__) = rNULL__
if(ar__2(72:76) .EQ. '') pmRA(i__) = rNULL__
if(ar__2(78:81) .EQ. '') e_pmRA(i__) = rNULL__
if(ar__2(83:87) .EQ. '') pmDE(i__) = rNULL__
if(ar__2(89:92) .EQ. '') e_pmDE(i__) = rNULL__
if(ar__2(111:113) .EQ. '') Sep2(i__) = rNULL__
if(ar__2(115:119) .EQ. '') Jmag(i__) = rNULL__
if(ar__2(121:124) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__2(126:130) .EQ. '') Hmag(i__) = rNULL__
if(ar__2(132:135) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__2(137:141) .EQ. '') Ksmag(i__) = rNULL__
if(ar__2(143:146) .EQ. '') e_Ksmag(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Fld_2(i__),Src_2(i__),CXOU_2(i__),NOMAD(i__),Sep1(i__),
+ v_MMB2004b_(i__),Vmag(i__),e_Vmag(i__),Rmag(i__),e_Rmag(i__),
+ pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__),v2MASS(i__),
+ Sep2(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),
+ Ksmag(i__),e_Ksmag(i__),Flags(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end