Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/160/68 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-13
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/AJ/160/68 GLIMPSE360: stellar objects photometry and classes (Winston+, 2020)
*================================================================================
*A Census of Star Formation in the Outer Galaxy. II. The GLIMPSE360 Field.
* Winston E., Hora J.L., Tolls V.
* <Astron. J., 160, 68 (2020)>
* =2020AJ....160...68W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Photometry of YSOs in Star formation Field in
the Outer Galaxy (SFOG)
integer*4 nr__
parameter (nr__=47405) ! Number of records
character*1383 ar__ ! Full-size record
C J2000 position composed of: RA2Mdeg DE2Mdeg
character*26 Source (nr__) ! Source designation (designation) (1)
character*16 v2MASS (nr__) ! 2MASS-PSC source identifier
* (HHMMSSss+DDMMSSs; J2000) (II/246)
* (tmass_designation)
integer*8 v2MASS_cntr(nr__) ! [68911969/1326367653]? 2MASS-PSC
* unique source identifier (II/246)
* (tmass_cntr)
real*8 GLON_2M (nr__) ! (deg) ? 2MASS-PSC galactic longitude (II/246)
* (l)
real*8 GLAT_2M (nr__) ! (deg) [-3.1/3.4]? 2MASS-PSC galactic
* latitude (II/246) (b)
real*4 e_GLON_2M (nr__) ! (arcsec) [0.3]? Uncertainty in GLON_2M (dl)
real*4 e_GLAT_2M (nr__) ! (arcsec) [0.3]? Uncertainty in GLAT_2M (db)
real*8 RA2Mdeg (nr__) ! (deg) ? 2MASS-PSC right ascension (J2000) in
* degrees (II/246) (ra)
real*8 DE2Mdeg (nr__) ! (deg) [-47/66]? 2MASS-PSC declination
* (J2000) in degrees (II/246) (dec)
real*4 e_RA2Mdeg (nr__) ! (arcsec) [0.3]? Uncertainty in RA2Mdeg (dra)
real*4 e_DE2Mdeg (nr__) ! (arcsec) [0.3]? Uncertainty in DE2Mdeg (ddec)
integer*4 csf (nr__) ! [0/5]? Close source flag (csf) (2)
real*4 Jmag (nr__) ! (mag) [6.6/18.2]? 2MASS-PSC J-band (1.25um)
* magnitude (II/246) (mag_J)
real*4 e_Jmag (nr__) ! (mag) [0.01/1.1]? 1{sigma} uncertainty in
* Jmag (dJ_m)
real*4 Hmag (nr__) ! (mag) [6.6/16.7]? 2MASS-PSC H-band (1.65um)
* magnitude (II/246) (mag_H)
real*4 e_Hmag (nr__) ! (mag) [0.01/0.3]? 1{sigma} uncertainty in
* Hmag (dH_m)
real*4 Kmag (nr__) ! (mag) [6.2/16.4]? 2MASS-PSC Ks-band (2.17um)
* magnitude (II/246) (mag_Ks)
real*4 e_Kmag (nr__) ! (mag) [0.01/10]? 1{sigma} uncertainty in
* Kmag (dKs_m)
real*4 v3_6mag (nr__) ! (mag) [3.6/17.5]? Spitzer/IRAC 3.6um (Band
* 1) magnitude (J/ApJ/880/9, II/293)
* (mag3_6)
real*4 e_3_6mag (nr__) ! (mag) [0.01/0.22]? 1{sigma} uncertainty in
* 3.6mag (d3_6m)
real*4 v4_5mag (nr__) ! (mag) [3.6/16.7]? Spitzer/IRAC 4.5um (Band
* 2) magnitude (J/ApJ/880/9, II/293)
* (mag4_5)
real*4 e_4_5mag (nr__) ! (mag) [0.01/0.21]? 1{sigma} uncertainty in
* 4.5mag (d4_5m)
real*4 v5_8mag (nr__) ! (mag) [3.3/15]? Spitzer/IRAC 5.8um (Band 3)
* magnitude (J/ApJ/880/9, II/293)
* (mag5_8)
real*4 e_5_8mag (nr__) ! (mag) [0.01/0.21]? 1{sigma} uncertainty in
* 5.8mag (d5_8m)
real*4 v8_0mag (nr__) ! (mag) [3.1/14.1]? Spitzer/IRAC 8.0um (Band
* 4) magnitude (J/ApJ/880/9, II/293)
* (mag8_0)
real*4 e_8_0mag (nr__) ! (mag) [0.01/0.22]? 1{sigma} uncertainty in
* 8.0mag (d8_0m)
real*8 F_J (nr__) ! (mJy) [0.08/3346]? 2MASS-PSC J-band (1.25um)
* flux (II/246, II/293) (f_J)
real*8 e_F_J (nr__) ! (mJy) [0.02/71]? 1{sigma} uncertainty in
* F(J) (df_J)
real*8 F_H (nr__) ! (mJy) [0.2/2206]? 2MASS-PSC H-band (1.65um)
* flux (II/246, II/293) (f_H)
real*8 e_F_H (nr__) ! (mJy) [0.03/75.3]? 1{sigma} uncertainty in
* F(H) (df_H)
real*8 F_K (nr__) ! (mJy) [0.19/2130]? 2MASS-PSC Ks-band
* (2.17um) flux (II/246, II/293) (f_Ks)
real*8 e_F_K (nr__) ! (mJy) [0.04/5181]? 1{sigma} uncertainty in
* F(K) (df_Ks)
real*8 F_3_6 (nr__) ! (mJy) [0.029/10050]? Spitzer/IRAC 3.6um
* (Band 1) flux (II/293) (f3_6)
real*8 e_F_3_6 (nr__) ! (mJy) [3e-3/1319]? 1{sigma} uncertainty in
* F(3.6) (df3_6)
real*8 F_4_5 (nr__) ! (mJy) [0.03/6195]? Spitzer/IRAC 4.5um
* (Band 2) flux (II/293) (f4_5)
real*8 e_F_4_5 (nr__) ! (mJy) [3e-3/279]? 1{sigma} uncertainty in
* F(4.5) (df4_5)
real*8 F_5_8 (nr__) ! (mJy) [0.1/5421]? Spitzer/IRAC 5.8um
* (Band 3) flux (II/293) (f5_8)
real*8 e_F_5_8 (nr__) ! (mJy) [0.01/125]? 1{sigma} uncertainty in
* F(5.8) (df5_8)
real*8 F_8_0 (nr__) ! (mJy) [0.1/3580]? Spitzer/IRAC 8.0um
* (Band 4) flux (II/293) (f8_0)
real*8 e_F_8_0 (nr__) ! (mJy) [0.01/103]? 1{sigma} uncertainty in
* F(8.0) (df8_0)
real*4 rms_3_6 (nr__) ! (mJy) [0/860]? Rms deviation of individual
* detections from final Spitzer/IRAC
* 3.6um (Band 1) flux (II/293)
* (rms\_f3\_6)
integer*4 rms_4_5 (nr__) ! (mJy) [0/509]? Rms deviation of individual
* detections from final Spitzer/IRAC
* 4.5um (Band 2) flux (II/293)
* (rms\_f4\_5)
real*4 rms_5_8 (nr__) ! (mJy) [0/233]? Rms deviation of individual
* detections from final Spitzer/IRAC
* 5.8um (Band 3) flux (II/293)
* (rms\_f5\_8)
real*4 rms_8_0 (nr__) ! (mJy) [0/96.4]? Rms deviation of individual
* detections from final Spitzer/IRAC
* 8.0um (Band 4) flux (II/293)
* (rms\_f8\_0)
real*8 sky_3_6 (nr__) ! (MJy/sr) [-16/55.8]? Local sky background for
* Spitzer/IRAC 3.6um (Band 1) (II/293)
* (sky\_3\_6)
real*4 sky_4_5 (nr__) ! (MJy/sr) [-1.3/54.1]? Local sky background for
* Spitzer/IRAC 4.5um (Band 2) (II/293)
* (sky\_4\_5)
real*8 sky_5_8 (nr__) ! (MJy/sr) [-0.7/81.1]? Local sky background for
* Spitzer/IRAC 5.8um (Band 3) (II/293)
* (sky\_5\_8)
real*8 sky_8_0 (nr__) ! (MJy/sr) [2.3/134]? Local sky background for
* Spitzer/IRAC 8.0um (Band 4) (II/293)
* (sky\_8\_0)
real*4 Jsn (nr__) ! [1/63.9]? 2MASS-PSC J-band (1.25um)
* signal-to-noise ratio (II/246) (sn_J)
real*4 Hsn (nr__) ! [3.7/77.6]? 2MASS-PSC H-band (1.65um)
* signal-to-noise ratio (II/246) (sn_H)
real*4 Ksn (nr__) ! [0.1/83.6]? 2MASS-PSC Ks-band (2.17um)
* signal-to-noise ratio (II/246) (sn_Ks)
real*4 v3_6sn (nr__) ! [5/71.8]? Spitzer/IRAC 3.6um (Band 1)
* signal-to-noise ratio (II/293)
* (sn\_3\_6)
real*4 v4_5sn (nr__) ! [5.3/87.9]? Spitzer/IRAC 4.5um (Band
* 2) signal-to-noise ratio (II/293)
* (sn\_4\_5)
real*4 v5_8sn (nr__) ! [5.2/76.2]? Spitzer/IRAC 5.8um (Band
* 3) signal-to-noise ratio (II/293)
* (sn\_5\_8)
real*4 v8_0sn (nr__) ! [5.1/86]? Spitzer/IRAC 8.0um (Band 4)
* signal-to-noise ratio (II/293)
* (sn\_8\_0)
real*4 d_3_6 (nr__) ! (arcmin-2) [4.1/142]? Local source density for
* Spitzer/IRAC 3.6um (Band 1) (II/293)
* (dens\_3\_6)
real*4 d_4_5 (nr__) ! (arcmin-2) [3.4/115]? Local source density for
* Spitzer/IRAC 4.5um (Band 2) (II/293)
* (dens\_4\_5)
real*4 d_5_8 (nr__) ! (arcmin-2) [3.4/32.3]? Local source density for
* Spitzer/IRAC 5.8um (Band 3) (II/293)
* (dens\_5\_8)
real*4 d_8_0 (nr__) ! (arcmin-2) [2.1/25.1]? Local source density for
* Spitzer/IRAC 8.0um (Band 4) (II/293)
* (dens\_8\_0)
integer*4 N_3_6 (nr__) ! [0/17]? Number of detections for
* Spitzer/IRAC 3.6um (Band 1) (II/293)
* (m3_6)
integer*4 N_4_5 (nr__) ! [0/15]? Number of detections for
* Spitzer/IRAC 4.5um (Band 2) (II/293)
* (m4_5)
integer*4 N_5_8 (nr__) ! [0/16]? Number of detections for
* Spitzer/IRAC 5.8um (Band 3) (II/293)
* (m5_8)
integer*4 N_8_0 (nr__) ! [0/20]? Number of detections for
* Spitzer/IRAC 8.0um (Band 4) (II/293)
* (m8_0)
integer*4 n_3_6_1 (nr__) ! [0/24]? Number of possible detections
* for Spitzer/IRAC 3.6um (Band 1)
* (II/293) (n3_6)
integer*4 n_4_5_1 (nr__) ! [1/16]? Number of possible detections
* for Spitzer/IRAC 4.5um (Band 2)
* (II/293) (n4_5)
integer*4 n_5_8_1 (nr__) ! [0/18]? Number of possible detections
* for Spitzer/IRAC 5.8um (Band 3)
* (II/293) (n5_8)
integer*4 n_8_0_1 (nr__) ! [0/20]? Number of possible detections
* for Spitzer/IRAC 8.0um (Band 4)
* (II/293) (n8_0)
integer*4 q_Jmag (nr__) ! ? 2MASS-PSC J-band (1.25um) source
* quality flag (II/246, II/293) (sqf_J)
* (3)
integer*4 q_Hmag (nr__) ! ? 2MASS-PSC H-band (1.65um) source
* quality flag (II/246, II/293) (sqf_H)
* (3)
integer*4 q_Kmag (nr__) ! ? 2MASS-PSC Ks-band (2.17um) source
* quality flag (II/246, II/293) (sqf_Ks)
* (3)
integer*8 q_3_6mag (nr__) ! ? Spitzer/IRAC 3.6um (Band 1) source
* quality flag (II/293) (sqf\_3\_6) (3)
integer*8 q_4_5mag (nr__) ! ? Spitzer/IRAC 4.5um (Band 2) source
* quality flag (II/293) (sqf\_4\_5) (3)
integer*8 q_5_8mag (nr__) ! ? Spitzer/IRAC 5.8um (Band 3) source
* quality flag (II/293) (sqf\_5\_8) (3)
integer*8 q_8_0mag (nr__) ! ? Spitzer/IRAC 8.0um (Band 4) source
* quality flag (II/293) (sqf\_8\_0) (3)
integer*8 m_3_6 (nr__) ! ? Spitzer/IRAC 3.6um (Band 1) flux
* calculation method flag (II/293)
* (mf3_6)
integer*8 m_4_5 (nr__) ! ? Spitzer/IRAC 4.5um (Band 2) flux
* calculation method flag (II/293)
* (mf4_5)
integer*8 m_5_8 (nr__) ! ? Spitzer/IRAC 5.8um (Band 3) flux
* calculation method flag (II/293)
* (mf5_8)
integer*8 m_8_0 (nr__) ! ? Spitzer/IRAC 8.0um (Band 4) flux
* calculation method flag (II/293)
* (mf8_0)
real*8 RAG360deg (nr__) ! (deg) ? GLIMPSE360 right ascension in
* degrees (J2000) (ra_1)
real*8 DEG360deg (nr__) ! (deg) [-46.5/65.5]? GLIMPSE360 declination
* in degrees (J2000) (dec_1)
character*19 WISEA (nr__) ! AllWISE source identifier
* (JHHMMSS.ss+DDMMSS.s; J2000) (II/328)
* (designation_2) (12)
real*8 RAAWdeg (nr__) ! (deg) ? AllWISE right ascension in degrees
* (J2000) (II/328) (ra_2)
real*8 DEAWdeg (nr__) ! (deg) [-46.5/65.5]? AllWISE declination in
* degrees (J2000) (II/328) (dec_2)
real*4 RAsig (nr__) ! (arcsec) [1e-3/0.6]? Uncertainty in RA (sigra)
real*4 DEsig (nr__) ! (arcsec) [1e-3/0.6215]? Uncertainty in DE
* (sigdec)
real*8 RADEsig (nr__) ! (arcsec) [-0.3/0.2]? Cross-term of RA and DE
* uncertainties (sigradec)
real*8 GLON_AW (nr__) ! (deg) ? AllWISE galactic longitude (II/328)
* (glon)
real*8 GLAT_AW (nr__) ! (deg) [-3.2/3.6]? AllWISE galactic latitude
* (II/328) (glat)
real*8 ELON_AW (nr__) ! (deg) ? AllWISE ecliptic longitude (II/328)
* (elon)
real*8 ELAT_AW (nr__) ! (deg) [-61.4/63.7]? AllWISE ecliptic
* latitude (II/328) (elat)
real*8 wx (nr__) ! (pix) [46.5/4058.4]? X-pixel coordinate on
* the Atlas Image from AllWISE, all
* bands (II/328)
real*8 wy (nr__) ! (pix) [38.6/4058.1]? Y-pixel coordinate on
* the Atlas Image from AllWISE, all
* bands (II/328)
real*4 W1mag (nr__) ! (mag) [5/16]? WISE W1 band (3.35um)
* magnitude (II/328) (w1mpro)
real*4 e_W1mag (nr__) ! (mag) [7e-3/0.52]? Mean error in W1mag
* (w1sigmpro)
real*4 snr1 (nr__) ! [2.1/157]? Signal to noise ratio for
* W1 filter (II/328) (w1snr)
real*8 chi2W1 (nr__) ! [0.02/2583]? Reduced {chi}^2^ of the
* W1 profile fit (w1rchi2)
real*4 W2mag (nr__) ! (mag) [2.9/17]? WISE W2 band (4.6um)
* magnitude (II/328) (w2mpro)
real*4 e_W2mag (nr__) ! (mag) [7e-3/0.55]? Mean error in W2mag
* (w2sigmpro)
real*4 snr2 (nr__) ! [2/156]? Signal to noise ratio for W2
* filter (II/328) (w2snr)
real*8 chi2W2 (nr__) ! [0.04/1345]? Reduced {chi}^2^ of the
* W2 profile fit (II/328) (w2rchi2)
real*8 W3mag (nr__) ! (mag) [-0.3/10]? WISE W3 band (11.6um)
* magnitude (II/328) (w3mpro)
real*4 e_W3mag (nr__) ! (mag) [4e-3/0.6]? Mean error in W3mag
* (w3sigmpro)
real*4 snr3 (nr__) ! [2/274]? Signal to noise ratio for W3
* filter (II/328) (w3snr)
real*8 chi2W3 (nr__) ! [0.01/3133]? Reduced {chi}^2^ of the
* W3 profile fit (II/328) (w3rchi2)
real*4 W4mag (nr__) ! (mag) [-4/9.7]? WISE W4 (22.1um) band
* magnitude (II/328) (w4mpro)
real*4 e_W4mag (nr__) ! (mag) [1e-3/0.6]? Mean error in W4mag
* (w4sigmpro)
real*4 snr4 (nr__) ! [-120/1307]? Signal to noise ratio for
* W4 filter (II/328) (w4snr)
real*8 chi2W4 (nr__) ! [0.02/10340]? Reduced {chi}^2^ of the
* W4 profile fit (II/328) (w4rchi2)
real*8 chi2 (nr__) ! [0.1/1908]? All-band combined reduced
* {chi}^2^ (II/328) (rchi2)
integer*4 nb (nr__) ! [1/4]? Number of PSF components used
* simultaneously in the profile-fitting
* (II/328)
integer*4 na (nr__) ! [0/1]? Active deblending flag (1 if a
* detection was split into multiple
* sources) (II/328)
real*4 sat1 (nr__) ! [0/0.67]? Saturated pixel fraction in
* W1 (II/328) (w1sat)
real*4 sat2 (nr__) ! [0/0.71]? Saturated pixel fraction in
* W2 (II/328) (w2sat)
real*4 sat3 (nr__) ! [0/0.32]? Saturated pixel fraction in
* W3 (II/328) (w3sat)
real*4 sat4 (nr__) ! [0/0.48]? Saturated pixel fraction in
* W4 (II/328) (w4sat)
character*4 satnum (nr__) ! Minimum sample at which saturation
* occurs in each band, 1 character per
* band (II/328)
real*8 RA_pm (nr__) ! (deg) ? Right ascension at epoch MJD=55400.0
* (2010.5589) from the profile-fitting
* measurement model that includes
* motion (II/328) (ra_pm)
real*8 DE_pm (nr__) ! (deg) [-47/66]? Declination at epoch
* MJD=55400.0 (2010.5589) from the
* profile-fitting measurement model that
* includes motion (II/328) (dec_pm)
real*4 e_RA_pm (nr__) ! (arcsec) [2e-3/4.5]? 1{sigma} uncertainty on
* RA_pm (sigra_pm)
real*4 e_DE_pm (nr__) ! (arcsec) [2e-3/2.8]? 1{sigma} uncertainty on
* DE_pm (sigdec_pm)
real*8 cosig_pm (nr__) ! (arcsec) [-1.9/2.3]? Co-sigma of the equatorial
* position uncertainties (II/328)
* (sigradec_pm)
integer*4 pmRA (nr__) ! (mas/yr) [-98190/130293]? Apparent motion in RA
* (II/328) (pmra) (4)
integer*4 e_pmRA (nr__) ! (mas/yr) [20/11196]? Mean error on pmRA
* (sigpmra)
integer*4 pmDE (nr__) ! (mas/yr) [-198842/169468]? Apparent motion in
* DE (II/328) (pmdec) (4)
integer*4 e_pmDE (nr__) ! (mas/yr) [20/12961]? Mean error on pmDE
* (sigpmdec)
real*8 chi2W1_pm (nr__) ! [0.04/2535]? Reduced {chi}^2^ of the
* W1 profile fit including motion
* estimation (II/328) (w1rchi2_pm)
real*8 chi2W2_pm (nr__) ! [0.05/1335]? Reduced {chi}^2^ of the
* W2 profile fit including motion
* estimation (II/328) (w2rchi2_pm)
real*8 chi2W3_pm (nr__) ! [0.01/3107]? Reduced {chi}^2^ of the
* W3 profile fit including motion
* estimation (II/328) (w3rchi2_pm)
real*8 chi2W4_pm (nr__) ! [0.02/10150]? Reduced {chi}^2^ of the
* W4 profile fit including motion
* estimation (II/328) (w4rchi2_pm)
real*8 chi2pm (nr__) ! [0.1/1905]? Combined reduced {chi}^2^
* profile fit including motion
* estimation (II/328) (rchi2_pm)
character*5 qpm (nr__) ! Motion estimation quality (II/328)
* (pmcode) (5)
character*4 ccf (nr__) ! Contamination and confusion flag, one
* per band (II/328) (cc_flags) (6)
integer*4 ex (nr__) ! [0/5]? Extended source flag (II/328)
* (ext_flg) (7)
character*4 var (nr__) ! Variability flag, one per band (II/328)
* (var_flg) (8)
character*4 qph (nr__) ! Photometric quality flag (II/328)
* (ph_qual) (9)
integer*4 fdet (nr__) ! [7/15]? Bit-encoded integer indicating
* bands in which a source has a snr?>2
* detection (II/328) (det_bit) (13)
character*4 fmoon (nr__) ! Scattered moonlight contamination
* flag, 1 per band (II/328) (moon_lev)
integer*4 nW1 (nr__) ! [2/81]? Frame detection count in W1
* (II/328) (w1nm) (10)
integer*4 mW1 (nr__) ! [16/81]? Integer frame coverage in W1
* (II/328) (w1m) (11)
integer*4 nW2 (nr__) ! [0/81]? Frame detection count in W2
* (II/328) (w2nm) (10)
integer*4 mW2 (nr__) ! [12/81]? Integer frame coverage in W2
* (II/328) (w2m) (11)
integer*4 nW3 (nr__) ! [0/55]? Frame detection count in W3
* (II/328) (w3nm) (10)
integer*4 mW3 (nr__) ! [0/57]? Integer frame coverage in W3
* (II/328) (w3m) (11)
integer*4 nW4 (nr__) ! [0/54]? Frame detection count in W4
* (II/328) (w4nm) (10)
integer*4 mW4 (nr__) ! [0/57]? Integer frame coverage in W4
* (II/328) (w4m) (11)
real*4 covW1 (nr__) ! [5.5/83.1]? Mean pixel coverage in W1
* (II/328) (w1cov)
real*4 covW2 (nr__) ! [2/83.1]? Mean pixel coverage in W2
* (II/328) (w2cov)
real*4 covW3 (nr__) ! [0.4/60.6]? Mean pixel coverage in W3
* (II/328) (w3cov)
real*4 covW4 (nr__) ! [1.4/60.3]? Mean pixel coverage in W4
* (II/328) (w4cov)
integer*8 v2Mkey (nr__) ! [68916543/1315293861]? 2MASS-PSC
* closest associated key
* (not identification) (II/328, II/246)
* (tmass_key)
real*4 Offset (nr__) ! (arcsec) [3e-3/3]? Radial offset between WISE
* source (II/328) and associated 2MASS
* source (II/246) (r_2mass)
integer*4 PA (nr__) ! (deg) [-180/180]? Position angle
* (North to East) of vector from WISE
* source (II/328) to associated 2MASS
* source (II/246) (pa_2mass)
integer*4 Num (nr__) ! [0/3]? Number of 2MASS-PSC entries
* found within search radius of WISE
* source (n_2mass)
real*4 W_Jmag (nr__) ! (mag) [5/19]? J-band (1.25 um) magnitude of
* the associated 2MASS-PSC source
* (II/246) (j\_m\_2mass)
real*4 e_W_Jmag (nr__) ! (mag) [0.01/0.7]? Mean error on W.Jmag
* (j_msig_2mass)
real*4 W_Hmag (nr__) ! (mag) [5.1/17.8]? H-band (1.65 um) magnitude
* of the associated 2MASS-PSC source
* (II/246) (h\_m\_2mass)
real*4 e_W_Hmag (nr__) ! (mag) [0.01/10]? Mean error on W.Hmag
* (h_msig_2mass)
real*4 W_Kmag (nr__) ! (mag) [5/16.5]? Ks-band (2.17 um) magnitude
* of the associated 2MASS-PSC source
* (II/246) (k\_m\_2mass)
real*4 e_W_Kmag (nr__) ! (mag) [0.01/10]? Mean error on W.Kmag
* (k_msig_2mass)
real*4 J2mag (nr__) ! (mag) [5.6/17.3]? 2MASS J-band (1.25um)
* magnitude from SMOG catalog/archive
* (J/ApJ/880/9, II/246) (mag_j2)
real*4 e_J2mag (nr__) ! (mag) [0.01/0.6]? Mean error on J2mag (dj2_m)
real*4 H2mag (nr__) ! (mag) [4.5/16.2]? 2MASS H-band (1.65um)
* magnitude from SMOG catalog/archive
* (J/ApJ/880/9, II/246) (mag_h2)
real*4 e_H2mag (nr__) ! (mag) [0.01/0.3]? Mean error on H2mag (dh2_m)
real*4 K2mag (nr__) ! (mag) [3.9/15.8]? 2MASS Ks-band (2.17um)
* magnitude from SMOG catalog/archive
* (J/ApJ/880/9, II/246) (mag_ks2)
real*4 e_K2mag (nr__) ! (mag) [0.01/10]? Mean error on K2mag (dks2_m)
real*8 J2Flux (nr__) ! (mJy) [0.1/8882]? 2MASS-PSC J-band (1.25um)
* flux from SMOG catalog/archive (f_j2)
real*8 e_J2Flux (nr__) ! (mJy) [0.03/207.1]? 1{sigma} uncertainty on
* J2Flux (df_j2)
real*8 H2Flux (nr__) ! (mJy) [0.3/14990]? 2MASS-PSC H-band (1.65um)
* flux from SMOG catalog/archive (f_h2)
real*8 e_H2Flux (nr__) ! (mJy) [0.05/3148]? 1{sigma} uncertainty on
* H2Flux (df_h2)
real*8 K2Flux (nr__) ! (mJy) [0.3/16790]? 2MASS-PSC Ks-band
* (2.17um) flux from SMOG
* catalog/archive (f_ks2)
real*8 e_K2Flux (nr__) ! (mJy) [0.06/19200]? 1{sigma} uncertainty on
* K2mag (df_ks2)
real*4 J2snr (nr__) ! [2/60.4]? 2MASS-PSC J-band (1.25um)
* signal-to-noise ratio from SMOG
* catalog/archive (sn_j2)
real*4 H2snr (nr__) ! [4/63.9]? 2MASS-PSC H-band (1.65um)
* signal-to-noise ratio from SMOG
* catalog/archive (sn_h2)
real*4 K2snr (nr__) ! [0.1/72.4]? 2MASS-PSC Ks-band (2.17um)
* signal-to-noise ratio from SMOG
* catalog/archive (sn_ks2)
integer*8 q_J2mag (nr__) ! ? 2MASS-PSC J-Band (1.25um) source
* quality flag from SMOG catalog/archive
* (sqf_j2) (3)
integer*8 q_H2mag (nr__) ! ? 2MASS-PSC H-Band (1.65um) source
* quality flag from SMOG catalog/archive
* (sqf_h2) (3)
integer*8 q_K2mag (nr__) ! ? 2MASS-PSC Ks-Band (2.17um) source
* quality flag from SMOG catalog/archive
* (sqf_ks2) (3)
integer*4 Counter (nr__) ! [23/2836315]? Entry counter (key)
* number (unique within table) (cntr)
real*8 v24mag (nr__) ! (mag) [-100/11.4]? Spitzer/MIPS 24um (Band
* 1) magnitude from SMOG catalog/archive
* (mag_24)
real*8 e_24mag (nr__) ! (mag) [4e-4/10]? Uncertainty in 24mag (d24_m)
*Note (1): This column is a merge of the original "designation" and
* "designation_1" columns. The two separated columns were an artifact
* of the way the catalog was assembled using different databases, and
* the merge was done with the authors' authorization.
* The object identifiers can start with the acronym SSTGLMA (GLIMPSE360
* database) or SSTSMOGA (SMOG database). GLIMPSE stands for Galactic Legacy
* Infrared Mid-Plane Survey Extraordinaire and SMOG stands for Spitzer
* Mapping of the Outer Galaxy.
* The objects with an identifier starting with "J" are WISE sources with no
* IRAC data.
*Note (2): The close source flag is set if there are Archive sources that are
* within 3" of the source. Flags are as follows:
* 0 = no Archive source within 3.0" of source
* 1 = Archive sources between 2.5" and 3.0" of source
* 2 = Archive sources between 2.0" and 2.5" of source
* 3 = Archive sources between 1.5" and 2.0" of source
* 4 = Archive sources between 1.0" and 1.5" of source
* 5 = Archive sources between 0.5" and 1.0" of source
* 6 = Archive sources within 0.5" of source
*Note (3): The value of the SQF is {Sigma}2^(bit-1)^. For example, a source with
* bits 1 and 4 set will have SQF=2^0^+2^3^=9. If the SQF is 0, the source
* has no detected problems. On the other hand, a source with SQF=29360128
* has a corresponding binary number of 000001110000000000000000000000,
* meaning the source has bits 23, 24, and 25.
* The source quality flags bits along with their origin are:
* ----------------------------------------------------------------------------
* Bit Explanation Origin
* ----------------------------------------------------------------------------
* 1 poor pixels in dark current SSC pmask
* 2 flat field questionable SSC dmask
* 3 latent image SSC dmask
* 3 persistence (p) 2MASS
* 4 photometric confusion (c) 2MASS
* 7 muxbleed correction applied GLIMPSE
* (for SSTSMOGA sources only)
* 8 hot, dead or otherwise unacceptable pixel SSC pmask,dmask,GLIMPSE
* 9 muxbleed correction applied is >3sig above bkg GLIMPSE
* (for SSTSMOGA sources only)
* 9 electronic stripe (s) 2MASS
* 10 DAOPHOT tweak positive GLIMPSE
* 11 DAOPHOT tweak negative GLIMPSE
* 13 confusion in in-band merge GLIMPSE
* 14 confusion in cross-band merge (IRAC) GLIMPSE
* 14 confusion in cross-band merge (2MASS) GLIMPSE
* 15 column pulldown corrected GLIMPSE
* 16 banding corrected GLIMPSE
* (for SSTSMOGA sources only)
* 19 data predicted to saturate GLIMPSE
* 20 saturated star wing region GLIMPSE
* 20 diffraction spike (d) 2MASS
* 21 pre-lumping in in-band merge GLIMPSE
* 22 post-lumping in cross-band merge (IRAC) GLIMPSE
* 22 post-lumping in cross-band merge (2MASS) GLIMPSE
* 23 photometry quality flag (*) 2MASS
* 24 photometry quality flag (*) 2MASS
* 25 photometry quality flag (*) 2MASS
* 30 within three pixels of edge of frame GLIMPSE
* ----------------------------------------------------------------------------
* (*) the 2MASS photometric quality is stored here on 3 bits, with:
*
* 2MASS SQF bits
* "ph" Flag => 23, 24, 25 value
* ------------- --------------- ----------
* X 0 0 0 0
* U 1 0 0 4194304
* F 0 1 0 8388608
* E 1 1 0 12582912
* D 0 0 1 16777216
* C 1 0 1 20971520
* B 0 1 1 25165824
* A 1 1 1 29360128
* ------------ --------------- ----------
*
* where
*
* X - There is a detection at this location, but no valid brightness
* estimate can be extracted using any algorithm.
* U - Upper limit on magnitude. Source is not detected in this band or it is
* detected, but not resolved in a consistent fashion with other bands.
* F - This category includes sources where a reliable estimate of the
* photometric error could not be determined.
* E - This category includes detections where the goodness-of-fit quality of
* the profile-fit photometry was very poor, or detections where psf fit
* photometry did not converge and an aperture magnitude is reported, or
* detections where the number of frames was too small in relation to the
* number of frames in which a detection was geometrically possible.
* D - Detections in any brightness regime where valid measurements were made
* with no [JHK]snr or [JHK]cmsig requirement. (**)
* C - Detections in any brightness regime where valid measurements were made
* with [JHK]snr > 5 AND [JHK]cmsig < 0.21714. (**)
* B - Detections in any brightness regime where valid measurements were made
* with [JHK]snr > 7 AND [JHK]cmsig < 0.15510. (**)
* A - Detections in any brightness regime where valid measurements were made
* with [JHK]snr > 10 AND [JHK]cmsig < 0.10857. (**)
*
* (**) See 2MASS-PSC (II/246) for [JHK]cmsig.
* ----------------------------------------------------------------------------
* Therefore, in the previous example, the source with SQF=29360128 has a
* photometric quality A flag from 2MASS-PSC.
*
* For more information, please refer to the following documentation:
* https://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/doc/\
* glimpse360_dataprod_v1.5.pdf (Table 5)
* https://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/doc/\
* smog_dataprod_v1.1.pdf (Table 4)
*Note (4): CAUTION - Note concerning "proper" motions: according to AllWISE
* documentation, the motions measured are not proper motions because
* they are affected by parallax motions. More information is given in:
* http://wise2.ipac.caltech.edu/docs/release/allwise/expsup/sec2_6.html\
* #How_To_Interpret
*Note (5): Five character string that encodes information about factors that
* impact the accuracy of the motion estimation. These include the
* original blend count, whether a blend-swap took place, and the
* distance in hundredths of an arcsecond between the non-motion
* position and the motion mean-observation-epoch position. The format
* is NQDDD where N is the original blend count, Q is either "Y" or "N"
* for "yes" or "no" a blend-swap occurred (i.e., the original primary
* component was not the final primary component), and DDD is the radial
* distance between the non-motion and motion at mean-observation epoch
* positions in units of 0.01 arcsec, clipped at 999 (almost 10 arcsec).
*Note (6): One character per band (W1/W2/W3/W4) that indicates that the
* photometry and/or position measurements of a source may be
* contaminated or biased due to proximity to an image artifact:
* D,d = Diffraction spike. Source may be a spurious detection of or
* contaminated by a diffraction spike from a nearby bright star on the
* same image
* P,p = Persistence. Source may be a spurious detection of or contaminated by
* a latent image left by a bright star
* H,h = Halo. Source photometry may be a spurious detection of or contaminated
* by the scattered light halo surrounding a nearby bright source
* O,o = (letter "o") Optical ghost. Source may be a spurious detection of or
* contaminated by an optical ghost image caused by a nearby bright source
* 0 = (number zero) Source is unaffected by known artifacts.
*Note (7): Extended source flag as follows:
* 0 = The source shape is consistent with a point-source and the source is
* not associated with or superimposed on a 2MASS XSC source 1.
* 1 = The profile-fit photometry goodness-of-fit, chi2W?, is >3.0 in one or
* more bands.
* 2 = The source falls within the extrapolated isophotal footprint of a
* 2MASS XSC source.
* 3 = The profile-fit photometry goodness-of-fit, chi2W?, is >3.0 in one or
* more bands, and The source falls within the extrapolated isophotal
* footprint of a 2MASS XSC source.
* 4 = The source position falls within 5" of a 2MASS XSC source.
* 5 = The profile-fit photometry goodness-of-fit, chi2W?, is >3.0 in one or
* more bands, and the source position falls within 5" of a 2MASS XSC
* source.
*Note (8): The variability flag is a 4-character string, one character per band,
* containing a value related to the probability that the source flux
* measured on the individual WISE exposures is variable.
* - value "n" indicates insufficient or inadequate data to make a
* determination (<6 exposures)
* - values 0 thru 9 indicate increasing probabilities of variation
* - values 0-5 are most likely not variable
* - values 6-7 are likely variables (but susceptible of false-positive
* variability)
* - values >7 have the highest probability of being true variables
*Note (9): Four character flag, one character per band (W1/W2/W3/W4), that
* provides a shorthand summary of the quality of the profile-fit
* photometry measurement in each band, as derived from the measurement
* signal-to-noise ratio:
* A = Source is detected in this band with a flux signal-to-noise ratio
* snr?>10.
* B = Source is detected in this band with a flux signal-to-noise ratio
* 3<snr?<10.
* C = Source is detected in this band with a flux signal-to-noise ratio
* 2<snr?<3.
* U = Upper limit on magnitude. Source measurement has snr?<2. The
* profile-fit magnitude W?mag is a 95% confidence upper limit.
* X = A profile-fit measurement was not possible at this location in this
* band. The value of W?mag and e_W?mag will be "null" in this band.
* Z = A profile-fit source flux measurement was made at this location, but
* the flux uncertainty could not be measured. The value of e_W?mag will
* be "null" in this band. The value of W?mag will be "null" if the
* measured flux, W?Flux, is negative, but will not be "null" if the flux
* is positive. If a non-null magnitude is present, it corresponds to the
* true flux, and not the 95% confidence upper limit. This occurs for a
* small number of sources found in a narrow range of ecliptic longitude
* which were covered by a large number of saturated pixels from 3-Band
* Cryo single-exposures.
*Note (10): Integer frame detection count. This column gives the number of
* individual 7.7s exposures on which this source was detected with
* SNR>3 in the W? profile-fit measurement. This number can be zero for
* sources that are well-detected on the coadded Atlas Image, but too
* faint for detection on the single exposures.
*Note (11): Number of individual 7.7s W? exposures on which a profile-fit
* measurement of this source was possible. This number can differ
* between the four bands because band-dependent criteria are used to
* select individual frames for inclusion in the coadded Atlas Images.
* This column is null if there is no frame coverage in this band at
* the position of this source.
*Note (12): These objects are from a crossmatch of the IRAC and AllWISE
* databases, with their AllWISE identifiers given here and their
* other source designations given in the "Source" column.
*Note (13): For example, a source detected in W1 only has fdet=1 (binary 0001).
* A source detected in W4 only has fdet=8 (binary 1000). A source
* detected in all four bands has fdet=15 (binary 1111).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! SFOG field YSOs: star classification table
integer*4 nr__1
parameter (nr__1=47405) ! Number of records
character*95 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*8 Field (nr__1) ! SMOG Field ID (SRC0 to SRC47404) (ID)
character*26 Source_1 (nr__1) ! Source designation (IR_name)
character*16 v2MASS_1 (nr__1) ! 2MASS identifier (2M name)
integer*4 RAh (nr__1) ! (h) Hour of right ascension (J2000)
integer*4 RAm (nr__1) ! (min) Minute of right ascension (J2000)
real*8 RAs (nr__1) ! (s) Second of right ascension (J2000)
character*1 DE_ (nr__1) ! Sign of declination (J2000)
integer*4 DEd (nr__1) ! (deg) Degree of declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Arcminute of declination (J2000)
real*8 DEs (nr__1) ! (arcsec) Arcsecond of declination (J2000)
integer*4 dIRAC (nr__1) ! [0/1] Identified as a YSO based on IRAC
* data set (1=yes) (dIRAC)
integer*4 dWISE (nr__1) ! [0/1] Identified as a YSO based on the WISE
* data set II/328 (1=yes) (dWISE)
integer*4 dIRACWISE (nr__1) ! [0/1] Identified as a YSO based on the
* combined IRAC+WISE data set (1=yes)
* (dIRACWISE) (1)
integer*4 dSMOG (nr__1) ! [0/1] Identified as a YSO based on SMOG
* data set (1=yes) (dSMOG)
integer*4 Seq (nr__1) ! [-1/620]? Cluster number (Cluster) (2)
integer*4 Class (nr__1) ! [1/3] Object class (Class) (3)
*Note (1): Paper I (Winston+2019, J/ApJ/880/9) compared the photometry in IRAC
* bands 1 and 2 and WISE bands 1 and 2 for sources matched to within
* 1.5" and found that the photometric magnitudes matched well to within
* the uncertainties. The combined IRAC+WISE data set uses photometry
* from the two available IRAC bands to replace the appropriate WISE
* band fluxes where the sources matched to within 1.5".
*Note (2): The cluster number of the YSO is indicated in this column.
* See Table 4.
* A "-1" indicates that there was no cluster affiliation identified.
*Note (3): The numerical YSO class, as described in Section 3.4.1., which
* depends on the slope, {alpha}, of the spectral energy distribution
* across the mid-IR bandpasses. Classes as follows:
* 1 = Class 0/I: protostellar objects with a rising slope, {alpha}>0;
* (10476 occurrences)
* 2 = Class II: characterized by decreasing slopes in the range -1.6<{alpha}<0;
* (29604 occurrences)
* 3 = Class IIa/III: class IIa sources lack optically thick emission from
* a disk and possess decreasing slopes {alpha}<-1.6, consistent with a
* weak emission above a stellar photosphere. Truly diskless class III
* objects show slopes in the range -2.7<{alpha}<-2.0;
* (7325 occurrences)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! SFOG Field YSOs: SIMBAD Matches within 2"
integer*4 nr__2
parameter (nr__2=7039) ! Number of records
character*816 ar__2 ! Full-size record
C J2000 position composed of: RASIMh RASIMm RASIMs DESIM- DESIMd DESIMm DESIMs
real*8 RAdeg_1 (nr__2) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__2) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*8 Field_1 (nr__2) ! SMOG Field ID (SRC9 to SRC47353) (SRC)
character*26 IRAC (nr__2) ! IRAC identifier (IRAC_ID)
real*8 RAIRACdeg (nr__2) ! (deg) Right ascension from IRAC (J2000) (IRAC_RA)
real*8 DEIRACdeg (nr__2) ! (deg) Declination from IRAC (J2000) (IRAC_Dec)
character*34 Simbad (nr__2) ! SIMBAD identifier (SIMBAD_ID)
character*19 Type (nr__2) ! Object type (OTYPE)
integer*4 RASIMh (nr__2) ! (h) Hour of right ascension from SIMBAD (J2000)
* (SIMBAD_RA)
integer*4 RASIMm (nr__2) ! (min) Minute of right ascension from Simbad
* (J2000)
real*8 RASIMs (nr__2) ! (s) Second of right ascension from Simbad
* (J2000)
character*1 DESIM_ (nr__2) ! Sign of declination from SIMBAD (J2000)
* (SIMBAD_Dec)
integer*4 DESIMd (nr__2) ! (deg) Degree of declination from SIMBAD (J2000)
real*4 DESIMm (nr__2) ! (arcmin) Arcminute of declination from SIMBAD (J2000)
real*4 DESIMs (nr__2) ! (arcsec) Arcsecond of declination from SIMBAD (J2000)
character*675 IDs (nr__2) ! List of SIMBAD identifiers for this object
* (IDS)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! SFOG YSO cluster descriptions
integer*4 nr__3
parameter (nr__3=621) ! Number of records
character*51 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq_1 (nr__3) ! [0/620] Cluster number (Cluster_#)
character*12 GLIMPSE (nr__3) ! GLIMPSE identifier (GLLL.llll+BB.bbbb)
* (Cluster_Name)
integer*4 YSO (nr__3) ! [5/1177] Number of YSO in the cluster (#_YSO)
real*8 RAdeg_2 (nr__3) ! (deg) Right ascension of cluster center (J2000)
* (Central_RA)
real*8 DEdeg_2 (nr__3) ! (deg) Declination of cluster center (J2000)
* (Central_Dec)
real*8 Rad (nr__3) ! (deg) [0.013/1.21] Circular radius (Circular_Radius)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 't4wise.dat' ! SFOG YSO clusters: WISE HII counterparts
(reformatted from Table 4 by CDS)
integer*4 nr__4
parameter (nr__4=870) ! Number of records
character*32 ar__4 ! Full-size record
integer*4 Seq_2 (nr__4) ! [0/620] Cluster number
character*12 GLIMPSE_1 (nr__4) ! GLIMPSE identifier (GLLL.llll+BB.bbbb)
* (Cluster_name)
character*15 WISE (nr__4) ! WISE H II counterparts
* (Anderson+2014, J/ApJS/212/1) (G1)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 't4skiff.dat' ! SFOG YSO clusters: O & B objects from Skiff+2014
(reformatted from Table 4 by CDS)
integer*4 nr__5
parameter (nr__5=6768) ! Number of records
character*43 ar__5 ! Full-size record
integer*4 Seq_3 (nr__5) ! [0/620] Cluster number
character*12 GLIMPSE_2 (nr__5) ! GLIMPSE identifier (GLLL.ll+B.bb)
* (Cluster_name)
character*26 Skiff (nr__5) ! Skiff identifier (Skiff+2009-, B/mk) (G1)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 't4simbad.dat' ! SFOG YSO clusters: objects from SIMBAD
(reformatted from Table 4 by CDS)
integer*4 nr__6
parameter (nr__6=38445) ! Number of records
character*51 ar__6 ! Full-size record
integer*4 Seq_4 (nr__6) ! [0/620] Cluster number
character*12 GLIMPSE_3 (nr__6) ! GLIMPSE identifier (G.LLL.ll+B.bb)
* (Cluster_name)
character*34 Simbad_1 (nr__6) ! SIMBAD objects located within the convex hull
* of the cluster (G1)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! SFOG Field YSOs: SEDFitter Results
integer*4 nr__7
parameter (nr__7=6234) ! Number of records
character*420 ar__7 ! Full-size record
character*8 Field_2 (nr__7) ! SMOG Field ID (SRC588 to SRC47400) (SRC_ID)
integer*4 Seq_5 (nr__7) ! [5/593] Cluster number
real*4 Dist (nr__7) ! (kpc) [0.7/13.3] Distance
real*4 RadGal (nr__7) ! (kpc) [7.8/18.1] Galactic radius
integer*4 Ndata (nr__7) ! [3/7] Number of SED data points (N_data)
real*4 Chi2_1 (nr__7) ! [7.2e-7/3] {chi}^2^
real*4 e_Chi2 (nr__7) ! [7.2e-7/3] Lower limit on {chi}^2^
* (Chi_2_l)
real*4 E_Chi2_1 (nr__7) ! [0.67/24] Upper limit on {chi}^2^ (Chi_2_h)
real*4 Av (nr__7) ! (mag) [0/40] Extinction A_V_
real*4 e_Av (nr__7) ! (mag) [0/24.3] Lower limit on A_V_ (Av_l)
real*4 E_Av_1 (nr__7) ! (mag) [0.3/40] Upper limit on A_V_ (Av_h)
real*4 Av_avg (nr__7) ! (mag) [0/40] Average A_V_ (Av_Ave) (1)
real*4 s_Av (nr__7) ! (mag) [0/16] Standard deviation on A_V_ (Av_SD)
real*4 Mass (nr__7) ! (Msun) [0.1/31.2] Object mass (Mass)
real*4 e_Mass (nr__7) ! (Msun) [0.1/17.6] Lower limit on Mass (Mass_l)
real*4 E_Mass_1 (nr__7) ! (Msun) [1.3/47.3] Upper limit on Mass (Mass_h)
real*4 Mass_avg (nr__7) ! (Msun) [0.1/20.4] Average of Mass (Mass_Ave) (1)
real*4 s_Mass (nr__7) ! (Msun) [0/10.4] Standard deviation on Mass
* (Mass_SD)
real*4 Age (nr__7) ! (yr) [1015/1e+7] Age of the object
real*4 e_Age (nr__7) ! (yr) [1008/5.5e+6] Lower limit on Age (Age_l)
real*4 E_Age_1 (nr__7) ! (yr) [2210/1e+7] Upper limit on Age (Age_h)
real*4 Age_avg (nr__7) ! (yr) [1015/1e+7] Average of Age (Age_Ave) (1)
real*4 s_Age (nr__7) ! (yr) [0/4.8e+06] Standard deviation on Age
* (Age_SD)
real*4 Mdisk (nr__7) ! (Msun) [0/0.96] Disk mass M_disk_ (Disk_Mass)
real*4 e_Mdisk (nr__7) ! (Msun) [0/0.21] Lower limit on M_disk_
* (Disk-Mass-l)
real*4 E_Mdisk_1 (nr__7) ! (Msun) [0/1.9] Upper limit on M_disk_
* (Disk-mass-h)
real*4 Mdisk_avg (nr__7) ! (Msun) [0/0.9] Average of M_disk_ (Disk-Mass-Ave)
* (1)
real*4 s_Mdisk (nr__7) ! (Msun) [0/0.45] Standard deviation on M_disk_
* (Disk-Mass-SD)
real*4 Menv (nr__7) ! (Msun) [6.5e-10/3773] Envelope mass M_env_
* (Envl-Mass)
real*4 e_Menv (nr__7) ! (Msun) [1.7e-10/482] Lower limit on M_env_
* (Envl-Mass-l)
real*4 E_Menv_1 (nr__7) ! (Msun) [4.3e-07/4351] Upper limit on M_env_
* (Envl-mass-h)
real*4 Menv_avg (nr__7) ! (Msun) [1.4e-09/3773] Average of M_env_
* (Envl-Mass-Ave) (1)
real*4 s_Menv (nr__7) ! (Msun) [0/1009] Standard deviation on M_env_
* (Envl-Mass-SD)
real*4 Mdot (nr__7) ! (Msun/yr) [3.5e-10/7e-3]? Disk accretion rate Mdot
* (Mdot_star)
real*4 e_Mdot (nr__7) ! (Msun/yr) [0/1.4e-3] Lower limit on Mdot
* (Mdot-star-l)
real*4 E_Mdot_1 (nr__7) ! (Msun/yr) [0/1e-2] Upper limit on Mdot (Mdot-star-h)
real*4 Mdot_avg (nr__7) ! (Msun/yr) [0/7e-3] Average of Mdot (Mdot-star-Ave)
* (1)
real*4 s_Mdot (nr__7) ! (Msun/yr) [0/2.4e-3] Standard deviation on Mdot
* (Mdot-star-SD)
integer*4 Temp (nr__7) ! (K) [2550/40780] Central temperature T_*_
* (Temp)
real*4 e_Temp (nr__7) ! (K) [2536/26150] Lower limit on T_*_ (Temp_l)
real*4 E_Temp_1 (nr__7) ! (K) [4024/45160] Upper limit on T_*_ (Temp_h)
real*4 Temp_avg (nr__7) ! (K) [2651/35340] Average of T_*_ (Temp_Ave) (1)
real*4 s_Temp (nr__7) ! (K) [0/13230] Standard deviation on T_*_
* (Temp_SD)
real*4 Lstar (nr__7) ! (Lsun) [0.05/138900] Total luminosity L_*_
* (Tot_Lum)
real*4 e_Lstar (nr__7) ! (Lsun) [0.04/21500] Lower limit on L_*_
* (Tot-Lum-l)
real*4 E_Lstar_1 (nr__7) ! (Lsun) [1.5/327300] Upper limit on L_*_
* (Tot-Lum-h)
real*4 Lstar_avg (nr__7) ! (Lsun) [0.09/47920] Average of L_*_ (Tot-Lum-Ave)
* (1)
real*4 s_Lstar (nr__7) ! (Lsun) [0/36870] Standard deviation on L_*_
* (Tot-Lum-SD)
*Note (1): Weighted average value of the parameter for fits with {chi}^2^<3.
C=============================================================================
C Loading file 'table1.dat' ! Photometry of YSOs in Star formation Field in
* the Outer Galaxy (SFOG)
C Format for file interpretation
1 format(
+ A26,1X,A16,1X,I11,1X,F10.6,1X,F9.6,1X,F3.1,1X,F3.1,1X,F10.6,
+ 1X,F10.6,1X,F3.1,1X,F3.1,1X,I1,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,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F6.3,1X,F5.3,1X,F10.5,1X,F8.5,1X,F8.3,1X,F8.5,1X,
+ F9.4,1X,F9.4,1X,F11.5,1X,F11.6,1X,F10.5,1X,F10.6,1X,F9.4,1X,
+ F9.5,1X,F9.4,1X,F9.5,1X,F5.1,1X,I3,1X,F5.1,1X,F5.2,1X,F7.3,1X,
+ F6.3,1X,F7.3,1X,F7.3,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F5.2,1X,F5.1,1X,F5.1,1X,F4.1,1X,F4.1,1X,I2,1X,I2,
+ 1X,I2,1X,I2,1X,I2,1X,I2,1X,I2,1X,I2,1X,I8,1X,I8,1X,I8,1X,I11,
+ 1X,I11,1X,I11,1X,I11,1X,I11,1X,I11,1X,I11,1X,I11,1X,F10.6,1X,
+ F10.6,1X,A19,1X,F10.6,1X,F8.4,1X,F6.4,1X,F6.4,1X,F7.4,1X,F8.4,
+ 1X,F8.5,1X,F11.7,1X,F8.4,1X,F8.3,1X,F8.3,1X,F6.3,1X,F5.3,1X,
+ F5.1,1X,F10.5,1X,F6.3,1X,F5.3,1X,F5.1,1X,F10.5,1X,F7.3,1X,
+ F5.3,1X,F5.1,1X,F10.5,1X,F6.3,1X,F5.3,1X,F6.1,1X,F11.5,1X,
+ F9.4,1X,I1,1X,I1,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,A4,1X,
+ F10.6,1X,F8.4,1X,F6.4,1X,F6.4,1X,F7.4,1X,I6,1X,I5,1X,I7,1X,I5,
+ 1X,F10.5,1X,F10.5,1X,F10.5,1X,F11.5,1X,F9.4,1X,A5,1X,A4,1X,I1,
+ 1X,A4,1X,A4,1X,I2,1X,A4,1X,I2,1X,I2,1X,I2,1X,I2,1X,I2,1X,I2,
+ 1X,I2,1X,I2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,I10,1X,F5.3,1X,
+ I4,1X,I1,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,F5.3,1X,F6.3,1X,F5.3,1X,F9.4,1X,F9.5,
+ 1X,F10.4,1X,F10.5,1X,F10.4,1X,F11.5,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,I11,1X,I11,1X,I11,1X,I7,1X,F11.6,1X,F12.9)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,47405
read(1,'(A1383)')ar__
read(ar__,1)
+ Source(i__),v2MASS(i__),v2MASS_cntr(i__),GLON_2M(i__),
+ GLAT_2M(i__),e_GLON_2M(i__),e_GLAT_2M(i__),RA2Mdeg(i__),
+ DE2Mdeg(i__),e_RA2Mdeg(i__),e_DE2Mdeg(i__),csf(i__),Jmag(i__),
+ e_Jmag(i__),Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(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__),
+ F_J(i__),e_F_J(i__),F_H(i__),e_F_H(i__),F_K(i__),e_F_K(i__),
+ F_3_6(i__),e_F_3_6(i__),F_4_5(i__),e_F_4_5(i__),F_5_8(i__),
+ e_F_5_8(i__),F_8_0(i__),e_F_8_0(i__),rms_3_6(i__),
+ rms_4_5(i__),rms_5_8(i__),rms_8_0(i__),sky_3_6(i__),
+ sky_4_5(i__),sky_5_8(i__),sky_8_0(i__),Jsn(i__),Hsn(i__),
+ Ksn(i__),v3_6sn(i__),v4_5sn(i__),v5_8sn(i__),v8_0sn(i__),
+ d_3_6(i__),d_4_5(i__),d_5_8(i__),d_8_0(i__),N_3_6(i__),
+ N_4_5(i__),N_5_8(i__),N_8_0(i__),n_3_6_1(i__),n_4_5_1(i__),
+ n_5_8_1(i__),n_8_0_1(i__),q_Jmag(i__),q_Hmag(i__),q_Kmag(i__),
+ q_3_6mag(i__),q_4_5mag(i__),q_5_8mag(i__),q_8_0mag(i__),
+ m_3_6(i__),m_4_5(i__),m_5_8(i__),m_8_0(i__),RAG360deg(i__),
+ DEG360deg(i__),WISEA(i__),RAAWdeg(i__),DEAWdeg(i__),
+ RAsig(i__),DEsig(i__),RADEsig(i__),GLON_AW(i__),GLAT_AW(i__),
+ ELON_AW(i__),ELAT_AW(i__),wx(i__),wy(i__),W1mag(i__),
+ e_W1mag(i__),snr1(i__),chi2W1(i__),W2mag(i__),e_W2mag(i__),
+ snr2(i__),chi2W2(i__),W3mag(i__),e_W3mag(i__),snr3(i__),
+ chi2W3(i__),W4mag(i__),e_W4mag(i__),snr4(i__),chi2W4(i__),
+ chi2(i__),nb(i__),na(i__),sat1(i__),sat2(i__),sat3(i__),
+ sat4(i__),satnum(i__),RA_pm(i__),DE_pm(i__),e_RA_pm(i__),
+ e_DE_pm(i__),cosig_pm(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),
+ e_pmDE(i__),chi2W1_pm(i__),chi2W2_pm(i__),chi2W3_pm(i__),
+ chi2W4_pm(i__),chi2pm(i__),qpm(i__),ccf(i__),ex(i__),var(i__),
+ qph(i__),fdet(i__),fmoon(i__),nW1(i__),mW1(i__),nW2(i__),
+ mW2(i__),nW3(i__),mW3(i__),nW4(i__),mW4(i__),covW1(i__),
+ covW2(i__),covW3(i__),covW4(i__),v2Mkey(i__),Offset(i__),
+ PA(i__),Num(i__),W_Jmag(i__),e_W_Jmag(i__),W_Hmag(i__),
+ e_W_Hmag(i__),W_Kmag(i__),e_W_Kmag(i__),J2mag(i__),
+ e_J2mag(i__),H2mag(i__),e_H2mag(i__),K2mag(i__),e_K2mag(i__),
+ J2Flux(i__),e_J2Flux(i__),H2Flux(i__),e_H2Flux(i__),
+ K2Flux(i__),e_K2Flux(i__),J2snr(i__),H2snr(i__),K2snr(i__),
+ q_J2mag(i__),q_H2mag(i__),q_K2mag(i__),Counter(i__),
+ v24mag(i__),e_24mag(i__)
if(ar__(45:55) .EQ. '') v2MASS_cntr(i__) = iNULL__
if(ar__(57:66) .EQ. '') GLON_2M(i__) = rNULL__
if(ar__(68:76) .EQ. '') GLAT_2M(i__) = rNULL__
if(ar__(78:80) .EQ. '') e_GLON_2M(i__) = rNULL__
if(ar__(82:84) .EQ. '') e_GLAT_2M(i__) = rNULL__
if(ar__(86:95) .EQ. '') RA2Mdeg(i__) = rNULL__
if(ar__(97:106) .EQ. '') DE2Mdeg(i__) = rNULL__
if(ar__(108:110) .EQ. '') e_RA2Mdeg(i__) = rNULL__
if(ar__(112:114) .EQ. '') e_DE2Mdeg(i__) = rNULL__
if(ar__(116:116) .EQ. '') csf(i__) = iNULL__
if(ar__(118:123) .EQ. '') Jmag(i__) = rNULL__
if(ar__(125:129) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(131:136) .EQ. '') Hmag(i__) = rNULL__
if(ar__(138:142) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__(144:149) .EQ. '') Kmag(i__) = rNULL__
if(ar__(151:155) .EQ. '') e_Kmag(i__) = rNULL__
if(ar__(157:162) .EQ. '') v3_6mag(i__) = rNULL__
if(ar__(164:168) .EQ. '') e_3_6mag(i__) = rNULL__
if(ar__(170:175) .EQ. '') v4_5mag(i__) = rNULL__
if(ar__(177:181) .EQ. '') e_4_5mag(i__) = rNULL__
if(ar__(183:188) .EQ. '') v5_8mag(i__) = rNULL__
if(ar__(190:194) .EQ. '') e_5_8mag(i__) = rNULL__
if(ar__(196:201) .EQ. '') v8_0mag(i__) = rNULL__
if(ar__(203:207) .EQ. '') e_8_0mag(i__) = rNULL__
if(ar__(209:218) .EQ. '') F_J(i__) = rNULL__
if(ar__(220:227) .EQ. '') e_F_J(i__) = rNULL__
if(ar__(229:236) .EQ. '') F_H(i__) = rNULL__
if(ar__(238:245) .EQ. '') e_F_H(i__) = rNULL__
if(ar__(247:255) .EQ. '') F_K(i__) = rNULL__
if(ar__(257:265) .EQ. '') e_F_K(i__) = rNULL__
if(ar__(267:277) .EQ. '') F_3_6(i__) = rNULL__
if(ar__(279:289) .EQ. '') e_F_3_6(i__) = rNULL__
if(ar__(291:300) .EQ. '') F_4_5(i__) = rNULL__
if(ar__(302:311) .EQ. '') e_F_4_5(i__) = rNULL__
if(ar__(313:321) .EQ. '') F_5_8(i__) = rNULL__
if(ar__(323:331) .EQ. '') e_F_5_8(i__) = rNULL__
if(ar__(333:341) .EQ. '') F_8_0(i__) = rNULL__
if(ar__(343:351) .EQ. '') e_F_8_0(i__) = rNULL__
if(ar__(353:357) .EQ. '') rms_3_6(i__) = rNULL__
if(ar__(359:361) .EQ. '') rms_4_5(i__) = iNULL__
if(ar__(363:367) .EQ. '') rms_5_8(i__) = rNULL__
if(ar__(369:373) .EQ. '') rms_8_0(i__) = rNULL__
if(ar__(375:381) .EQ. '') sky_3_6(i__) = rNULL__
if(ar__(383:388) .EQ. '') sky_4_5(i__) = rNULL__
if(ar__(390:396) .EQ. '') sky_5_8(i__) = rNULL__
if(ar__(398:404) .EQ. '') sky_8_0(i__) = rNULL__
if(ar__(406:410) .EQ. '') Jsn(i__) = rNULL__
if(ar__(412:416) .EQ. '') Hsn(i__) = rNULL__
if(ar__(418:422) .EQ. '') Ksn(i__) = rNULL__
if(ar__(424:428) .EQ. '') v3_6sn(i__) = rNULL__
if(ar__(430:434) .EQ. '') v4_5sn(i__) = rNULL__
if(ar__(436:440) .EQ. '') v5_8sn(i__) = rNULL__
if(ar__(442:446) .EQ. '') v8_0sn(i__) = rNULL__
if(ar__(448:452) .EQ. '') d_3_6(i__) = rNULL__
if(ar__(454:458) .EQ. '') d_4_5(i__) = rNULL__
if(ar__(460:463) .EQ. '') d_5_8(i__) = rNULL__
if(ar__(465:468) .EQ. '') d_8_0(i__) = rNULL__
if(ar__(470:471) .EQ. '') N_3_6(i__) = iNULL__
if(ar__(473:474) .EQ. '') N_4_5(i__) = iNULL__
if(ar__(476:477) .EQ. '') N_5_8(i__) = iNULL__
if(ar__(479:480) .EQ. '') N_8_0(i__) = iNULL__
if(ar__(482:483) .EQ. '') n_3_6_1(i__) = iNULL__
if(ar__(485:486) .EQ. '') n_4_5_1(i__) = iNULL__
if(ar__(488:489) .EQ. '') n_5_8_1(i__) = iNULL__
if(ar__(491:492) .EQ. '') n_8_0_1(i__) = iNULL__
if(ar__(494:501) .EQ. '') q_Jmag(i__) = iNULL__
if(ar__(503:510) .EQ. '') q_Hmag(i__) = iNULL__
if(ar__(512:519) .EQ. '') q_Kmag(i__) = iNULL__
if(ar__(521:531) .EQ. '') q_3_6mag(i__) = iNULL__
if(ar__(533:543) .EQ. '') q_4_5mag(i__) = iNULL__
if(ar__(545:555) .EQ. '') q_5_8mag(i__) = iNULL__
if(ar__(557:567) .EQ. '') q_8_0mag(i__) = iNULL__
if(ar__(569:579) .EQ. '') m_3_6(i__) = iNULL__
if(ar__(581:591) .EQ. '') m_4_5(i__) = iNULL__
if(ar__(593:603) .EQ. '') m_5_8(i__) = iNULL__
if(ar__(605:615) .EQ. '') m_8_0(i__) = iNULL__
if(ar__(617:626) .EQ. '') RAG360deg(i__) = rNULL__
if(ar__(628:637) .EQ. '') DEG360deg(i__) = rNULL__
if(ar__(659:668) .EQ. '') RAAWdeg(i__) = rNULL__
if(ar__(670:677) .EQ. '') DEAWdeg(i__) = rNULL__
if(ar__(679:684) .EQ. '') RAsig(i__) = rNULL__
if(ar__(686:691) .EQ. '') DEsig(i__) = rNULL__
if(ar__(693:699) .EQ. '') RADEsig(i__) = rNULL__
if(ar__(701:708) .EQ. '') GLON_AW(i__) = rNULL__
if(ar__(710:717) .EQ. '') GLAT_AW(i__) = rNULL__
if(ar__(719:729) .EQ. '') ELON_AW(i__) = rNULL__
if(ar__(731:738) .EQ. '') ELAT_AW(i__) = rNULL__
if(ar__(740:747) .EQ. '') wx(i__) = rNULL__
if(ar__(749:756) .EQ. '') wy(i__) = rNULL__
if(ar__(758:763) .EQ. '') W1mag(i__) = rNULL__
if(ar__(765:769) .EQ. '') e_W1mag(i__) = rNULL__
if(ar__(771:775) .EQ. '') snr1(i__) = rNULL__
if(ar__(777:786) .EQ. '') chi2W1(i__) = rNULL__
if(ar__(788:793) .EQ. '') W2mag(i__) = rNULL__
if(ar__(795:799) .EQ. '') e_W2mag(i__) = rNULL__
if(ar__(801:805) .EQ. '') snr2(i__) = rNULL__
if(ar__(807:816) .EQ. '') chi2W2(i__) = rNULL__
if(ar__(818:824) .EQ. '') W3mag(i__) = rNULL__
if(ar__(826:830) .EQ. '') e_W3mag(i__) = rNULL__
if(ar__(832:836) .EQ. '') snr3(i__) = rNULL__
if(ar__(838:847) .EQ. '') chi2W3(i__) = rNULL__
if(ar__(849:854) .EQ. '') W4mag(i__) = rNULL__
if(ar__(856:860) .EQ. '') e_W4mag(i__) = rNULL__
if(ar__(862:867) .EQ. '') snr4(i__) = rNULL__
if(ar__(869:879) .EQ. '') chi2W4(i__) = rNULL__
if(ar__(881:889) .EQ. '') chi2(i__) = rNULL__
if(ar__(891:891) .EQ. '') nb(i__) = iNULL__
if(ar__(893:893) .EQ. '') na(i__) = iNULL__
if(ar__(895:899) .EQ. '') sat1(i__) = rNULL__
if(ar__(901:905) .EQ. '') sat2(i__) = rNULL__
if(ar__(907:911) .EQ. '') sat3(i__) = rNULL__
if(ar__(913:917) .EQ. '') sat4(i__) = rNULL__
if(ar__(924:933) .EQ. '') RA_pm(i__) = rNULL__
if(ar__(935:942) .EQ. '') DE_pm(i__) = rNULL__
if(ar__(944:949) .EQ. '') e_RA_pm(i__) = rNULL__
if(ar__(951:956) .EQ. '') e_DE_pm(i__) = rNULL__
if(ar__(958:964) .EQ. '') cosig_pm(i__) = rNULL__
if(ar__(966:971) .EQ. '') pmRA(i__) = iNULL__
if(ar__(973:977) .EQ. '') e_pmRA(i__) = iNULL__
if(ar__(979:985) .EQ. '') pmDE(i__) = iNULL__
if(ar__(987:991) .EQ. '') e_pmDE(i__) = iNULL__
if(ar__(993:1002) .EQ. '') chi2W1_pm(i__) = rNULL__
if(ar__(1004:1013) .EQ. '') chi2W2_pm(i__) = rNULL__
if(ar__(1015:1024) .EQ. '') chi2W3_pm(i__) = rNULL__
if(ar__(1026:1036) .EQ. '') chi2W4_pm(i__) = rNULL__
if(ar__(1038:1046) .EQ. '') chi2pm(i__) = rNULL__
if(ar__(1059:1059) .EQ. '') ex(i__) = iNULL__
if(ar__(1071:1072) .EQ. '') fdet(i__) = iNULL__
if(ar__(1079:1080) .EQ. '') nW1(i__) = iNULL__
if(ar__(1082:1083) .EQ. '') mW1(i__) = iNULL__
if(ar__(1085:1086) .EQ. '') nW2(i__) = iNULL__
if(ar__(1088:1089) .EQ. '') mW2(i__) = iNULL__
if(ar__(1091:1092) .EQ. '') nW3(i__) = iNULL__
if(ar__(1094:1095) .EQ. '') mW3(i__) = iNULL__
if(ar__(1097:1098) .EQ. '') nW4(i__) = iNULL__
if(ar__(1100:1101) .EQ. '') mW4(i__) = iNULL__
if(ar__(1103:1108) .EQ. '') covW1(i__) = rNULL__
if(ar__(1110:1115) .EQ. '') covW2(i__) = rNULL__
if(ar__(1117:1122) .EQ. '') covW3(i__) = rNULL__
if(ar__(1124:1129) .EQ. '') covW4(i__) = rNULL__
if(ar__(1131:1140) .EQ. '') v2Mkey(i__) = iNULL__
if(ar__(1142:1146) .EQ. '') Offset(i__) = rNULL__
if(ar__(1148:1151) .EQ. '') PA(i__) = iNULL__
if(ar__(1153:1153) .EQ. '') Num(i__) = iNULL__
if(ar__(1155:1160) .EQ. '') W_Jmag(i__) = rNULL__
if(ar__(1162:1166) .EQ. '') e_W_Jmag(i__) = rNULL__
if(ar__(1168:1173) .EQ. '') W_Hmag(i__) = rNULL__
if(ar__(1175:1179) .EQ. '') e_W_Hmag(i__) = rNULL__
if(ar__(1181:1186) .EQ. '') W_Kmag(i__) = rNULL__
if(ar__(1188:1192) .EQ. '') e_W_Kmag(i__) = rNULL__
if(ar__(1194:1199) .EQ. '') J2mag(i__) = rNULL__
if(ar__(1201:1205) .EQ. '') e_J2mag(i__) = rNULL__
if(ar__(1207:1212) .EQ. '') H2mag(i__) = rNULL__
if(ar__(1214:1218) .EQ. '') e_H2mag(i__) = rNULL__
if(ar__(1220:1225) .EQ. '') K2mag(i__) = rNULL__
if(ar__(1227:1231) .EQ. '') e_K2mag(i__) = rNULL__
if(ar__(1233:1241) .EQ. '') J2Flux(i__) = rNULL__
if(ar__(1243:1251) .EQ. '') e_J2Flux(i__) = rNULL__
if(ar__(1253:1262) .EQ. '') H2Flux(i__) = rNULL__
if(ar__(1264:1273) .EQ. '') e_H2Flux(i__) = rNULL__
if(ar__(1275:1284) .EQ. '') K2Flux(i__) = rNULL__
if(ar__(1286:1296) .EQ. '') e_K2Flux(i__) = rNULL__
if(ar__(1298:1302) .EQ. '') J2snr(i__) = rNULL__
if(ar__(1304:1308) .EQ. '') H2snr(i__) = rNULL__
if(ar__(1310:1314) .EQ. '') K2snr(i__) = rNULL__
if(ar__(1316:1326) .EQ. '') q_J2mag(i__) = iNULL__
if(ar__(1328:1338) .EQ. '') q_H2mag(i__) = iNULL__
if(ar__(1340:1350) .EQ. '') q_K2mag(i__) = iNULL__
if(ar__(1352:1358) .EQ. '') Counter(i__) = iNULL__
if(ar__(1360:1370) .EQ. '') v24mag(i__) = rNULL__
if(ar__(1372:1383) .EQ. '') e_24mag(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Source(i__),v2MASS(i__),v2MASS_cntr(i__),GLON_2M(i__),
+ GLAT_2M(i__),e_GLON_2M(i__),e_GLAT_2M(i__),RA2Mdeg(i__),
+ DE2Mdeg(i__),e_RA2Mdeg(i__),e_DE2Mdeg(i__),csf(i__),Jmag(i__),
+ e_Jmag(i__),Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(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__),
+ F_J(i__),e_F_J(i__),F_H(i__),e_F_H(i__),F_K(i__),e_F_K(i__),
+ F_3_6(i__),e_F_3_6(i__),F_4_5(i__),e_F_4_5(i__),F_5_8(i__),
+ e_F_5_8(i__),F_8_0(i__),e_F_8_0(i__),rms_3_6(i__),
+ rms_4_5(i__),rms_5_8(i__),rms_8_0(i__),sky_3_6(i__),
+ sky_4_5(i__),sky_5_8(i__),sky_8_0(i__),Jsn(i__),Hsn(i__),
+ Ksn(i__),v3_6sn(i__),v4_5sn(i__),v5_8sn(i__),v8_0sn(i__),
+ d_3_6(i__),d_4_5(i__),d_5_8(i__),d_8_0(i__),N_3_6(i__),
+ N_4_5(i__),N_5_8(i__),N_8_0(i__),n_3_6_1(i__),n_4_5_1(i__),
+ n_5_8_1(i__),n_8_0_1(i__),q_Jmag(i__),q_Hmag(i__),q_Kmag(i__),
+ q_3_6mag(i__),q_4_5mag(i__),q_5_8mag(i__),q_8_0mag(i__),
+ m_3_6(i__),m_4_5(i__),m_5_8(i__),m_8_0(i__),RAG360deg(i__),
+ DEG360deg(i__),WISEA(i__),RAAWdeg(i__),DEAWdeg(i__),
+ RAsig(i__),DEsig(i__),RADEsig(i__),GLON_AW(i__),GLAT_AW(i__),
+ ELON_AW(i__),ELAT_AW(i__),wx(i__),wy(i__),W1mag(i__),
+ e_W1mag(i__),snr1(i__),chi2W1(i__),W2mag(i__),e_W2mag(i__),
+ snr2(i__),chi2W2(i__),W3mag(i__),e_W3mag(i__),snr3(i__),
+ chi2W3(i__),W4mag(i__),e_W4mag(i__),snr4(i__),chi2W4(i__),
+ chi2(i__),nb(i__),na(i__),sat1(i__),sat2(i__),sat3(i__),
+ sat4(i__),satnum(i__),RA_pm(i__),DE_pm(i__),e_RA_pm(i__),
+ e_DE_pm(i__),cosig_pm(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),
+ e_pmDE(i__),chi2W1_pm(i__),chi2W2_pm(i__),chi2W3_pm(i__),
+ chi2W4_pm(i__),chi2pm(i__),qpm(i__),ccf(i__),ex(i__),var(i__),
+ qph(i__),fdet(i__),fmoon(i__),nW1(i__),mW1(i__),nW2(i__),
+ mW2(i__),nW3(i__),mW3(i__),nW4(i__),mW4(i__),covW1(i__),
+ covW2(i__),covW3(i__),covW4(i__),v2Mkey(i__),Offset(i__),
+ PA(i__),Num(i__),W_Jmag(i__),e_W_Jmag(i__),W_Hmag(i__),
+ e_W_Hmag(i__),W_Kmag(i__),e_W_Kmag(i__),J2mag(i__),
+ e_J2mag(i__),H2mag(i__),e_H2mag(i__),K2mag(i__),e_K2mag(i__),
+ J2Flux(i__),e_J2Flux(i__),H2Flux(i__),e_H2Flux(i__),
+ K2Flux(i__),e_K2Flux(i__),J2snr(i__),H2snr(i__),K2snr(i__),
+ q_J2mag(i__),q_H2mag(i__),q_K2mag(i__),Counter(i__),
+ v24mag(i__),e_24mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! SFOG field YSOs: star classification table
C Format for file interpretation
2 format(
+ A8,1X,A26,1X,A16,1X,I2,1X,I2,1X,F7.4,1X,A1,I2,1X,I2,1X,F7.4,
+ 1X,I1,1X,I1,1X,I1,1X,I1,1X,I3,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,47405
read(1,'(A95)')ar__1
read(ar__1,2)
+ Field(i__),Source_1(i__),v2MASS_1(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),dIRAC(i__),
+ dWISE(i__),dIRACWISE(i__),dSMOG(i__),Seq(i__),Class(i__)
if(ar__1(91:93) .EQ. '') Seq(i__) = iNULL__
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,2)
+ Field(i__),Source_1(i__),v2MASS_1(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),dIRAC(i__),
+ dWISE(i__),dIRACWISE(i__),dSMOG(i__),Seq(i__),Class(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 'table3.dat' ! SFOG Field YSOs: SIMBAD Matches within 2"
C Format for file interpretation
3 format(
+ A8,1X,A26,1X,F10.6,1X,F8.4,1X,A34,1X,A19,1X,I2,1X,I2,1X,F7.4,
+ 1X,A1,I2,1X,F4.1,1X,F6.3,1X,A675)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,7039
read(1,'(A816)')ar__2
read(ar__2,3)
+ Field_1(i__),IRAC(i__),RAIRACdeg(i__),DEIRACdeg(i__),
+ Simbad(i__),Type(i__),RASIMh(i__),RASIMm(i__),RASIMs(i__),
+ DESIM_(i__),DESIMd(i__),DESIMm(i__),DESIMs(i__),IDs(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(RASIMh(i__) .GT. -180) RAdeg_1(i__)=RASIMh(i__)*15.
if(RASIMm(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RASIMm(i__)/4.
if(RASIMs(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RASIMs(i__)/240.
if(DESIMd(i__) .GE. 0) DEdeg_1(i__)=DESIMd(i__)
if(DESIMm(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DESIMm(i__)/60.
if(DESIMs(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DESIMs(i__)/3600.
if(DESIM_(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,3)
+ Field_1(i__),IRAC(i__),RAIRACdeg(i__),DEIRACdeg(i__),
+ Simbad(i__),Type(i__),RASIMh(i__),RASIMm(i__),RASIMs(i__),
+ DESIM_(i__),DESIMd(i__),DESIMm(i__),DESIMs(i__),IDs(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 'table4.dat' ! SFOG YSO cluster descriptions
C Format for file interpretation
4 format(I3,1X,A12,1X,I4,1X,F10.6,1X,F8.4,1X,F9.7)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,621
read(1,'(A51)')ar__3
read(ar__3,4)
+ Seq_1(i__),GLIMPSE(i__),YSO(i__),RAdeg_2(i__),DEdeg_2(i__),
+ Rad(i__)
c ..............Just test output...........
write(6,4)
+ Seq_1(i__),GLIMPSE(i__),YSO(i__),RAdeg_2(i__),DEdeg_2(i__),
+ Rad(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 't4wise.dat' ! SFOG YSO clusters: WISE HII counterparts
* (reformatted from Table 4 by CDS)
C Format for file interpretation
5 format(I3,1X,A12,1X,A15)
C Effective file loading
open(unit=1,status='old',file=
+'t4wise.dat')
write(6,*) '....Loading file: t4wise.dat'
do i__=1,870
read(1,'(A32)')ar__4
read(ar__4,5)Seq_2(i__),GLIMPSE_1(i__),WISE(i__)
c ..............Just test output...........
write(6,5)Seq_2(i__),GLIMPSE_1(i__),WISE(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 't4skiff.dat' ! SFOG YSO clusters: O & B objects from Skiff+2014
* (reformatted from Table 4 by CDS)
C Format for file interpretation
6 format(I3,1X,A12,1X,A26)
C Effective file loading
open(unit=1,status='old',file=
+'t4skiff.dat')
write(6,*) '....Loading file: t4skiff.dat'
do i__=1,6768
read(1,'(A43)')ar__5
read(ar__5,6)Seq_3(i__),GLIMPSE_2(i__),Skiff(i__)
c ..............Just test output...........
write(6,6)Seq_3(i__),GLIMPSE_2(i__),Skiff(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 't4simbad.dat' ! SFOG YSO clusters: objects from SIMBAD
* (reformatted from Table 4 by CDS)
C Format for file interpretation
7 format(I3,1X,A12,1X,A34)
C Effective file loading
open(unit=1,status='old',file=
+'t4simbad.dat')
write(6,*) '....Loading file: t4simbad.dat'
do i__=1,38445
read(1,'(A51)')ar__6
read(ar__6,7)Seq_4(i__),GLIMPSE_3(i__),Simbad_1(i__)
c ..............Just test output...........
write(6,7)Seq_4(i__),GLIMPSE_3(i__),Simbad_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! SFOG Field YSOs: SEDFitter Results
C Format for file interpretation
8 format(
+ A8,1X,I3,1X,F4.1,1X,F4.1,1X,I1,1X,E9.4,1X,E9.4,1X,E9.4,1X,
+ F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,
+ 1X,F6.3,1X,F6.3,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,
+ E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,
+ 1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,I5,
+ 1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,E9.4,1X,
+ E9.4,1X,E9.4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,6234
read(1,'(A420)')ar__7
read(ar__7,8)
+ Field_2(i__),Seq_5(i__),Dist(i__),RadGal(i__),Ndata(i__),
+ Chi2_1(i__),e_Chi2(i__),E_Chi2_1(i__),Av(i__),e_Av(i__),
+ E_Av_1(i__),Av_avg(i__),s_Av(i__),Mass(i__),e_Mass(i__),
+ E_Mass_1(i__),Mass_avg(i__),s_Mass(i__),Age(i__),e_Age(i__),
+ E_Age_1(i__),Age_avg(i__),s_Age(i__),Mdisk(i__),e_Mdisk(i__),
+ E_Mdisk_1(i__),Mdisk_avg(i__),s_Mdisk(i__),Menv(i__),
+ e_Menv(i__),E_Menv_1(i__),Menv_avg(i__),s_Menv(i__),Mdot(i__),
+ e_Mdot(i__),E_Mdot_1(i__),Mdot_avg(i__),s_Mdot(i__),Temp(i__),
+ e_Temp(i__),E_Temp_1(i__),Temp_avg(i__),s_Temp(i__),
+ Lstar(i__),e_Lstar(i__),E_Lstar_1(i__),Lstar_avg(i__),
+ s_Lstar(i__)
if(ar__7(276:284) .EQ. '') Mdot(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ Field_2(i__),Seq_5(i__),Dist(i__),RadGal(i__),Ndata(i__),
+ Chi2_1(i__),e_Chi2(i__),E_Chi2_1(i__),Av(i__),e_Av(i__),
+ E_Av_1(i__),Av_avg(i__),s_Av(i__),Mass(i__),e_Mass(i__),
+ E_Mass_1(i__),Mass_avg(i__),s_Mass(i__),Age(i__),e_Age(i__),
+ E_Age_1(i__),Age_avg(i__),s_Age(i__),Mdisk(i__),e_Mdisk(i__),
+ E_Mdisk_1(i__),Mdisk_avg(i__),s_Mdisk(i__),Menv(i__),
+ e_Menv(i__),E_Menv_1(i__),Menv_avg(i__),s_Menv(i__),Mdot(i__),
+ e_Mdot(i__),E_Mdot_1(i__),Mdot_avg(i__),s_Mdot(i__),Temp(i__),
+ e_Temp(i__),E_Temp_1(i__),Temp_avg(i__),s_Temp(i__),
+ Lstar(i__),e_Lstar(i__),E_Lstar_1(i__),Lstar_avg(i__),
+ s_Lstar(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end