Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/131/2986 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 2024-Apr-18
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/131/2986 Multiplicity of contact binaries (Pribulla+, 2006)
*================================================================================
*Contact binaries with additional components. I. The extant data.
* Pribulla T., Rucinski S.M.
* <Astron. J., 131, 2986-3007 (2006)>
* =2006AJ....131.2986P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Multiplicity of contact binary stars: summary data
integer*4 nr__
parameter (nr__=181) ! Number of records
character*111 ar__ ! Full-size record
character*9 Name (nr__) ! The General Catalog of Variable Stars name
* (Cat. II/250) (1)
character*1 n_Name (nr__) ! [*] Indicates a note detailed in table3
real*4 Vmag (nr__) ! (mag) Johnson or transformed Hipparcos H_p_ maximum
* magnitude
character*10 WDS (nr__) ! Washington Double Star catalogue designation
* (Cat. B/wds) (2)
integer*4 PA (nr__) ! (deg) ? Position angle of the secondary/fainter
* component
real*4 Sep (nr__) ! (arcsec) ? Angular separation of the components (3)
integer*4 d (nr__) ! (AU) ? Projected separation (4)
character*1 f_d (nr__) ! [*] Flag for large uncertainty on d (5)
real*4 DelV (nr__) ! (mag) ? Visual brightness difference of
* components (6)
character*1 H59 (nr__) ! [CGX] Hipparcos multiple system flag H59 (7)
character*1 H61 (nr__) ! [ABCS] Hipparcos solution quality flag H61 (8)
real*4 e_pi (nr__) ! (mas) ? Hipparcos parallax error
real*4 e_mu (nr__) ! (mas/yr) ? Proper-motion error in TYCHO 2 catalog
* (Cat. I/259) (9)
real*4 beta (nr__) ! ? Ratio of apparent luminosities (10)
integer*4 dRV (nr__) ! (km/s) ? Difference of third component velocity and
* systemic velocity of the underlying binary
character*1 f_dRV (nr__) ! [vs] Flag on dRV variability (11)
integer*4 P3 (nr__) ! ? Orbital period of the third component as
* determined from LITE
real*4 asin_i (nr__) ! (AU) ? Projected semi-major axis of the eclipsing
* pair (12)
real*4 dB_V (nr__) ! (mag) ? Difference between TYCHO-2 B-V color and
* the short-period, blue envelope (SPBE)
* (see the text)
real*4 logFX_Fbol (nr__) ! ([-]) ? Log of the RASS X-ray flux to the
* bolometric flux (see Section 4.3)
character*1 Occ (nr__) ! [M] Indicates a Moon occultation (13)
character*1 F1 (nr__) ! [DPN] Visual binary present in the WDS catalog
* flag (14)
character*1 F2 (nr__) ! [DPN-] CFHT Adaptative Optics observations
* flag (14)
character*1 F3 (nr__) ! [DSN-] Hipparcos Catalog flags (15)
character*1 F4 (nr__) ! [SN-] {sigma}_{mu}_ flag (16)
character*1 F5 (nr__) ! [DSN-] Spectroscopic detection (17)
character*1 F6 (nr__) ! [DSN-] Light-time effect (LITE) indication (18)
character*1 F7 (nr__) ! [SN-] Short-period, blue envelope (SPBE)
* indication (19)
character*1 F8 (nr__) ! [LEN] RASS (<IX/10>,<IX/29>) X-ray to visual
* flux ratio flag (20)
real*4 W (nr__) ! Combined weight (21)
character*1 T3 (nr__) ! [Y] Final evaluation of triplicity (22)
*Note (1): Systems without GCVS designations:
* V2 Mon = TYC 4824-00153-1,
* V3 Mon = TYC 4835-01947-1,
* V2 Peg = TYC 1720-00658-1.
*Note (2): Designation of the system (if it is visual binary or common
* proper-motion pair) in the most recent on-line version of the catalog
* (http://ad.usno.navy.mil/wds/wds.html) or CFHT indicating a new
* detection, as given in Table 4.
*Note (3): Angular separation of components of LL Vel is <0.03".
*Note (4): Computed from angular separation and Hipparcos parallax.
*Note (5): A "*" indicates a relative error of parallax larger than 0.2 or a
* negative parallax, the distance estimated using absolute magnitude
* calibration of Rucinski & Duerbeck (1997PASP..109.1340R) or estimated
* otherwise (see text).
*Note (6): Visual brightness difference of components of ES Lib is not given
* in the WDS catalogue.
*Note (7): H59 Hipparcos flag represents the type of astrometric solution,
* defined as follows:
* C = linear motion,
* G = acceleration term necessary (i.e. very probably an astrometric binary)
* X = stochastic solution (i.e. very probably an astrometric binary with
* P<3years);
*Note (8): H61 Hipparcos flag represents the quality of the orbital solution:
* grades ABCD (A=best) or suspected non-single flag S.
*Note (9): e_mu = sqrt[({sigma}({mu}_{alpha})_cos({delta})^2^
* +{sigma}^2^({mu}_{delta})_].
*Note (10): l_3_/(l_1_+l_2_) determined from broadening functions (see the text).
*Note (11): Flag on DelRV as follows:
* v = variable velocity;
* s = similar velocity.
*Note (12): Around the common center of gravity determined using the LITE.
*Note (13): Occultation within 5{deg}9' from the ecliptic.
*Note (14): Flags as follows:
* D = for systems with projected separation less than 2000AU,
* P = for systems with separation between 2000 and 20000AU,
* N = system not present in WDS or with separation larger than 20000AU.
* - = system not observed.
*Note (15): Hipparcos Catalog flags as follows:
* N = no detection
* - = not observed or insufficient data
* D = if X or G in H59,
* S = if S or C in H61. (Three systems not recognized before as the "S" cases,
* but having parallax errors larger than 3{sigma} relative to the error
* median, have been added here: UX Eri, V1363 Eri, DY Cet.)
*Note (16): The proper-motion error is defined as follows:
* S = if proper-motion error higher than 3-{sigma} from the median.
* N = no detection
* - = not observed or insufficient data
*Note (17): Spectroscopic detection as follows:
* D = detection in broadening functions or through spectral features
* N = no detection in broadening functions or through spectral features
* - = not observed or insufficient data
* S = large changes of systemic velocity
*Note (18): LITE indication, defined as follows:
* D = a stable orbit with more than 1.5 cycles covered,
* S = a somewhat unstable solution or less than 1.5 cycles covered,
* N = no detection
* - = fewer than 25 pe or CCD minima or the interval of observations <10yr.
*Note (19): SPBE indication, defined as follows:
* S = if the system is bluer than the corresponding blue envelope color.
* N = no detection
* - = not observed or insufficient data
*Note (20): RASS X-ray to visual flux ratio, defined as follows:
* L = when a late-type companion is indicated by a large logX ratio in
* an early-type contact binary,
* E = when the ratio is below the lower envelope of detection limit
* indicating an early-type third component.
* N = no detection
*Note (21): The sum of weights of individual detection techniques (any
* detection counts with weight 1.0, any possibility or suspicion with 0.5)
*Note (22): A "Y" (=Yes) indicates the sum of weights is >= 1.0
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Notes on selected contact binaries
integer*4 nr__1
parameter (nr__1=72) ! Number of records
character*80 ar__1 ! Full-size record
character*9 Name_1 (nr__1) ! The General Catalog of Variable Stars name
* (Cat. II/250)
character*70 Note (nr__1) ! Text of note
C=============================================================================
C Loading file 'table1.dat' ! Multiplicity of contact binary stars: summary data
C Format for file interpretation
1 format(
+ A9,1X,A1,1X,F5.2,1X,A10,1X,I3,1X,F5.2,1X,I5,1X,A1,1X,F4.2,1X,
+ A1,A1,1X,F4.2,1X,F4.2,1X,F4.2,1X,I3,1X,A1,1X,I2,1X,F3.1,1X,
+ F5.2,1X,F5.2,1X,A1,1X,A1,A1,A1,A1,A1,A1,A1,A1,1X,F3.1,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,181
read(1,'(A111)')ar__
read(ar__,1)
+ Name(i__),n_Name(i__),Vmag(i__),WDS(i__),PA(i__),Sep(i__),
+ d(i__),f_d(i__),DelV(i__),H59(i__),H61(i__),e_pi(i__),
+ e_mu(i__),beta(i__),dRV(i__),f_dRV(i__),P3(i__),asin_i(i__),
+ dB_V(i__),logFX_Fbol(i__),Occ(i__),F1(i__),F2(i__),F3(i__),
+ F4(i__),F5(i__),F6(i__),F7(i__),F8(i__),W(i__),T3(i__)
if(ar__(30:32) .EQ. '') PA(i__) = iNULL__
if(ar__(34:38) .EQ. '') Sep(i__) = rNULL__
if(ar__(40:44) .EQ. '') d(i__) = iNULL__
if(ar__(48:51) .EQ. '') DelV(i__) = rNULL__
if(ar__(56:59) .EQ. '') e_pi(i__) = rNULL__
if(ar__(61:64) .EQ. '') e_mu(i__) = rNULL__
if(ar__(66:69) .EQ. '') beta(i__) = rNULL__
if(ar__(71:73) .EQ. '') dRV(i__) = iNULL__
if(ar__(77:78) .EQ. '') P3(i__) = iNULL__
if(ar__(80:82) .EQ. '') asin_i(i__) = rNULL__
if(ar__(84:88) .EQ. '') dB_V(i__) = rNULL__
if(ar__(90:94) .EQ. '') logFX_Fbol(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),n_Name(i__),Vmag(i__),WDS(i__),PA(i__),Sep(i__),
+ d(i__),f_d(i__),DelV(i__),H59(i__),H61(i__),e_pi(i__),
+ e_mu(i__),beta(i__),dRV(i__),f_dRV(i__),P3(i__),asin_i(i__),
+ dB_V(i__),logFX_Fbol(i__),Occ(i__),F1(i__),F2(i__),F3(i__),
+ F4(i__),F5(i__),F6(i__),F7(i__),F8(i__),W(i__),T3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Notes on selected contact binaries
C Format for file interpretation
2 format(A9,1X,A70)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,72
read(1,'(A80)')ar__1
read(ar__1,2)Name_1(i__),Note(i__)
c ..............Just test output...........
write(6,2)Name_1(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end