Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/229/28 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/ApJS/229/28 SFiNCs: X-ray, IR and membership catalogs (Getman+, 2017)
*================================================================================
*Star Formation in Nearby Clouds (SFiNCs): X-ray and infrared source catalogs
*and membership.
* Getman K.V., Broos P.S., Kuhn M.A., Feigelson E.D., Richert A.J.W., Ota Y.,
* Bate M.R., Garmire G.P.
* <Astrophys. J. Suppl. Ser., 229, 28-28 (2017)>
* =2017ApJS..229...28G (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Star Formation in Nearby Clouds (SFiNCs)
star-forming regions
integer*4 nr__
parameter (nr__=22) ! Number of records
character*99 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*15 Reg (nr__) ! SFiNCs star forming region name
character*5 OName (nr__) ! Other name
integer*4 RAh (nr__) ! (h) Hour of right ascension (J2000)
integer*4 RAm (nr__) ! (min) Minute of right ascension (J2000)
real*4 RAs (nr__) ! (s) Second of right ascension (J2000)
character*1 DE_ (nr__) ! Sign of declination (J2000)
integer*4 DEd (nr__) ! (deg) Degree of declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Arcminute of declination (J2000)
real*4 DEs (nr__) ! (arcsec) Arcsecond of declination (J2000)
real*8 GLON (nr__) ! (deg) Galactic longitude
real*8 GLAT (nr__) ! (deg) Galactic latitude
integer*4 Dist (nr__) ! (pc) [235/1400] Distance from the Sun
character*9 r_Dist (nr__) ! Reference for Dist (see refs.dat file)
character*5 OB_ (nr__) ! Earliest OB star (1)
character*4 RefIR (nr__) ! Spitzer reference (see refs.dat file)
character*9 RefX (nr__) ! Chandra reference(s) (see refs.dat file)
*Note (1): The earliest OB star within the Chandra field from the catalogs of
* Skiff (Cat. B/mk) and SIMBAD.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Log of SFiNCs Chandra-ACIS observations
integer*4 nr__1
parameter (nr__1=65) ! Number of records
character*99 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*15 Reg_1 (nr__1) ! SFiNCs star forming region name
integer*4 ObsID (nr__1) ! [606/17643] Chandra observation ID
integer*4 Seq (nr__1) ! Chandra sequence number
character*16 Date (nr__1) ! UT start time of Chandra observation
integer*4 Exp (nr__1) ! (s) [1841/97485] Exposure time (1)
integer*4 RAh_1 (nr__1) ! (h) Hour of right ascension (J2000) (2)
integer*4 RAm_1 (nr__1) ! (min) Minute of right ascension (J2000)
real*4 RAs_1 (nr__1) ! (s) Second of right ascension (J2000)
character*1 DE__1 (nr__1) ! Sign of declination (J2000) (2)
integer*4 DEd_1 (nr__1) ! (deg) Degree of declination (J2000) (2)
integer*4 DEm_1 (nr__1) ! (arcmin) Arcminute of declination (J2000)
real*4 DEs_1 (nr__1) ! (arcsec) Arcsecond of declination (J2000)
integer*4 Roll (nr__1) ! (deg) [7/355] Roll angle of Chandra observation
character*6 Mode (nr__1) ! Chandra-ACIS observing mode (3)
character*12 PI (nr__1) ! Principal investigator of a Chandra observation
*Note (1): Exposure time is the net usable time after various filtering steps
* are applied in the data reduction process.
*Note (2): The aimpoint of a Chandra observation is obtained from the satellite
* aspect solution before astrometric correction is applied.
*Note (3): The ACIS modes are described in Section 6.12 of Chandra
* Proposers' Observatory Guide, http://asc.harvard.edu/proposer/POG/.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! SFiNCs X-ray sources and basic properties
integer*4 nr__2
parameter (nr__2=15364) ! Number of records
character*125 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*15 Reg_2 (nr__2) ! Star Formation Region identifier
integer*4 Seq_1 (nr__2) ! [1/2196] X-ray catalog sequence number
* within Reg; RA sorted
character*18 CXOU (nr__2) ! CXOU identifier (HHMMSS.ss+DDMMSS.s; J2000)
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension in Right Ascension (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination in Right Ascension (J2000)
real*4 PosErr (nr__2) ! (arcsec) [0/3.5] Estimated standard deviation
* in position (1)
real*4 theta (nr__2) ! (arcmin) [0/21] Off-axis angle
real*8 NCts (nr__2) ! (ct) [-13/21921] Total band (0.5-8keV) net counts
real*4 e_NCts (nr__2) ! (ct) [1/149]? The 1{sigma} error in NCts (2)
real*4 BCts (nr__2) ! (ct) [0/715] Total band background counts
real*8 HCts (nr__2) ! (ct) [0/10189] Hard band (2-8keV) net counts
real*4 PSF (nr__2) ! [0.3/1] Fraction of PSF (3)
real*4 SNR (nr__2) ! [-0.8/146.5] Photometric significance (4)
character*1 l_logPB (nr__2) ! Limit flag on logPB
real*4 logPB (nr__2) ! ([-]) [-5/-1.3] Log probability (5)
character*3 Anom (nr__2) ! Source anomalie(s) (6)
character*1 Var (nr__2) ! [a-c] Variability characterization (7)
real*4 EMed (nr__2) ! (keV) [0.5/7.7]? Background-corrected total band
* median photon energy
*Note (1): Estimated standard deviation of the random component of the position
* error, ({sigma}_x_^2^+{sigma}_y_^2^)^0.5^. 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 (2): Average of the upper and lower 1{sigma} errors.
*Note (3): (at 1.497keV) enclosed within the extraction region. A reduced PSF
* fraction (significantly below 90%) may indicate that the source
* is in a crowded region.
*Note (4): Computed as net counts divided by the upper error on net counts.
*Note (5): that extracted counts (total band) are solely from background. Some
* sources have PB values above the 1% threshold that defines the
* catalog because local background estimates can rise during the final
* extraction iteration after sources are removed from the catalog.
*Note (6): Anomalies as follows:
* g = fractional time that source was on a detector
* (FRACEXPO from mkarf) is <0.9;
* e = source on field edge;
* p = X-ray properties may be biased due to photon pile-up;
* s = source on readout streak;
* a = photometry and spectrum may contain >10% afterglow events.
*Note (7): Based on K-S statistic (total band) from the single ObsId showing
* the most variability. Variabilitiy 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
* with fewer than four counts or for sources in chip gaps or on
* field edges.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! SFiNCs X-ray fluxes
integer*4 nr__3
parameter (nr__3=15364) ! Number of records
character*138 ar__3 ! Full-size record
character*15 Reg_3 (nr__3) ! Star Formation Region identifier
integer*4 Seq_2 (nr__3) ! [1/2196] X-ray catalog sequence number
* within Reg; RA sorted
character*18 CXOU_1 (nr__3) ! CXOU identifier (HHMMSS.ss+DDMMSS.s; J2000)
real*8 logPFh (nr__3) ! ([ct/cm2/s]) [-11/-3]? AE log hard band (2-8keV)
* photon flux (1)
real*4 logPFt (nr__3) ! ([ct/cm2/s]) [-8/-2]? AE log total band (0.5-8keV)
* photon flux (1)
real*8 logFh (nr__3) ! ([mW/m2]) [-16/-11]? XPHOT log apparent hard band
* flux (2)
real*4 e_logFh (nr__3) ! ([mW/m2]) [0.01/2.2]? The 1{sigma} error in logFh
real*8 logFt (nr__3) ! ([mW/m2]) [-16/-11]? XPHOT log apparent total band
* flux (2)
real*4 e_logFt (nr__3) ! ([mW/m2]) [0.007/1.4]? The 1{sigma} error in logFt
real*4 logNH (nr__3) ! ([cm-2]) [20/24]? XPHOT log Hydrogen column density (2)
real*4 e_logNH (nr__3) ! ([cm-2]) [0.01/1.1]? logNH 1{sigma} statistical error
real*4 s_logNH (nr__3) ! ([cm-2]) [0.02/1.2]? logNH 1{sigma} systematic error
real*8 logFhc (nr__3) ! ([mW/m2]) [-16/-11]? XPHOT log corrected hard band
* flux (2)
real*4 e_logFhc (nr__3) ! ([mW/m2]) [0.01/2.2]? logFhc 1{sigma} statistical error
real*4 s_logFhc (nr__3) ! ([mW/m2]) [0/2.2]? logFhc 1{sigma} systematic error
real*8 logFtc (nr__3) ! ([mW/m2]) [-16/-11]? XPHOT log corrected hard band
* flux (2)
real*4 e_logFtc (nr__3) ! ([mW/m2]) [0.007/1.4]? logFtc 1{sigma} statistical error
real*4 s_logFtc (nr__3) ! ([mW/m2]) [0.001/5.1]? logFtc 1{sigma} systematic error
*Note (1): Produced by AE (Broos et al. 2010ApJ...714.1582B).
*Note (2): Produced by XPHOT (Getman+ 2010, J/ApJ/708/1760). XPHOT
* assumes X-ray spectral shapes of young, low-mass stars. Intrinsic
* XPHOT quantities will be unreliable for high-mass stellar members of
* the SFiNCs SFRs as well as for non-members, such as Galactic field
* stars and extragalactic objects.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Log of SFiNCs Spitzer-IRAC observations
integer*4 nr__4
parameter (nr__4=437) ! Number of records
character*136 ar__4 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*15 Reg_4 (nr__4) ! Star Formation Region identifier
character*11 Reg2 (nr__4) ! Overlapping Star Formation Region identifier
integer*4 AOR (nr__4) ! Spitzer Astronomcail Object Request number
integer*4 PID (nr__4) ! Spitzer program identification number
character*10 st_date (nr__4) ! ("Y/M/D") UT date of observation start
character*8 st_time (nr__4) ! ("h:m:s") Time of observation start
character*10 end_date (nr__4) ! ("Y/M/D") UT date of observation end
character*8 end_time (nr__4) ! ("h:m:s") Time of observation end
real*8 RAdeg_3 (nr__4) ! (deg) Right Ascension in Right Ascension (J2000) (1)
real*8 DEdeg_3 (nr__4) ! (deg) Declination in Right Ascension (J2000) (1)
character*11 Mode_1 (nr__4) ! IRAC mode (2)
character*20 PI_1 (nr__4) ! Principal investigator
*Note (1): Approximate center of the observation.
*Note (2): The PC (Post-Cryo) mode has been introduced during the warm mission
* of the Spitzer observatory with only two shortest-wavelength IRAC
* modules in operation.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! SFiNCs IRAC sources and photometry
integer*4 nr__5
parameter (nr__5=1638654) ! Number of records
character*123 ar__5 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*15 Reg_5 (nr__5) ! Star Formation Region identifier
character*11 Reg2_1 (nr__5) ! Overlapping Star Formation Region identifier
character*17 IRAC (nr__5) ! Source identifier (GLLL.llll+BB.bbbb)
real*8 RAdeg_4 (nr__5) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_4 (nr__5) ! (deg) Declination in decimal degrees (J2000)
real*4 v3_6mag (nr__5) ! (mag) [4.7/20.4] Spitzer/IRAC 3.6um band magnitude
real*4 e_3_6mag (nr__5) ! (mag) [0.001/0.3] The 1{sigma} uncertainty in 3.6mag
real*4 v4_5mag (nr__5) ! (mag) [4.2/20] Spitzer/IRAC 4.5 micron band magnitude
real*4 e_4_5mag (nr__5) ! (mag) [0.001/0.3] The 1{sigma} uncertainty in 4.5mag
real*4 v5_8mag (nr__5) ! (mag) [2.3/17.5]? Spitzer/IRAC 5.8um band magnitude
real*4 e_5_8mag (nr__5) ! (mag) [0.001/0.3]? The 1{sigma} uncertainty in 5.8mag
real*4 v8_0mag (nr__5) ! (mag) [2.4/17.5]? Spitzer/IRAC 8.0um band magnitude
real*4 e_8_0mag (nr__5) ! (mag) [0.001/0.3]? The 1{sigma} uncertainty in 8.0mag
integer*4 Ap (nr__5) ! [-1/4] Aperture size flag (1)
*Note (1): 2, 3, and 4 = 2, 3, and 4-pixel aperture with contaminating
* flux <10% of source flux; -1 = 2-pixel aperture with contaminating
* flux >10% of source flux.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! SFiNCs probable cluster members: IR photometry
integer*4 nr__6
parameter (nr__6=8492) ! Number of records
character*176 ar__6 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*15 Reg_6 (nr__6) ! Star Formation Region identifier
character*18 Source (nr__6) ! Source identifier (HHMMSS.ss+DDMMSS.s; J2000)
real*8 RAdeg_5 (nr__6) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_5 (nr__6) ! (deg) Declination in decimal degrees (J2000)
character*17 v2MASS (nr__6) ! 2MASS identifier (HHMMSSss+DDMMSSs; J2000)
real*4 Jmag (nr__6) ! (mag) [5.3/19]? 2MASS J band magnitude
real*4 e_Jmag (nr__6) ! (mag) [0.02/0.8]? Uncertainty in Jmag
real*4 Hmag (nr__6) ! (mag) [5/18]? 2MASS H band magnitude
real*4 e_Hmag (nr__6) ! (mag) [0.01/0.7]? Uncertainty in Hmag
real*4 Ksmag (nr__6) ! (mag) [3/16.5]? 2MASS Ks band magnitude
real*4 e_Ksmag (nr__6) ! (mag) [0.01/9.99]? Uncertainty in Kmag
character*6 Flag (nr__6) ! 2MASS flags (1)
character*17 IRAC_1 (nr__6) ! IRAC source identifier (GLLL.llll+BB.bbbb)
real*4 v3_6mag_1 (nr__6) ! (mag) [4.9/19]? Spitzer/IRAC 3.6 micron band magnitude
real*4 e_3_6mag_1 (nr__6) ! (mag) [0/0.3]? The 1{sigma} uncertainty in 3.6mag
real*4 v4_5mag_1 (nr__6) ! (mag) [4/17.1]? Spitzer/IRAC 4.5 micron band magnitude
real*4 e_4_5mag_1 (nr__6) ! (mag) [0/0.3]? The 1{sigma} uncertainty in 4.5mag
real*4 v5_8mag_1 (nr__6) ! (mag) [2.5/15.7]? Spitzer/IRAC 5.8um band magnitude
real*4 e_5_8mag_1 (nr__6) ! (mag) [0/0.3]? The 1{sigma} uncertainty in 5.8mag
real*4 v8_0mag_1 (nr__6) ! (mag) [2.4/15]? Spitzer/IRAC 8.0 micron band magnitude
real*4 e_8_0mag_1 (nr__6) ! (mag) [0.01/0.3]? The 1{sigma} uncertainty in 8.0mag
*Note (1): The first 3 bytes are the JHK photometry quality codes. See
* http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec1_6b.html#phqual
* for details. The next 3 bytes are the JHK confusion/contamination codes.
* See http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec4_7.html
* for details.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! SFiNCs probable cluster members: main X-ray and
other properties
integer*4 nr__7
parameter (nr__7=8492) ! Number of records
character*262 ar__7 ! Full-size record
character*15 Reg_7 (nr__7) ! Star Formation Region identifier
character*18 Source_1 (nr__7) ! Source identifier
* (HHMMSS.ss+DDMMSS.s; J2000)
character*5 ACIS (nr__7) ! ACIS source label
character*1 XFOV (nr__7) ! [0/1] XFOV flag (1)
real*8 NCts_1 (nr__7) ! (ct) [1/21921]? Total band (0.5-8keV) net counts
real*4 EMed_1 (nr__7) ! (keV) [0.6/6.7]? Background-corrected total band
* median photon energy
real*4 logPFt_1 (nr__7) ! ([ct/cm2/s]) [-7.6/-2.9]? AE log total band photon flux
real*4 logNH_1 (nr__7) ! ([cm-2]) [20/24]? XPHOT log Hydrogen column density
real*4 logLtc (nr__7) ! ([10-7W]) [27.6/32.6]? log total band luminosity
real*4 alpha (nr__7) ! [-5.2/5.5]? Apparent SED slope (2)
real*4 e_alpha (nr__7) ! [0.01/1.2]? Error in alpha
character*7 Class (nr__7) ! YSO class(es) (3)
real*4 Av (nr__7) ! (mag) [0/34.5]? The V band extinction (4)
real*4 AgeJX (nr__7) ! [0.1/5]? Stellar age (4)
character*44 Note (nr__7) ! Know OB-type star note (5)
character*116 X_IR (nr__7) ! X-ray-IR source match note (6)
*Note (1): A flag indicating whether the source is located inside (XFOV= 1)
* or outside (XFOV=0) the Chandra ACIS-I field of view.
*Note (2): Measured in the IRAC wavelength range from 3.6 to 8.0 microns.
*Note (3): Class as follows:
* NOD = diskless (2532 occurrences);
* DSK = disky (5628 occurrences);
* PMB = possible member without a disk class (303 occurrences);
* FRG = source could be a foreground star or a YSO member of the region
* (29 occurrences).
*Note (4): Using the methods of Getman et al. (2014ApJ...787..108G).
*Note (5): The related information on the OB star' name, spectral type, and
* primary catalog of origin.
*Note (6): A note indicates that the X-ray-IR source match was rejected by the
* tool match xy but was reinstated as a legitimate match upon the
* visual inspection. Such cases are often associated with source'
* binarity and multiplicity. Related comments are provided.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! Literature references to the previous Spitzer
and Chandra studies of the region
integer*4 nr__8
parameter (nr__8=38) ! Number of records
character*69 ar__8 ! Full-size record
character*4 Ref (nr__8) ! Reference code
character*19 BibCode (nr__8) ! Bibcode of the reference
character*22 Auth (nr__8) ! First author's name
character*21 Comm (nr__8) ! VizieR catalog reference
C=============================================================================
C Loading file 'table1.dat' ! Star Formation in Nearby Clouds (SFiNCs)
* star-forming regions
C Format for file interpretation
1 format(
+ A15,1X,A5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F8.4,
+ 1X,F8.4,1X,I4,1X,A9,1X,A5,1X,A4,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,22
read(1,'(A99)')ar__
read(ar__,1)
+ Reg(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Dist(i__),
+ r_Dist(i__),OB_(i__),RefIR(i__),RefX(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Reg(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Dist(i__),
+ r_Dist(i__),OB_(i__),RefIR(i__),RefX(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Log of SFiNCs Chandra-ACIS observations
C Format for file interpretation
2 format(
+ A15,1X,I5,1X,I6,1X,A16,1X,I5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,
+ I2,1X,F4.1,1X,I3,1X,A6,1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,65
read(1,'(A99)')ar__1
read(ar__1,2)
+ Reg_1(i__),ObsID(i__),Seq(i__),Date(i__),Exp(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),Roll(i__),Mode(i__),PI(i__)
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ Reg_1(i__),ObsID(i__),Seq(i__),Date(i__),Exp(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),Roll(i__),Mode(i__),PI(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! SFiNCs X-ray sources and basic properties
C Format for file interpretation
3 format(
+ A15,1X,I4,1X,A18,1X,F10.6,1X,F10.6,1X,F3.1,1X,F4.1,1X,F7.1,1X,
+ F5.1,1X,F5.1,1X,F7.1,1X,F4.2,1X,F5.1,1X,A1,F4.1,1X,A3,1X,A1,
+ 1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,15364
read(1,'(A125)')ar__2
read(ar__2,3)
+ Reg_2(i__),Seq_1(i__),CXOU(i__),RAdeg_2(i__),DEdeg_2(i__),
+ PosErr(i__),theta(i__),NCts(i__),e_NCts(i__),BCts(i__),
+ HCts(i__),PSF(i__),SNR(i__),l_logPB(i__),logPB(i__),Anom(i__),
+ Var(i__),EMed(i__)
if(ar__2(80:84) .EQ. '') e_NCts(i__) = rNULL__
if(ar__2(123:125) .EQ. '') EMed(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Reg_2(i__),Seq_1(i__),CXOU(i__),RAdeg_2(i__),DEdeg_2(i__),
+ PosErr(i__),theta(i__),NCts(i__),e_NCts(i__),BCts(i__),
+ HCts(i__),PSF(i__),SNR(i__),l_logPB(i__),logPB(i__),Anom(i__),
+ Var(i__),EMed(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! SFiNCs X-ray fluxes
C Format for file interpretation
4 format(
+ A15,1X,I4,1X,A18,1X,F7.3,1X,F6.3,1X,F7.3,1X,F5.3,1X,F7.3,1X,
+ F5.3,1X,F5.2,1X,F4.2,1X,F4.2,1X,F7.3,1X,F5.3,1X,F5.3,1X,F7.3,
+ 1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,15364
read(1,'(A138)')ar__3
read(ar__3,4)
+ Reg_3(i__),Seq_2(i__),CXOU_1(i__),logPFh(i__),logPFt(i__),
+ logFh(i__),e_logFh(i__),logFt(i__),e_logFt(i__),logNH(i__),
+ e_logNH(i__),s_logNH(i__),logFhc(i__),e_logFhc(i__),
+ s_logFhc(i__),logFtc(i__),e_logFtc(i__),s_logFtc(i__)
if(ar__3(41:47) .EQ. '') logPFh(i__) = rNULL__
if(ar__3(49:54) .EQ. '') logPFt(i__) = rNULL__
if(ar__3(56:62) .EQ. '') logFh(i__) = rNULL__
if(ar__3(64:68) .EQ. '') e_logFh(i__) = rNULL__
if(ar__3(70:76) .EQ. '') logFt(i__) = rNULL__
if(ar__3(78:82) .EQ. '') e_logFt(i__) = rNULL__
if(ar__3(84:88) .EQ. '') logNH(i__) = rNULL__
if(ar__3(90:93) .EQ. '') e_logNH(i__) = rNULL__
if(ar__3(95:98) .EQ. '') s_logNH(i__) = rNULL__
if(ar__3(100:106) .EQ. '') logFhc(i__) = rNULL__
if(ar__3(108:112) .EQ. '') e_logFhc(i__) = rNULL__
if(ar__3(114:118) .EQ. '') s_logFhc(i__) = rNULL__
if(ar__3(120:126) .EQ. '') logFtc(i__) = rNULL__
if(ar__3(128:132) .EQ. '') e_logFtc(i__) = rNULL__
if(ar__3(134:138) .EQ. '') s_logFtc(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Reg_3(i__),Seq_2(i__),CXOU_1(i__),logPFh(i__),logPFt(i__),
+ logFh(i__),e_logFh(i__),logFt(i__),e_logFt(i__),logNH(i__),
+ e_logNH(i__),s_logNH(i__),logFhc(i__),e_logFhc(i__),
+ s_logFhc(i__),logFtc(i__),e_logFtc(i__),s_logFtc(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Log of SFiNCs Spitzer-IRAC observations
C Format for file interpretation
5 format(
+ A15,A11,1X,I8,1X,I5,1X,A10,1X,A8,1X,A10,1X,A8,1X,F10.6,1X,
+ F10.6,1X,A11,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,437
read(1,'(A136)')ar__4
read(ar__4,5)
+ Reg_4(i__),Reg2(i__),AOR(i__),PID(i__),st_date(i__),
+ st_time(i__),end_date(i__),end_time(i__),RAdeg_3(i__),
+ DEdeg_3(i__),Mode_1(i__),PI_1(i__)
c ..............Just test output...........
write(6,5)
+ Reg_4(i__),Reg2(i__),AOR(i__),PID(i__),st_date(i__),
+ st_time(i__),end_date(i__),end_time(i__),RAdeg_3(i__),
+ DEdeg_3(i__),Mode_1(i__),PI_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! SFiNCs IRAC sources and photometry
C Format for file interpretation
6 format(
+ A15,A11,1X,A17,1X,F11.7,1X,F11.7,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,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,1638654
read(1,'(A123)')ar__5
read(ar__5,6)
+ Reg_5(i__),Reg2_1(i__),IRAC(i__),RAdeg_4(i__),DEdeg_4(i__),
+ v3_6mag(i__),e_3_6mag(i__),v4_5mag(i__),e_4_5mag(i__),
+ v5_8mag(i__),e_5_8mag(i__),v8_0mag(i__),e_8_0mag(i__),Ap(i__)
if(ar__5(96:101) .EQ. '') v5_8mag(i__) = rNULL__
if(ar__5(103:107) .EQ. '') e_5_8mag(i__) = rNULL__
if(ar__5(109:114) .EQ. '') v8_0mag(i__) = rNULL__
if(ar__5(116:120) .EQ. '') e_8_0mag(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ Reg_5(i__),Reg2_1(i__),IRAC(i__),RAdeg_4(i__),DEdeg_4(i__),
+ v3_6mag(i__),e_3_6mag(i__),v4_5mag(i__),e_4_5mag(i__),
+ v5_8mag(i__),e_5_8mag(i__),v8_0mag(i__),e_8_0mag(i__),Ap(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! SFiNCs probable cluster members: IR photometry
C Format for file interpretation
7 format(
+ A15,1X,A18,1X,F10.6,1X,F10.6,1X,A17,1X,F5.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.2,1X,F4.2,1X,A6,1X,A17,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=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,8492
read(1,'(A176)')ar__6
read(ar__6,7)
+ Reg_6(i__),Source(i__),RAdeg_5(i__),DEdeg_5(i__),v2MASS(i__),
+ Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),Ksmag(i__),
+ e_Ksmag(i__),Flag(i__),IRAC_1(i__),v3_6mag_1(i__),
+ e_3_6mag_1(i__),v4_5mag_1(i__),e_4_5mag_1(i__),v5_8mag_1(i__),
+ e_5_8mag_1(i__),v8_0mag_1(i__),e_8_0mag_1(i__)
if(ar__6(76:80) .EQ. '') Jmag(i__) = rNULL__
if(ar__6(82:85) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__6(87:91) .EQ. '') Hmag(i__) = rNULL__
if(ar__6(93:96) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__6(98:102) .EQ. '') Ksmag(i__) = rNULL__
if(ar__6(104:107) .EQ. '') e_Ksmag(i__) = rNULL__
if(ar__6(134:138) .EQ. '') v3_6mag_1(i__) = rNULL__
if(ar__6(140:143) .EQ. '') e_3_6mag_1(i__) = rNULL__
if(ar__6(145:149) .EQ. '') v4_5mag_1(i__) = rNULL__
if(ar__6(151:154) .EQ. '') e_4_5mag_1(i__) = rNULL__
if(ar__6(156:160) .EQ. '') v5_8mag_1(i__) = rNULL__
if(ar__6(162:165) .EQ. '') e_5_8mag_1(i__) = rNULL__
if(ar__6(167:171) .EQ. '') v8_0mag_1(i__) = rNULL__
if(ar__6(173:176) .EQ. '') e_8_0mag_1(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ Reg_6(i__),Source(i__),RAdeg_5(i__),DEdeg_5(i__),v2MASS(i__),
+ Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),Ksmag(i__),
+ e_Ksmag(i__),Flag(i__),IRAC_1(i__),v3_6mag_1(i__),
+ e_3_6mag_1(i__),v4_5mag_1(i__),e_4_5mag_1(i__),v5_8mag_1(i__),
+ e_5_8mag_1(i__),v8_0mag_1(i__),e_8_0mag_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! SFiNCs probable cluster members: main X-ray and
* other properties
C Format for file interpretation
8 format(
+ A15,1X,A18,1X,A5,1X,A1,1X,F7.1,1X,F4.2,1X,F5.2,1X,F4.1,1X,
+ F5.2,1X,F5.2,1X,F4.2,1X,A7,1X,F4.1,1X,F3.1,1X,A44,1X,A116)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,8492
read(1,'(A262)')ar__7
read(ar__7,8)
+ Reg_7(i__),Source_1(i__),ACIS(i__),XFOV(i__),NCts_1(i__),
+ EMed_1(i__),logPFt_1(i__),logNH_1(i__),logLtc(i__),alpha(i__),
+ e_alpha(i__),Class(i__),Av(i__),AgeJX(i__),Note(i__),X_IR(i__)
if(ar__7(44:50) .EQ. '') NCts_1(i__) = rNULL__
if(ar__7(52:55) .EQ. '') EMed_1(i__) = rNULL__
if(ar__7(57:61) .EQ. '') logPFt_1(i__) = rNULL__
if(ar__7(63:66) .EQ. '') logNH_1(i__) = rNULL__
if(ar__7(68:72) .EQ. '') logLtc(i__) = rNULL__
if(ar__7(74:78) .EQ. '') alpha(i__) = rNULL__
if(ar__7(80:83) .EQ. '') e_alpha(i__) = rNULL__
if(ar__7(93:96) .EQ. '') Av(i__) = rNULL__
if(ar__7(98:100) .EQ. '') AgeJX(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ Reg_7(i__),Source_1(i__),ACIS(i__),XFOV(i__),NCts_1(i__),
+ EMed_1(i__),logPFt_1(i__),logNH_1(i__),logLtc(i__),alpha(i__),
+ e_alpha(i__),Class(i__),Av(i__),AgeJX(i__),Note(i__),X_IR(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! Literature references to the previous Spitzer
* and Chandra studies of the region
C Format for file interpretation
9 format(A4,1X,A19,1X,A22,1X,A21)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,38
read(1,'(A69)')ar__8
read(ar__8,9)Ref(i__),BibCode(i__),Auth(i__),Comm(i__)
c ..............Just test output...........
write(6,9)Ref(i__),BibCode(i__),Auth(i__),Comm(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end