Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/164/80 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-Aug-14
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/164/80 Periods and supporting data for UCL/LCC non/members (Rebull+, 2022)
*================================================================================
*Rotation of Low-mass Stars in Upper Centaurus-Lupus and Lower Centaurus-Crux
*with TESS.
* Rebull L.M., Stauffer J.R., Hillenbrand L.A., Cody A.M., Kruse E.,
* Powell B.P.
* <Astron. J., 164, 80 (2022)>
* =2022AJ....164...80R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Periods and supporting data for UCL/LCC members
with viable light curves
integer*4 nr__
parameter (nr__=3714) ! Number of records
character*243 ar__ ! Full-size record
integer*8 TIC ! TESS Input Catalog identifier
character*18 Name ! Coordinate-based (RA/DE, J2000) name for target
character*23 OName ! Alternate name for target
integer*8 Gaia ! ? Gaia DR2 ID
real*8 Dist ! (pc) ?=-9.00 Distance from Bailer-Jones+, 2018, I/347
character*6 Bin ! Membership sample (gold, silver, or bronze)
real*4 Kmag ! (mag) ? Ks magnitude, Vega, if observed
real*4 V_Ks ! (mag) ? V-Ks color used, in Vega mags, observed or
* inferred; see text
real*4 EV_Ks ! (mag) ? E(V-Ks) adopted for this star, Sec. 3.3
real*4 Kmag0 ! (mag) ? Dereddened Ks magnitude, Vega, as inferred,
* Sec. 3.3
real*4 V_K0 ! (mag) ? Dereddened V-Ks color, Vega, as inferred,
* Sec. 3.3 (1)
integer*4 r_V_K0 ! Source of vmk-used, Sec. 3.1 (2)
integer*4 r_EV_Ks ! Source of evmk, Sec. 3.3 (3)
real*8 Per1 ! (d) ? Primary period (4)
real*4 Per2 ! (d) ? Secondary period
real*4 Per3 ! (d) ? Tertiary period
real*4 Per4 ! (d) ? Quaternary period
character*3 Code_P ! uncertainty code for period (5)
character*5 IRex_ ! Whether an IR excess is present or not, Sec. 3.4
character*10 IRex_wave ! Minimum wavelength of detected IR excess,
* Sec. 3.4
real*4 SED ! ? best-fit slope to all detections between 2 and
* 25 microns
character*3 SEDclass ! SED class (6)
character*3 dipper ! LC matches dipper characteristics, Sec. 4.5
character*3 burster ! LC matches burster characteristics, Sec. 4.5
character*6 multi_Per ! Single or multi-period star
character*3 dd ! LC and power spectrum matches double-dip
* characteristics, Sec. 4.5
character*3 ddmoving ! LC and power spectrum matches moving double-dip
* characteristics, Sec. 4.5
character*3 shape ! LC matches shape changer characteristics,
* Sec. 4.5
character*3 beater ! LC has beating visible, Sec. 4.5
character*3 peak ! power spectrum has a complex, structured peak
* and/or has a wide peak, Sec. 4.5
character*3 peak_r ! power spectrum has resolved close peaks,
* Sec. 4.5
character*3 Dist_r ! power spectrum has resolved distant peaks,
* Sec. 4.5
character*3 pulsator ! power spectrum and LC match pulsator
* characteristics, Sec. 4.5
character*5 scallop ! LC matches scallop or flux dip characteristics,
* Sec. 4.5, Appendix C
character*5 EB ! LC has characteristics of eclipsing binary,
* Sec. 4.5, Appendix C
*Note (1): Rounded to nearest 0.1 to emphasize the relatively low accuracy
*Note (2): References as follows:
* 1 = V measured directly from the literature (including SIMBAD) and Ks from
* 2MASS;
* 2 = V from the literature (see Sec. 3.1) and Ks from 2MASS;
* 3 = (V - Ks) inferred from Gaia DR1 G and Ks from 2MASS (see Sec. 3.1);
* 4 = (V -Ks) inferred from Pan-STARRS1 g and Ks from 2MASS (see Sec. 3.1);
* 6 = V inferred from well-populated optical SED and Ks from 2MASS
* (see Sec. 3.1);
* 7 = (V - Ks) inferred from Gaia DR2 G and Ks from 2MASS (see Sec. 3.1);
* -9 = no measure of (V-Ks) [Not used].
*Note (3): References as follows:
* 1 = dereddening from JHKs diagram (see Sec. 3.3);
* 2 = dereddening back to (V - Ks)0 expected for spectral type;
* 3 = dereddening from SED fits;
* 4 = used median E(V-Ks)=0 (see Sec. 3.3);
* -9 = no measure of E(V-Ks) [Not used].
*Note (4): Taken to be rotation period in cases where there is >1 period
*Note (5): Codes as follows:
* n = no worry, full confidence; by far the most common value;
* (n) = (no period);
* n? = progressively less confident periods
* y? = progressively less confident periods
*Note (6): Based on the SED slope between 2 and 25 microns.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Periods and supporting data for discarded/rejected
sources
integer*4 nr__1
parameter (nr__1=387) ! Number of records
character*243 ar__1 ! Full-size record
integer*8 TIC_1 ! TESS Input Catalog identifier
character*18 Name_1 ! Coordinate-based (RA/DE, J2000) name for target
character*23 OName_1 ! Alternate name for target
integer*8 Gaia_1 ! ? Gaia DR2 ID
real*8 Dist_1 ! (pc) ?=-9.00 Distance from Bailer-Jones+, 2018, I/347
character*6 Bin_1 ! Membership sample (gold, silver, or bronze)
real*4 Kmag_1 ! (mag) ? Ks magnitude, Vega, if observed
real*4 V_Ks_1 ! (mag) ? V-Ks color used, in Vega mags, observed or
* inferred; see text
real*4 EV_Ks_1 ! (mag) ? E(V-Ks) adopted for this star, Sec. 3.3
real*4 Kmag0_1 ! (mag) ? Dereddened Ks magnitude, Vega, as inferred,
* Sec. 3.3
real*4 V_K0_1 ! (mag) ? Dereddened V-Ks color, Vega, as inferred,
* Sec. 3.3 (1)
integer*4 r_V_K0_1 ! Source of vmk-used, Sec. 3.1 (2)
integer*4 r_EV_Ks_1 ! Source of evmk, Sec. 3.3 (3)
real*8 Per1_1 ! (d) ? Primary period (4)
real*4 Per2_1 ! (d) ? Secondary period
real*4 Per3_1 ! (d) ? Tertiary period
real*4 Per4_1 ! (d) ? Quaternary period
character*3 Code_P_1 ! uncertainty code for period (5)
character*5 IRex__1 ! Whether an IR excess is present or not, Sec. 3.4
character*10 IRex_wave_1 ! Minimum wavelength of detected IR excess,
* Sec. 3.4
real*4 SED_1 ! ? best-fit slope to all detections between 2 and
* 25 microns
character*3 SEDclass_1 ! SED class (6)
character*3 dipper_1 ! LC matches dipper characteristics, Sec. 4.5
character*3 burster_1 ! LC matches burster characteristics, Sec. 4.5
character*6 multi_Per_1 ! Single or multi-period star
character*3 dd_1 ! LC and power spectrum matches double-dip
* characteristics, Sec. 4.5
character*3 ddmoving_1 ! LC and power spectrum matches moving double-dip
* characteristics, Sec. 4.5
character*3 shape_1 ! LC matches shape changer characteristics,
* Sec. 4.5
character*3 beater_1 ! LC has beating visible, Sec. 4.5
character*3 peak_1 ! power spectrum has a complex, structured peak
* and/or has a wide peak, Sec. 4.5
character*3 peak_r_1 ! power spectrum has resolved close peaks,
* Sec. 4.5
character*3 Dist_r_1 ! power spectrum has resolved distant peaks,
* Sec. 4.5
character*3 pulsator_1 ! power spectrum and LC match pulsator
* characteristics, Sec. 4.5
character*5 scallop_1 ! LC matches scallop or flux dip characteristics,
* Sec. 4.5, Appendix C
character*5 EB_1 ! LC has characteristics of eclipsing binary,
* Sec. 4.5, Appendix C
*Note (1): Rounded to nearest 0.1 to emphasize the relatively low accuracy
*Note (2): References as follows:
* 1 = V measured directly from the literature (including SIMBAD) and Ks from
* 2MASS;
* 2 = V from the literature (see Sec. 3.1) and Ks from 2MASS;
* 3 = (V - Ks) inferred from Gaia DR1 G and Ks from 2MASS (see Sec. 3.1);
* 4 = (V -Ks) inferred from Pan-STARRS1 g and Ks from 2MASS (see Sec. 3.1);
* 6 = V inferred from well-populated optical SED and Ks from 2MASS
* (see Sec. 3.1);
* 7 = (V - Ks) inferred from Gaia DR2 G and Ks from 2MASS (see Sec. 3.1);
* -9 = no measure of (V-Ks) [Not used].
*Note (3): References as follows:
* 1 = dereddening from JHKs diagram (see Sec. 3.3);
* 2 = dereddening back to (V - Ks)0 expected for spectral type;
* 3 = dereddening from SED fits;
* 4 = used median E(V-Ks)=0 (see Sec. 3.3);
* -9 = no measure of E(V-Ks) [Not used].
*Note (4): Taken to be rotation period in cases where there is >1 period
*Note (5): Codes as follows:
* n = no worry, full confidence; by far the most common value;
* (n) = (no period);
* n? = progressively less confident periods
* y? = progressively less confident periods
*Note (6): Based on the SED slope between 2 and 25 microns.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Comparison to literature periods
integer*4 nr__2
parameter (nr__2=596) ! Number of records
character*317 ar__2 ! Full-size record
integer*8 TIC_2 ! TESS Input Catalog identifier
character*18 Name_2 ! Coordinate-based (RA/DE, J2000) name for target
real*8 Per_K ! (d) [0.31/18.4]? Period, Kiraga, 2012, J/AcA/62/67
real*4 Per_M ! (d) [0.31/8.58]? Period Mellon+, 2017, J/ApJ/844/66
real*4 Per_Gaia ! (d) [0.46/31.5]? Period Gaia DR2
real*8 Per_ASAS ! (d) [0.27/815]? Period ASAS-SN
real*4 Per_Lit ! (d) [0.44/11.3]? Period other literature
character*1 q_Per_Lit ! [:] Quality flag on plit
character*23 r_Per_Lit ! Citation for literature source of plit (1)
real*4 Per1_2 ! (d) [0.23/16.2]? first period reported here
real*4 Per2_2 ! (d) [0.17/6.52]? second period reported here
real*4 Per3_2 ! (d) [0.53/6.17]? third period reported here
real*4 Per4_2 ! (d) [4.33]? fourth period reported here
character*8 match ! Is this a match?
character*192 Note ! Notes
*Note (1): References as follows:
* Alfonso+12 = Alfonso-Garzon+, 2012, J/A+A/548/A79
* Batalha+98 = Batalha+, 1998, J/A+AS/128/561
* Christiansen+08 = Christiansen+, 2008, J/MNRAS/385/1749
* Desidera+15 = Desidera+, 2015, J/A+A/573/A126
* Distefano+16 = Distefano+, 2016, J/A+A/591/A43
* Donati+12 = Donati+, 2012MNRAS.425.2948D
* Drake+17 = Drake+, 2017, J/MNRAS/469/3688
* Fruth+13 = Fruth+, 2013, J/AJ/146/136
* Kospal+14 = Kospal+, 2014A&A...561A..61K
* Kiraga12 = Kiraga+, 2012, J/AcA/62/67
* Messina+10 = Messina+, 2010, J/A+A/520/A15
* Messina+10;Samus+17 = Messina+, 2010, J/A+A/520/A15 and Samus+, 2017, B/gcvs
* Messina+11 = Messina+, 2011, J/A+A/532/A10
* Samus+17 = Samus+, 2017, B/gcvs
* Wichmann+98 = Wichmann+, 1998A&A...330..521W
* Zuniga-Fernandez+21 = Zuniga-Fernandez+, 2021, J/A+A/645/A30
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! List of unusual objects
integer*4 nr__3
parameter (nr__3=151) ! Number of records
character*239 ar__3 ! Full-size record
integer*8 TIC_3 ! TESS Input Catalog identifier
character*18 Name_3 ! Coordinate-based (RA/DE, J2000) name for target
character*6 Bin_2 ! Membership sample (gold, silver, or bronze)
character*21 Conf ! Sample confidence (1)
character*180 Notes ! Notes
*Note (1): Sample confidence means how confident we are in both the placement
* in the category (e.g., 'scallop' or 'EB') and in the measured period
* being rotation.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! List of objects with timescales
integer*4 nr__4
parameter (nr__4=69) ! Number of records
character*74 ar__4 ! Full-size record
integer*8 TIC_4 ! TESS Input Catalog identifier
character*18 Name_4 ! Coordinate-based (RA/DE, J2000) name for target
character*44 Timescale ! (d) Timescale and/or Notes
C=============================================================================
C Loading file 'table2.dat' ! Periods and supporting data for UCL/LCC members
* with viable light curves
C Format for file interpretation
1 format(
+ I10,1X,A18,1X,A23,1X,I19,1X,F8.3,1X,A6,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F5.2,1X,F4.2,1X,I2,1X,I2,1X,F7.4,1X,F6.4,1X,F6.4,1X,
+ F6.4,1X,A3,1X,A5,1X,A10,1X,F6.2,1X,A3,1X,A3,1X,A3,1X,A6,1X,A3,
+ 1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A5,1X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,3714
read(1,'(A243)')ar__
read(ar__,1)
+ TIC,Name,OName,Gaia,Dist,Bin,Kmag,V_Ks,EV_Ks,Kmag0,V_K0,
+ r_V_K0,r_EV_Ks,Per1,Per2,Per3,Per4,Code_P,IRex_,IRex_wave,SED,
+ SEDclass,dipper,burster,multi_Per,dd,ddmoving,shape,beater,
+ peak,peak_r,Dist_r,pulsator,scallop,EB
if(ar__(55:73) .EQ. '') Gaia = iNULL__
if(ar__(91:95) .EQ. '') Kmag = rNULL__
if(ar__(97:100) .EQ. '') V_Ks = rNULL__
if(ar__(102:106) .EQ. '') EV_Ks = rNULL__
if(ar__(108:112) .EQ. '') Kmag0 = rNULL__
if(ar__(114:117) .EQ. '') V_K0 = rNULL__
if(ar__(125:131) .EQ. '') Per1 = rNULL__
if(ar__(133:138) .EQ. '') Per2 = rNULL__
if(ar__(140:145) .EQ. '') Per3 = rNULL__
if(ar__(147:152) .EQ. '') Per4 = rNULL__
if(ar__(175:180) .EQ. '') SED = rNULL__
c ..............Just test output...........
write(6,1)
+ TIC,Name,OName,Gaia,Dist,Bin,Kmag,V_Ks,EV_Ks,Kmag0,V_K0,
+ r_V_K0,r_EV_Ks,Per1,Per2,Per3,Per4,Code_P,IRex_,IRex_wave,SED,
+ SEDclass,dipper,burster,multi_Per,dd,ddmoving,shape,beater,
+ peak,peak_r,Dist_r,pulsator,scallop,EB
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Periods and supporting data for discarded/rejected
* sources
C Format for file interpretation
2 format(
+ I10,1X,A18,1X,A23,1X,I19,1X,F8.3,1X,A6,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F5.2,1X,F4.2,1X,I2,1X,I2,1X,F7.4,1X,F6.4,1X,F6.4,1X,
+ F6.4,1X,A3,1X,A5,1X,A10,1X,F6.2,1X,A3,1X,A3,1X,A3,1X,A6,1X,A3,
+ 1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A5,1X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,387
read(1,'(A243)')ar__1
read(ar__1,2)
+ TIC_1,Name_1,OName_1,Gaia_1,Dist_1,Bin_1,Kmag_1,V_Ks_1,
+ EV_Ks_1,Kmag0_1,V_K0_1,r_V_K0_1,r_EV_Ks_1,Per1_1,Per2_1,
+ Per3_1,Per4_1,Code_P_1,IRex__1,IRex_wave_1,SED_1,SEDclass_1,
+ dipper_1,burster_1,multi_Per_1,dd_1,ddmoving_1,shape_1,
+ beater_1,peak_1,peak_r_1,Dist_r_1,pulsator_1,scallop_1,EB_1
if(ar__1(55:73) .EQ. '') Gaia_1 = iNULL__
if(ar__1(91:95) .EQ. '') Kmag_1 = rNULL__
if(ar__1(97:100) .EQ. '') V_Ks_1 = rNULL__
if(ar__1(102:106) .EQ. '') EV_Ks_1 = rNULL__
if(ar__1(108:112) .EQ. '') Kmag0_1 = rNULL__
if(ar__1(114:117) .EQ. '') V_K0_1 = rNULL__
if(ar__1(125:131) .EQ. '') Per1_1 = rNULL__
if(ar__1(133:138) .EQ. '') Per2_1 = rNULL__
if(ar__1(140:145) .EQ. '') Per3_1 = rNULL__
if(ar__1(147:152) .EQ. '') Per4_1 = rNULL__
if(ar__1(175:180) .EQ. '') SED_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ TIC_1,Name_1,OName_1,Gaia_1,Dist_1,Bin_1,Kmag_1,V_Ks_1,
+ EV_Ks_1,Kmag0_1,V_K0_1,r_V_K0_1,r_EV_Ks_1,Per1_1,Per2_1,
+ Per3_1,Per4_1,Code_P_1,IRex__1,IRex_wave_1,SED_1,SEDclass_1,
+ dipper_1,burster_1,multi_Per_1,dd_1,ddmoving_1,shape_1,
+ beater_1,peak_1,peak_r_1,Dist_r_1,pulsator_1,scallop_1,EB_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Comparison to literature periods
C Format for file interpretation
3 format(
+ I10,1X,A18,1X,F7.4,1X,F5.3,1X,F6.3,1X,F7.3,1X,F6.3,1X,A1,1X,
+ A23,1X,F6.3,1X,F5.3,1X,F5.3,1X,F4.2,1X,A8,1X,A192)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,596
read(1,'(A317)')ar__2
read(ar__2,3)
+ TIC_2,Name_2,Per_K,Per_M,Per_Gaia,Per_ASAS,Per_Lit,q_Per_Lit,
+ r_Per_Lit,Per1_2,Per2_2,Per3_2,Per4_2,match,Note
if(ar__2(31:37) .EQ. '') Per_K = rNULL__
if(ar__2(39:43) .EQ. '') Per_M = rNULL__
if(ar__2(45:50) .EQ. '') Per_Gaia = rNULL__
if(ar__2(52:58) .EQ. '') Per_ASAS = rNULL__
if(ar__2(60:65) .EQ. '') Per_Lit = rNULL__
if(ar__2(93:98) .EQ. '') Per1_2 = rNULL__
if(ar__2(100:104) .EQ. '') Per2_2 = rNULL__
if(ar__2(106:110) .EQ. '') Per3_2 = rNULL__
if(ar__2(112:115) .EQ. '') Per4_2 = rNULL__
c ..............Just test output...........
write(6,3)
+ TIC_2,Name_2,Per_K,Per_M,Per_Gaia,Per_ASAS,Per_Lit,q_Per_Lit,
+ r_Per_Lit,Per1_2,Per2_2,Per3_2,Per4_2,match,Note
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! List of unusual objects
C Format for file interpretation
4 format(I10,1X,A18,1X,A6,1X,A21,1X,A180)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,151
read(1,'(A239)')ar__3
read(ar__3,4)TIC_3,Name_3,Bin_2,Conf,Notes
c ..............Just test output...........
write(6,4)TIC_3,Name_3,Bin_2,Conf,Notes
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! List of objects with timescales
C Format for file interpretation
5 format(I10,1X,A18,1X,A44)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,69
read(1,'(A74)')ar__4
read(ar__4,5)TIC_4,Name_4,Timescale
c ..............Just test output...........
write(6,5)TIC_4,Name_4,Timescale
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end