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