Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/209/32 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 2026-Oct-10
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/ApJS/209/32 Probable young stars in the MYStIX project (Broos+, 2013)
*================================================================================
*Identifying young stars in massive star-forming regions for the MYStIX project.
* Broos P.S., Getman K.V., Povich M.S., Feigelson E.D., Townsley L.K.,
* Naylor T., Kuhn M.A., King R.R., Busk H.A.
* <Astrophys. J. Suppl. Ser., 209, 32 (2013)>
* =2013ApJS..209...32B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'reg.dat' ! List of the 20 nearby Galactic massive
star-forming regions; table added by CDS
integer*4 nr__
parameter (nr__=20) ! Number of records
character*22 ar__ ! Full-size record
character*16 Name (nr__) ! MYStIX Galactic massive star-forming region
* (MSFR) name
character*5 Reg (nr__) ! Abbreviated region name
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mpcm.dat' ! MPCM (MYStIX Probable Complex Members) sources
and properties (table 2)
integer*4 nr__1
parameter (nr__1=31549) ! Number of records
character*971 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 Reg_1 (nr__1) ! MYStIX Galactic massive star-forming region
* (MSFR) name
character*18 MPCM (nr__1) ! IAU source name (HHMMSS.ss+DDMMSS.s; J2000)
character*1 m_MPCM (nr__1) ! [ab] Multiple MPCM name component (1)
real*8 RAdeg (nr__1) ! (deg) Right Ascension (J2000)
real*8 DEdeg (nr__1) ! (deg) Declination (J2000)
real*4 ePos (nr__1) ! (arcsec) [0/2.1]?=- 1{sigma} error around (RA,DE)
character*11 rPos (nr__1) ! Origin of position (G1)
character*18 CXO (nr__1) ! X-ray source name in IAU format
* (HHMMSS.s+DDMMSS; J2000)
character*11 Xlab (nr__1) ! X-ray source name used within the MYStIX
* project (G2)
character*18 NName (nr__1) ! Name in NIR catalog (G3)
character*16 Nlab (nr__1) ! Label in NIR catalog
character*18 MName (nr__1) ! Name in MIR catalog
* (HHMMSS.ss+DDMMSS.s; J2000 or NNNN) (G4)
character*29 Mlab (nr__1) ! Label in MIR catalog
character*30 Olab (nr__1) ! Label in OB catalog (2)
integer*4 Xind (nr__1) ! ?=-99 0-based index in X-ray catalog (Sec. 4)
integer*4 Iind (nr__1) ! ?=-99 0-based index in IR excess catalog
* (Povich et al. 2013, Cat. J/ApJS/209/31)
character*30 SpT (nr__1) ! MK spectral type
character*14 r_SpT (nr__1) ! Reference for spectral type (3)
real*4 OBmag (nr__1) ! (mag) [4.4/24.5]? Visual photometry in OB Catalog
character*3 OBband (nr__1) ! Visual band for OBmag (J, K or V)
real*4 pr_0 (nr__1) ! [0/1]? p-value for no-source hypothesis (4)(5)
real*4 pr_1 (nr__1) ! [0/1]? p-value for single observation (4)(6)
real*4 pr_m (nr__1) ! [0/1]? p-value for merged observations (4)(6)
real*8 t_exp (nr__1) ! (s) ? Total exposure time in merged
* observations (4)
real*4 f_exp (nr__1) ! [0/1]? Fraction of "Exp" that source was
* observed (4)(7)
integer*4 Nobs (nr__1) ! [1/9]?=-99 Total number of observations
* extracted (4)
integer*4 Nph (nr__1) ! [1/9]?=-99 Number of observations merged to
* estimate photometry properties (4)
real*4 b_theta (nr__1) ! (arcmin) ? Smallest off-axis angle for merged
* observations (4)
real*4 theta (nr__1) ! (arcmin) ? Average off-axis angle {theta} for merged
* observations (4)
real*4 B_theta_1 (nr__1) ! (arcmin) ? Largest off-axis angle for merged obs. (4)
real*4 PSFf (nr__1) ! [0/1]? Average PSF fraction (at 1.5 keV) for
* merged observations (4)
real*4 aglow (nr__1) ! [0/1]? Suspected afterglow fraction (4)(8)
integer*4 Oct_t (nr__1) ! (ct) ?=-99 Extracted total (0.5-8keV) counts in
* merged apertures (4)
real*8 Nct_t (nr__1) ! (ct) [1/89438]? Net total (0.5-8keV) counts in
* merged apertures (4)
real*8 Nct_s (nr__1) ! (ct) ? Net soft (0.5-2 keV) counts in merged
* apertures (4)
real*8 Nct_h (nr__1) ! (ct) ? Net hard (2-8 keV) counts in merged
* apertures (4)
real*8 b_Nct_t (nr__1) ! (ct) ? 1{sigma} lower bound on Nct-t (4)(9)
real*8 B_Nct_t_1 (nr__1) ! (ct) ? 1{sigma} upper bound on Nct-t (4)(9)
real*8 b_Nct_s (nr__1) ! (ct) ? 1{sigma} lower bound on Nct-s (4)(9)
real*8 B_Nct_s_1 (nr__1) ! (ct) ? 1{sigma} upper bound on Nct-s (4)(9)
real*8 b_Nct_h (nr__1) ! (ct) ? 1{sigma} lower bound on Nct-h (4)(9)
real*8 B_Nct_h_1 (nr__1) ! (ct) ? 1{sigma} upper bound on Nct-h (4)(9)
real*4 E_t (nr__1) ! (keV) ? Median 0.5-8keV energy, observed
* spectrum (4)(10)
real*8 XFlx_t (nr__1) ! ([ph/cm2/s]) ? Log of incident photon flux in total
* 0.5-8 keV band (4)(11)
real*8 XFlx_s (nr__1) ! ([ph/cm2/s]) ? Log of incident photon flux in soft
* 0.5-2 keV band (4)(11)
real*8 XFlx_h (nr__1) ! ([ph/cm2/s]) ? Log of incident photon flux in hard
* 2-8 keV band (4)(11)
real*4 LX_h (nr__1) ! (10-7W) ? X-ray luminosity in erg/s, 2-8 keV (12)
real*4 LXc_h (nr__1) ! (10-7W) ? Absorption-corrected X-ray luminosity,
* in erg/s; 2-8 keV (12)
real*4 e_LXc_h (nr__1) ! (10-7W) ? LXhc 1{sigma} statistical uncertainty (12)
real*4 s_LXc_h (nr__1) ! (10-7W) ? LXhc 1{sigma} systematic uncertainty (12)
real*4 LX_t (nr__1) ! (10-7W) ? X-ray luminosity in erg/s, 0.5-8 keV (12)
real*4 LXc_t (nr__1) ! (10-7W) ? Absorption-corrected X-ray luminosity,
* in erg/s; 0.5-8 kev (12)
real*4 e_LXc_t (nr__1) ! (10-7W) ? LXtc 1{sigma} statistical uncertainty (12)
real*4 s_LXc_t (nr__1) ! (10-7W) ? LXtc 1{sigma} systematic uncertainty (12)
real*4 logNH (nr__1) ! (cm-2) ? Gas column density (12)
real*4 e_logNH (nr__1) ! (cm-2) ? logNH 1{sigma} statistical error (12)
real*4 s_logNH (nr__1) ! (cm-2) ? logNH 1{sigma} systematic uncertainty (12)
real*4 pr_N (nr__1) ! [0/1]? Counterpart probability, X-ray/NIR (G5)
real*4 pr_M_1 (nr__1) ! [0/1]? counterpart probability, X-ray/MIR (G5)
real*4 Jmag (nr__1) ! (mag) [4.7/25.2]? J-band magnitude
real*8 e_Jmag (nr__1) ! (mag) [-25.1/3.3]? Jmag 1{sigma} uncertainty
real*4 Hmag (nr__1) ! (mag) [4.5/29.9]? H-band magnitude
real*8 e_Hmag (nr__1) ! (mag) [-22.2/6.6]? Hmag 1{sigma} uncertainty
real*4 Kmag (nr__1) ! (mag) [4.3/21.5]? K-band magnitude
real*4 e_Kmag (nr__1) ! (mag) ? Kmag 1-{sigma} uncertainty
real*4 v_3_6_ (nr__1) ! (mag) ? Spitzer/IRAC 3.6um band magnitude
real*4 e__3_6_ (nr__1) ! (mag) ? [3.6] 1-{sigma} uncertainty
real*4 v_4_5_ (nr__1) ! (mag) ? Spitzer/IRAC 4.5um band magnitude
real*4 e__4_5_ (nr__1) ! (mag) ? [4.5] 1-{sigma} uncertainty
real*4 v_5_8_ (nr__1) ! (mag) [2/17.5]? Spitzer/IRAC 5.8um band magnitude
real*4 e__5_8_ (nr__1) ! (mag) [0/10]? [5.8] 1-{sigma} uncertainty
real*4 v_8_0_ (nr__1) ! (mag) [1.5/17]? Spitzer/IRAC 8.0um band magnitude
real*4 e__8_0_ (nr__1) ! (mag) [0/15]? [8.0] 1-{sigma} uncertainty
character*3 f_Jmag (nr__1) ! UKIRT photometry flag on Jmag (G6)
character*3 f_Hmag (nr__1) ! UKIRT photometry flag on Hmag (G6)
character*3 f_Kmag (nr__1) ! UKIRT photometry flag on Kmag (G6)
character*3 v2M (nr__1) ! 2MASS photometry flag
character*3 q_2M (nr__1) ! 2MASS quality photometry flag
integer*4 n_Jmag (nr__1) ! GLIMPSE photometry flag on Jmag
integer*4 n_Hmag (nr__1) ! GLIMPSE photometry flag on Hmag
integer*4 n_Kmag (nr__1) ! GLIMPSE photometry flag on Kmag
integer*4 f__3_6_ (nr__1) ! GLIMPSE photometry flag on [3.6]
integer*4 f__4_5_ (nr__1) ! GLIMPSE photometry flag on [4.5]
integer*4 f__5_8_ (nr__1) ! GLIMPSE photometry flag on [5.8]
integer*4 f__8_0_ (nr__1) ! GLIMPSE photometry flag on [8.0]
integer*4 LS (nr__1) ! [-2/4]?=-99 Local Spitzer photometry flag (G7)
character*19 r_Jmag (nr__1) ! Origin of Jmag photometry
character*19 r_Hmag (nr__1) ! Origin of Hmag photometry
character*19 r_Kmag (nr__1) ! Origin of Kmag photometry
character*19 r__3_6_ (nr__1) ! Origin of [3.6] photometry
character*19 r__4_5_ (nr__1) ! Origin of [4.5] photometry
character*19 r__5_8_ (nr__1) ! Origin of [5.8] photometry
character*19 r__8_0_ (nr__1) ! Origin of [8.0] photometry
integer*4 Cl (nr__1) ! [-2/2]?=-99 Classification (G8)
real*4 Av (nr__1) ! (mag) [0/81]? Visual extinction A_V_
integer*4 St (nr__1) ! [-1/2]?=-99 YSO stage (G9)
real*4 H1pr (nr__1) ! [0/1]? H1 class prior probability
* (position-dependent)
real*4 H2pr (nr__1) ! [0/1]? H2 class prior probability
* (position-dependent)
real*4 H3pr (nr__1) ! [0/1]? H3 class prior probability
* (position-dependent)
real*4 H4pr (nr__1) ! [0/1]? H4 class prior probability
* (position-dependent)
real*4 H1ps (nr__1) ! [0/1]? H1 class posterior probability
real*4 H2ps (nr__1) ! [0/1]? H2 class posterior probability
real*4 H3ps (nr__1) ! [0/1]? H3 class posterior probability
real*4 H4ps (nr__1) ! [0/1]? H4 class posterior probability
integer*4 H2f (nr__1) ! [1/7]?=-99 H2 dominant classification
* term (G10)
integer*4 Xcl (nr__1) ! [1/2]?=-99 X-ray classification (G11)
*Note (1): The MYStIX source name (MPCM) is the sexagesimal representation
* of the MYStIX source position. In rare cases two MPCM rows will share
* the same source position, because two X-ray sources have been assigned
* the same NIR counterpart by our counterpart identification algorithm.
* In such cases we append "a" or "b" to the source names to make them
* unique. See section 2.
*Note (2): When a source has been associated with a known massive star
* (column "Olab") all XPHOT quantities (see Note (12)) are unreliable.
*Note (3): References for SpT (number 0-10 are cited in Gagne et al., 2011,
* Cat. J/ApJS/194/5) - see refs.dat file.
*Note (4): These X-ray column labels were previously published by Broos et
* al. (2011, Cat. J/ApJS/194/4; 2011, Cat. J/ApJS/194/2) and are
* produced by the ACIS Extract (AE) software package (Broos et al.
* 2010ApJ...714.1582B, 2012ascl.soft03001B). The AE software and User's
* Guide are available at
* http://www.astro.psu.edu/xray/acis/acis_analysis.html
*Note (5): In statistical hypothesis testing, the p-value is the probability
* of obtaining a test statistic at least as extreme as the one that was
* actually observed when the null hypothesis is true.
*Note (6): Smallest p-value for the one-sample Kolmogorov-Smirnov statistic
* under the no-variability null hypothesis within a single-observation
* or over merged observations. See Broos et al. (2010ApJ...714.1582B,
* Section 7.6) for a description of the variability metrics and caveats
* regarding possible spurious indications of variability using the
* ProbKS_merge metric.
*Note (7): Due to dithering over inactive portions of the focal plane, a
* Chandra source is often not observed during some fraction of the
* nominal exposure time (See http://cxc.harvard.edu/ciao/why/dither.html).
* The reported quantity is FRACEXPO produced by the CIAO tool mkarf.
*Note (8): Some background events arising from an effect known as
* "afterglow" (http://cxc.harvard.edu/ciao/why/afterglow.html) may
* contaminate source extractions, despite careful procedures to identify
* and remove them during data preparation (Broos et al.
* 2010ApJ...714.1582B, Section 3). After extraction, we attempt to
* identify afterglow events using the tool ae_afterglow_report, and
* report the fraction of extracted events attributed to afterglow; see
* the ACIS Extract manual
* (http://www.astro.psu.edu/xray/acis/acis_analysis.html).
*Note (9): Confidence intervals (68%) for NetCounts quantities are estimated
* by the CIAO tool aprates:
* http://asc.harvard.edu/ciao/ahelp/aprates.html
*Note (10): MedianEnergy ("MedEt") is the ACIS Extract quantity
* "ENERG_PCT50_OBSERVED", the median energy of extracted events,
* corrected for background (Broos et al. 2010ApJ...714.1582B, #7.3).
*Note (11): PhotonFlux
* ("PFlx")=(NetCounts/MeanEffectiveArea/ExposureTimeNominal) (Broos et
* al. 2010ApJ...714.1582B, Section 7.4).
*Note (12): X-ray Spectral Model (XPHOT) from Getman et al. 2010, Cat.
* J/ApJ/708/1760.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Counterparts to X-ray sources and classification
probabilities
integer*4 nr__2
parameter (nr__2=25277) ! Number of records
character*491 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 Reg_2 (nr__2) ! MYStIX Galactic massive star-forming region
* (MSFR) name
character*18 MPCM_1 (nr__2) ! IAU source name (HHMMSS.ss+DDMMSS.s; J2000)
real*8 RAdeg_1 (nr__2) ! (deg) Right ascension (J2000)
real*8 DEdeg_1 (nr__2) ! (deg) Declination (J2000)
real*4 ePos_1 (nr__2) ! (arcsec) [0/2.5] 1{sigma} error around (RAdeg,DEdeg)
character*7 rPos_1 (nr__2) ! Origin of position (G1)
character*18 CXO_1 (nr__2) ! X-ray source name in IAU format
* (HHMMSS.ss+DDMMSS.s; J2000)
character*11 Xlab_1 (nr__2) ! X-ray source name used within the MYStIX
* project (G2)
character*18 NName_1 (nr__2) ! Name in NIR catalog (G3)
* (HHMMSS.ss+DDMMSS.s; J2000)
character*16 Nlab_1 (nr__2) ! Label in NIR catalog
character*18 MName_1 (nr__2) ! Name in MIR catalog
* (HHMMSS.ss+DDMMSS.s; J2000) (G4)
character*29 Mlab_1 (nr__2) ! Label in MIR catalog
real*4 pr_N_1 (nr__2) ! [0/1]? Counterpart probability, X-ray/NIR (G5)
real*4 pr_M_2 (nr__2) ! [0/1]? Counterpart probability, X-ray/MIR (G5)
real*4 Jmag_1 (nr__2) ! (mag) [4.7/25.7]? J-band magnitude
real*8 e_Jmag_1 (nr__2) ! (mag) [-25.1/3.6]? Jmag 1{sigma} uncertainty
real*4 Hmag_1 (nr__2) ! (mag) [4.7/23.9]? H-band magnitude
real*8 e_Hmag_1 (nr__2) ! (mag) [-22.2/2.9]? Hmag 1{sigma} uncertainty
real*4 Kmag_1 (nr__2) ! (mag) [4.4/19.5]? K-band magnitude
real*4 e_Kmag_1 (nr__2) ! (mag) ? Kmag 1{sigma} uncertainty
real*4 v_3_6__1 (nr__2) ! (mag) ? Spitzer/IRAC 3.6um band magnitude
real*4 e__3_6__1 (nr__2) ! (mag) ? [3.6] 1-{sigma} uncertainty
real*4 v_4_5__1 (nr__2) ! (mag) ? Spitzer/IRAC 4.5um band magnitude
real*4 e__4_5__1 (nr__2) ! (mag) ? [4.5] 1-{sigma} uncertainty
real*4 v_5_8__1 (nr__2) ! (mag) [2.3/17.7]? Spitzer/IRAC 5.8um band magnitude
real*4 e__5_8__1 (nr__2) ! (mag) [0/9.8]? [5.8] 1-{sigma} uncertainty
real*4 v_8_0__1 (nr__2) ! (mag) [1.5/17]? Spitzer/IRAC 8.0um band magnitude
real*8 e__8_0__1 (nr__2) ! (mag) [0/15.1]? [8.0] 1-{sigma} uncertainty
character*3 f_Jmag_1 (nr__2) ! UKIRT photometry flag on Jmag (G6)
character*3 f_Hmag_1 (nr__2) ! UKIRT photometry flag on Hmag (G6)
character*3 f_Kmag_1 (nr__2) ! UKIRT photometry flag on Kmag (G6)
character*3 v2M_1 (nr__2) ! 2MASS photometry flag
character*3 q_2M_1 (nr__2) ! 2MASS photometry quality flag
integer*4 n_Jmag_1 (nr__2) ! [0] GLIMPSE photometry flag on Jmag
integer*4 n_Hmag_1 (nr__2) ! [0] GLIMPSE photometry flag on Hmag
integer*4 n_Kmag_1 (nr__2) ! [0] GLIMPSE photometry flag on Kmag
integer*4 f__3_6__1 (nr__2) ! GLIMPSE photometry flag on [3.6]
integer*4 f__4_5__1 (nr__2) ! GLIMPSE photometry flag on [4.5]
integer*4 f__5_8__1 (nr__2) ! GLIMPSE photometry flag on [5.8]
integer*4 f__8_0__1 (nr__2) ! GLIMPSE photometry flag on [8.0]
integer*4 LS_1 (nr__2) ! [-2/4]?=-99 Local Spitzer photometry flag (G7)
character*7 r_Jmag_1 (nr__2) ! Origin of Jmag photometry
character*7 r_Hmag_1 (nr__2) ! Origin of Hmag photometry
character*7 r_Kmag_1 (nr__2) ! Origin of Kmag photometry
character*13 r__3_6__1 (nr__2) ! Origin of [3.6] photometry
character*13 r__4_5__1 (nr__2) ! Origin of [4.5] photometry
character*13 r__5_8__1 (nr__2) ! Origin of [5.8] photometry
character*13 r__8_0__1 (nr__2) ! Origin of [8.0] photometry
integer*4 Cl_1 (nr__2) ! [-2/2]?=-99 Classification (G8)
real*4 Av_1 (nr__2) ! (mag) [0/61]? Visual extinction A_V_
integer*4 St_1 (nr__2) ! [-1/2]?=-99 YSO stage (G9)
real*4 H1pr_1 (nr__2) ! [0/1] H1 class prior probability
* (position-dependent)
real*4 H2pr_1 (nr__2) ! [0/1] H2 class prior probability
* (position-dependent)
real*4 H3pr_1 (nr__2) ! [0/1] H3 class prior probability
* (position-dependent)
real*4 H4pr_1 (nr__2) ! [0/1] H4 class prior probability
* (position-dependent)
real*4 H1ps_1 (nr__2) ! [0/1]? H1 class posterior probability
real*4 H2ps_1 (nr__2) ! [0/1]? H2 class posterior probability
real*4 H3ps_1 (nr__2) ! [0/1]? H3 class posterior probability
real*4 H4ps_1 (nr__2) ! [0/1]? H4 class posterior probability
integer*4 H2f_1 (nr__2) ! [1/7] H2 dominant classification term (G10)
integer*4 Xcl_1 (nr__2) ! [0/4]? X-ray classification (G11)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__3
parameter (nr__3=42) ! Number of records
character*97 ar__3 ! Full-size record
character*19 Ref (nr__3) ! Reference code
character*19 BibCode (nr__3) ! Bibcode
character*22 Aut (nr__3) ! Author's name(s)
character*34 Comm (nr__3) ! Comment
C=============================================================================
C Loading file 'reg.dat' ! List of the 20 nearby Galactic massive
* star-forming regions; table added by CDS
C Format for file interpretation
1 format(A16,1X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'reg.dat')
write(6,*) '....Loading file: reg.dat'
do i__=1,20
read(1,'(A22)')ar__
read(ar__,1)Name(i__),Reg(i__)
c ..............Just test output...........
write(6,1)Name(i__),Reg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'mpcm.dat' ! MPCM (MYStIX Probable Complex Members) sources
* and properties (table 2)
C Format for file interpretation
2 format(
+ A16,1X,A18,A1,1X,F10.6,1X,F10.6,1X,F5.3,1X,A11,1X,A18,1X,A11,
+ 1X,A18,1X,A16,1X,A18,1X,A29,1X,A30,1X,I5,1X,I4,1X,A30,1X,A14,
+ 1X,F6.3,1X,A3,1X,F4.2,1X,F4.2,1X,F4.2,1X,F9.2,1X,F4.2,1X,I3,
+ 1X,I3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F4.2,1X,F4.2,1X,I5,1X,F8.2,
+ 1X,F8.2,1X,F8.2,1X,F8.2,1X,F8.2,1X,F8.2,1X,F8.2,1X,F8.2,1X,
+ F8.2,1X,F5.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.3,1X,F7.3,1X,F6.3,1X,F7.3,
+ 1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F6.3,1X,F6.3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,I4,1X,
+ I4,1X,I4,1X,I5,1X,I5,1X,I6,1X,I6,1X,I3,1X,A19,1X,A19,1X,A19,
+ 1X,A19,1X,A19,1X,A19,1X,A19,1X,I3,1X,F6.3,1X,I3,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,I3,1X,
+ I3)
C Effective file loading
open(unit=1,status='old',file=
+'mpcm.dat')
write(6,*) '....Loading file: mpcm.dat'
do i__=1,31549
read(1,'(A971)')ar__1
read(ar__1,2)
+ Reg_1(i__),MPCM(i__),m_MPCM(i__),RAdeg(i__),DEdeg(i__),
+ ePos(i__),rPos(i__),CXO(i__),Xlab(i__),NName(i__),Nlab(i__),
+ MName(i__),Mlab(i__),Olab(i__),Xind(i__),Iind(i__),SpT(i__),
+ r_SpT(i__),OBmag(i__),OBband(i__),pr_0(i__),pr_1(i__),
+ pr_m(i__),t_exp(i__),f_exp(i__),Nobs(i__),Nph(i__),
+ b_theta(i__),theta(i__),B_theta_1(i__),PSFf(i__),aglow(i__),
+ Oct_t(i__),Nct_t(i__),Nct_s(i__),Nct_h(i__),b_Nct_t(i__),
+ B_Nct_t_1(i__),b_Nct_s(i__),B_Nct_s_1(i__),b_Nct_h(i__),
+ B_Nct_h_1(i__),E_t(i__),XFlx_t(i__),XFlx_s(i__),XFlx_h(i__),
+ LX_h(i__),LXc_h(i__),e_LXc_h(i__),s_LXc_h(i__),LX_t(i__),
+ LXc_t(i__),e_LXc_t(i__),s_LXc_t(i__),logNH(i__),e_logNH(i__),
+ s_logNH(i__),pr_N(i__),pr_M_1(i__),Jmag(i__),e_Jmag(i__),
+ Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),v_3_6_(i__),
+ e__3_6_(i__),v_4_5_(i__),e__4_5_(i__),v_5_8_(i__),
+ e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),f_Jmag(i__),f_Hmag(i__),
+ f_Kmag(i__),v2M(i__),q_2M(i__),n_Jmag(i__),n_Hmag(i__),
+ n_Kmag(i__),f__3_6_(i__),f__4_5_(i__),f__5_8_(i__),
+ f__8_0_(i__),LS(i__),r_Jmag(i__),r_Hmag(i__),r_Kmag(i__),
+ r__3_6_(i__),r__4_5_(i__),r__5_8_(i__),r__8_0_(i__),Cl(i__),
+ Av(i__),St(i__),H1pr(i__),H2pr(i__),H3pr(i__),H4pr(i__),
+ H1ps(i__),H2ps(i__),H3ps(i__),H4ps(i__),H2f(i__),Xcl(i__)
if (idig(ar__1(60:64)).EQ.0) ePos(i__) = rNULL__
if(ar__1(282:287) .EQ. '') OBmag(i__) = rNULL__
if(ar__1(293:296) .EQ. '') pr_0(i__) = rNULL__
if(ar__1(298:301) .EQ. '') pr_1(i__) = rNULL__
if(ar__1(303:306) .EQ. '') pr_m(i__) = rNULL__
if(ar__1(308:316) .EQ. '') t_exp(i__) = rNULL__
if(ar__1(318:321) .EQ. '') f_exp(i__) = rNULL__
if(ar__1(331:336) .EQ. '') b_theta(i__) = rNULL__
if(ar__1(338:343) .EQ. '') theta(i__) = rNULL__
if(ar__1(345:350) .EQ. '') B_theta_1(i__) = rNULL__
if(ar__1(352:355) .EQ. '') PSFf(i__) = rNULL__
if(ar__1(357:360) .EQ. '') aglow(i__) = rNULL__
if(ar__1(368:375) .EQ. '') Nct_t(i__) = rNULL__
if(ar__1(377:384) .EQ. '') Nct_s(i__) = rNULL__
if(ar__1(386:393) .EQ. '') Nct_h(i__) = rNULL__
if(ar__1(395:402) .EQ. '') b_Nct_t(i__) = rNULL__
if(ar__1(404:411) .EQ. '') B_Nct_t_1(i__) = rNULL__
if(ar__1(413:420) .EQ. '') b_Nct_s(i__) = rNULL__
if(ar__1(422:429) .EQ. '') B_Nct_s_1(i__) = rNULL__
if(ar__1(431:438) .EQ. '') b_Nct_h(i__) = rNULL__
if(ar__1(440:447) .EQ. '') B_Nct_h_1(i__) = rNULL__
if(ar__1(449:453) .EQ. '') E_t(i__) = rNULL__
if(ar__1(455:461) .EQ. '') XFlx_t(i__) = rNULL__
if(ar__1(463:469) .EQ. '') XFlx_s(i__) = rNULL__
if(ar__1(471:477) .EQ. '') XFlx_h(i__) = rNULL__
if(ar__1(479:490) .EQ. '') LX_h(i__) = rNULL__
if(ar__1(492:503) .EQ. '') LXc_h(i__) = rNULL__
if(ar__1(505:516) .EQ. '') e_LXc_h(i__) = rNULL__
if(ar__1(518:529) .EQ. '') s_LXc_h(i__) = rNULL__
if(ar__1(531:542) .EQ. '') LX_t(i__) = rNULL__
if(ar__1(544:555) .EQ. '') LXc_t(i__) = rNULL__
if(ar__1(557:568) .EQ. '') e_LXc_t(i__) = rNULL__
if(ar__1(570:581) .EQ. '') s_LXc_t(i__) = rNULL__
if(ar__1(583:587) .EQ. '') logNH(i__) = rNULL__
if(ar__1(589:592) .EQ. '') e_logNH(i__) = rNULL__
if(ar__1(594:597) .EQ. '') s_logNH(i__) = rNULL__
if(ar__1(599:602) .EQ. '') pr_N(i__) = rNULL__
if(ar__1(604:607) .EQ. '') pr_M_1(i__) = rNULL__
if(ar__1(609:614) .EQ. '') Jmag(i__) = rNULL__
if(ar__1(616:622) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__1(624:629) .EQ. '') Hmag(i__) = rNULL__
if(ar__1(631:637) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__1(639:644) .EQ. '') Kmag(i__) = rNULL__
if(ar__1(646:650) .EQ. '') e_Kmag(i__) = rNULL__
if(ar__1(652:657) .EQ. '') v_3_6_(i__) = rNULL__
if(ar__1(659:663) .EQ. '') e__3_6_(i__) = rNULL__
if(ar__1(665:670) .EQ. '') v_4_5_(i__) = rNULL__
if(ar__1(672:676) .EQ. '') e__4_5_(i__) = rNULL__
if(ar__1(678:683) .EQ. '') v_5_8_(i__) = rNULL__
if(ar__1(685:689) .EQ. '') e__5_8_(i__) = rNULL__
if(ar__1(691:696) .EQ. '') v_8_0_(i__) = rNULL__
if(ar__1(698:703) .EQ. '') e__8_0_(i__) = rNULL__
if(ar__1(914:919) .EQ. '') Av(i__) = rNULL__
if(ar__1(925:928) .EQ. '') H1pr(i__) = rNULL__
if(ar__1(930:933) .EQ. '') H2pr(i__) = rNULL__
if(ar__1(935:938) .EQ. '') H3pr(i__) = rNULL__
if(ar__1(940:943) .EQ. '') H4pr(i__) = rNULL__
if(ar__1(945:948) .EQ. '') H1ps(i__) = rNULL__
if(ar__1(950:953) .EQ. '') H2ps(i__) = rNULL__
if(ar__1(955:958) .EQ. '') H3ps(i__) = rNULL__
if(ar__1(960:963) .EQ. '') H4ps(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Reg_1(i__),MPCM(i__),m_MPCM(i__),RAdeg(i__),DEdeg(i__),
+ ePos(i__),rPos(i__),CXO(i__),Xlab(i__),NName(i__),Nlab(i__),
+ MName(i__),Mlab(i__),Olab(i__),Xind(i__),Iind(i__),SpT(i__),
+ r_SpT(i__),OBmag(i__),OBband(i__),pr_0(i__),pr_1(i__),
+ pr_m(i__),t_exp(i__),f_exp(i__),Nobs(i__),Nph(i__),
+ b_theta(i__),theta(i__),B_theta_1(i__),PSFf(i__),aglow(i__),
+ Oct_t(i__),Nct_t(i__),Nct_s(i__),Nct_h(i__),b_Nct_t(i__),
+ B_Nct_t_1(i__),b_Nct_s(i__),B_Nct_s_1(i__),b_Nct_h(i__),
+ B_Nct_h_1(i__),E_t(i__),XFlx_t(i__),XFlx_s(i__),XFlx_h(i__),
+ LX_h(i__),LXc_h(i__),e_LXc_h(i__),s_LXc_h(i__),LX_t(i__),
+ LXc_t(i__),e_LXc_t(i__),s_LXc_t(i__),logNH(i__),e_logNH(i__),
+ s_logNH(i__),pr_N(i__),pr_M_1(i__),Jmag(i__),e_Jmag(i__),
+ Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),v_3_6_(i__),
+ e__3_6_(i__),v_4_5_(i__),e__4_5_(i__),v_5_8_(i__),
+ e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),f_Jmag(i__),f_Hmag(i__),
+ f_Kmag(i__),v2M(i__),q_2M(i__),n_Jmag(i__),n_Hmag(i__),
+ n_Kmag(i__),f__3_6_(i__),f__4_5_(i__),f__5_8_(i__),
+ f__8_0_(i__),LS(i__),r_Jmag(i__),r_Hmag(i__),r_Kmag(i__),
+ r__3_6_(i__),r__4_5_(i__),r__5_8_(i__),r__8_0_(i__),Cl(i__),
+ Av(i__),St(i__),H1pr(i__),H2pr(i__),H3pr(i__),H4pr(i__),
+ H1ps(i__),H2ps(i__),H3ps(i__),H4ps(i__),H2f(i__),Xcl(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Counterparts to X-ray sources and classification
* probabilities
C Format for file interpretation
3 format(
+ A16,1X,A18,1X,F10.6,1X,F10.6,1X,F5.3,1X,A7,1X,A18,1X,A11,1X,
+ A18,1X,A16,1X,A18,1X,A29,1X,F4.2,1X,F4.2,1X,F6.3,1X,F7.3,1X,
+ F6.3,1X,F7.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
+ 1X,F6.3,1X,F5.3,1X,F6.3,1X,F8.3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,
+ 1X,I1,1X,I1,1X,I1,1X,I5,1X,I6,1X,I6,1X,I6,1X,I3,1X,A7,1X,A7,
+ 1X,A7,1X,A13,1X,A13,1X,A13,1X,A13,1X,I3,1X,F6.3,1X,I3,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,I1,
+ 1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,25277
read(1,'(A491)')ar__2
read(ar__2,3)
+ Reg_2(i__),MPCM_1(i__),RAdeg_1(i__),DEdeg_1(i__),ePos_1(i__),
+ rPos_1(i__),CXO_1(i__),Xlab_1(i__),NName_1(i__),Nlab_1(i__),
+ MName_1(i__),Mlab_1(i__),pr_N_1(i__),pr_M_2(i__),Jmag_1(i__),
+ e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),Kmag_1(i__),
+ e_Kmag_1(i__),v_3_6__1(i__),e__3_6__1(i__),v_4_5__1(i__),
+ e__4_5__1(i__),v_5_8__1(i__),e__5_8__1(i__),v_8_0__1(i__),
+ e__8_0__1(i__),f_Jmag_1(i__),f_Hmag_1(i__),f_Kmag_1(i__),
+ v2M_1(i__),q_2M_1(i__),n_Jmag_1(i__),n_Hmag_1(i__),
+ n_Kmag_1(i__),f__3_6__1(i__),f__4_5__1(i__),f__5_8__1(i__),
+ f__8_0__1(i__),LS_1(i__),r_Jmag_1(i__),r_Hmag_1(i__),
+ r_Kmag_1(i__),r__3_6__1(i__),r__4_5__1(i__),r__5_8__1(i__),
+ r__8_0__1(i__),Cl_1(i__),Av_1(i__),St_1(i__),H1pr_1(i__),
+ H2pr_1(i__),H3pr_1(i__),H4pr_1(i__),H1ps_1(i__),H2ps_1(i__),
+ H3ps_1(i__),H4ps_1(i__),H2f_1(i__),Xcl_1(i__)
if(ar__2(189:192) .EQ. '') pr_N_1(i__) = rNULL__
if(ar__2(194:197) .EQ. '') pr_M_2(i__) = rNULL__
if(ar__2(199:204) .EQ. '') Jmag_1(i__) = rNULL__
if(ar__2(206:212) .EQ. '') e_Jmag_1(i__) = rNULL__
if(ar__2(214:219) .EQ. '') Hmag_1(i__) = rNULL__
if(ar__2(221:227) .EQ. '') e_Hmag_1(i__) = rNULL__
if(ar__2(229:234) .EQ. '') Kmag_1(i__) = rNULL__
if(ar__2(236:240) .EQ. '') e_Kmag_1(i__) = rNULL__
if(ar__2(242:247) .EQ. '') v_3_6__1(i__) = rNULL__
if(ar__2(249:253) .EQ. '') e__3_6__1(i__) = rNULL__
if(ar__2(255:260) .EQ. '') v_4_5__1(i__) = rNULL__
if(ar__2(262:266) .EQ. '') e__4_5__1(i__) = rNULL__
if(ar__2(268:273) .EQ. '') v_5_8__1(i__) = rNULL__
if(ar__2(275:279) .EQ. '') e__5_8__1(i__) = rNULL__
if(ar__2(281:286) .EQ. '') v_8_0__1(i__) = rNULL__
if(ar__2(288:295) .EQ. '') e__8_0__1(i__) = rNULL__
if(ar__2(438:443) .EQ. '') Av_1(i__) = rNULL__
if(ar__2(469:472) .EQ. '') H1ps_1(i__) = rNULL__
if(ar__2(474:477) .EQ. '') H2ps_1(i__) = rNULL__
if(ar__2(479:482) .EQ. '') H3ps_1(i__) = rNULL__
if(ar__2(484:487) .EQ. '') H4ps_1(i__) = rNULL__
if(ar__2(491:491) .EQ. '') Xcl_1(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Reg_2(i__),MPCM_1(i__),RAdeg_1(i__),DEdeg_1(i__),ePos_1(i__),
+ rPos_1(i__),CXO_1(i__),Xlab_1(i__),NName_1(i__),Nlab_1(i__),
+ MName_1(i__),Mlab_1(i__),pr_N_1(i__),pr_M_2(i__),Jmag_1(i__),
+ e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),Kmag_1(i__),
+ e_Kmag_1(i__),v_3_6__1(i__),e__3_6__1(i__),v_4_5__1(i__),
+ e__4_5__1(i__),v_5_8__1(i__),e__5_8__1(i__),v_8_0__1(i__),
+ e__8_0__1(i__),f_Jmag_1(i__),f_Hmag_1(i__),f_Kmag_1(i__),
+ v2M_1(i__),q_2M_1(i__),n_Jmag_1(i__),n_Hmag_1(i__),
+ n_Kmag_1(i__),f__3_6__1(i__),f__4_5__1(i__),f__5_8__1(i__),
+ f__8_0__1(i__),LS_1(i__),r_Jmag_1(i__),r_Hmag_1(i__),
+ r_Kmag_1(i__),r__3_6__1(i__),r__4_5__1(i__),r__5_8__1(i__),
+ r__8_0__1(i__),Cl_1(i__),Av_1(i__),St_1(i__),H1pr_1(i__),
+ H2pr_1(i__),H3pr_1(i__),H4pr_1(i__),H1ps_1(i__),H2ps_1(i__),
+ H3ps_1(i__),H4ps_1(i__),H2f_1(i__),Xcl_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
4 format(A19,1X,A19,1X,A22,1X,A34)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,42
read(1,'(A97)')ar__3
read(ar__3,4)Ref(i__),BibCode(i__),Aut(i__),Comm(i__)
c ..............Just test output...........
write(6,4)Ref(i__),BibCode(i__),Aut(i__),Comm(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end