FORTRAN Generation
(/./ftp/cats/V/119)

Conversion of standardized ReadMe file for file /./ftp/cats/V/119 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-Mar-19
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. V/119          Catalogue of field contact binary stars         (Pribulla+, 2003)
*================================================================================
*Catalogue of the field contact binary stars.
*    Pribulla T., Kreiner J.M., Tremko J.
*   <Contr. Astron. Obs. Skalnate Pleso, 33, 38 (2003)>
*   =2003CoSka..33...38P
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'catalog.dat'	! *Catalogue of the field binary stars

      integer*4 nr__
      parameter (nr__=361)	! Number of records
      character*189 ar__   	! Full-size record

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg      (nr__) ! (deg) Right Ascension J2000
      real*8        DEdeg      (nr__) ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*9   Star       (nr__) ! Star name
      character*1   n_Star     (nr__) ! [*] *: note detailed in "notes.dat" file
      character*3   Type       (nr__) ! Type of the system (AWEB+PT) (1)
      integer*4     G          (nr__) ! [1/5] Group number (G1 to G5)  (2)
      real*8        HJD0_2400000(nr__) ! (d) Ephemeris of minimum, minus 2400000
      real*8        Per        (nr__) ! (d) Period
      integer*4     N          (nr__) ! Number of photoelectric and CCD minima
*                                   available
      character*3   Note1      (nr__) ! Note on orbital period change (3)
      real*4        qph        (nr__) ! ? Photometric mass ratio, m2/m1
      real*4        f          (nr__) ! ? Fill-out factor
      real*4        i          (nr__) ! (deg) ? Inclination
      integer*4     T1         (nr__) ! (K) ? Temperature of the primary
      integer*4     T2         (nr__) ! (K) ? Temperature of the secondary
      real*4        DeltaM     (nr__) ! (mag) ? Magnitude of the O'Connell effect (usually
*                                    explained by spots on the stellar surface)
      character*3   Code       (nr__) ! Code used for the determination of the
*                                   photometric elements (4)
      real*4        qsp        (nr__) ! ? Spectroscopic mass ratio (m2/m1)
      real*4        Mass       (nr__) ! (solMass) ? Mass, (m1+m2)*(sin)i^3^
      character*6   SpType     (nr__) ! MK Spectral type of the primary component
      integer*4     RAh        (nr__) ! (h) Right ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Right ascension (J2000)
      integer*4     RAs        (nr__) ! (s) Right ascension (J2000)
      character*1   DE_        (nr__) ! Declination sign (J2000)
      integer*4     DEd        (nr__) ! (deg) Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Declination (J2000)
      integer*4     DEs        (nr__) ! (arcsec) Declination (J2000)
      real*4        plx        (nr__) ! (mas) ? Parallax
      real*4        e_plx      (nr__) ! (mas) ? rms uncertainty on plx
      real*4        Vmax       (nr__) ! (mag) ? Maximum visual magnitude
      real*4        Vmin       (nr__) ! (mag) ? Minimum visual magnitude
      character*7   r_Star     (nr__) ! References, in refs.dat file
      character*5   Notes      (nr__) ! Additional informations (5)
*Note (1): The systems are classified according to the Binnendijk
*    (1970VA.....12..217B) scheme: the systems with the primary minimum
*    being the transit are of the A type, systems with the primary minimum
*    being the occultation are of the W type. Contact binaries with hot
*    components of the OB spectral type are assigned as E type. Systems in
*    which the components are in physical but not in thermal contact
*    (represented by W Crv) are denoted by B (reminding {beta} Lyrae LC).
*    The type of the eclipses is also given: P = partial, T = total.
*Note (2):
*    The contact binary stars are divided into five groups depending
*    on available data:
*    1: systems with both spectroscopic and photometric orbits
*    2: systems without reliable photometry but with spectroscopic orbit
*    3: systems with photometric  orbit
*    4: systems with unresolved photoelectric or CCD light curve
*    5: systems without good CCD or photoelectric light curves but with
*       frequent observations of minima which allowed to determine
*       reliable ephemeris.
*Note (3): Observed orbital period chance are coded as follows:
*     v: orbital period is variable
*    v!: fast or complicated variation of the orbital period
*    3b: orbital period is modulated by the light-time effect
*   3b?: orbital period is probably modulated by the light-time effect
*     ":" uncertain  or spurious period
*    "::" very uncertain or unreliable period
*Note (4): The following code is used for the determination of the
*          photometric elements:
*    WD: various versions of the W&D program
*        (see Wilson & Devinney, 1971ApJ...166..605W)
*    L2: LIGHT2 software (Hill & Rucinski, 1993, in Light Curve Modeling
*        of Eclipsing Binary Stars, E.F. Milone (ed.), Springer-Verlag, 135)
*    LT: LIGHT software (Hill, 1979PDAO...15..298H)
*    BM: Binary Maker 2.0 (Bradstreet, 1993, Binary Maker 2.0, Dept. of
*        Physics, Eastern College, St. Davis)
*    W3: WUMA3 code (Rucinski, 1976PASP...88..244R)
*   WDp: Adapted W&D code using the Price algorithm (see Barone et al.,
*        1993ApJ...407..237B), grid search or genetic algorithms
*   WDr: W&D code taking into account radiation pressure (see Drechsel
*        et al., 1995A&A...294..723D)
*    BS: BINSYN code (see Vinko et al., 1996MNRAS.280..489V) using the
*        Price algorithm
*    ME: Maximum entropy mapping (see Hendry & Mochnacki, 2000ApJ...531..467H)
*    DJ: Djurasevic's (1992Ap&SS.196..241D, 1992Ap&SS.197...17D) code
*    EL: Wood's (1972, A Computer Program for Modeling Non-Spherical
*        Eclipsing Binary Systems, Greenbelt, U.S.A., Goddard Space
*        Flight Center) ellipsoidal model
*    KO: Kopal's (1982Ap&SS..81..123K, 1982Ap&SS..81..411K) Fourier
*        analysis method
*    An additional parameter "s" denotes objects were a simultaneous fit
*        of the radial velocities and light curves has been performed
*Note (5): Additional informations has the following meaning:
*    H: Hipparcos discovery
*    R: ROTSE discovery
*   mv: member of the visual pair (at least triple system)
*   ms: third or multiple components to the system found in the spectra
*    A: unsure equatorial coordinates of the system (uncertainty >0.1')
*    C: star in cluster
*    7: For BH CMi only, not explicited in the paper

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'notes.dat'	! Individual notes

      integer*4 nr__1
      parameter (nr__1=88)	! Number of records
      character*80 ar__1  	! Full-size record

      character*9   Star_1     (nr__1) ! Star name
      character*70  Note       (nr__1) ! Text of the note

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'refs.dat'	! References

      integer*4 nr__2
      parameter (nr__2=312)	! Number of records
      character*183 ar__2  	! Full-size record

      integer*4     Ref        (nr__2) ! Reference number
      character*19  BibCode    (nr__2) ! BibCode
      character*34  Aut        (nr__2) ! Author's name
      character*124 Com        (nr__2) ! Comments

C=============================================================================

C  Loading file 'catalog.dat'	! *Catalogue of the field binary stars

C  Format for file interpretation

    1 format(
     +  A9,A1,1X,A3,2X,I1,2X,F10.4,2X,F10.8,2X,I3,2X,A3,2X,F6.4,2X,
     +  F4.2,2X,F5.2,1X,I5,2X,I5,2X,F5.2,2X,A3,2X,F5.3,3X,F6.3,3X,A6,
     +  4X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,3X,F5.2,3X,F4.2,3X,
     +  F6.3,1X,F5.2,2X,A7,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,361
        read(1,'(A189)')ar__
        read(ar__,1)
     +  Star(i__),n_Star(i__),Type(i__),G(i__),HJD0_2400000(i__),
     +  Per(i__),N(i__),Note1(i__),qph(i__),f(i__),i(i__),T1(i__),
     +  T2(i__),DeltaM(i__),Code(i__),qsp(i__),Mass(i__),SpType(i__),
     +  RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
     +  DEs(i__),plx(i__),e_plx(i__),Vmax(i__),Vmin(i__),r_Star(i__),
     +  Notes(i__)
        if(ar__(54:59) .EQ. '') qph(i__) = rNULL__
        if(ar__(62:65) .EQ. '') f(i__) = rNULL__
        if(ar__(68:72) .EQ. '') i(i__) = rNULL__
        if(ar__(74:78) .EQ. '') T1(i__) = iNULL__
        if(ar__(81:85) .EQ. '') T2(i__) = iNULL__
        if(ar__(88:92) .EQ. '') DeltaM(i__) = rNULL__
        if(ar__(100:104) .EQ. '') qsp(i__) = rNULL__
        if(ar__(108:113) .EQ. '') Mass(i__) = rNULL__
        if(ar__(148:152) .EQ. '') plx(i__) = rNULL__
        if(ar__(156:159) .EQ. '') e_plx(i__) = rNULL__
        if(ar__(163:168) .EQ. '') Vmax(i__) = rNULL__
        if(ar__(170:174) .EQ. '') Vmin(i__) = rNULL__
        RAdeg(i__) = rNULL__
        DEdeg(i__) = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
        if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
        if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
        if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
        if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
        if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
        if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c    ..............Just test output...........
        write(6,1)
     +  Star(i__),n_Star(i__),Type(i__),G(i__),HJD0_2400000(i__),
     +  Per(i__),N(i__),Note1(i__),qph(i__),f(i__),i(i__),T1(i__),
     +  T2(i__),DeltaM(i__),Code(i__),qsp(i__),Mass(i__),SpType(i__),
     +  RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
     +  DEs(i__),plx(i__),e_plx(i__),Vmax(i__),Vmin(i__),r_Star(i__),
     +  Notes(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'notes.dat'	! Individual notes

C  Format for file interpretation

    2 format(A9,1X,A70)

C  Effective file loading

      open(unit=1,status='old',file=
     +'notes.dat')
      write(6,*) '....Loading file: notes.dat'
      do i__=1,88
        read(1,'(A80)')ar__1
        read(ar__1,2)Star_1(i__),Note(i__)
c    ..............Just test output...........
        write(6,2)Star_1(i__),Note(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    3 format(I3,1X,A19,1X,A34,1X,A124)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,312
        read(1,'(A183)')ar__2
        read(ar__2,3)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c    ..............Just test output...........
        write(6,3)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end