Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/209/32 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 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 ! MYStIX Galactic massive star-forming region
* (MSFR) name
character*5 Reg ! 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 ! MYStIX Galactic massive star-forming region
* (MSFR) name
character*18 MPCM ! IAU source name (HHMMSS.ss+DDMMSS.s; J2000)
character*1 m_MPCM ! [ab] Multiple MPCM name component (1)
real*8 RAdeg ! (deg) Right Ascension (J2000)
real*8 DEdeg ! (deg) Declination (J2000)
real*4 ePos ! (arcsec) [0/2.1]?=- 1{sigma} error around (RA,DE)
character*11 rPos ! Origin of position (G1)
character*18 CXO ! X-ray source name in IAU format
* (HHMMSS.s+DDMMSS; J2000)
character*11 Xlab ! X-ray source name used within the MYStIX
* project (G2)
character*18 NName ! Name in NIR catalog (G3)
character*16 Nlab ! Label in NIR catalog
character*18 MName ! Name in MIR catalog
* (HHMMSS.ss+DDMMSS.s; J2000 or NNNN) (G4)
character*29 Mlab ! Label in MIR catalog
character*30 Olab ! Label in OB catalog (2)
integer*4 Xind ! ?=-99 0-based index in X-ray catalog (Sec. 4)
integer*4 Iind ! ?=-99 0-based index in IR excess catalog
* (Povich et al. 2013, Cat. J/ApJS/209/31)
character*30 SpT ! MK spectral type
character*14 r_SpT ! Reference for spectral type (3)
real*4 OBmag ! (mag) [4.4/24.5]? Visual photometry in OB Catalog
character*3 OBband ! Visual band for OBmag (J, K or V)
real*4 pr_0 ! [0/1]? p-value for no-source hypothesis (4)(5)
real*4 pr_1 ! [0/1]? p-value for single observation (4)(6)
real*4 pr_m ! [0/1]? p-value for merged observations (4)(6)
real*8 t_exp ! (s) ? Total exposure time in merged
* observations (4)
real*4 f_exp ! [0/1]? Fraction of "Exp" that source was
* observed (4)(7)
integer*4 Nobs ! [1/9]?=-99 Total number of observations
* extracted (4)
integer*4 Nph ! [1/9]?=-99 Number of observations merged to
* estimate photometry properties (4)
real*4 b_theta ! (arcmin) ? Smallest off-axis angle for merged
* observations (4)
real*4 theta ! (arcmin) ? Average off-axis angle {theta} for merged
* observations (4)
real*4 B_theta_1 ! (arcmin) ? Largest off-axis angle for merged obs. (4)
real*4 PSFf ! [0/1]? Average PSF fraction (at 1.5 keV) for
* merged observations (4)
real*4 aglow ! [0/1]? Suspected afterglow fraction (4)(8)
integer*4 Oct_t ! (ct) ?=-99 Extracted total (0.5-8keV) counts in
* merged apertures (4)
real*8 Nct_t ! (ct) [1/89438]? Net total (0.5-8keV) counts in
* merged apertures (4)
real*8 Nct_s ! (ct) ? Net soft (0.5-2 keV) counts in merged
* apertures (4)
real*8 Nct_h ! (ct) ? Net hard (2-8 keV) counts in merged
* apertures (4)
real*8 b_Nct_t ! (ct) ? 1{sigma} lower bound on Nct-t (4)(9)
real*8 B_Nct_t_1 ! (ct) ? 1{sigma} upper bound on Nct-t (4)(9)
real*8 b_Nct_s ! (ct) ? 1{sigma} lower bound on Nct-s (4)(9)
real*8 B_Nct_s_1 ! (ct) ? 1{sigma} upper bound on Nct-s (4)(9)
real*8 b_Nct_h ! (ct) ? 1{sigma} lower bound on Nct-h (4)(9)
real*8 B_Nct_h_1 ! (ct) ? 1{sigma} upper bound on Nct-h (4)(9)
real*4 E_t ! (keV) ? Median 0.5-8keV energy, observed
* spectrum (4)(10)
real*8 XFlx_t ! ([ph/cm2/s]) ? Log of incident photon flux in total
* 0.5-8 keV band (4)(11)
real*8 XFlx_s ! ([ph/cm2/s]) ? Log of incident photon flux in soft
* 0.5-2 keV band (4)(11)
real*8 XFlx_h ! ([ph/cm2/s]) ? Log of incident photon flux in hard
* 2-8 keV band (4)(11)
real*4 LX_h ! (10-7W) ? X-ray luminosity in erg/s, 2-8 keV (12)
real*4 LXc_h ! (10-7W) ? Absorption-corrected X-ray luminosity,
* in erg/s; 2-8 keV (12)
real*4 e_LXc_h ! (10-7W) ? LXhc 1{sigma} statistical uncertainty (12)
real*4 s_LXc_h ! (10-7W) ? LXhc 1{sigma} systematic uncertainty (12)
real*4 LX_t ! (10-7W) ? X-ray luminosity in erg/s, 0.5-8 keV (12)
real*4 LXc_t ! (10-7W) ? Absorption-corrected X-ray luminosity,
* in erg/s; 0.5-8 kev (12)
real*4 e_LXc_t ! (10-7W) ? LXtc 1{sigma} statistical uncertainty (12)
real*4 s_LXc_t ! (10-7W) ? LXtc 1{sigma} systematic uncertainty (12)
real*4 logNH ! (cm-2) ? Gas column density (12)
real*4 e_logNH ! (cm-2) ? logNH 1{sigma} statistical error (12)
real*4 s_logNH ! (cm-2) ? logNH 1{sigma} systematic uncertainty (12)
real*4 pr_N ! [0/1]? Counterpart probability, X-ray/NIR (G5)
real*4 pr_M_1 ! [0/1]? counterpart probability, X-ray/MIR (G5)
real*4 Jmag ! (mag) [4.7/25.2]? J-band magnitude
real*8 e_Jmag ! (mag) [-25.1/3.3]? Jmag 1{sigma} uncertainty
real*4 Hmag ! (mag) [4.5/29.9]? H-band magnitude
real*8 e_Hmag ! (mag) [-22.2/6.6]? Hmag 1{sigma} uncertainty
real*4 Kmag ! (mag) [4.3/21.5]? K-band magnitude
real*4 e_Kmag ! (mag) ? Kmag 1-{sigma} uncertainty
real*4 v_3_6_ ! (mag) ? Spitzer/IRAC 3.6um band magnitude
real*4 e__3_6_ ! (mag) ? [3.6] 1-{sigma} uncertainty
real*4 v_4_5_ ! (mag) ? Spitzer/IRAC 4.5um band magnitude
real*4 e__4_5_ ! (mag) ? [4.5] 1-{sigma} uncertainty
real*4 v_5_8_ ! (mag) [2/17.5]? Spitzer/IRAC 5.8um band magnitude
real*4 e__5_8_ ! (mag) [0/10]? [5.8] 1-{sigma} uncertainty
real*4 v_8_0_ ! (mag) [1.5/17]? Spitzer/IRAC 8.0um band magnitude
real*4 e__8_0_ ! (mag) [0/15]? [8.0] 1-{sigma} uncertainty
character*3 f_Jmag ! UKIRT photometry flag on Jmag (G6)
character*3 f_Hmag ! UKIRT photometry flag on Hmag (G6)
character*3 f_Kmag ! UKIRT photometry flag on Kmag (G6)
character*3 v2M ! 2MASS photometry flag
character*3 q_2M ! 2MASS quality photometry flag
integer*4 n_Jmag ! GLIMPSE photometry flag on Jmag
integer*4 n_Hmag ! GLIMPSE photometry flag on Hmag
integer*4 n_Kmag ! GLIMPSE photometry flag on Kmag
integer*4 f__3_6_ ! GLIMPSE photometry flag on [3.6]
integer*4 f__4_5_ ! GLIMPSE photometry flag on [4.5]
integer*4 f__5_8_ ! GLIMPSE photometry flag on [5.8]
integer*4 f__8_0_ ! GLIMPSE photometry flag on [8.0]
integer*4 LS ! [-2/4]?=-99 Local Spitzer photometry flag (G7)
character*19 r_Jmag ! Origin of Jmag photometry
character*19 r_Hmag ! Origin of Hmag photometry
character*19 r_Kmag ! Origin of Kmag photometry
character*19 r__3_6_ ! Origin of [3.6] photometry
character*19 r__4_5_ ! Origin of [4.5] photometry
character*19 r__5_8_ ! Origin of [5.8] photometry
character*19 r__8_0_ ! Origin of [8.0] photometry
integer*4 Cl ! [-2/2]?=-99 Classification (G8)
real*4 Av ! (mag) [0/81]? Visual extinction A_V_
integer*4 St ! [-1/2]?=-99 YSO stage (G9)
real*4 H1pr ! [0/1]? H1 class prior probability
* (position-dependent)
real*4 H2pr ! [0/1]? H2 class prior probability
* (position-dependent)
real*4 H3pr ! [0/1]? H3 class prior probability
* (position-dependent)
real*4 H4pr ! [0/1]? H4 class prior probability
* (position-dependent)
real*4 H1ps ! [0/1]? H1 class posterior probability
real*4 H2ps ! [0/1]? H2 class posterior probability
real*4 H3ps ! [0/1]? H3 class posterior probability
real*4 H4ps ! [0/1]? H4 class posterior probability
integer*4 H2f ! [1/7]?=-99 H2 dominant classification
* term (G10)
integer*4 Xcl ! [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 ! MYStIX Galactic massive star-forming region
* (MSFR) name
character*18 MPCM_1 ! IAU source name (HHMMSS.ss+DDMMSS.s; J2000)
real*8 RAdeg_1 ! (deg) Right ascension (J2000)
real*8 DEdeg_1 ! (deg) Declination (J2000)
real*4 ePos_1 ! (arcsec) [0/2.5] 1{sigma} error around (RAdeg,DEdeg)
character*7 rPos_1 ! Origin of position (G1)
character*18 CXO_1 ! X-ray source name in IAU format
* (HHMMSS.ss+DDMMSS.s; J2000)
character*11 Xlab_1 ! X-ray source name used within the MYStIX
* project (G2)
character*18 NName_1 ! Name in NIR catalog (G3)
* (HHMMSS.ss+DDMMSS.s; J2000)
character*16 Nlab_1 ! Label in NIR catalog
character*18 MName_1 ! Name in MIR catalog
* (HHMMSS.ss+DDMMSS.s; J2000) (G4)
character*29 Mlab_1 ! Label in MIR catalog
real*4 pr_N_1 ! [0/1]? Counterpart probability, X-ray/NIR (G5)
real*4 pr_M_2 ! [0/1]? Counterpart probability, X-ray/MIR (G5)
real*4 Jmag_1 ! (mag) [4.7/25.7]? J-band magnitude
real*8 e_Jmag_1 ! (mag) [-25.1/3.6]? Jmag 1{sigma} uncertainty
real*4 Hmag_1 ! (mag) [4.7/23.9]? H-band magnitude
real*8 e_Hmag_1 ! (mag) [-22.2/2.9]? Hmag 1{sigma} uncertainty
real*4 Kmag_1 ! (mag) [4.4/19.5]? K-band magnitude
real*4 e_Kmag_1 ! (mag) ? Kmag 1{sigma} uncertainty
real*4 v_3_6__1 ! (mag) ? Spitzer/IRAC 3.6um band magnitude
real*4 e__3_6__1 ! (mag) ? [3.6] 1-{sigma} uncertainty
real*4 v_4_5__1 ! (mag) ? Spitzer/IRAC 4.5um band magnitude
real*4 e__4_5__1 ! (mag) ? [4.5] 1-{sigma} uncertainty
real*4 v_5_8__1 ! (mag) [2.3/17.7]? Spitzer/IRAC 5.8um band magnitude
real*4 e__5_8__1 ! (mag) [0/9.8]? [5.8] 1-{sigma} uncertainty
real*4 v_8_0__1 ! (mag) [1.5/17]? Spitzer/IRAC 8.0um band magnitude
real*8 e__8_0__1 ! (mag) [0/15.1]? [8.0] 1-{sigma} uncertainty
character*3 f_Jmag_1 ! UKIRT photometry flag on Jmag (G6)
character*3 f_Hmag_1 ! UKIRT photometry flag on Hmag (G6)
character*3 f_Kmag_1 ! UKIRT photometry flag on Kmag (G6)
character*3 v2M_1 ! 2MASS photometry flag
character*3 q_2M_1 ! 2MASS photometry quality flag
integer*4 n_Jmag_1 ! [0] GLIMPSE photometry flag on Jmag
integer*4 n_Hmag_1 ! [0] GLIMPSE photometry flag on Hmag
integer*4 n_Kmag_1 ! [0] GLIMPSE photometry flag on Kmag
integer*4 f__3_6__1 ! GLIMPSE photometry flag on [3.6]
integer*4 f__4_5__1 ! GLIMPSE photometry flag on [4.5]
integer*4 f__5_8__1 ! GLIMPSE photometry flag on [5.8]
integer*4 f__8_0__1 ! GLIMPSE photometry flag on [8.0]
integer*4 LS_1 ! [-2/4]?=-99 Local Spitzer photometry flag (G7)
character*7 r_Jmag_1 ! Origin of Jmag photometry
character*7 r_Hmag_1 ! Origin of Hmag photometry
character*7 r_Kmag_1 ! Origin of Kmag photometry
character*13 r__3_6__1 ! Origin of [3.6] photometry
character*13 r__4_5__1 ! Origin of [4.5] photometry
character*13 r__5_8__1 ! Origin of [5.8] photometry
character*13 r__8_0__1 ! Origin of [8.0] photometry
integer*4 Cl_1 ! [-2/2]?=-99 Classification (G8)
real*4 Av_1 ! (mag) [0/61]? Visual extinction A_V_
integer*4 St_1 ! [-1/2]?=-99 YSO stage (G9)
real*4 H1pr_1 ! [0/1] H1 class prior probability
* (position-dependent)
real*4 H2pr_1 ! [0/1] H2 class prior probability
* (position-dependent)
real*4 H3pr_1 ! [0/1] H3 class prior probability
* (position-dependent)
real*4 H4pr_1 ! [0/1] H4 class prior probability
* (position-dependent)
real*4 H1ps_1 ! [0/1]? H1 class posterior probability
real*4 H2ps_1 ! [0/1]? H2 class posterior probability
real*4 H3ps_1 ! [0/1]? H3 class posterior probability
real*4 H4ps_1 ! [0/1]? H4 class posterior probability
integer*4 H2f_1 ! [1/7] H2 dominant classification term (G10)
integer*4 Xcl_1 ! [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 ! Reference code
character*19 BibCode ! Bibcode
character*22 Aut ! Author's name(s)
character*34 Comm ! 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,Reg
c ..............Just test output...........
write(6,1)Name,Reg
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,MPCM,m_MPCM,RAdeg,DEdeg,ePos,rPos,CXO,Xlab,NName,Nlab,
+ MName,Mlab,Olab,Xind,Iind,SpT,r_SpT,OBmag,OBband,pr_0,pr_1,
+ pr_m,t_exp,f_exp,Nobs,Nph,b_theta,theta,B_theta_1,PSFf,aglow,
+ Oct_t,Nct_t,Nct_s,Nct_h,b_Nct_t,B_Nct_t_1,b_Nct_s,B_Nct_s_1,
+ b_Nct_h,B_Nct_h_1,E_t,XFlx_t,XFlx_s,XFlx_h,LX_h,LXc_h,e_LXc_h,
+ s_LXc_h,LX_t,LXc_t,e_LXc_t,s_LXc_t,logNH,e_logNH,s_logNH,pr_N,
+ pr_M_1,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,v_3_6_,e__3_6_,
+ v_4_5_,e__4_5_,v_5_8_,e__5_8_,v_8_0_,e__8_0_,f_Jmag,f_Hmag,
+ f_Kmag,v2M,q_2M,n_Jmag,n_Hmag,n_Kmag,f__3_6_,f__4_5_,f__5_8_,
+ f__8_0_,LS,r_Jmag,r_Hmag,r_Kmag,r__3_6_,r__4_5_,r__5_8_,
+ r__8_0_,Cl,Av,St,H1pr,H2pr,H3pr,H4pr,H1ps,H2ps,H3ps,H4ps,H2f,
+ Xcl
if (idig(ar__1(60:64)).EQ.0) ePos = rNULL__
if(ar__1(282:287) .EQ. '') OBmag = rNULL__
if(ar__1(293:296) .EQ. '') pr_0 = rNULL__
if(ar__1(298:301) .EQ. '') pr_1 = rNULL__
if(ar__1(303:306) .EQ. '') pr_m = rNULL__
if(ar__1(308:316) .EQ. '') t_exp = rNULL__
if(ar__1(318:321) .EQ. '') f_exp = rNULL__
if(ar__1(331:336) .EQ. '') b_theta = rNULL__
if(ar__1(338:343) .EQ. '') theta = rNULL__
if(ar__1(345:350) .EQ. '') B_theta_1 = rNULL__
if(ar__1(352:355) .EQ. '') PSFf = rNULL__
if(ar__1(357:360) .EQ. '') aglow = rNULL__
if(ar__1(368:375) .EQ. '') Nct_t = rNULL__
if(ar__1(377:384) .EQ. '') Nct_s = rNULL__
if(ar__1(386:393) .EQ. '') Nct_h = rNULL__
if(ar__1(395:402) .EQ. '') b_Nct_t = rNULL__
if(ar__1(404:411) .EQ. '') B_Nct_t_1 = rNULL__
if(ar__1(413:420) .EQ. '') b_Nct_s = rNULL__
if(ar__1(422:429) .EQ. '') B_Nct_s_1 = rNULL__
if(ar__1(431:438) .EQ. '') b_Nct_h = rNULL__
if(ar__1(440:447) .EQ. '') B_Nct_h_1 = rNULL__
if(ar__1(449:453) .EQ. '') E_t = rNULL__
if(ar__1(455:461) .EQ. '') XFlx_t = rNULL__
if(ar__1(463:469) .EQ. '') XFlx_s = rNULL__
if(ar__1(471:477) .EQ. '') XFlx_h = rNULL__
if(ar__1(479:490) .EQ. '') LX_h = rNULL__
if(ar__1(492:503) .EQ. '') LXc_h = rNULL__
if(ar__1(505:516) .EQ. '') e_LXc_h = rNULL__
if(ar__1(518:529) .EQ. '') s_LXc_h = rNULL__
if(ar__1(531:542) .EQ. '') LX_t = rNULL__
if(ar__1(544:555) .EQ. '') LXc_t = rNULL__
if(ar__1(557:568) .EQ. '') e_LXc_t = rNULL__
if(ar__1(570:581) .EQ. '') s_LXc_t = rNULL__
if(ar__1(583:587) .EQ. '') logNH = rNULL__
if(ar__1(589:592) .EQ. '') e_logNH = rNULL__
if(ar__1(594:597) .EQ. '') s_logNH = rNULL__
if(ar__1(599:602) .EQ. '') pr_N = rNULL__
if(ar__1(604:607) .EQ. '') pr_M_1 = rNULL__
if(ar__1(609:614) .EQ. '') Jmag = rNULL__
if(ar__1(616:622) .EQ. '') e_Jmag = rNULL__
if(ar__1(624:629) .EQ. '') Hmag = rNULL__
if(ar__1(631:637) .EQ. '') e_Hmag = rNULL__
if(ar__1(639:644) .EQ. '') Kmag = rNULL__
if(ar__1(646:650) .EQ. '') e_Kmag = rNULL__
if(ar__1(652:657) .EQ. '') v_3_6_ = rNULL__
if(ar__1(659:663) .EQ. '') e__3_6_ = rNULL__
if(ar__1(665:670) .EQ. '') v_4_5_ = rNULL__
if(ar__1(672:676) .EQ. '') e__4_5_ = rNULL__
if(ar__1(678:683) .EQ. '') v_5_8_ = rNULL__
if(ar__1(685:689) .EQ. '') e__5_8_ = rNULL__
if(ar__1(691:696) .EQ. '') v_8_0_ = rNULL__
if(ar__1(698:703) .EQ. '') e__8_0_ = rNULL__
if(ar__1(914:919) .EQ. '') Av = rNULL__
if(ar__1(925:928) .EQ. '') H1pr = rNULL__
if(ar__1(930:933) .EQ. '') H2pr = rNULL__
if(ar__1(935:938) .EQ. '') H3pr = rNULL__
if(ar__1(940:943) .EQ. '') H4pr = rNULL__
if(ar__1(945:948) .EQ. '') H1ps = rNULL__
if(ar__1(950:953) .EQ. '') H2ps = rNULL__
if(ar__1(955:958) .EQ. '') H3ps = rNULL__
if(ar__1(960:963) .EQ. '') H4ps = rNULL__
c ..............Just test output...........
write(6,2)
+ Reg_1,MPCM,m_MPCM,RAdeg,DEdeg,ePos,rPos,CXO,Xlab,NName,Nlab,
+ MName,Mlab,Olab,Xind,Iind,SpT,r_SpT,OBmag,OBband,pr_0,pr_1,
+ pr_m,t_exp,f_exp,Nobs,Nph,b_theta,theta,B_theta_1,PSFf,aglow,
+ Oct_t,Nct_t,Nct_s,Nct_h,b_Nct_t,B_Nct_t_1,b_Nct_s,B_Nct_s_1,
+ b_Nct_h,B_Nct_h_1,E_t,XFlx_t,XFlx_s,XFlx_h,LX_h,LXc_h,e_LXc_h,
+ s_LXc_h,LX_t,LXc_t,e_LXc_t,s_LXc_t,logNH,e_logNH,s_logNH,pr_N,
+ pr_M_1,Jmag,e_Jmag,Hmag,e_Hmag,Kmag,e_Kmag,v_3_6_,e__3_6_,
+ v_4_5_,e__4_5_,v_5_8_,e__5_8_,v_8_0_,e__8_0_,f_Jmag,f_Hmag,
+ f_Kmag,v2M,q_2M,n_Jmag,n_Hmag,n_Kmag,f__3_6_,f__4_5_,f__5_8_,
+ f__8_0_,LS,r_Jmag,r_Hmag,r_Kmag,r__3_6_,r__4_5_,r__5_8_,
+ r__8_0_,Cl,Av,St,H1pr,H2pr,H3pr,H4pr,H1ps,H2ps,H3ps,H4ps,H2f,
+ Xcl
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,MPCM_1,RAdeg_1,DEdeg_1,ePos_1,rPos_1,CXO_1,Xlab_1,
+ NName_1,Nlab_1,MName_1,Mlab_1,pr_N_1,pr_M_2,Jmag_1,e_Jmag_1,
+ Hmag_1,e_Hmag_1,Kmag_1,e_Kmag_1,v_3_6__1,e__3_6__1,v_4_5__1,
+ e__4_5__1,v_5_8__1,e__5_8__1,v_8_0__1,e__8_0__1,f_Jmag_1,
+ f_Hmag_1,f_Kmag_1,v2M_1,q_2M_1,n_Jmag_1,n_Hmag_1,n_Kmag_1,
+ f__3_6__1,f__4_5__1,f__5_8__1,f__8_0__1,LS_1,r_Jmag_1,
+ r_Hmag_1,r_Kmag_1,r__3_6__1,r__4_5__1,r__5_8__1,r__8_0__1,
+ Cl_1,Av_1,St_1,H1pr_1,H2pr_1,H3pr_1,H4pr_1,H1ps_1,H2ps_1,
+ H3ps_1,H4ps_1,H2f_1,Xcl_1
if(ar__2(189:192) .EQ. '') pr_N_1 = rNULL__
if(ar__2(194:197) .EQ. '') pr_M_2 = rNULL__
if(ar__2(199:204) .EQ. '') Jmag_1 = rNULL__
if(ar__2(206:212) .EQ. '') e_Jmag_1 = rNULL__
if(ar__2(214:219) .EQ. '') Hmag_1 = rNULL__
if(ar__2(221:227) .EQ. '') e_Hmag_1 = rNULL__
if(ar__2(229:234) .EQ. '') Kmag_1 = rNULL__
if(ar__2(236:240) .EQ. '') e_Kmag_1 = rNULL__
if(ar__2(242:247) .EQ. '') v_3_6__1 = rNULL__
if(ar__2(249:253) .EQ. '') e__3_6__1 = rNULL__
if(ar__2(255:260) .EQ. '') v_4_5__1 = rNULL__
if(ar__2(262:266) .EQ. '') e__4_5__1 = rNULL__
if(ar__2(268:273) .EQ. '') v_5_8__1 = rNULL__
if(ar__2(275:279) .EQ. '') e__5_8__1 = rNULL__
if(ar__2(281:286) .EQ. '') v_8_0__1 = rNULL__
if(ar__2(288:295) .EQ. '') e__8_0__1 = rNULL__
if(ar__2(438:443) .EQ. '') Av_1 = rNULL__
if(ar__2(469:472) .EQ. '') H1ps_1 = rNULL__
if(ar__2(474:477) .EQ. '') H2ps_1 = rNULL__
if(ar__2(479:482) .EQ. '') H3ps_1 = rNULL__
if(ar__2(484:487) .EQ. '') H4ps_1 = rNULL__
if(ar__2(491:491) .EQ. '') Xcl_1 = iNULL__
c ..............Just test output...........
write(6,3)
+ Reg_2,MPCM_1,RAdeg_1,DEdeg_1,ePos_1,rPos_1,CXO_1,Xlab_1,
+ NName_1,Nlab_1,MName_1,Mlab_1,pr_N_1,pr_M_2,Jmag_1,e_Jmag_1,
+ Hmag_1,e_Hmag_1,Kmag_1,e_Kmag_1,v_3_6__1,e__3_6__1,v_4_5__1,
+ e__4_5__1,v_5_8__1,e__5_8__1,v_8_0__1,e__8_0__1,f_Jmag_1,
+ f_Hmag_1,f_Kmag_1,v2M_1,q_2M_1,n_Jmag_1,n_Hmag_1,n_Kmag_1,
+ f__3_6__1,f__4_5__1,f__5_8__1,f__8_0__1,LS_1,r_Jmag_1,
+ r_Hmag_1,r_Kmag_1,r__3_6__1,r__4_5__1,r__5_8__1,r__8_0__1,
+ Cl_1,Av_1,St_1,H1pr_1,H2pr_1,H3pr_1,H4pr_1,H1ps_1,H2ps_1,
+ H3ps_1,H4ps_1,H2f_1,Xcl_1
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,BibCode,Aut,Comm
c ..............Just test output...........
write(6,4)Ref,BibCode,Aut,Comm
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