Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/640/A128 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-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 (nr__) ! [1/155618] Internal identifier
integer*4 GDW15 (nr__) ! [119/5699]? Internal identifier in
* Guarcello et al., 2015 (2015arXiv150103761G)
integer*4 GDW13 (nr__) ! [229/326542]? Internal identifier in
* Guarcello et al., 2013 (J/ApJ/773/135)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000.0)
real*8 DEdeg (nr__) ! (deg) Declination (J2000.0)
real*4 Stet (nr__) ! [-0.11/15.64]? Stetson Variability Index
integer*4 Sx (nr__) ! [0/2]? Stetson Index selection flag (1)
integer*4 Type (nr__) ! [1/3]? Light-curve morphology type (2)
character*1 Asymm (nr__) ! [+/-] Light Curve Symetry (3)
integer*4 Disk (nr__) ! [0/1]? Disk flag (4)
real*4 Per (nr__) ! (d) ? Detected Period (5)
integer*4 Compound (nr__) ! [1/5]? Compound flag (6)
integer*4 Slope (nr__) ! [1]? Slope Flag (7)
integer*4 nJ (nr__) ! [0/107] Number of J observations
integer*4 nH (nr__) ! [0/104] Number of H observations
integer*4 nK (nr__) ! [0/104] Number of K observations
real*4 Jmag (nr__) ! (mag) [10.45/18.83]? Median UKIRT/WFCAM J magnitude
real*4 Hmag (nr__) ! (mag) [9.65/18.17]? Median UKIRT/WFCAM H magnitude
real*4 Kmag (nr__) ! (mag) [9.15/17.47]? Median UKIRT/WFCAM K magnitude
real*4 e_Jmag (nr__) ! (mag) [0.02/0.1]? Representative uncertainty in J
real*4 e_Hmag (nr__) ! (mag) [0.02/0.38]? Representative uncertainty in H
real*4 e_Kmag (nr__) ! (mag) [0.02/0.09]? Representative uncertainty in K
real*4 J_H (nr__) ! (mag) [0.23/3.54]? Median J-H color
real*4 J_K (nr__) ! (mag) [0.33/5.94]? Median J-K color
real*4 H_K (nr__) ! (mag) [0.06/2.4]? Median H-K color
real*4 e_J_H (nr__) ! (mag) [0.03/0.12]? Representative J-H uncertainty
real*4 e_J_K (nr__) ! (mag) [0.03/0.13]? Representative J-K uncertainty
real*4 e_H_K (nr__) ! (mag) [0.03/0.13]? Representative H-K uncertainty
real*4 Jamp (nr__) ! (mag) [0.02/1.79]? Peak-to-peak J amplitude
real*4 Hamp (nr__) ! (mag) [0.02/1.84]? Peak-to-peak H amplitude
real*4 Kamp (nr__) ! (mag) [0.02/1.61]? Peak-to-peak K amplitude
real*4 J_Hamp (nr__) ! (mag) [0.02/0.72]? Peak-to-peak J-H amplitude
real*4 J_Kamp (nr__) ! (mag) [0.01/1.08]? Peak-to-peak J-K amplitude
real*4 H_Kamp (nr__) ! (mag) [0.01/0.97]? Peak-to-peak H-K amplitude
real*4 rmsJamp (nr__) ! (mag) [0.02/1.29]? rms J amplitude
real*4 rmsHamp (nr__) ! (mag) [0.02/1.32]? rms H amplitude
real*4 rmsKamp (nr__) ! (mag) [0.02/1.14]? rms K amplitude
real*4 rmsJ_Hamp (nr__) ! (mag) [0.01/0.27]? rms J-H amplitude
real*4 rmsJ_Kamp (nr__) ! (mag) [0.01/0.39]? rms J-K amplitude
real*4 rmsH_Kamp (nr__) ! (mag) [0.01/0.35]? rms H-K amplitude
real*4 IncKHK (nr__) ! (deg) [-79.5/74.2]? The K,(H-K) slope
real*4 IncJJH (nr__) ! (deg) [-58.3/80.5]? The J,(J-H) slope
real*4 IncJHHK (nr__) ! (deg) [21/57.5]? The (J-H),(H-K) slope
real*4 e_IncJHHK (nr__) ! (deg) [0./0.8]? Uncertainty in IncJHHK
real*4 e_IncKHK (nr__) ! (deg) [0.5/9.4]? Uncertainty in IncKHK
real*4 e_IncJJH (nr__) ! (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 (nr__1) ! (d) ? J-band Julian date (HJD-2400000)
real*4 Jmag_1 (nr__1) ! (mag) ? UKIRT/WFCAM J magnitude
real*4 e_Jmag_1 (nr__1) ! (mag) ? rms uncertainty in Jmag
real*8 RAJdeg (nr__1) ! (deg) ? J-band right ascension (J2000)
real*8 DEJdeg (nr__1) ! (deg) ? J-band declination (J2000)
real*8 HJD_H (nr__1) ! (d) ? H-band Julian date (HJD-2400000)
real*4 Hmag_1 (nr__1) ! (mag) ? UKIRT/WFCAM H magnitude
real*4 e_Hmag_1 (nr__1) ! (mag) ? rms uncertainty in Hmag
real*8 RAHdeg (nr__1) ! (deg) ? H-band right ascension (J2000)
real*8 DEHdeg (nr__1) ! (deg) ? H-band declination (J2000)
real*8 HJD_K (nr__1) ! (d) ? K-band Julian date (HJD-2400000)
real*4 Kmag_1 (nr__1) ! (mag) ? UKIRT/WFCAM K magnitude
real*4 e_Kmag_1 (nr__1) ! (mag) ? rms uncertainty in Kmag
real*8 RAKdeg (nr__1) ! (deg) ? K-band right ascension (J2000)
real*8 DEKdeg (nr__1) ! (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(i__),GDW15(i__),GDW13(i__),RAdeg(i__),DEdeg(i__),Stet(i__),
+ Sx(i__),Type(i__),Asymm(i__),Disk(i__),Per(i__),Compound(i__),
+ Slope(i__),nJ(i__),nH(i__),nK(i__),Jmag(i__),Hmag(i__),
+ Kmag(i__),e_Jmag(i__),e_Hmag(i__),e_Kmag(i__),J_H(i__),
+ J_K(i__),H_K(i__),e_J_H(i__),e_J_K(i__),e_H_K(i__),Jamp(i__),
+ Hamp(i__),Kamp(i__),J_Hamp(i__),J_Kamp(i__),H_Kamp(i__),
+ rmsJamp(i__),rmsHamp(i__),rmsKamp(i__),rmsJ_Hamp(i__),
+ rmsJ_Kamp(i__),rmsH_Kamp(i__),IncKHK(i__),IncJJH(i__),
+ IncJHHK(i__),e_IncJHHK(i__),e_IncKHK(i__),e_IncJJH(i__)
if(ar__(8:11) .EQ. '') GDW15(i__) = iNULL__
if(ar__(13:18) .EQ. '') GDW13(i__) = iNULL__
if(ar__(39:43) .EQ. '') Stet(i__) = rNULL__
if(ar__(45:45) .EQ. '') Sx(i__) = iNULL__
if(ar__(47:47) .EQ. '') Type(i__) = iNULL__
if(ar__(51:51) .EQ. '') Disk(i__) = iNULL__
if(ar__(53:58) .EQ. '') Per(i__) = rNULL__
if(ar__(60:60) .EQ. '') Compound(i__) = iNULL__
if(ar__(62:62) .EQ. '') Slope(i__) = iNULL__
if(ar__(76:80) .EQ. '') Jmag(i__) = rNULL__
if(ar__(82:86) .EQ. '') Hmag(i__) = rNULL__
if(ar__(88:92) .EQ. '') Kmag(i__) = rNULL__
if(ar__(94:97) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(99:102) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__(104:107) .EQ. '') e_Kmag(i__) = rNULL__
if(ar__(109:112) .EQ. '') J_H(i__) = rNULL__
if(ar__(114:117) .EQ. '') J_K(i__) = rNULL__
if(ar__(119:122) .EQ. '') H_K(i__) = rNULL__
if(ar__(124:127) .EQ. '') e_J_H(i__) = rNULL__
if(ar__(129:132) .EQ. '') e_J_K(i__) = rNULL__
if(ar__(134:137) .EQ. '') e_H_K(i__) = rNULL__
if(ar__(139:142) .EQ. '') Jamp(i__) = rNULL__
if(ar__(144:147) .EQ. '') Hamp(i__) = rNULL__
if(ar__(149:152) .EQ. '') Kamp(i__) = rNULL__
if(ar__(154:157) .EQ. '') J_Hamp(i__) = rNULL__
if(ar__(159:162) .EQ. '') J_Kamp(i__) = rNULL__
if(ar__(164:167) .EQ. '') H_Kamp(i__) = rNULL__
if(ar__(169:172) .EQ. '') rmsJamp(i__) = rNULL__
if(ar__(174:177) .EQ. '') rmsHamp(i__) = rNULL__
if(ar__(179:182) .EQ. '') rmsKamp(i__) = rNULL__
if(ar__(184:187) .EQ. '') rmsJ_Hamp(i__) = rNULL__
if(ar__(189:192) .EQ. '') rmsJ_Kamp(i__) = rNULL__
if(ar__(194:197) .EQ. '') rmsH_Kamp(i__) = rNULL__
if(ar__(199:203) .EQ. '') IncKHK(i__) = rNULL__
if(ar__(205:209) .EQ. '') IncJJH(i__) = rNULL__
if(ar__(211:214) .EQ. '') IncJHHK(i__) = rNULL__
if(ar__(216:218) .EQ. '') e_IncJHHK(i__) = rNULL__
if(ar__(220:222) .EQ. '') e_IncKHK(i__) = rNULL__
if(ar__(224:226) .EQ. '') e_IncJJH(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),GDW15(i__),GDW13(i__),RAdeg(i__),DEdeg(i__),Stet(i__),
+ Sx(i__),Type(i__),Asymm(i__),Disk(i__),Per(i__),Compound(i__),
+ Slope(i__),nJ(i__),nH(i__),nK(i__),Jmag(i__),Hmag(i__),
+ Kmag(i__),e_Jmag(i__),e_Hmag(i__),e_Kmag(i__),J_H(i__),
+ J_K(i__),H_K(i__),e_J_H(i__),e_J_K(i__),e_H_K(i__),Jamp(i__),
+ Hamp(i__),Kamp(i__),J_Hamp(i__),J_Kamp(i__),H_Kamp(i__),
+ rmsJamp(i__),rmsHamp(i__),rmsKamp(i__),rmsJ_Hamp(i__),
+ rmsJ_Kamp(i__),rmsH_Kamp(i__),IncKHK(i__),IncJJH(i__),
+ IncJHHK(i__),e_IncJHHK(i__),e_IncKHK(i__),e_IncJJH(i__)
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(i__),Jmag_1(i__),e_Jmag_1(i__),RAJdeg(i__),DEJdeg(i__),
+ HJD_H(i__),Hmag_1(i__),e_Hmag_1(i__),RAHdeg(i__),DEHdeg(i__),
+ HJD_K(i__),Kmag_1(i__),e_Kmag_1(i__),RAKdeg(i__),DEKdeg(i__)
if(ar__1(1:11) .EQ. '') HJD_J(i__) = rNULL__
if(ar__1(14:18) .EQ. '') Jmag_1(i__) = rNULL__
if(ar__1(21:24) .EQ. '') e_Jmag_1(i__) = rNULL__
if(ar__1(27:35) .EQ. '') RAJdeg(i__) = rNULL__
if(ar__1(38:45) .EQ. '') DEJdeg(i__) = rNULL__
if(ar__1(48:58) .EQ. '') HJD_H(i__) = rNULL__
if(ar__1(61:65) .EQ. '') Hmag_1(i__) = rNULL__
if(ar__1(68:71) .EQ. '') e_Hmag_1(i__) = rNULL__
if(ar__1(74:82) .EQ. '') RAHdeg(i__) = rNULL__
if(ar__1(85:92) .EQ. '') DEHdeg(i__) = rNULL__
if(ar__1(95:105) .EQ. '') HJD_K(i__) = rNULL__
if(ar__1(108:112) .EQ. '') Kmag_1(i__) = rNULL__
if(ar__1(115:118) .EQ. '') e_Kmag_1(i__) = rNULL__
if(ar__1(121:129) .EQ. '') RAKdeg(i__) = rNULL__
if(ar__1(132:139) .EQ. '') DEKdeg(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ HJD_J(i__),Jmag_1(i__),e_Jmag_1(i__),RAJdeg(i__),DEJdeg(i__),
+ HJD_H(i__),Hmag_1(i__),e_Hmag_1(i__),RAHdeg(i__),DEHdeg(i__),
+ HJD_K(i__),Kmag_1(i__),e_Kmag_1(i__),RAKdeg(i__),DEKdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end