Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/640/A128 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-12
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/A+A/640/A128 Cygnus OB2 association NIR light curves (Roquette+, 2020)
*================================================================================
*Near-infrared time-series photometry in the field of Cygnus OB2 association.
*II. Mapping the variability of candidate members.
* Roquette J., Alencar S.H.P., Bouvier J., Guarcello M.G., Reipurth B.
* <Astron. Astrophys. 640, A128 (2020)>
* =2020A&A...640A.128R (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1-5.dat' ! Paper's tables 1 to 5
integer*4 nr__
parameter (nr__=5083) ! Number of records
character*226 ar__ ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 ID ! [1/155618] Internal identifier
integer*4 GDW15 ! [119/5699]? Internal identifier in
* Guarcello et al., 2015 (2015arXiv150103761G)
integer*4 GDW13 ! [229/326542]? Internal identifier in
* Guarcello et al., 2013 (J/ApJ/773/135)
real*8 RAdeg ! (deg) Right ascension (J2000.0)
real*8 DEdeg ! (deg) Declination (J2000.0)
real*4 Stet ! [-0.11/15.64]? Stetson Variability Index
integer*4 Sx ! [0/2]? Stetson Index selection flag (1)
integer*4 Type ! [1/3]? Light-curve morphology type (2)
character*1 Asymm ! [+/-] Light Curve Symetry (3)
integer*4 Disk ! [0/1]? Disk flag (4)
real*4 Per ! (d) ? Detected Period (5)
integer*4 Compound ! [1/5]? Compound flag (6)
integer*4 Slope ! [1]? Slope Flag (7)
integer*4 nJ ! [0/107] Number of J observations
integer*4 nH ! [0/104] Number of H observations
integer*4 nK ! [0/104] Number of K observations
real*4 Jmag ! (mag) [10.45/18.83]? Median UKIRT/WFCAM J magnitude
real*4 Hmag ! (mag) [9.65/18.17]? Median UKIRT/WFCAM H magnitude
real*4 Kmag ! (mag) [9.15/17.47]? Median UKIRT/WFCAM K magnitude
real*4 e_Jmag ! (mag) [0.02/0.1]? Representative uncertainty in J
real*4 e_Hmag ! (mag) [0.02/0.38]? Representative uncertainty in H
real*4 e_Kmag ! (mag) [0.02/0.09]? Representative uncertainty in K
real*4 J_H ! (mag) [0.23/3.54]? Median J-H color
real*4 J_K ! (mag) [0.33/5.94]? Median J-K color
real*4 H_K ! (mag) [0.06/2.4]? Median H-K color
real*4 e_J_H ! (mag) [0.03/0.12]? Representative J-H uncertainty
real*4 e_J_K ! (mag) [0.03/0.13]? Representative J-K uncertainty
real*4 e_H_K ! (mag) [0.03/0.13]? Representative H-K uncertainty
real*4 Jamp ! (mag) [0.02/1.79]? Peak-to-peak J amplitude
real*4 Hamp ! (mag) [0.02/1.84]? Peak-to-peak H amplitude
real*4 Kamp ! (mag) [0.02/1.61]? Peak-to-peak K amplitude
real*4 J_Hamp ! (mag) [0.02/0.72]? Peak-to-peak J-H amplitude
real*4 J_Kamp ! (mag) [0.01/1.08]? Peak-to-peak J-K amplitude
real*4 H_Kamp ! (mag) [0.01/0.97]? Peak-to-peak H-K amplitude
real*4 rmsJamp ! (mag) [0.02/1.29]? rms J amplitude
real*4 rmsHamp ! (mag) [0.02/1.32]? rms H amplitude
real*4 rmsKamp ! (mag) [0.02/1.14]? rms K amplitude
real*4 rmsJ_Hamp ! (mag) [0.01/0.27]? rms J-H amplitude
real*4 rmsJ_Kamp ! (mag) [0.01/0.39]? rms J-K amplitude
real*4 rmsH_Kamp ! (mag) [0.01/0.35]? rms H-K amplitude
real*4 IncKHK ! (deg) [-79.5/74.2]? The K,(H-K) slope
real*4 IncJJH ! (deg) [-58.3/80.5]? The J,(J-H) slope
real*4 IncJHHK ! (deg) [21/57.5]? The (J-H),(H-K) slope
real*4 e_IncJHHK ! (deg) [0./0.8]? Uncertainty in IncJHHK
real*4 e_IncKHK ! (deg) [0.5/9.4]? Uncertainty in IncKHK
real*4 e_IncJJH ! (deg) [0.2/7.5]? Uncertainty in IncJJH
*Note (1): Stetson Index selection Flag as follows:
* 1 = Variable stars in the S1 sample with S>1
* 2 = Variable stars in the S2 sample with 0.25<S<1
* - = Star was not variable in this survey
*Note (2): Light-curve morphology type as follows:
* 1 = Periodic stars
* 2 = Other types of variability
* 3 = Occulation variable
*Note (3): Light Curve Symetry as follows:
* + = If light-curve had burst asymmetry
* - = If light-curve had dimming asymmetry
*Note (4): Disk flag as follows:
* 0 = diskless
* 1 = disk-bearing
*Note (5): Period measured in Roquette et al., 2017, J/A+A/603/A106.
*Note (6): Compound flag as follows:
* 1 = Compound variable with a feature of distinct color
* 2 = Compound variable with variable IR-excess
* 3 = Compound variable with changing color-slope
* 5 = both 2 and 3
*Note (7): Slope flag as follows:
* 1 = if the slope of the star's trajectories in the color-magnitude
* and/or color-color diagram measured
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lc/*' ! Individual light curves
integer*4 nr__1
parameter (nr__1=5083) ! Number of records
character*139 ar__1 ! Full-size record
C J2000 position composed of: RAJdeg DEJdeg
real*8 HJD_J ! (d) ? J-band Julian date (HJD-2400000)
real*4 Jmag_1 ! (mag) ? UKIRT/WFCAM J magnitude
real*4 e_Jmag_1 ! (mag) ? rms uncertainty in Jmag
real*8 RAJdeg ! (deg) ? J-band right ascension (J2000)
real*8 DEJdeg ! (deg) ? J-band declination (J2000)
real*8 HJD_H ! (d) ? H-band Julian date (HJD-2400000)
real*4 Hmag_1 ! (mag) ? UKIRT/WFCAM H magnitude
real*4 e_Hmag_1 ! (mag) ? rms uncertainty in Hmag
real*8 RAHdeg ! (deg) ? H-band right ascension (J2000)
real*8 DEHdeg ! (deg) ? H-band declination (J2000)
real*8 HJD_K ! (d) ? K-band Julian date (HJD-2400000)
real*4 Kmag_1 ! (mag) ? UKIRT/WFCAM K magnitude
real*4 e_Kmag_1 ! (mag) ? rms uncertainty in Kmag
real*8 RAKdeg ! (deg) ? K-band right ascension (J2000)
real*8 DEKdeg ! (deg) ? K-band declination (J2000)
C=============================================================================
C Loading file 'table1-5.dat' ! Paper's tables 1 to 5
C Format for file interpretation
1 format(
+ I6,1X,I4,1X,I6,1X,F9.5,1X,F8.5,1X,F5.2,1X,I1,1X,I1,1X,A1,1X,
+ I1,1X,F6.3,1X,I1,1X,I1,1X,I3,1X,I3,1X,I3,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.1,
+ 1X,F5.1,1X,F4.1,1X,F3.1,1X,F3.1,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1-5.dat')
write(6,*) '....Loading file: table1-5.dat'
do i__=1,5083
read(1,'(A226)')ar__
read(ar__,1)
+ ID,GDW15,GDW13,RAdeg,DEdeg,Stet,Sx,Type,Asymm,Disk,Per,
+ Compound,Slope,nJ,nH,nK,Jmag,Hmag,Kmag,e_Jmag,e_Hmag,e_Kmag,
+ J_H,J_K,H_K,e_J_H,e_J_K,e_H_K,Jamp,Hamp,Kamp,J_Hamp,J_Kamp,
+ H_Kamp,rmsJamp,rmsHamp,rmsKamp,rmsJ_Hamp,rmsJ_Kamp,rmsH_Kamp,
+ IncKHK,IncJJH,IncJHHK,e_IncJHHK,e_IncKHK,e_IncJJH
if(ar__(8:11) .EQ. '') GDW15 = iNULL__
if(ar__(13:18) .EQ. '') GDW13 = iNULL__
if(ar__(39:43) .EQ. '') Stet = rNULL__
if(ar__(45:45) .EQ. '') Sx = iNULL__
if(ar__(47:47) .EQ. '') Type = iNULL__
if(ar__(51:51) .EQ. '') Disk = iNULL__
if(ar__(53:58) .EQ. '') Per = rNULL__
if(ar__(60:60) .EQ. '') Compound = iNULL__
if(ar__(62:62) .EQ. '') Slope = iNULL__
if(ar__(76:80) .EQ. '') Jmag = rNULL__
if(ar__(82:86) .EQ. '') Hmag = rNULL__
if(ar__(88:92) .EQ. '') Kmag = rNULL__
if(ar__(94:97) .EQ. '') e_Jmag = rNULL__
if(ar__(99:102) .EQ. '') e_Hmag = rNULL__
if(ar__(104:107) .EQ. '') e_Kmag = rNULL__
if(ar__(109:112) .EQ. '') J_H = rNULL__
if(ar__(114:117) .EQ. '') J_K = rNULL__
if(ar__(119:122) .EQ. '') H_K = rNULL__
if(ar__(124:127) .EQ. '') e_J_H = rNULL__
if(ar__(129:132) .EQ. '') e_J_K = rNULL__
if(ar__(134:137) .EQ. '') e_H_K = rNULL__
if(ar__(139:142) .EQ. '') Jamp = rNULL__
if(ar__(144:147) .EQ. '') Hamp = rNULL__
if(ar__(149:152) .EQ. '') Kamp = rNULL__
if(ar__(154:157) .EQ. '') J_Hamp = rNULL__
if(ar__(159:162) .EQ. '') J_Kamp = rNULL__
if(ar__(164:167) .EQ. '') H_Kamp = rNULL__
if(ar__(169:172) .EQ. '') rmsJamp = rNULL__
if(ar__(174:177) .EQ. '') rmsHamp = rNULL__
if(ar__(179:182) .EQ. '') rmsKamp = rNULL__
if(ar__(184:187) .EQ. '') rmsJ_Hamp = rNULL__
if(ar__(189:192) .EQ. '') rmsJ_Kamp = rNULL__
if(ar__(194:197) .EQ. '') rmsH_Kamp = rNULL__
if(ar__(199:203) .EQ. '') IncKHK = rNULL__
if(ar__(205:209) .EQ. '') IncJJH = rNULL__
if(ar__(211:214) .EQ. '') IncJHHK = rNULL__
if(ar__(216:218) .EQ. '') e_IncJHHK = rNULL__
if(ar__(220:222) .EQ. '') e_IncKHK = rNULL__
if(ar__(224:226) .EQ. '') e_IncJJH = rNULL__
c ..............Just test output...........
write(6,1)
+ ID,GDW15,GDW13,RAdeg,DEdeg,Stet,Sx,Type,Asymm,Disk,Per,
+ Compound,Slope,nJ,nH,nK,Jmag,Hmag,Kmag,e_Jmag,e_Hmag,e_Kmag,
+ J_H,J_K,H_K,e_J_H,e_J_K,e_H_K,Jamp,Hamp,Kamp,J_Hamp,J_Kamp,
+ H_Kamp,rmsJamp,rmsHamp,rmsKamp,rmsJ_Hamp,rmsJ_Kamp,rmsH_Kamp,
+ IncKHK,IncJJH,IncJHHK,e_IncJHHK,e_IncKHK,e_IncJJH
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'lc/*' ! Individual light curves
C Format for file interpretation
2 format(
+ F11.5,2X,F5.2,2X,F4.2,2X,F9.5,2X,F8.5,2X,F11.5,2X,F5.2,2X,
+ F4.2,2X,F9.5,2X,F8.5,2X,F11.5,2X,F5.2,2X,F4.2,2X,F9.5,2X,F8.5)
C Effective file loading
open(unit=1,status='old',file=
+'lc/*')
write(6,*) '....Loading file: lc/*'
do i__=1,5083
read(1,'(A139)')ar__1
read(ar__1,2)
+ HJD_J,Jmag_1,e_Jmag_1,RAJdeg,DEJdeg,HJD_H,Hmag_1,e_Hmag_1,
+ RAHdeg,DEHdeg,HJD_K,Kmag_1,e_Kmag_1,RAKdeg,DEKdeg
if(ar__1(1:11) .EQ. '') HJD_J = rNULL__
if(ar__1(14:18) .EQ. '') Jmag_1 = rNULL__
if(ar__1(21:24) .EQ. '') e_Jmag_1 = rNULL__
if(ar__1(27:35) .EQ. '') RAJdeg = rNULL__
if(ar__1(38:45) .EQ. '') DEJdeg = rNULL__
if(ar__1(48:58) .EQ. '') HJD_H = rNULL__
if(ar__1(61:65) .EQ. '') Hmag_1 = rNULL__
if(ar__1(68:71) .EQ. '') e_Hmag_1 = rNULL__
if(ar__1(74:82) .EQ. '') RAHdeg = rNULL__
if(ar__1(85:92) .EQ. '') DEHdeg = rNULL__
if(ar__1(95:105) .EQ. '') HJD_K = rNULL__
if(ar__1(108:112) .EQ. '') Kmag_1 = rNULL__
if(ar__1(115:118) .EQ. '') e_Kmag_1 = rNULL__
if(ar__1(121:129) .EQ. '') RAKdeg = rNULL__
if(ar__1(132:139) .EQ. '') DEKdeg = rNULL__
c ..............Just test output...........
write(6,2)
+ HJD_J,Jmag_1,e_Jmag_1,RAJdeg,DEJdeg,HJD_H,Hmag_1,e_Hmag_1,
+ RAHdeg,DEHdeg,HJD_K,Kmag_1,e_Kmag_1,RAKdeg,DEKdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end