Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/160/68 into FORTRAN code for reading data files line by line.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Oct-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 ! Source designation (designation) (1)
character*16 v2MASS ! 2MASS-PSC source identifier
* (HHMMSSss+DDMMSSs; J2000) (II/246)
* (tmass_designation)
integer*8 v2MASS_cntr ! [68911969/1326367653]? 2MASS-PSC
* unique source identifier (II/246)
* (tmass_cntr)
real*8 GLON_2M ! (deg) ? 2MASS-PSC galactic longitude (II/246)
* (l)
real*8 GLAT_2M ! (deg) [-3.1/3.4]? 2MASS-PSC galactic
* latitude (II/246) (b)
real*4 e_GLON_2M ! (arcsec) [0.3]? Uncertainty in GLON_2M (dl)
real*4 e_GLAT_2M ! (arcsec) [0.3]? Uncertainty in GLAT_2M (db)
real*8 RA2Mdeg ! (deg) ? 2MASS-PSC right ascension (J2000) in
* degrees (II/246) (ra)
real*8 DE2Mdeg ! (deg) [-47/66]? 2MASS-PSC declination
* (J2000) in degrees (II/246) (dec)
real*4 e_RA2Mdeg ! (arcsec) [0.3]? Uncertainty in RA2Mdeg (dra)
real*4 e_DE2Mdeg ! (arcsec) [0.3]? Uncertainty in DE2Mdeg (ddec)
integer*4 csf ! [0/5]? Close source flag (csf) (2)
real*4 Jmag ! (mag) [6.6/18.2]? 2MASS-PSC J-band (1.25um)
* magnitude (II/246) (mag_J)
real*4 e_Jmag ! (mag) [0.01/1.1]? 1{sigma} uncertainty in
* Jmag (dJ_m)
real*4 Hmag ! (mag) [6.6/16.7]? 2MASS-PSC H-band (1.65um)
* magnitude (II/246) (mag_H)
real*4 e_Hmag ! (mag) [0.01/0.3]? 1{sigma} uncertainty in
* Hmag (dH_m)
real*4 Kmag ! (mag) [6.2/16.4]? 2MASS-PSC Ks-band (2.17um)
* magnitude (II/246) (mag_Ks)
real*4 e_Kmag ! (mag) [0.01/10]? 1{sigma} uncertainty in
* Kmag (dKs_m)
real*4 v3_6mag ! (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 ! (mag) [0.01/0.22]? 1{sigma} uncertainty in
* 3.6mag (d3_6m)
real*4 v4_5mag ! (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 ! (mag) [0.01/0.21]? 1{sigma} uncertainty in
* 4.5mag (d4_5m)
real*4 v5_8mag ! (mag) [3.3/15]? Spitzer/IRAC 5.8um (Band 3)
* magnitude (J/ApJ/880/9, II/293)
* (mag5_8)
real*4 e_5_8mag ! (mag) [0.01/0.21]? 1{sigma} uncertainty in
* 5.8mag (d5_8m)
real*4 v8_0mag ! (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 ! (mag) [0.01/0.22]? 1{sigma} uncertainty in
* 8.0mag (d8_0m)
real*8 F_J ! (mJy) [0.08/3346]? 2MASS-PSC J-band (1.25um)
* flux (II/246, II/293) (f_J)
real*8 e_F_J ! (mJy) [0.02/71]? 1{sigma} uncertainty in
* F(J) (df_J)
real*8 F_H ! (mJy) [0.2/2206]? 2MASS-PSC H-band (1.65um)
* flux (II/246, II/293) (f_H)
real*8 e_F_H ! (mJy) [0.03/75.3]? 1{sigma} uncertainty in
* F(H) (df_H)
real*8 F_K ! (mJy) [0.19/2130]? 2MASS-PSC Ks-band
* (2.17um) flux (II/246, II/293) (f_Ks)
real*8 e_F_K ! (mJy) [0.04/5181]? 1{sigma} uncertainty in
* F(K) (df_Ks)
real*8 F_3_6 ! (mJy) [0.029/10050]? Spitzer/IRAC 3.6um
* (Band 1) flux (II/293) (f3_6)
real*8 e_F_3_6 ! (mJy) [3e-3/1319]? 1{sigma} uncertainty in
* F(3.6) (df3_6)
real*8 F_4_5 ! (mJy) [0.03/6195]? Spitzer/IRAC 4.5um
* (Band 2) flux (II/293) (f4_5)
real*8 e_F_4_5 ! (mJy) [3e-3/279]? 1{sigma} uncertainty in
* F(4.5) (df4_5)
real*8 F_5_8 ! (mJy) [0.1/5421]? Spitzer/IRAC 5.8um
* (Band 3) flux (II/293) (f5_8)
real*8 e_F_5_8 ! (mJy) [0.01/125]? 1{sigma} uncertainty in
* F(5.8) (df5_8)
real*8 F_8_0 ! (mJy) [0.1/3580]? Spitzer/IRAC 8.0um
* (Band 4) flux (II/293) (f8_0)
real*8 e_F_8_0 ! (mJy) [0.01/103]? 1{sigma} uncertainty in
* F(8.0) (df8_0)
real*4 rms_3_6 ! (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 ! (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 ! (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 ! (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 ! (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 ! (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 ! (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 ! (MJy/sr) [2.3/134]? Local sky background for
* Spitzer/IRAC 8.0um (Band 4) (II/293)
* (sky\_8\_0)
real*4 Jsn ! [1/63.9]? 2MASS-PSC J-band (1.25um)
* signal-to-noise ratio (II/246) (sn_J)
real*4 Hsn ! [3.7/77.6]? 2MASS-PSC H-band (1.65um)
* signal-to-noise ratio (II/246) (sn_H)
real*4 Ksn ! [0.1/83.6]? 2MASS-PSC Ks-band (2.17um)
* signal-to-noise ratio (II/246) (sn_Ks)
real*4 v3_6sn ! [5/71.8]? Spitzer/IRAC 3.6um (Band 1)
* signal-to-noise ratio (II/293)
* (sn\_3\_6)
real*4 v4_5sn ! [5.3/87.9]? Spitzer/IRAC 4.5um (Band
* 2) signal-to-noise ratio (II/293)
* (sn\_4\_5)
real*4 v5_8sn ! [5.2/76.2]? Spitzer/IRAC 5.8um (Band
* 3) signal-to-noise ratio (II/293)
* (sn\_5\_8)
real*4 v8_0sn ! [5.1/86]? Spitzer/IRAC 8.0um (Band 4)
* signal-to-noise ratio (II/293)
* (sn\_8\_0)
real*4 d_3_6 ! (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 ! (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 ! (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 ! (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 ! [0/17]? Number of detections for
* Spitzer/IRAC 3.6um (Band 1) (II/293)
* (m3_6)
integer*4 N_4_5 ! [0/15]? Number of detections for
* Spitzer/IRAC 4.5um (Band 2) (II/293)
* (m4_5)
integer*4 N_5_8 ! [0/16]? Number of detections for
* Spitzer/IRAC 5.8um (Band 3) (II/293)
* (m5_8)
integer*4 N_8_0 ! [0/20]? Number of detections for
* Spitzer/IRAC 8.0um (Band 4) (II/293)
* (m8_0)
integer*4 n_3_6_1 ! [0/24]? Number of possible detections
* for Spitzer/IRAC 3.6um (Band 1)
* (II/293) (n3_6)
integer*4 n_4_5_1 ! [1/16]? Number of possible detections
* for Spitzer/IRAC 4.5um (Band 2)
* (II/293) (n4_5)
integer*4 n_5_8_1 ! [0/18]? Number of possible detections
* for Spitzer/IRAC 5.8um (Band 3)
* (II/293) (n5_8)
integer*4 n_8_0_1 ! [0/20]? Number of possible detections
* for Spitzer/IRAC 8.0um (Band 4)
* (II/293) (n8_0)
integer*4 q_Jmag ! ? 2MASS-PSC J-band (1.25um) source
* quality flag (II/246, II/293) (sqf_J)
* (3)
integer*4 q_Hmag ! ? 2MASS-PSC H-band (1.65um) source
* quality flag (II/246, II/293) (sqf_H)
* (3)
integer*4 q_Kmag ! ? 2MASS-PSC Ks-band (2.17um) source
* quality flag (II/246, II/293) (sqf_Ks)
* (3)
integer*8 q_3_6mag ! ? Spitzer/IRAC 3.6um (Band 1) source
* quality flag (II/293) (sqf\_3\_6) (3)
integer*8 q_4_5mag ! ? Spitzer/IRAC 4.5um (Band 2) source
* quality flag (II/293) (sqf\_4\_5) (3)
integer*8 q_5_8mag ! ? Spitzer/IRAC 5.8um (Band 3) source
* quality flag (II/293) (sqf\_5\_8) (3)
integer*8 q_8_0mag ! ? Spitzer/IRAC 8.0um (Band 4) source
* quality flag (II/293) (sqf\_8\_0) (3)
integer*8 m_3_6 ! ? Spitzer/IRAC 3.6um (Band 1) flux
* calculation method flag (II/293)
* (mf3_6)
integer*8 m_4_5 ! ? Spitzer/IRAC 4.5um (Band 2) flux
* calculation method flag (II/293)
* (mf4_5)
integer*8 m_5_8 ! ? Spitzer/IRAC 5.8um (Band 3) flux
* calculation method flag (II/293)
* (mf5_8)
integer*8 m_8_0 ! ? Spitzer/IRAC 8.0um (Band 4) flux
* calculation method flag (II/293)
* (mf8_0)
real*8 RAG360deg ! (deg) ? GLIMPSE360 right ascension in
* degrees (J2000) (ra_1)
real*8 DEG360deg ! (deg) [-46.5/65.5]? GLIMPSE360 declination
* in degrees (J2000) (dec_1)
character*19 WISEA ! AllWISE source identifier
* (JHHMMSS.ss+DDMMSS.s; J2000) (II/328)
* (designation_2) (12)
real*8 RAAWdeg ! (deg) ? AllWISE right ascension in degrees
* (J2000) (II/328) (ra_2)
real*8 DEAWdeg ! (deg) [-46.5/65.5]? AllWISE declination in
* degrees (J2000) (II/328) (dec_2)
real*4 RAsig ! (arcsec) [1e-3/0.6]? Uncertainty in RA (sigra)
real*4 DEsig ! (arcsec) [1e-3/0.6215]? Uncertainty in DE
* (sigdec)
real*8 RADEsig ! (arcsec) [-0.3/0.2]? Cross-term of RA and DE
* uncertainties (sigradec)
real*8 GLON_AW ! (deg) ? AllWISE galactic longitude (II/328)
* (glon)
real*8 GLAT_AW ! (deg) [-3.2/3.6]? AllWISE galactic latitude
* (II/328) (glat)
real*8 ELON_AW ! (deg) ? AllWISE ecliptic longitude (II/328)
* (elon)
real*8 ELAT_AW ! (deg) [-61.4/63.7]? AllWISE ecliptic
* latitude (II/328) (elat)
real*8 wx ! (pix) [46.5/4058.4]? X-pixel coordinate on
* the Atlas Image from AllWISE, all
* bands (II/328)
real*8 wy ! (pix) [38.6/4058.1]? Y-pixel coordinate on
* the Atlas Image from AllWISE, all
* bands (II/328)
real*4 W1mag ! (mag) [5/16]? WISE W1 band (3.35um)
* magnitude (II/328) (w1mpro)
real*4 e_W1mag ! (mag) [7e-3/0.52]? Mean error in W1mag
* (w1sigmpro)
real*4 snr1 ! [2.1/157]? Signal to noise ratio for
* W1 filter (II/328) (w1snr)
real*8 chi2W1 ! [0.02/2583]? Reduced {chi}^2^ of the
* W1 profile fit (w1rchi2)
real*4 W2mag ! (mag) [2.9/17]? WISE W2 band (4.6um)
* magnitude (II/328) (w2mpro)
real*4 e_W2mag ! (mag) [7e-3/0.55]? Mean error in W2mag
* (w2sigmpro)
real*4 snr2 ! [2/156]? Signal to noise ratio for W2
* filter (II/328) (w2snr)
real*8 chi2W2 ! [0.04/1345]? Reduced {chi}^2^ of the
* W2 profile fit (II/328) (w2rchi2)
real*8 W3mag ! (mag) [-0.3/10]? WISE W3 band (11.6um)
* magnitude (II/328) (w3mpro)
real*4 e_W3mag ! (mag) [4e-3/0.6]? Mean error in W3mag
* (w3sigmpro)
real*4 snr3 ! [2/274]? Signal to noise ratio for W3
* filter (II/328) (w3snr)
real*8 chi2W3 ! [0.01/3133]? Reduced {chi}^2^ of the
* W3 profile fit (II/328) (w3rchi2)
real*4 W4mag ! (mag) [-4/9.7]? WISE W4 (22.1um) band
* magnitude (II/328) (w4mpro)
real*4 e_W4mag ! (mag) [1e-3/0.6]? Mean error in W4mag
* (w4sigmpro)
real*4 snr4 ! [-120/1307]? Signal to noise ratio for
* W4 filter (II/328) (w4snr)
real*8 chi2W4 ! [0.02/10340]? Reduced {chi}^2^ of the
* W4 profile fit (II/328) (w4rchi2)
real*8 chi2 ! [0.1/1908]? All-band combined reduced
* {chi}^2^ (II/328) (rchi2)
integer*4 nb ! [1/4]? Number of PSF components used
* simultaneously in the profile-fitting
* (II/328)
integer*4 na ! [0/1]? Active deblending flag (1 if a
* detection was split into multiple
* sources) (II/328)
real*4 sat1 ! [0/0.67]? Saturated pixel fraction in
* W1 (II/328) (w1sat)
real*4 sat2 ! [0/0.71]? Saturated pixel fraction in
* W2 (II/328) (w2sat)
real*4 sat3 ! [0/0.32]? Saturated pixel fraction in
* W3 (II/328) (w3sat)
real*4 sat4 ! [0/0.48]? Saturated pixel fraction in
* W4 (II/328) (w4sat)
character*4 satnum ! Minimum sample at which saturation
* occurs in each band, 1 character per
* band (II/328)
real*8 RA_pm ! (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 ! (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 ! (arcsec) [2e-3/4.5]? 1{sigma} uncertainty on
* RA_pm (sigra_pm)
real*4 e_DE_pm ! (arcsec) [2e-3/2.8]? 1{sigma} uncertainty on
* DE_pm (sigdec_pm)
real*8 cosig_pm ! (arcsec) [-1.9/2.3]? Co-sigma of the equatorial
* position uncertainties (II/328)
* (sigradec_pm)
integer*4 pmRA ! (mas/yr) [-98190/130293]? Apparent motion in RA
* (II/328) (pmra) (4)
integer*4 e_pmRA ! (mas/yr) [20/11196]? Mean error on pmRA
* (sigpmra)
integer*4 pmDE ! (mas/yr) [-198842/169468]? Apparent motion in
* DE (II/328) (pmdec) (4)
integer*4 e_pmDE ! (mas/yr) [20/12961]? Mean error on pmDE
* (sigpmdec)
real*8 chi2W1_pm ! [0.04/2535]? Reduced {chi}^2^ of the
* W1 profile fit including motion
* estimation (II/328) (w1rchi2_pm)
real*8 chi2W2_pm ! [0.05/1335]? Reduced {chi}^2^ of the
* W2 profile fit including motion
* estimation (II/328) (w2rchi2_pm)
real*8 chi2W3_pm ! [0.01/3107]? Reduced {chi}^2^ of the
* W3 profile fit including motion
* estimation (II/328) (w3rchi2_pm)
real*8 chi2W4_pm ! [0.02/10150]? Reduced {chi}^2^ of the
* W4 profile fit including motion
* estimation (II/328) (w4rchi2_pm)
real*8 chi2pm ! [0.1/1905]? Combined reduced {chi}^2^
* profile fit including motion
* estimation (II/328) (rchi2_pm)
character*5 qpm ! Motion estimation quality (II/328)
* (pmcode) (5)
character*4 ccf ! Contamination and confusion flag, one
* per band (II/328) (cc_flags) (6)
integer*4 ex ! [0/5]? Extended source flag (II/328)
* (ext_flg) (7)
character*4 var ! Variability flag, one per band (II/328)
* (var_flg) (8)
character*4 qph ! Photometric quality flag (II/328)
* (ph_qual) (9)
integer*4 fdet ! [7/15]? Bit-encoded integer indicating
* bands in which a source has a snr?>2
* detection (II/328) (det_bit) (13)
character*4 fmoon ! Scattered moonlight contamination
* flag, 1 per band (II/328) (moon_lev)
integer*4 nW1 ! [2/81]? Frame detection count in W1
* (II/328) (w1nm) (10)
integer*4 mW1 ! [16/81]? Integer frame coverage in W1
* (II/328) (w1m) (11)
integer*4 nW2 ! [0/81]? Frame detection count in W2
* (II/328) (w2nm) (10)
integer*4 mW2 ! [12/81]? Integer frame coverage in W2
* (II/328) (w2m) (11)
integer*4 nW3 ! [0/55]? Frame detection count in W3
* (II/328) (w3nm) (10)
integer*4 mW3 ! [0/57]? Integer frame coverage in W3
* (II/328) (w3m) (11)
integer*4 nW4 ! [0/54]? Frame detection count in W4
* (II/328) (w4nm) (10)
integer*4 mW4 ! [0/57]? Integer frame coverage in W4
* (II/328) (w4m) (11)
real*4 covW1 ! [5.5/83.1]? Mean pixel coverage in W1
* (II/328) (w1cov)
real*4 covW2 ! [2/83.1]? Mean pixel coverage in W2
* (II/328) (w2cov)
real*4 covW3 ! [0.4/60.6]? Mean pixel coverage in W3
* (II/328) (w3cov)
real*4 covW4 ! [1.4/60.3]? Mean pixel coverage in W4
* (II/328) (w4cov)
integer*8 v2Mkey ! [68916543/1315293861]? 2MASS-PSC
* closest associated key
* (not identification) (II/328, II/246)
* (tmass_key)
real*4 Offset ! (arcsec) [3e-3/3]? Radial offset between WISE
* source (II/328) and associated 2MASS
* source (II/246) (r_2mass)
integer*4 PA ! (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 ! [0/3]? Number of 2MASS-PSC entries
* found within search radius of WISE
* source (n_2mass)
real*4 W_Jmag ! (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 ! (mag) [0.01/0.7]? Mean error on W.Jmag
* (j_msig_2mass)
real*4 W_Hmag ! (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 ! (mag) [0.01/10]? Mean error on W.Hmag
* (h_msig_2mass)
real*4 W_Kmag ! (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 ! (mag) [0.01/10]? Mean error on W.Kmag
* (k_msig_2mass)
real*4 J2mag ! (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 ! (mag) [0.01/0.6]? Mean error on J2mag (dj2_m)
real*4 H2mag ! (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 ! (mag) [0.01/0.3]? Mean error on H2mag (dh2_m)
real*4 K2mag ! (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 ! (mag) [0.01/10]? Mean error on K2mag (dks2_m)
real*8 J2Flux ! (mJy) [0.1/8882]? 2MASS-PSC J-band (1.25um)
* flux from SMOG catalog/archive (f_j2)
real*8 e_J2Flux ! (mJy) [0.03/207.1]? 1{sigma} uncertainty on
* J2Flux (df_j2)
real*8 H2Flux ! (mJy) [0.3/14990]? 2MASS-PSC H-band (1.65um)
* flux from SMOG catalog/archive (f_h2)
real*8 e_H2Flux ! (mJy) [0.05/3148]? 1{sigma} uncertainty on
* H2Flux (df_h2)
real*8 K2Flux ! (mJy) [0.3/16790]? 2MASS-PSC Ks-band
* (2.17um) flux from SMOG
* catalog/archive (f_ks2)
real*8 e_K2Flux ! (mJy) [0.06/19200]? 1{sigma} uncertainty on
* K2mag (df_ks2)
real*4 J2snr ! [2/60.4]? 2MASS-PSC J-band (1.25um)
* signal-to-noise ratio from SMOG
* catalog/archive (sn_j2)
real*4 H2snr ! [4/63.9]? 2MASS-PSC H-band (1.65um)
* signal-to-noise ratio from SMOG
* catalog/archive (sn_h2)
real*4 K2snr ! [0.1/72.4]? 2MASS-PSC Ks-band (2.17um)
* signal-to-noise ratio from SMOG
* catalog/archive (sn_ks2)
integer*8 q_J2mag ! ? 2MASS-PSC J-Band (1.25um) source
* quality flag from SMOG catalog/archive
* (sqf_j2) (3)
integer*8 q_H2mag ! ? 2MASS-PSC H-Band (1.65um) source
* quality flag from SMOG catalog/archive
* (sqf_h2) (3)
integer*8 q_K2mag ! ? 2MASS-PSC Ks-Band (2.17um) source
* quality flag from SMOG catalog/archive
* (sqf_ks2) (3)
integer*4 Counter ! [23/2836315]? Entry counter (key)
* number (unique within table) (cntr)
real*8 v24mag ! (mag) [-100/11.4]? Spitzer/MIPS 24um (Band
* 1) magnitude from SMOG catalog/archive
* (mag_24)
real*8 e_24mag ! (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 ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*8 Field ! SMOG Field ID (SRC0 to SRC47404) (ID)
character*26 Source_1 ! Source designation (IR_name)
character*16 v2MASS_1 ! 2MASS identifier (2M name)
integer*4 RAh ! (h) Hour of right ascension (J2000)
integer*4 RAm ! (min) Minute of right ascension (J2000)
real*8 RAs ! (s) Second of right ascension (J2000)
character*1 DE_ ! Sign of declination (J2000)
integer*4 DEd ! (deg) Degree of declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of declination (J2000)
real*8 DEs ! (arcsec) Arcsecond of declination (J2000)
integer*4 dIRAC ! [0/1] Identified as a YSO based on IRAC
* data set (1=yes) (dIRAC)
integer*4 dWISE ! [0/1] Identified as a YSO based on the WISE
* data set II/328 (1=yes) (dWISE)
integer*4 dIRACWISE ! [0/1] Identified as a YSO based on the
* combined IRAC+WISE data set (1=yes)
* (dIRACWISE) (1)
integer*4 dSMOG ! [0/1] Identified as a YSO based on SMOG
* data set (1=yes) (dSMOG)
integer*4 Seq ! [-1/620]? Cluster number (Cluster) (2)
integer*4 Class ! [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 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*8 Field_1 ! SMOG Field ID (SRC9 to SRC47353) (SRC)
character*26 IRAC ! IRAC identifier (IRAC_ID)
real*8 RAIRACdeg ! (deg) Right ascension from IRAC (J2000) (IRAC_RA)
real*8 DEIRACdeg ! (deg) Declination from IRAC (J2000) (IRAC_Dec)
character*34 Simbad ! SIMBAD identifier (SIMBAD_ID)
character*19 Type ! Object type (OTYPE)
integer*4 RASIMh ! (h) Hour of right ascension from SIMBAD (J2000)
* (SIMBAD_RA)
integer*4 RASIMm ! (min) Minute of right ascension from Simbad
* (J2000)
real*8 RASIMs ! (s) Second of right ascension from Simbad
* (J2000)
character*1 DESIM_ ! Sign of declination from SIMBAD (J2000)
* (SIMBAD_Dec)
integer*4 DESIMd ! (deg) Degree of declination from SIMBAD (J2000)
real*4 DESIMm ! (arcmin) Arcminute of declination from SIMBAD (J2000)
real*4 DESIMs ! (arcsec) Arcsecond of declination from SIMBAD (J2000)
character*675 IDs ! 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 ! [0/620] Cluster number (Cluster_#)
character*12 GLIMPSE ! GLIMPSE identifier (GLLL.llll+BB.bbbb)
* (Cluster_Name)
integer*4 YSO ! [5/1177] Number of YSO in the cluster (#_YSO)
real*8 RAdeg_2 ! (deg) Right ascension of cluster center (J2000)
* (Central_RA)
real*8 DEdeg_2 ! (deg) Declination of cluster center (J2000)
* (Central_Dec)
real*8 Rad ! (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 ! [0/620] Cluster number
character*12 GLIMPSE_1 ! GLIMPSE identifier (GLLL.llll+BB.bbbb)
* (Cluster_name)
character*15 WISE ! 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 ! [0/620] Cluster number
character*12 GLIMPSE_2 ! GLIMPSE identifier (GLLL.ll+B.bb)
* (Cluster_name)
character*26 Skiff ! 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 ! [0/620] Cluster number
character*12 GLIMPSE_3 ! GLIMPSE identifier (G.LLL.ll+B.bb)
* (Cluster_name)
character*34 Simbad_1 ! 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 ! SMOG Field ID (SRC588 to SRC47400) (SRC_ID)
integer*4 Seq_5 ! [5/593] Cluster number
real*4 Dist ! (kpc) [0.7/13.3] Distance
real*4 RadGal ! (kpc) [7.8/18.1] Galactic radius
integer*4 Ndata ! [3/7] Number of SED data points (N_data)
real*4 Chi2_1 ! [7.2e-7/3] {chi}^2^
real*4 e_Chi2 ! [7.2e-7/3] Lower limit on {chi}^2^
* (Chi_2_l)
real*4 E_Chi2_1 ! [0.67/24] Upper limit on {chi}^2^ (Chi_2_h)
real*4 Av ! (mag) [0/40] Extinction A_V_
real*4 e_Av ! (mag) [0/24.3] Lower limit on A_V_ (Av_l)
real*4 E_Av_1 ! (mag) [0.3/40] Upper limit on A_V_ (Av_h)
real*4 Av_avg ! (mag) [0/40] Average A_V_ (Av_Ave) (1)
real*4 s_Av ! (mag) [0/16] Standard deviation on A_V_ (Av_SD)
real*4 Mass ! (Msun) [0.1/31.2] Object mass (Mass)
real*4 e_Mass ! (Msun) [0.1/17.6] Lower limit on Mass (Mass_l)
real*4 E_Mass_1 ! (Msun) [1.3/47.3] Upper limit on Mass (Mass_h)
real*4 Mass_avg ! (Msun) [0.1/20.4] Average of Mass (Mass_Ave) (1)
real*4 s_Mass ! (Msun) [0/10.4] Standard deviation on Mass
* (Mass_SD)
real*4 Age ! (yr) [1015/1e+7] Age of the object
real*4 e_Age ! (yr) [1008/5.5e+6] Lower limit on Age (Age_l)
real*4 E_Age_1 ! (yr) [2210/1e+7] Upper limit on Age (Age_h)
real*4 Age_avg ! (yr) [1015/1e+7] Average of Age (Age_Ave) (1)
real*4 s_Age ! (yr) [0/4.8e+06] Standard deviation on Age
* (Age_SD)
real*4 Mdisk ! (Msun) [0/0.96] Disk mass M_disk_ (Disk_Mass)
real*4 e_Mdisk ! (Msun) [0/0.21] Lower limit on M_disk_
* (Disk-Mass-l)
real*4 E_Mdisk_1 ! (Msun) [0/1.9] Upper limit on M_disk_
* (Disk-mass-h)
real*4 Mdisk_avg ! (Msun) [0/0.9] Average of M_disk_ (Disk-Mass-Ave)
* (1)
real*4 s_Mdisk ! (Msun) [0/0.45] Standard deviation on M_disk_
* (Disk-Mass-SD)
real*4 Menv ! (Msun) [6.5e-10/3773] Envelope mass M_env_
* (Envl-Mass)
real*4 e_Menv ! (Msun) [1.7e-10/482] Lower limit on M_env_
* (Envl-Mass-l)
real*4 E_Menv_1 ! (Msun) [4.3e-07/4351] Upper limit on M_env_
* (Envl-mass-h)
real*4 Menv_avg ! (Msun) [1.4e-09/3773] Average of M_env_
* (Envl-Mass-Ave) (1)
real*4 s_Menv ! (Msun) [0/1009] Standard deviation on M_env_
* (Envl-Mass-SD)
real*4 Mdot ! (Msun/yr) [3.5e-10/7e-3]? Disk accretion rate Mdot
* (Mdot_star)
real*4 e_Mdot ! (Msun/yr) [0/1.4e-3] Lower limit on Mdot
* (Mdot-star-l)
real*4 E_Mdot_1 ! (Msun/yr) [0/1e-2] Upper limit on Mdot (Mdot-star-h)
real*4 Mdot_avg ! (Msun/yr) [0/7e-3] Average of Mdot (Mdot-star-Ave)
* (1)
real*4 s_Mdot ! (Msun/yr) [0/2.4e-3] Standard deviation on Mdot
* (Mdot-star-SD)
integer*4 Temp ! (K) [2550/40780] Central temperature T_*_
* (Temp)
real*4 e_Temp ! (K) [2536/26150] Lower limit on T_*_ (Temp_l)
real*4 E_Temp_1 ! (K) [4024/45160] Upper limit on T_*_ (Temp_h)
real*4 Temp_avg ! (K) [2651/35340] Average of T_*_ (Temp_Ave) (1)
real*4 s_Temp ! (K) [0/13230] Standard deviation on T_*_
* (Temp_SD)
real*4 Lstar ! (Lsun) [0.05/138900] Total luminosity L_*_
* (Tot_Lum)
real*4 e_Lstar ! (Lsun) [0.04/21500] Lower limit on L_*_
* (Tot-Lum-l)
real*4 E_Lstar_1 ! (Lsun) [1.5/327300] Upper limit on L_*_
* (Tot-Lum-h)
real*4 Lstar_avg ! (Lsun) [0.09/47920] Average of L_*_ (Tot-Lum-Ave)
* (1)
real*4 s_Lstar ! (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,v2MASS,v2MASS_cntr,GLON_2M,GLAT_2M,e_GLON_2M,e_GLAT_2M,
+ RA2Mdeg,DE2Mdeg,e_RA2Mdeg,e_DE2Mdeg,csf,Jmag,e_Jmag,Hmag,
+ e_Hmag,Kmag,e_Kmag,v3_6mag,e_3_6mag,v4_5mag,e_4_5mag,v5_8mag,
+ e_5_8mag,v8_0mag,e_8_0mag,F_J,e_F_J,F_H,e_F_H,F_K,e_F_K,F_3_6,
+ e_F_3_6,F_4_5,e_F_4_5,F_5_8,e_F_5_8,F_8_0,e_F_8_0,rms_3_6,
+ rms_4_5,rms_5_8,rms_8_0,sky_3_6,sky_4_5,sky_5_8,sky_8_0,Jsn,
+ Hsn,Ksn,v3_6sn,v4_5sn,v5_8sn,v8_0sn,d_3_6,d_4_5,d_5_8,d_8_0,
+ N_3_6,N_4_5,N_5_8,N_8_0,n_3_6_1,n_4_5_1,n_5_8_1,n_8_0_1,
+ q_Jmag,q_Hmag,q_Kmag,q_3_6mag,q_4_5mag,q_5_8mag,q_8_0mag,
+ m_3_6,m_4_5,m_5_8,m_8_0,RAG360deg,DEG360deg,WISEA,RAAWdeg,
+ DEAWdeg,RAsig,DEsig,RADEsig,GLON_AW,GLAT_AW,ELON_AW,ELAT_AW,
+ wx,wy,W1mag,e_W1mag,snr1,chi2W1,W2mag,e_W2mag,snr2,chi2W2,
+ W3mag,e_W3mag,snr3,chi2W3,W4mag,e_W4mag,snr4,chi2W4,chi2,nb,
+ na,sat1,sat2,sat3,sat4,satnum,RA_pm,DE_pm,e_RA_pm,e_DE_pm,
+ cosig_pm,pmRA,e_pmRA,pmDE,e_pmDE,chi2W1_pm,chi2W2_pm,
+ chi2W3_pm,chi2W4_pm,chi2pm,qpm,ccf,ex,var,qph,fdet,fmoon,nW1,
+ mW1,nW2,mW2,nW3,mW3,nW4,mW4,covW1,covW2,covW3,covW4,v2Mkey,
+ Offset,PA,Num,W_Jmag,e_W_Jmag,W_Hmag,e_W_Hmag,W_Kmag,e_W_Kmag,
+ J2mag,e_J2mag,H2mag,e_H2mag,K2mag,e_K2mag,J2Flux,e_J2Flux,
+ H2Flux,e_H2Flux,K2Flux,e_K2Flux,J2snr,H2snr,K2snr,q_J2mag,
+ q_H2mag,q_K2mag,Counter,v24mag,e_24mag
if(ar__(45:55) .EQ. '') v2MASS_cntr = iNULL__
if(ar__(57:66) .EQ. '') GLON_2M = rNULL__
if(ar__(68:76) .EQ. '') GLAT_2M = rNULL__
if(ar__(78:80) .EQ. '') e_GLON_2M = rNULL__
if(ar__(82:84) .EQ. '') e_GLAT_2M = rNULL__
if(ar__(86:95) .EQ. '') RA2Mdeg = rNULL__
if(ar__(97:106) .EQ. '') DE2Mdeg = rNULL__
if(ar__(108:110) .EQ. '') e_RA2Mdeg = rNULL__
if(ar__(112:114) .EQ. '') e_DE2Mdeg = rNULL__
if(ar__(116:116) .EQ. '') csf = iNULL__
if(ar__(118:123) .EQ. '') Jmag = rNULL__
if(ar__(125:129) .EQ. '') e_Jmag = rNULL__
if(ar__(131:136) .EQ. '') Hmag = rNULL__
if(ar__(138:142) .EQ. '') e_Hmag = rNULL__
if(ar__(144:149) .EQ. '') Kmag = rNULL__
if(ar__(151:155) .EQ. '') e_Kmag = rNULL__
if(ar__(157:162) .EQ. '') v3_6mag = rNULL__
if(ar__(164:168) .EQ. '') e_3_6mag = rNULL__
if(ar__(170:175) .EQ. '') v4_5mag = rNULL__
if(ar__(177:181) .EQ. '') e_4_5mag = rNULL__
if(ar__(183:188) .EQ. '') v5_8mag = rNULL__
if(ar__(190:194) .EQ. '') e_5_8mag = rNULL__
if(ar__(196:201) .EQ. '') v8_0mag = rNULL__
if(ar__(203:207) .EQ. '') e_8_0mag = rNULL__
if(ar__(209:218) .EQ. '') F_J = rNULL__
if(ar__(220:227) .EQ. '') e_F_J = rNULL__
if(ar__(229:236) .EQ. '') F_H = rNULL__
if(ar__(238:245) .EQ. '') e_F_H = rNULL__
if(ar__(247:255) .EQ. '') F_K = rNULL__
if(ar__(257:265) .EQ. '') e_F_K = rNULL__
if(ar__(267:277) .EQ. '') F_3_6 = rNULL__
if(ar__(279:289) .EQ. '') e_F_3_6 = rNULL__
if(ar__(291:300) .EQ. '') F_4_5 = rNULL__
if(ar__(302:311) .EQ. '') e_F_4_5 = rNULL__
if(ar__(313:321) .EQ. '') F_5_8 = rNULL__
if(ar__(323:331) .EQ. '') e_F_5_8 = rNULL__
if(ar__(333:341) .EQ. '') F_8_0 = rNULL__
if(ar__(343:351) .EQ. '') e_F_8_0 = rNULL__
if(ar__(353:357) .EQ. '') rms_3_6 = rNULL__
if(ar__(359:361) .EQ. '') rms_4_5 = iNULL__
if(ar__(363:367) .EQ. '') rms_5_8 = rNULL__
if(ar__(369:373) .EQ. '') rms_8_0 = rNULL__
if(ar__(375:381) .EQ. '') sky_3_6 = rNULL__
if(ar__(383:388) .EQ. '') sky_4_5 = rNULL__
if(ar__(390:396) .EQ. '') sky_5_8 = rNULL__
if(ar__(398:404) .EQ. '') sky_8_0 = rNULL__
if(ar__(406:410) .EQ. '') Jsn = rNULL__
if(ar__(412:416) .EQ. '') Hsn = rNULL__
if(ar__(418:422) .EQ. '') Ksn = rNULL__
if(ar__(424:428) .EQ. '') v3_6sn = rNULL__
if(ar__(430:434) .EQ. '') v4_5sn = rNULL__
if(ar__(436:440) .EQ. '') v5_8sn = rNULL__
if(ar__(442:446) .EQ. '') v8_0sn = rNULL__
if(ar__(448:452) .EQ. '') d_3_6 = rNULL__
if(ar__(454:458) .EQ. '') d_4_5 = rNULL__
if(ar__(460:463) .EQ. '') d_5_8 = rNULL__
if(ar__(465:468) .EQ. '') d_8_0 = rNULL__
if(ar__(470:471) .EQ. '') N_3_6 = iNULL__
if(ar__(473:474) .EQ. '') N_4_5 = iNULL__
if(ar__(476:477) .EQ. '') N_5_8 = iNULL__
if(ar__(479:480) .EQ. '') N_8_0 = iNULL__
if(ar__(482:483) .EQ. '') n_3_6_1 = iNULL__
if(ar__(485:486) .EQ. '') n_4_5_1 = iNULL__
if(ar__(488:489) .EQ. '') n_5_8_1 = iNULL__
if(ar__(491:492) .EQ. '') n_8_0_1 = iNULL__
if(ar__(494:501) .EQ. '') q_Jmag = iNULL__
if(ar__(503:510) .EQ. '') q_Hmag = iNULL__
if(ar__(512:519) .EQ. '') q_Kmag = iNULL__
if(ar__(521:531) .EQ. '') q_3_6mag = iNULL__
if(ar__(533:543) .EQ. '') q_4_5mag = iNULL__
if(ar__(545:555) .EQ. '') q_5_8mag = iNULL__
if(ar__(557:567) .EQ. '') q_8_0mag = iNULL__
if(ar__(569:579) .EQ. '') m_3_6 = iNULL__
if(ar__(581:591) .EQ. '') m_4_5 = iNULL__
if(ar__(593:603) .EQ. '') m_5_8 = iNULL__
if(ar__(605:615) .EQ. '') m_8_0 = iNULL__
if(ar__(617:626) .EQ. '') RAG360deg = rNULL__
if(ar__(628:637) .EQ. '') DEG360deg = rNULL__
if(ar__(659:668) .EQ. '') RAAWdeg = rNULL__
if(ar__(670:677) .EQ. '') DEAWdeg = rNULL__
if(ar__(679:684) .EQ. '') RAsig = rNULL__
if(ar__(686:691) .EQ. '') DEsig = rNULL__
if(ar__(693:699) .EQ. '') RADEsig = rNULL__
if(ar__(701:708) .EQ. '') GLON_AW = rNULL__
if(ar__(710:717) .EQ. '') GLAT_AW = rNULL__
if(ar__(719:729) .EQ. '') ELON_AW = rNULL__
if(ar__(731:738) .EQ. '') ELAT_AW = rNULL__
if(ar__(740:747) .EQ. '') wx = rNULL__
if(ar__(749:756) .EQ. '') wy = rNULL__
if(ar__(758:763) .EQ. '') W1mag = rNULL__
if(ar__(765:769) .EQ. '') e_W1mag = rNULL__
if(ar__(771:775) .EQ. '') snr1 = rNULL__
if(ar__(777:786) .EQ. '') chi2W1 = rNULL__
if(ar__(788:793) .EQ. '') W2mag = rNULL__
if(ar__(795:799) .EQ. '') e_W2mag = rNULL__
if(ar__(801:805) .EQ. '') snr2 = rNULL__
if(ar__(807:816) .EQ. '') chi2W2 = rNULL__
if(ar__(818:824) .EQ. '') W3mag = rNULL__
if(ar__(826:830) .EQ. '') e_W3mag = rNULL__
if(ar__(832:836) .EQ. '') snr3 = rNULL__
if(ar__(838:847) .EQ. '') chi2W3 = rNULL__
if(ar__(849:854) .EQ. '') W4mag = rNULL__
if(ar__(856:860) .EQ. '') e_W4mag = rNULL__
if(ar__(862:867) .EQ. '') snr4 = rNULL__
if(ar__(869:879) .EQ. '') chi2W4 = rNULL__
if(ar__(881:889) .EQ. '') chi2 = rNULL__
if(ar__(891:891) .EQ. '') nb = iNULL__
if(ar__(893:893) .EQ. '') na = iNULL__
if(ar__(895:899) .EQ. '') sat1 = rNULL__
if(ar__(901:905) .EQ. '') sat2 = rNULL__
if(ar__(907:911) .EQ. '') sat3 = rNULL__
if(ar__(913:917) .EQ. '') sat4 = rNULL__
if(ar__(924:933) .EQ. '') RA_pm = rNULL__
if(ar__(935:942) .EQ. '') DE_pm = rNULL__
if(ar__(944:949) .EQ. '') e_RA_pm = rNULL__
if(ar__(951:956) .EQ. '') e_DE_pm = rNULL__
if(ar__(958:964) .EQ. '') cosig_pm = rNULL__
if(ar__(966:971) .EQ. '') pmRA = iNULL__
if(ar__(973:977) .EQ. '') e_pmRA = iNULL__
if(ar__(979:985) .EQ. '') pmDE = iNULL__
if(ar__(987:991) .EQ. '') e_pmDE = iNULL__
if(ar__(993:1002) .EQ. '') chi2W1_pm = rNULL__
if(ar__(1004:1013) .EQ. '') chi2W2_pm = rNULL__
if(ar__(1015:1024) .EQ. '') chi2W3_pm = rNULL__
if(ar__(1026:1036) .EQ. '') chi2W4_pm = rNULL__
if(ar__(1038:1046) .EQ. '') chi2pm = rNULL__
if(ar__(1059:1059) .EQ. '') ex = iNULL__
if(ar__(1071:1072) .EQ. '') fdet = iNULL__
if(ar__(1079:1080) .EQ. '') nW1 = iNULL__
if(ar__(1082:1083) .EQ. '') mW1 = iNULL__
if(ar__(1085:1086) .EQ. '') nW2 = iNULL__
if(ar__(1088:1089) .EQ. '') mW2 = iNULL__
if(ar__(1091:1092) .EQ. '') nW3 = iNULL__
if(ar__(1094:1095) .EQ. '') mW3 = iNULL__
if(ar__(1097:1098) .EQ. '') nW4 = iNULL__
if(ar__(1100:1101) .EQ. '') mW4 = iNULL__
if(ar__(1103:1108) .EQ. '') covW1 = rNULL__
if(ar__(1110:1115) .EQ. '') covW2 = rNULL__
if(ar__(1117:1122) .EQ. '') covW3 = rNULL__
if(ar__(1124:1129) .EQ. '') covW4 = rNULL__
if(ar__(1131:1140) .EQ. '') v2Mkey = iNULL__
if(ar__(1142:1146) .EQ. '') Offset = rNULL__
if(ar__(1148:1151) .EQ. '') PA = iNULL__
if(ar__(1153:1153) .EQ. '') Num = iNULL__
if(ar__(1155:1160) .EQ. '') W_Jmag = rNULL__
if(ar__(1162:1166) .EQ. '') e_W_Jmag = rNULL__
if(ar__(1168:1173) .EQ. '') W_Hmag = rNULL__
if(ar__(1175:1179) .EQ. '') e_W_Hmag = rNULL__
if(ar__(1181:1186) .EQ. '') W_Kmag = rNULL__
if(ar__(1188:1192) .EQ. '') e_W_Kmag = rNULL__
if(ar__(1194:1199) .EQ. '') J2mag = rNULL__
if(ar__(1201:1205) .EQ. '') e_J2mag = rNULL__
if(ar__(1207:1212) .EQ. '') H2mag = rNULL__
if(ar__(1214:1218) .EQ. '') e_H2mag = rNULL__
if(ar__(1220:1225) .EQ. '') K2mag = rNULL__
if(ar__(1227:1231) .EQ. '') e_K2mag = rNULL__
if(ar__(1233:1241) .EQ. '') J2Flux = rNULL__
if(ar__(1243:1251) .EQ. '') e_J2Flux = rNULL__
if(ar__(1253:1262) .EQ. '') H2Flux = rNULL__
if(ar__(1264:1273) .EQ. '') e_H2Flux = rNULL__
if(ar__(1275:1284) .EQ. '') K2Flux = rNULL__
if(ar__(1286:1296) .EQ. '') e_K2Flux = rNULL__
if(ar__(1298:1302) .EQ. '') J2snr = rNULL__
if(ar__(1304:1308) .EQ. '') H2snr = rNULL__
if(ar__(1310:1314) .EQ. '') K2snr = rNULL__
if(ar__(1316:1326) .EQ. '') q_J2mag = iNULL__
if(ar__(1328:1338) .EQ. '') q_H2mag = iNULL__
if(ar__(1340:1350) .EQ. '') q_K2mag = iNULL__
if(ar__(1352:1358) .EQ. '') Counter = iNULL__
if(ar__(1360:1370) .EQ. '') v24mag = rNULL__
if(ar__(1372:1383) .EQ. '') e_24mag = rNULL__
c ..............Just test output...........
write(6,1)
+ Source,v2MASS,v2MASS_cntr,GLON_2M,GLAT_2M,e_GLON_2M,e_GLAT_2M,
+ RA2Mdeg,DE2Mdeg,e_RA2Mdeg,e_DE2Mdeg,csf,Jmag,e_Jmag,Hmag,
+ e_Hmag,Kmag,e_Kmag,v3_6mag,e_3_6mag,v4_5mag,e_4_5mag,v5_8mag,
+ e_5_8mag,v8_0mag,e_8_0mag,F_J,e_F_J,F_H,e_F_H,F_K,e_F_K,F_3_6,
+ e_F_3_6,F_4_5,e_F_4_5,F_5_8,e_F_5_8,F_8_0,e_F_8_0,rms_3_6,
+ rms_4_5,rms_5_8,rms_8_0,sky_3_6,sky_4_5,sky_5_8,sky_8_0,Jsn,
+ Hsn,Ksn,v3_6sn,v4_5sn,v5_8sn,v8_0sn,d_3_6,d_4_5,d_5_8,d_8_0,
+ N_3_6,N_4_5,N_5_8,N_8_0,n_3_6_1,n_4_5_1,n_5_8_1,n_8_0_1,
+ q_Jmag,q_Hmag,q_Kmag,q_3_6mag,q_4_5mag,q_5_8mag,q_8_0mag,
+ m_3_6,m_4_5,m_5_8,m_8_0,RAG360deg,DEG360deg,WISEA,RAAWdeg,
+ DEAWdeg,RAsig,DEsig,RADEsig,GLON_AW,GLAT_AW,ELON_AW,ELAT_AW,
+ wx,wy,W1mag,e_W1mag,snr1,chi2W1,W2mag,e_W2mag,snr2,chi2W2,
+ W3mag,e_W3mag,snr3,chi2W3,W4mag,e_W4mag,snr4,chi2W4,chi2,nb,
+ na,sat1,sat2,sat3,sat4,satnum,RA_pm,DE_pm,e_RA_pm,e_DE_pm,
+ cosig_pm,pmRA,e_pmRA,pmDE,e_pmDE,chi2W1_pm,chi2W2_pm,
+ chi2W3_pm,chi2W4_pm,chi2pm,qpm,ccf,ex,var,qph,fdet,fmoon,nW1,
+ mW1,nW2,mW2,nW3,mW3,nW4,mW4,covW1,covW2,covW3,covW4,v2Mkey,
+ Offset,PA,Num,W_Jmag,e_W_Jmag,W_Hmag,e_W_Hmag,W_Kmag,e_W_Kmag,
+ J2mag,e_J2mag,H2mag,e_H2mag,K2mag,e_K2mag,J2Flux,e_J2Flux,
+ H2Flux,e_H2Flux,K2Flux,e_K2Flux,J2snr,H2snr,K2snr,q_J2mag,
+ q_H2mag,q_K2mag,Counter,v24mag,e_24mag
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,Source_1,v2MASS_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,dIRAC,
+ dWISE,dIRACWISE,dSMOG,Seq,Class
if(ar__1(91:93) .EQ. '') Seq = iNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,2)
+ Field,Source_1,v2MASS_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,dIRAC,
+ dWISE,dIRACWISE,dSMOG,Seq,Class
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
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,IRAC,RAIRACdeg,DEIRACdeg,Simbad,Type,RASIMh,RASIMm,
+ RASIMs,DESIM_,DESIMd,DESIMm,DESIMs,IDs
RAdeg_1 = rNULL__
DEdeg_1 = 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 .GT. -180) RAdeg_1=RASIMh*15.
if(RASIMm .GT. -180) RAdeg_1=RAdeg_1+RASIMm/4.
if(RASIMs .GT. -180) RAdeg_1=RAdeg_1+RASIMs/240.
if(DESIMd .GE. 0) DEdeg_1=DESIMd
if(DESIMm .GE. 0) DEdeg_1=DEdeg_1+DESIMm/60.
if(DESIMs .GE. 0) DEdeg_1=DEdeg_1+DESIMs/3600.
if(DESIM_.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,3)
+ Field_1,IRAC,RAIRACdeg,DEIRACdeg,Simbad,Type,RASIMh,RASIMm,
+ RASIMs,DESIM_,DESIMd,DESIMm,DESIMs,IDs
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
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,GLIMPSE,YSO,RAdeg_2,DEdeg_2,Rad
c ..............Just test output...........
write(6,4)Seq_1,GLIMPSE,YSO,RAdeg_2,DEdeg_2,Rad
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,GLIMPSE_1,WISE
c ..............Just test output...........
write(6,5)Seq_2,GLIMPSE_1,WISE
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,GLIMPSE_2,Skiff
c ..............Just test output...........
write(6,6)Seq_3,GLIMPSE_2,Skiff
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,GLIMPSE_3,Simbad_1
c ..............Just test output...........
write(6,7)Seq_4,GLIMPSE_3,Simbad_1
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,Seq_5,Dist,RadGal,Ndata,Chi2_1,e_Chi2,E_Chi2_1,Av,
+ e_Av,E_Av_1,Av_avg,s_Av,Mass,e_Mass,E_Mass_1,Mass_avg,s_Mass,
+ Age,e_Age,E_Age_1,Age_avg,s_Age,Mdisk,e_Mdisk,E_Mdisk_1,
+ Mdisk_avg,s_Mdisk,Menv,e_Menv,E_Menv_1,Menv_avg,s_Menv,Mdot,
+ e_Mdot,E_Mdot_1,Mdot_avg,s_Mdot,Temp,e_Temp,E_Temp_1,Temp_avg,
+ s_Temp,Lstar,e_Lstar,E_Lstar_1,Lstar_avg,s_Lstar
if(ar__7(276:284) .EQ. '') Mdot = rNULL__
c ..............Just test output...........
write(6,8)
+ Field_2,Seq_5,Dist,RadGal,Ndata,Chi2_1,e_Chi2,E_Chi2_1,Av,
+ e_Av,E_Av_1,Av_avg,s_Av,Mass,e_Mass,E_Mass_1,Mass_avg,s_Mass,
+ Age,e_Age,E_Age_1,Age_avg,s_Age,Mdisk,e_Mdisk,E_Mdisk_1,
+ Mdisk_avg,s_Mdisk,Menv,e_Menv,E_Menv_1,Menv_avg,s_Menv,Mdot,
+ e_Mdot,E_Mdot_1,Mdot_avg,s_Mdot,Temp,e_Temp,E_Temp_1,Temp_avg,
+ s_Temp,Lstar,e_Lstar,E_Lstar_1,Lstar_avg,s_Lstar
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end