Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/170/172 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-Sep-15
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/170/172 2MASS & Gaia DR3 phot. of YSOs in three clusters (Patel+, 2025)
*================================================================================
*Identification of outer galaxy cluster members using Gaia DR3 and
*multidimensional simulation.
* Patel V., Hora J.L., Ashby M.L.N., Vig S.
* <Astron. J., 170, 172 (2025)>
* =2025AJ....170..172P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Cluster table
integer*4 nr__
parameter (nr__=1860) ! Number of records
character*518 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__) ! [123/257] Cluster number (1)
character*8 ID (nr__) ! Unique source identifier; SMOG Field ID
character*25 IRname (nr__) ! Infrared source name
character*16 v2MASS (nr__) ! 2MASS (II/246) ID
real*4 Cluster (nr__) ! [-1/509]? Cluster designation from
* Winston+2020 (J/AJ/160/68) (4)
real*4 Class (nr__) ! [1/3]? YSO classification from
* Winston+2020 (J/AJ/160/68) (2)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 e_RAdeg (nr__) ! (mas) ? Right ascension uncertainty (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
real*8 e_DEdeg (nr__) ! (mas) ? Declination uncertainty (J2000)
real*4 Jmag (nr__) ! (mag) [6.7/18.3]? Apparent 2MASS J-band Vega
* magnitude
real*4 e_Jmag (nr__) ! (mag) [0.01/0.31]? Uncertainty in Jmag
real*8 FJ (nr__) ! (mJy) [0.16/1186]? J-band flux
real*8 e_FJ (nr__) ! (mJy) [0.03/29.5]? Uncertainty in FJ
real*4 Hmag (nr__) ! (mag) [5.8/17.1]? Apparent 2MASS H-band Vega
* magnitude
real*4 e_Hmag (nr__) ! (mag) [0.01/0.21]? Uncertainty in Hmag
real*8 FH (nr__) ! (mJy) [0.3/973.4]? H-band flux
real*8 e_FH (nr__) ! (mJy) [0.05/47.6]? Uncertainty in FH
real*4 Ksmag (nr__) ! (mag) [5.4/15.9]? Apparent 2MASS Ks-band Vega
* magnitude
real*4 e_Ksmag (nr__) ! (mag) [0.01/0.27]? Uncertainty in Ksmag
real*8 FKs (nr__) ! (mJy) [0.3/860.5]? Ks-band flux
real*8 e_FKs (nr__) ! (mJy) [0.06/21.4]? Uncertainty in FKs
real*4 v3_6mag (nr__) ! (mag) [6.6/15.5]? Apparent 3.6 micron IRAC band
* 1 Vega magnitude
real*4 e_3_6mag (nr__) ! (mag) [0.01/0.21]? Uncertainty in 3.6mag
real*8 F3_6 (nr__) ! (mJy) [0.1/614.4]? IRAC band 1 flux
real*8 e_F3_6 (nr__) ! (mJy) [0.01/24.5]? Uncertainty in F3.6
real*4 v4_5mag (nr__) ! (mag) [5.6/14.8]? Apparent 4.5 micron IRAC band
* 2 Vega magnitude
real*4 e_4_5mag (nr__) ! (mag) [0.01/0.2]? Uncertainty in 4.5mag
real*8 F4_5 (nr__) ! (mJy) [0.2/1020]? IRAC band 2 flux
real*8 e_F4_5 (nr__) ! (mJy) [0.01/27.3]? Uncertainty in F4.5
real*4 W3mag (nr__) ! (mag) [1.6/10]? Apparent WISE W3 Vega magnitude
* from AllWISE (II/328)
real*4 e_W3mag (nr__) ! (mag) [0.01/0.51]? Uncertainty in W3mag
real*4 W4mag (nr__) ! (mag) [0.07/9.2]? Apparent WISE W4 Vega
* magnitude from All WISE (II/328)
real*4 e_W4mag (nr__) ! (mag) [0.009/0.51]? Uncertainty in W4mag
character*28 Designation(nr__) ! Gaia DR3 (I/355) designation
integer*8 GaiaDR3 (nr__) ! ?=-1 Gaia DR3 (I/355) source identifier
real*4 plx (nr__) ! (mas) [-3.17/4.53]? Gaia DR3 parallax
real*4 e_plx (nr__) ! (mas) [0.01/2.87]? Uncertainty in plx
real*4 RPlx (nr__) ! (mas) [-3.11/108.3]? Parallax over error
* (plxOverErr)
real*4 pmRA (nr__) ! (mas/yr) [-12.46/10.58]? Proper motion in
* direction of right ascension
real*4 e_pmRA (nr__) ! (mas/yr) [0.01/2.68]? Uncertainty in pmRA
real*4 pmDE (nr__) ! (mas/yr) [-16.26/9.06]? Proper motion in direction
* of declination
real*4 e_pmDE (nr__) ! (mas/yr) [0.01/2.98]? Uncertainty in pmDE
real*4 RUWE (nr__) ! [0.88/1.4]? Renormalized unit weight
* error from Gaia DR3
real*4 Gmag (nr__) ! (mag) [9.1/20.92]? Apparent Gaia G-band mean
* Vega magnitude
real*4 e_Gmag (nr__) ! (mag) [0.003/0.035]? Uncertainty in Gmag
real*4 BPmag (nr__) ! (mag) [9.5/22.4]? Apparent G_BP_ magnitude
real*4 e_BPmag (nr__) ! (mag) [0.003/0.5]? Uncertainty in BPmag
real*4 RPmag (nr__) ! (mag) [8.5/19.8]? Apparent G_RP_ magnitude
real*4 e_RPmag (nr__) ! (mag) [0.004/0.14]? Uncertainty in RPmag
real*4 BP_RP (nr__) ! (mag) [0.8/5.9]? Gaia BP-RP color
real*4 BP_G (nr__) ! (mag) [0.1/4.3]? Gaia BP-G color
real*4 G_RP (nr__) ! (mag) [0.5/2.5]? Gaia G-RP color
real*8 Dist (nr__) ! (pc) [345.7/11431]? ICRS source distance
real*8 e_Dist (nr__) ! (pc) [269.9/6962.2]? Lower bound on Dist
real*8 E_Dist_1 (nr__) ! (pc) [579/13968]? Upper bound on Dist
real*8 Xgal (nr__) ! (pc) [-14967/-7844]? Galactocentric
* X-coordinate
real*8 Ygal (nr__) ! (pc) [234.6/11391]? Galactocentric
* Y-coordinate
real*8 Zgal (nr__) ! (pc) [0.7/443.6]? Galactocentric Z-coordinate
integer*4 Memb (nr__) ! [0/3] Membership flag for YSO cluster
* (memFlag) (3)
character*6 Prob (nr__) ! Membership probability
character*6 e_Prob (nr__) ! Uncertainty in Prob
*Note (1): Cluster occurrences as follows:
* 123 = 563 occurrences
* 163 = 262 occurrences
* 257 = 1035 occurrences
*Note (2): 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.0 = Class 0/I: protostellar objects with a rising slope, {alpha}>0;
* (178 occurrences)
* 2.0 = Class II: characterized by decreasing slopes in the range
* -1.6<{alpha}<0; (1119 occurrences)
* 3.0 = 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; (142 occurrences)
*Note (3): Membership flag as follows:
* 0 = field YSO; 466 occurrences
* 1 = cluster member; 947 occurrences
* 2 = SFOG-only YSO identified as cluster member in Winston+2020 (J/AJ/160/68);
* 356 occurrences
* 3 = additional SFOG-only cluster member identified in this work based on
* the 30th (70th for Cluster 163) percentile YSO density contour (see
* Section 6); 91 occurrences
*Note (4): A "-1" indicates that there was no cluster affiliation identified.
C=============================================================================
C Loading file 'table5.dat' ! Cluster table
C Format for file interpretation
1 format(
+ I3,1X,A8,1X,A25,1X,A16,1X,F5.1,1X,F3.1,1X,F13.9,1X,F7.4,1X,
+ F12.9,1X,F7.4,1X,F6.3,1X,F5.3,1X,F9.4,1X,F8.5,1X,F6.3,1X,F5.3,
+ 1X,F8.4,1X,F8.5,1X,F6.3,1X,F5.3,1X,F8.4,1X,F8.5,1X,F6.3,1X,
+ F5.3,1X,F8.4,1X,F8.5,1X,F6.3,1X,F5.3,1X,F9.4,1X,F8.5,1X,F6.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,A28,1X,I19,1X,F5.2,1X,F4.2,1X,F6.2,
+ 1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F5.3,1X,F6.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F9.3,
+ 1X,F8.3,1X,F9.3,1X,F10.3,1X,F9.3,1X,F7.3,1X,I1,1X,A6,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,1860
read(1,'(A518)')ar__
read(ar__,1)
+ Seq(i__),ID(i__),IRname(i__),v2MASS(i__),Cluster(i__),
+ Class(i__),RAdeg(i__),e_RAdeg(i__),DEdeg(i__),e_DEdeg(i__),
+ Jmag(i__),e_Jmag(i__),FJ(i__),e_FJ(i__),Hmag(i__),e_Hmag(i__),
+ FH(i__),e_FH(i__),Ksmag(i__),e_Ksmag(i__),FKs(i__),e_FKs(i__),
+ v3_6mag(i__),e_3_6mag(i__),F3_6(i__),e_F3_6(i__),v4_5mag(i__),
+ e_4_5mag(i__),F4_5(i__),e_F4_5(i__),W3mag(i__),e_W3mag(i__),
+ W4mag(i__),e_W4mag(i__),Designation(i__),GaiaDR3(i__),
+ plx(i__),e_plx(i__),RPlx(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),
+ e_pmDE(i__),RUWE(i__),Gmag(i__),e_Gmag(i__),BPmag(i__),
+ e_BPmag(i__),RPmag(i__),e_RPmag(i__),BP_RP(i__),BP_G(i__),
+ G_RP(i__),Dist(i__),e_Dist(i__),E_Dist_1(i__),Xgal(i__),
+ Ygal(i__),Zgal(i__),Memb(i__),Prob(i__),e_Prob(i__)
if(ar__(57:61) .EQ. '') Cluster(i__) = rNULL__
if(ar__(63:65) .EQ. '') Class(i__) = rNULL__
if(ar__(81:87) .EQ. '') e_RAdeg(i__) = rNULL__
if(ar__(102:108) .EQ. '') e_DEdeg(i__) = rNULL__
if(ar__(110:115) .EQ. '') Jmag(i__) = rNULL__
if(ar__(117:121) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(123:131) .EQ. '') FJ(i__) = rNULL__
if(ar__(133:140) .EQ. '') e_FJ(i__) = rNULL__
if(ar__(142:147) .EQ. '') Hmag(i__) = rNULL__
if(ar__(149:153) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__(155:162) .EQ. '') FH(i__) = rNULL__
if(ar__(164:171) .EQ. '') e_FH(i__) = rNULL__
if(ar__(173:178) .EQ. '') Ksmag(i__) = rNULL__
if(ar__(180:184) .EQ. '') e_Ksmag(i__) = rNULL__
if(ar__(186:193) .EQ. '') FKs(i__) = rNULL__
if(ar__(195:202) .EQ. '') e_FKs(i__) = rNULL__
if(ar__(204:209) .EQ. '') v3_6mag(i__) = rNULL__
if(ar__(211:215) .EQ. '') e_3_6mag(i__) = rNULL__
if(ar__(217:224) .EQ. '') F3_6(i__) = rNULL__
if(ar__(226:233) .EQ. '') e_F3_6(i__) = rNULL__
if(ar__(235:240) .EQ. '') v4_5mag(i__) = rNULL__
if(ar__(242:246) .EQ. '') e_4_5mag(i__) = rNULL__
if(ar__(248:256) .EQ. '') F4_5(i__) = rNULL__
if(ar__(258:265) .EQ. '') e_F4_5(i__) = rNULL__
if(ar__(267:272) .EQ. '') W3mag(i__) = rNULL__
if(ar__(274:278) .EQ. '') e_W3mag(i__) = rNULL__
if(ar__(280:284) .EQ. '') W4mag(i__) = rNULL__
if(ar__(286:290) .EQ. '') e_W4mag(i__) = rNULL__
if(ar__(341:345) .EQ. '') plx(i__) = rNULL__
if(ar__(347:350) .EQ. '') e_plx(i__) = rNULL__
if(ar__(352:357) .EQ. '') RPlx(i__) = rNULL__
if(ar__(359:364) .EQ. '') pmRA(i__) = rNULL__
if(ar__(366:369) .EQ. '') e_pmRA(i__) = rNULL__
if(ar__(371:376) .EQ. '') pmDE(i__) = rNULL__
if(ar__(378:381) .EQ. '') e_pmDE(i__) = rNULL__
if(ar__(383:387) .EQ. '') RUWE(i__) = rNULL__
if(ar__(389:394) .EQ. '') Gmag(i__) = rNULL__
if(ar__(396:400) .EQ. '') e_Gmag(i__) = rNULL__
if(ar__(402:407) .EQ. '') BPmag(i__) = rNULL__
if(ar__(409:413) .EQ. '') e_BPmag(i__) = rNULL__
if(ar__(415:420) .EQ. '') RPmag(i__) = rNULL__
if(ar__(422:426) .EQ. '') e_RPmag(i__) = rNULL__
if(ar__(428:432) .EQ. '') BP_RP(i__) = rNULL__
if(ar__(434:438) .EQ. '') BP_G(i__) = rNULL__
if(ar__(440:444) .EQ. '') G_RP(i__) = rNULL__
if(ar__(446:454) .EQ. '') Dist(i__) = rNULL__
if(ar__(456:463) .EQ. '') e_Dist(i__) = rNULL__
if(ar__(465:473) .EQ. '') E_Dist_1(i__) = rNULL__
if(ar__(475:484) .EQ. '') Xgal(i__) = rNULL__
if(ar__(486:494) .EQ. '') Ygal(i__) = rNULL__
if(ar__(496:502) .EQ. '') Zgal(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),ID(i__),IRname(i__),v2MASS(i__),Cluster(i__),
+ Class(i__),RAdeg(i__),e_RAdeg(i__),DEdeg(i__),e_DEdeg(i__),
+ Jmag(i__),e_Jmag(i__),FJ(i__),e_FJ(i__),Hmag(i__),e_Hmag(i__),
+ FH(i__),e_FH(i__),Ksmag(i__),e_Ksmag(i__),FKs(i__),e_FKs(i__),
+ v3_6mag(i__),e_3_6mag(i__),F3_6(i__),e_F3_6(i__),v4_5mag(i__),
+ e_4_5mag(i__),F4_5(i__),e_F4_5(i__),W3mag(i__),e_W3mag(i__),
+ W4mag(i__),e_W4mag(i__),Designation(i__),GaiaDR3(i__),
+ plx(i__),e_plx(i__),RPlx(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),
+ e_pmDE(i__),RUWE(i__),Gmag(i__),e_Gmag(i__),BPmag(i__),
+ e_BPmag(i__),RPmag(i__),e_RPmag(i__),BP_RP(i__),BP_G(i__),
+ G_RP(i__),Dist(i__),e_Dist(i__),E_Dist_1(i__),Xgal(i__),
+ Ygal(i__),Zgal(i__),Memb(i__),Prob(i__),e_Prob(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end