Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/170/172 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-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 ! [123/257] Cluster number (1)
character*8 ID ! Unique source identifier; SMOG Field ID
character*25 IRname ! Infrared source name
character*16 v2MASS ! 2MASS (II/246) ID
real*4 Cluster ! [-1/509]? Cluster designation from
* Winston+2020 (J/AJ/160/68) (4)
real*4 Class ! [1/3]? YSO classification from
* Winston+2020 (J/AJ/160/68) (2)
real*8 RAdeg ! (deg) Right ascension (J2000)
real*8 e_RAdeg ! (mas) ? Right ascension uncertainty (J2000)
real*8 DEdeg ! (deg) Declination (J2000)
real*8 e_DEdeg ! (mas) ? Declination uncertainty (J2000)
real*4 Jmag ! (mag) [6.7/18.3]? Apparent 2MASS J-band Vega
* magnitude
real*4 e_Jmag ! (mag) [0.01/0.31]? Uncertainty in Jmag
real*8 FJ ! (mJy) [0.16/1186]? J-band flux
real*8 e_FJ ! (mJy) [0.03/29.5]? Uncertainty in FJ
real*4 Hmag ! (mag) [5.8/17.1]? Apparent 2MASS H-band Vega
* magnitude
real*4 e_Hmag ! (mag) [0.01/0.21]? Uncertainty in Hmag
real*8 FH ! (mJy) [0.3/973.4]? H-band flux
real*8 e_FH ! (mJy) [0.05/47.6]? Uncertainty in FH
real*4 Ksmag ! (mag) [5.4/15.9]? Apparent 2MASS Ks-band Vega
* magnitude
real*4 e_Ksmag ! (mag) [0.01/0.27]? Uncertainty in Ksmag
real*8 FKs ! (mJy) [0.3/860.5]? Ks-band flux
real*8 e_FKs ! (mJy) [0.06/21.4]? Uncertainty in FKs
real*4 v3_6mag ! (mag) [6.6/15.5]? Apparent 3.6 micron IRAC band
* 1 Vega magnitude
real*4 e_3_6mag ! (mag) [0.01/0.21]? Uncertainty in 3.6mag
real*8 F3_6 ! (mJy) [0.1/614.4]? IRAC band 1 flux
real*8 e_F3_6 ! (mJy) [0.01/24.5]? Uncertainty in F3.6
real*4 v4_5mag ! (mag) [5.6/14.8]? Apparent 4.5 micron IRAC band
* 2 Vega magnitude
real*4 e_4_5mag ! (mag) [0.01/0.2]? Uncertainty in 4.5mag
real*8 F4_5 ! (mJy) [0.2/1020]? IRAC band 2 flux
real*8 e_F4_5 ! (mJy) [0.01/27.3]? Uncertainty in F4.5
real*4 W3mag ! (mag) [1.6/10]? Apparent WISE W3 Vega magnitude
* from AllWISE (II/328)
real*4 e_W3mag ! (mag) [0.01/0.51]? Uncertainty in W3mag
real*4 W4mag ! (mag) [0.07/9.2]? Apparent WISE W4 Vega
* magnitude from All WISE (II/328)
real*4 e_W4mag ! (mag) [0.009/0.51]? Uncertainty in W4mag
character*28 Designation ! Gaia DR3 (I/355) designation
integer*8 GaiaDR3 ! ?=-1 Gaia DR3 (I/355) source identifier
real*4 plx ! (mas) [-3.17/4.53]? Gaia DR3 parallax
real*4 e_plx ! (mas) [0.01/2.87]? Uncertainty in plx
real*4 RPlx ! (mas) [-3.11/108.3]? Parallax over error
* (plxOverErr)
real*4 pmRA ! (mas/yr) [-12.46/10.58]? Proper motion in
* direction of right ascension
real*4 e_pmRA ! (mas/yr) [0.01/2.68]? Uncertainty in pmRA
real*4 pmDE ! (mas/yr) [-16.26/9.06]? Proper motion in direction
* of declination
real*4 e_pmDE ! (mas/yr) [0.01/2.98]? Uncertainty in pmDE
real*4 RUWE ! [0.88/1.4]? Renormalized unit weight
* error from Gaia DR3
real*4 Gmag ! (mag) [9.1/20.92]? Apparent Gaia G-band mean
* Vega magnitude
real*4 e_Gmag ! (mag) [0.003/0.035]? Uncertainty in Gmag
real*4 BPmag ! (mag) [9.5/22.4]? Apparent G_BP_ magnitude
real*4 e_BPmag ! (mag) [0.003/0.5]? Uncertainty in BPmag
real*4 RPmag ! (mag) [8.5/19.8]? Apparent G_RP_ magnitude
real*4 e_RPmag ! (mag) [0.004/0.14]? Uncertainty in RPmag
real*4 BP_RP ! (mag) [0.8/5.9]? Gaia BP-RP color
real*4 BP_G ! (mag) [0.1/4.3]? Gaia BP-G color
real*4 G_RP ! (mag) [0.5/2.5]? Gaia G-RP color
real*8 Dist ! (pc) [345.7/11431]? ICRS source distance
real*8 e_Dist ! (pc) [269.9/6962.2]? Lower bound on Dist
real*8 E_Dist_1 ! (pc) [579/13968]? Upper bound on Dist
real*8 Xgal ! (pc) [-14967/-7844]? Galactocentric
* X-coordinate
real*8 Ygal ! (pc) [234.6/11391]? Galactocentric
* Y-coordinate
real*8 Zgal ! (pc) [0.7/443.6]? Galactocentric Z-coordinate
integer*4 Memb ! [0/3] Membership flag for YSO cluster
* (memFlag) (3)
character*6 Prob ! Membership probability
character*6 e_Prob ! 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,ID,IRname,v2MASS,Cluster,Class,RAdeg,e_RAdeg,DEdeg,
+ e_DEdeg,Jmag,e_Jmag,FJ,e_FJ,Hmag,e_Hmag,FH,e_FH,Ksmag,e_Ksmag,
+ FKs,e_FKs,v3_6mag,e_3_6mag,F3_6,e_F3_6,v4_5mag,e_4_5mag,F4_5,
+ e_F4_5,W3mag,e_W3mag,W4mag,e_W4mag,Designation,GaiaDR3,plx,
+ e_plx,RPlx,pmRA,e_pmRA,pmDE,e_pmDE,RUWE,Gmag,e_Gmag,BPmag,
+ e_BPmag,RPmag,e_RPmag,BP_RP,BP_G,G_RP,Dist,e_Dist,E_Dist_1,
+ Xgal,Ygal,Zgal,Memb,Prob,e_Prob
if(ar__(57:61) .EQ. '') Cluster = rNULL__
if(ar__(63:65) .EQ. '') Class = rNULL__
if(ar__(81:87) .EQ. '') e_RAdeg = rNULL__
if(ar__(102:108) .EQ. '') e_DEdeg = rNULL__
if(ar__(110:115) .EQ. '') Jmag = rNULL__
if(ar__(117:121) .EQ. '') e_Jmag = rNULL__
if(ar__(123:131) .EQ. '') FJ = rNULL__
if(ar__(133:140) .EQ. '') e_FJ = rNULL__
if(ar__(142:147) .EQ. '') Hmag = rNULL__
if(ar__(149:153) .EQ. '') e_Hmag = rNULL__
if(ar__(155:162) .EQ. '') FH = rNULL__
if(ar__(164:171) .EQ. '') e_FH = rNULL__
if(ar__(173:178) .EQ. '') Ksmag = rNULL__
if(ar__(180:184) .EQ. '') e_Ksmag = rNULL__
if(ar__(186:193) .EQ. '') FKs = rNULL__
if(ar__(195:202) .EQ. '') e_FKs = rNULL__
if(ar__(204:209) .EQ. '') v3_6mag = rNULL__
if(ar__(211:215) .EQ. '') e_3_6mag = rNULL__
if(ar__(217:224) .EQ. '') F3_6 = rNULL__
if(ar__(226:233) .EQ. '') e_F3_6 = rNULL__
if(ar__(235:240) .EQ. '') v4_5mag = rNULL__
if(ar__(242:246) .EQ. '') e_4_5mag = rNULL__
if(ar__(248:256) .EQ. '') F4_5 = rNULL__
if(ar__(258:265) .EQ. '') e_F4_5 = rNULL__
if(ar__(267:272) .EQ. '') W3mag = rNULL__
if(ar__(274:278) .EQ. '') e_W3mag = rNULL__
if(ar__(280:284) .EQ. '') W4mag = rNULL__
if(ar__(286:290) .EQ. '') e_W4mag = rNULL__
if(ar__(341:345) .EQ. '') plx = rNULL__
if(ar__(347:350) .EQ. '') e_plx = rNULL__
if(ar__(352:357) .EQ. '') RPlx = rNULL__
if(ar__(359:364) .EQ. '') pmRA = rNULL__
if(ar__(366:369) .EQ. '') e_pmRA = rNULL__
if(ar__(371:376) .EQ. '') pmDE = rNULL__
if(ar__(378:381) .EQ. '') e_pmDE = rNULL__
if(ar__(383:387) .EQ. '') RUWE = rNULL__
if(ar__(389:394) .EQ. '') Gmag = rNULL__
if(ar__(396:400) .EQ. '') e_Gmag = rNULL__
if(ar__(402:407) .EQ. '') BPmag = rNULL__
if(ar__(409:413) .EQ. '') e_BPmag = rNULL__
if(ar__(415:420) .EQ. '') RPmag = rNULL__
if(ar__(422:426) .EQ. '') e_RPmag = rNULL__
if(ar__(428:432) .EQ. '') BP_RP = rNULL__
if(ar__(434:438) .EQ. '') BP_G = rNULL__
if(ar__(440:444) .EQ. '') G_RP = rNULL__
if(ar__(446:454) .EQ. '') Dist = rNULL__
if(ar__(456:463) .EQ. '') e_Dist = rNULL__
if(ar__(465:473) .EQ. '') E_Dist_1 = rNULL__
if(ar__(475:484) .EQ. '') Xgal = rNULL__
if(ar__(486:494) .EQ. '') Ygal = rNULL__
if(ar__(496:502) .EQ. '') Zgal = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq,ID,IRname,v2MASS,Cluster,Class,RAdeg,e_RAdeg,DEdeg,
+ e_DEdeg,Jmag,e_Jmag,FJ,e_FJ,Hmag,e_Hmag,FH,e_FH,Ksmag,e_Ksmag,
+ FKs,e_FKs,v3_6mag,e_3_6mag,F3_6,e_F3_6,v4_5mag,e_4_5mag,F4_5,
+ e_F4_5,W3mag,e_W3mag,W4mag,e_W4mag,Designation,GaiaDR3,plx,
+ e_plx,RPlx,pmRA,e_pmRA,pmDE,e_pmDE,RUWE,Gmag,e_Gmag,BPmag,
+ e_BPmag,RPmag,e_RPmag,BP_RP,BP_G,G_RP,Dist,e_Dist,E_Dist_1,
+ Xgal,Ygal,Zgal,Memb,Prob,e_Prob
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end