Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/673/331 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-Mar-28 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 ! Chandra ACIS-I fields surveyed in the CG 12 * region (G1) integer*4 Src ! X-ray source number within the field, * [GFL2008] "Field" "Source" in Simbad character*18 CXOU ! CXOUJ designation (HHMMSS.dd+DDMMSS.s) real*8 RAdeg ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg ! (deg) Declination in decimal degrees (J2000) real*4 e_Pos ! (arcsec) Positional uncertainty real*4 Theta ! (arcmin) Off-axis angle real*4 Fnet ! (ct) Estimated net counts from extracted area in * total energy band 0.5-8.0keV real*4 e_Fnet ! (ct) rms uncertainty on ctNetF real*4 Fbg ! (ct) Estimated background counts from extracted * area in total energy band real*4 Hnet ! (ct) Estimated net counts from extracted area in * hard energy band 2-8keV real*4 Frac ! Fraction of the PSF at the fiducial energy of * 1.497keV enclosed within the extracted * area (2) real*4 Signi ! Source significance in {sigma} character*2 Anom ! [ge ] Source anomalies (3) character*1 Var ! [abc] Variability characterization based on * K-S statistics (4) real*4 Texp ! (ks) Source effective exposure real*4 Emed ! (keV) Background-corrected median photon energy in * the total 0.5-8keV energy band real*4 logPB ! ([-]) 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 ! Chandra ACIS-I fields surveyed in the CG 12 * region (G1) integer*4 Src_1 ! X-ray source number within the field character*18 CXOU_1 ! CXOUJ designation (HHMMSS.dd+DDMMSS.s) real*4 Fnet_1 ! (ct) Estimated net counts from extracted area in * total energy band 0.5-8.0keV real*4 Emed_1 ! (keV) Background-corrected median photon energy in * the total 0.5-8keV energy band real*4 logNH ! ([cm-2]) Best-fit value or the column density (1) real*4 E_logNH ! ([cm-2]) ? Error on logNH (upper limit) real*4 e_logNH_1 ! ([cm-2]) ? Error on logNH (lower limit) real*4 kT ! (keV) Best-fit value for plasma temperature (1) real*4 E_kT ! (keV) ? Error on kT (upper limit) real*4 e_kT_1 ! (keV) ? Error on kT (lower limit) real*4 logLh ! ([10-7W]) Hard band (2-8keV) X-ray luminosity (2) real*4 e_logLh ! ([10-7W]) rms uncertainty on logLh real*4 logLhc ! ([10-7W]) Hard band (2-8keV) absorption-corrected * X-ray luminosity (2) real*4 e_logLhc ! ([10-7W]) rms uncertainty on logLhc real*4 logLt ! ([10-7W]) Total band (-0.5-8keV) X-ray luminosity (2) real*4 e_logLt ! ([10-7W]) rms uncertainty on logLt real*4 logLtc ! ([10-7W]) Total band (-0.5-8keV) absorption-corrected * X-ray luminosity (2) real*4 e_logLtc ! ([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 ! Chandra ACIS-I fields surveyed in the CG 12 * region (G1) integer*4 Src_2 ! X-ray source number within the field character*18 CXOU_2 ! CXOUJ designation (HHMMSS.dd+DDMMSS.s) character*12 NOMAD ! NOMAD identifier (NNNN-NNNNNNN) real*4 Sep1 ! (arcsec) ? Chandra-NOMAD positional offset character*5 v_MMB2004b_ ! Source identifier from Maheswar et al. * (2004, Cat. <J/MNRAS/355/1272>) real*4 Vmag ! (mag) ? V magnitude (1) real*4 e_Vmag ! (mag) ? rms uncertainty on Vmag real*4 Rmag ! (mag) ? R magnitude (1) real*4 e_Rmag ! (mag) ? rms uncertainty on Rmag real*4 pmRA ! (mas/yr) ? Best proper motion along RA * (pmRA*cosDE) (2) real*4 e_pmRA ! (mas/yr) ? rms uncertainty on pmRA real*4 pmDE ! (mas/yr) ? Best proper motion along DE (2) real*4 e_pmDE ! (mas/yr) ? rms uncertainty on pmDE character*16 v2MASS ! 2MASS designation (HHMMSSss+DDMMSSs) real*4 Sep2 ! (arcsec) ? Chandra-2MASS positional offset real*4 Jmag ! (mag) ? 2MASS J magnitude real*4 e_Jmag ! (mag) ? rms uncertainty on Jmag real*4 Hmag ! (mag) ? 2MASS H magnitude real*4 e_Hmag ! (mag) ? rms uncertainty on Hmag real*4 Ksmag ! (mag) ? 2MASS Ks magnitude real*4 e_Ksmag ! (mag) ? rms uncertainty on Ksmag character*6 Flags ! 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,Src,CXOU,RAdeg,DEdeg,e_Pos,Theta,Fnet,e_Fnet,Fbg,Hnet, + Frac,Signi,Anom,Var,Texp,Emed,logPB c ..............Just test output........... write(6,1) + Fld,Src,CXOU,RAdeg,DEdeg,e_Pos,Theta,Fnet,e_Fnet,Fbg,Hnet, + Frac,Signi,Anom,Var,Texp,Emed,logPB 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,Src_1,CXOU_1,Fnet_1,Emed_1,logNH,E_logNH,e_logNH_1,kT, + E_kT,e_kT_1,logLh,e_logLh,logLhc,e_logLhc,logLt,e_logLt, + logLtc,e_logLtc if(ar__1(43:45) .EQ. '') E_logNH = rNULL__ if(ar__1(47:49) .EQ. '') e_logNH_1 = rNULL__ if(ar__1(56:58) .EQ. '') E_kT = rNULL__ if(ar__1(60:63) .EQ. '') e_kT_1 = rNULL__ c ..............Just test output........... write(6,2) + Fld_1,Src_1,CXOU_1,Fnet_1,Emed_1,logNH,E_logNH,e_logNH_1,kT, + E_kT,e_kT_1,logLh,e_logLh,logLhc,e_logLhc,logLt,e_logLt, + logLtc,e_logLtc 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,Src_2,CXOU_2,NOMAD,Sep1,v_MMB2004b_,Vmag,e_Vmag,Rmag, + e_Rmag,pmRA,e_pmRA,pmDE,e_pmDE,v2MASS,Sep2,Jmag,e_Jmag,Hmag, + e_Hmag,Ksmag,e_Ksmag,Flags if(ar__2(40:42) .EQ. '') Sep1 = rNULL__ if(ar__2(50:54) .EQ. '') Vmag = rNULL__ if(ar__2(56:59) .EQ. '') e_Vmag = rNULL__ if(ar__2(61:65) .EQ. '') Rmag = rNULL__ if(ar__2(67:70) .EQ. '') e_Rmag = rNULL__ if(ar__2(72:76) .EQ. '') pmRA = rNULL__ if(ar__2(78:81) .EQ. '') e_pmRA = rNULL__ if(ar__2(83:87) .EQ. '') pmDE = rNULL__ if(ar__2(89:92) .EQ. '') e_pmDE = rNULL__ if(ar__2(111:113) .EQ. '') Sep2 = rNULL__ if(ar__2(115:119) .EQ. '') Jmag = rNULL__ if(ar__2(121:124) .EQ. '') e_Jmag = rNULL__ if(ar__2(126:130) .EQ. '') Hmag = rNULL__ if(ar__2(132:135) .EQ. '') e_Hmag = rNULL__ if(ar__2(137:141) .EQ. '') Ksmag = rNULL__ if(ar__2(143:146) .EQ. '') e_Ksmag = rNULL__ c ..............Just test output........... write(6,3) + Fld_2,Src_2,CXOU_2,NOMAD,Sep1,v_MMB2004b_,Vmag,e_Vmag,Rmag, + e_Rmag,pmRA,e_pmRA,pmDE,e_pmDE,v2MASS,Sep2,Jmag,e_Jmag,Hmag, + e_Hmag,Ksmag,e_Ksmag,Flags c .......End.of.Just test output........... end do close(1) C============================================================================= stop end