FORTRAN Generation
(/./ftp/cats/J/A_A/460/695)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/460/695 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-29
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/460/695    Search for Associations Containing Young stars  (Torres+, 2006)
*================================================================================
*Search for Associations Containing Young stars (SACY).
*I. Sample & Searching Method.
*    Torres C.A.O., Quast G.R., da Silva L., de la Reza R., Melo C.H.F.,
*    Sterzik M.
*   <Astron. Astrophys. 460, 695 (2006)>
*   =2006A&A...460..695T
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Astrometric, spectroscopic and photometric data
                              of the SACY sample

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

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

      integer*4     Seq        (nr__) ! Running number
      character*16  v1RXS      (nr__) ! 1RXS (<IX/10>,<IX/29>) designation (1)
      character*16  Name       (nr__) ! Hipparcos (<I/239>), Tycho-2 (<II/259>) or
*                                 UCAC2 (<I/289>) identification (2)
      character*13  OName      (nr__) ! Other identification (HD, DM, other) (3)
      integer*4     RAh        (nr__) ! (h) Right ascension (J2000.0)
      integer*4     RAm        (nr__) ! (min) Right ascension (J2000.0)
      real*4        RAs        (nr__) ! (s) Right ascension (J2000.0)
      character*1   DE_        (nr__) ! Declination sign (J2000.0)
      integer*4     DEd        (nr__) ! (deg) Declination (J2000.0)
      integer*4     DEm        (nr__) ! (arcmin) Declination (J2000.0)
      integer*4     DEs        (nr__) ! (arcsec) Declination (J2000.0)
      real*8        pmRA       (nr__) ! (mas/yr) Proper motion in RA (pmRA*cosDE)
      real*8        pmDE       (nr__) ! (mas/yr) Proper motion in DE (pmDE)
      character*3   r_pmRA     (nr__) ! Source of proper motions (4)
      real*4        HRV        (nr__) ! (km/s) Heliocentric radial velocity (5)
      character*1   f_HRV      (nr__) ! [:] Quality flag for RV
      real*4        e_HRV      (nr__) ! (km/s) ? rms of the observed RV
      integer*4     Nobs       (nr__) ! ? Number of observations (6)
      character*1   f_Nobs     (nr__) ! [*] *: RV value comes from the literature (7)
      character*3   r_HRV      (nr__) ! ? Source of radial velocity (8)
      character*5   n_HRV      (nr__) ! [D,SB1234 ?] Multiplicity flag (9)
      real*4        Sep        (nr__) ! (arcsec) ? Angular separation of visual binaries (10)
      real*4        dV         (nr__) ! (mag) ? Magnitude difference of components (11)
      character*1   u_dV       (nr__) ! [?] Uncertainty flag on dV
      real*4        vsini      (nr__) ! (km/s) ? Rotational velocity
      real*4        e_vsini    (nr__) ! (km/s) ? rms uncertainty on vsini
      character*3   r_vsini    (nr__) ! ? Source of vsini (8)
      real*4        Plx        (nr__) ! (mas) ? Trigonometric parallax
      real*4        e_Plx      (nr__) ! (mas) ? Standard error in Plx
      character*3   r_Plx      (nr__) ! Source of Plx (8)
      real*4        Vmag       (nr__) ! (mag) Johnson V magnitude
      character*1   f_Vmag     (nr__) ! [V] Quality flag indicating variability (12)
      real*4        B_V        (nr__) ! (mag) ? Johnson B-V colour index
      real*4        V_I        (nr__) ! (mag) ? Cousins V-I colour index
      character*3   r_Vmag     (nr__) ! Source of photometry (13)
      character*1   n_Vmag     (nr__) ! [*] *: photometry corrected for companions (14)
      character*8   SpType     (nr__) ! MK spectral type
      integer*4     EWLi       (nr__) ! (0.1pm) ? LiI equivalent width in milli-Angstroms (15)
      character*1   f_EWLi     (nr__) ! [:?] Quality flag for EWLi
      real*4        EWHa       (nr__) ! (0.1nm) ? H{alpha} line emission in Angstroms
      character*5   r_EWLi     (nr__) ! Source of spectral information (16)
*Note (1):
*    ROSAT All-Sky Survey Catalogue source name, in general from the Bright
*    Source Catalogue (<IX/10>) but sometimes refers to the Faint Source
*    Catalogue (<IX/29>). A repeated entry means more than one observed
*    possible optical counterpart. The same source may have entries in both
*    Tables 3 and 4, but the best optical candidate may be in any of the
*    tables.
*Note (2):
*    Order of preference of the identification is Hipparcos, TYCHO-2 or
*    UCAC2, independently of the astrometric data used.
*Note (3):
*    Order of preference is HD, DM, and other. GSC is only used if there is
*    no TYCHO-2 identification for the star.
*Note (4):
*    Order of preference of the sources of the proper motions is the same
*    as of the Name, adding other sources, but in some cases we used the
*    best value from our judgment.
*Note (5):
*    Radial velocities from SB2 observed once or SB1 are mean values, and
*    may not represent the systemic velocity. When there is a systemic
*    velocity from the literature, it is the used value.
*Note (6):
*    Null refers to no published number of observations in the literature
*    cited in r_HRV. Nobs from SACY includes references N=25 and 26.
*Note (7):
*    The * means that this is the number of spectral observations in the
*    SACY project, but the systemic HRV value was taken from the literature.
*    For V1005 Ori the non variable HRV is from literature as our spectrum
*    was not properly calibrated.
*Note (8):
*    If there is no number, the radial velocity is from this work.
*Note (9):
*    D means visual binary; sometimes the components were observed
*    separately and there is another entry, designated by A, B, C...,
*    in Table 3 or Table 4. If the bright component is also a SB the D
*    is given in the entry of the faint one.
*
*    A ? means possible evidence for radial velocity variations, sometimes
*    by comparing with published data. In some cases it is difficult to
*    decide between high rotation (+ spots?) and SB2. For other SB2? with
*    low rotation and large e_HRV the star should be considered at least
*    a SB1. For close visual doubles observed together, one of the
*    spectroscopic components may be also the visual one, mainly in
*    SB3 cases. For SBn all other data of this line are for the bright
*    component. In the future we will publish the data for spectroscopic
*    binaries found in the SACY project.
*Note (10):
*    Separation between the visual components, even if in n_HRV there is a
*    SBn. This means there is D for a faint companion in n_HRV in another
*    entry.
*Note (11):
*    Magnitude difference between faint companions and the bright one. A
*    negative value means that the faint companions together are brighter
*    than the primary, or the so called primary is a fainter companion with
*    sharper lines.
*
*    For visual binaries the negative value means that primary is also a
*    SBn, and the visual difference is given in this entry or the primary
*    is in the other table.
*Note (12):
*    V: photometric variability.
*Note (13): Source of photometry:
*       H = origin of the photometry from Hipparcos (<I/239>)
*       T = from TYCHO2 (<I/259>) converted to Johnson
*   blank = estimated magnitude
*Note (14):
*    for the proposes of the SACY the photometric measures for close visual
*    doubles and SBs are corrected to take into account the secondaries. As
*    these values are model dependent they must be used with caution.
*Note (15):
*    For SB2 or SB3 the values are tentatively corrected for the continuum
*    flux ratios.
*Note (16): For observations from this paper:
*       F = observation with ESO/FEROS
*       d = observation with the OPD coude spectrograph

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

C  Declarations for 'table4.dat'	! Astrometric, spectroscopic and photometric data
                              of other observed stars

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

C  J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1    (nr__1) ! (deg) Right Ascension J2000.0
      real*8        DEdeg_1    (nr__1) ! (deg)     Declination J2000.0
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     Seq_1      (nr__1) ! Running number
      character*16  v1RXS_1    (nr__1) ! 1RXS (<IX/10>,<IX/29>) designation (1)
      character*16  Name_1     (nr__1) ! Hipparcos (<I/239>), Tycho-2 (<II/259>) or
*                                 UCAC2 (<I/289>) identification (2)
      character*13  OName_1    (nr__1) ! Other identification (HD, DM, other) (3)
      integer*4     RAh_1      (nr__1) ! (h) Right ascension (J2000.0)
      integer*4     RAm_1      (nr__1) ! (min) Right ascension (J2000.0)
      real*4        RAs_1      (nr__1) ! (s) Right ascension (J2000.0)
      character*1   DE__1      (nr__1) ! Declination sign (J2000.0)
      integer*4     DEd_1      (nr__1) ! (deg) Declination (J2000.0)
      integer*4     DEm_1      (nr__1) ! (arcmin) Declination (J2000.0)
      integer*4     DEs_1      (nr__1) ! (arcsec) Declination (J2000.0)
      real*8        pmRA_1     (nr__1) ! (mas/yr) Proper motion in RA (pmRA*cosDE)
      real*8        pmDE_1     (nr__1) ! (mas/yr) Proper motion in DE (pmDE)
      character*3   r_pmRA_1   (nr__1) ! Source of proper motions (4)
      real*4        HRV_1      (nr__1) ! (km/s) Heliocentric radial velocity (5)
      character*1   f_HRV_1    (nr__1) ! [:] Quality flag for RV
      real*4        e_HRV_1    (nr__1) ! (km/s) ? rms of the observed RV
      integer*4     Nobs_1     (nr__1) ! ? Number of observations (6)
      character*1   f_Nobs_1   (nr__1) ! [*] *: RV value comes from the literature (7)
      character*3   r_HRV_1    (nr__1) ! ? Source of radial velocity (8)
      character*5   n_HRV_1    (nr__1) ! [D,SB1234 ?] Multiplicity flag (9)
      real*4        Sep_1      (nr__1) ! (arcsec) ? Angular separation of visual binaries (10)
      real*4        dV_1       (nr__1) ! (mag) ? Magnitude difference of components (11)
      character*1   u_dV_1     (nr__1) ! [?] Uncertainty flag on dV
      real*4        vsini_1    (nr__1) ! (km/s) ? Rotational velocity
      real*4        e_vsini_1  (nr__1) ! (km/s) ? rms uncertainty on vsini
      character*3   r_vsini_1  (nr__1) ! ? Source of vsini (8)
      real*4        Plx_1      (nr__1) ! (mas) ? Trigonometric parallax
      real*4        e_Plx_1    (nr__1) ! (mas) ? Standard error in Plx
      character*3   r_Plx_1    (nr__1) ! Source of Plx (8)
      real*4        Vmag_1     (nr__1) ! (mag) Johnson V magnitude
      character*1   f_Vmag_1   (nr__1) ! [V] Quality flag indicating variability (12)
      real*4        B_V_1      (nr__1) ! (mag) ? Johnson B-V colour index
      real*4        V_I_1      (nr__1) ! (mag) ? Cousins V-I colour index
      character*3   r_Vmag_1   (nr__1) ! Source of photometry (13)
      character*1   n_Vmag_1   (nr__1) ! [*] *: photometry corrected for companions (14)
      character*8   SpType_1   (nr__1) ! MK spectral type
      integer*4     EWLi_1     (nr__1) ! (0.1pm) ? LiI equivalent width in milli-Angstroms (15)
      character*1   f_EWLi_1   (nr__1) ! [:?] Quality flag for EWLi
      real*4        EWHa_1     (nr__1) ! (0.1nm) ? H{alpha} line emission in Angstroms
      character*5   r_EWLi_1   (nr__1) ! Source of spectral information (16)
*Note (1):
*    ROSAT All-Sky Survey Catalogue source name, in general from the Bright
*    Source Catalogue (<IX/10>) but sometimes refers to the Faint Source
*    Catalogue (<IX/29>). A repeated entry means more than one observed
*    possible optical counterpart. The same source may have entries in both
*    Tables 3 and 4, but the best optical candidate may be in any of the
*    tables.
*Note (2):
*    Order of preference of the identification is Hipparcos, TYCHO-2 or
*    UCAC2, independently of the astrometric data used.
*Note (3):
*    Order of preference is HD, DM, and other. GSC is only used if there is
*    no TYCHO-2 identification for the star.
*Note (4):
*    Order of preference of the sources of the proper motions is the same
*    as of the Name, adding other sources, but in some cases we used the
*    best value from our judgment.
*Note (5):
*    Radial velocities from SB2 observed once or SB1 are mean values, and
*    may not represent the systemic velocity. When there is a systemic
*    velocity from the literature, it is the used value.
*Note (6):
*    Null refers to no published number of observations in the literature
*    cited in r_HRV. Nobs from SACY includes references N=25 and 26.
*Note (7):
*    The * means that this is the number of spectral observations in the
*    SACY project, but the systemic HRV value was taken from the literature.
*    For V1005 Ori the non variable HRV is from literature as our spectrum
*    was not properly calibrated.
*Note (8):
*    If there is no number, the radial velocity is from this work.
*Note (9):
*    D means visual binary; sometimes the components were observed
*    separately and there is another entry, designated by A, B, C...,
*    in Table 3 or Table 4. If the bright component is also a SB the D
*    is given in the entry of the faint one.
*
*    A ? means possible evidence for radial velocity variations, sometimes
*    by comparing with published data. In some cases it is difficult to
*    decide between high rotation (+ spots?) and SB2. For other SB2? with
*    low rotation and large e_HRV the star should be considered at least
*    a SB1. For close visual doubles observed together, one of the
*    spectroscopic components may be also the visual one, mainly in
*    SB3 cases. For SBn all other data of this line are for the bright
*    component. In the future we will publish the data for spectroscopic
*    binaries found in the SACY project.
*Note (10):
*    Separation between the visual components, even if in n_HRV there is a
*    SBn. This means there is D for a faint companion in n_HRV in another
*    entry.
*Note (11):
*    Magnitude difference between faint companions and the bright one. A
*    negative value means that the faint companions together are brighter
*    than the primary, or the so called primary is a fainter companion with
*    sharper lines.
*
*    For visual binaries the negative value means that primary is also a
*    SBn, and the visual difference is given in this entry or the primary
*    is in the other table.
*Note (12):
*    V: photometric variability.
*Note (13): Source of photometry:
*       H = origin of the photometry from Hipparcos (<I/239>)
*       T = from TYCHO2 (<I/259>) converted to Johnson
*   blank = estimated magnitude
*Note (14):
*    for the proposes of the SACY the photometric measures for close visual
*    doubles and SBs are corrected to take into account the secondaries. As
*    these values are model dependent they must be used with caution.
*Note (15):
*    For SB2 or SB3 the values are tentatively corrected for the continuum
*    flux ratios.
*Note (16): For observations from this paper:
*       F = observation with ESO/FEROS
*       d = observation with the OPD coude spectrograph

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

C  Declarations for 'refs.dat'	! References (table5 of the paper)

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

      character*3   Ref        (nr__2) ! Reference number or letter
      character*19  BibCode    (nr__2) ! BibCode
      character*24  Aut        (nr__2) ! Author's name
      character*48  Com        (nr__2) ! Comments

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

C  Loading file 'table3.dat'	! Astrometric, spectroscopic and photometric data
*                              of the SACY sample

C  Format for file interpretation

    1 format(
     +  I4,2X,A16,2X,A16,1X,A13,1X,I2,1X,I2,1X,F4.1,2X,A1,I2,1X,I2,1X,
     +  I2,2X,F7.1,1X,F7.1,2X,A3,1X,F6.1,A1,1X,F4.1,1X,I3,A1,1X,A3,2X,
     +  A5,1X,F4.1,2X,F4.1,A1,2X,F5.1,1X,F4.1,2X,A3,2X,F5.1,2X,F4.1,
     +  2X,A3,2X,F5.2,A1,F5.2,1X,F5.2,2X,A3,A1,3X,A8,1X,I3,A1,1X,F5.1,
     +  2X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,1626
        read(1,'(A228)')ar__
        read(ar__,1)
     +  Seq(i__),v1RXS(i__),Name(i__),OName(i__),RAh(i__),RAm(i__),
     +  RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),pmRA(i__),
     +  pmDE(i__),r_pmRA(i__),HRV(i__),f_HRV(i__),e_HRV(i__),
     +  Nobs(i__),f_Nobs(i__),r_HRV(i__),n_HRV(i__),Sep(i__),dV(i__),
     +  u_dV(i__),vsini(i__),e_vsini(i__),r_vsini(i__),Plx(i__),
     +  e_Plx(i__),r_Plx(i__),Vmag(i__),f_Vmag(i__),B_V(i__),V_I(i__),
     +  r_Vmag(i__),n_Vmag(i__),SpType(i__),EWLi(i__),f_EWLi(i__),
     +  EWHa(i__),r_EWLi(i__)
        if(ar__(108:111) .EQ. '') e_HRV(i__) = rNULL__
        if(ar__(113:115) .EQ. '') Nobs(i__) = iNULL__
        if(ar__(129:132) .EQ. '') Sep(i__) = rNULL__
        if(ar__(135:138) .EQ. '') dV(i__) = rNULL__
        if(ar__(142:146) .EQ. '') vsini(i__) = rNULL__
        if(ar__(148:151) .EQ. '') e_vsini(i__) = rNULL__
        if(ar__(159:163) .EQ. '') Plx(i__) = rNULL__
        if(ar__(166:169) .EQ. '') e_Plx(i__) = rNULL__
        if(ar__(183:187) .EQ. '') B_V(i__) = rNULL__
        if(ar__(189:193) .EQ. '') V_I(i__) = rNULL__
        if(ar__(212:214) .EQ. '') EWLi(i__) = iNULL__
        if(ar__(217:221) .EQ. '') EWHa(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)
     +  Seq(i__),v1RXS(i__),Name(i__),OName(i__),RAh(i__),RAm(i__),
     +  RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),pmRA(i__),
     +  pmDE(i__),r_pmRA(i__),HRV(i__),f_HRV(i__),e_HRV(i__),
     +  Nobs(i__),f_Nobs(i__),r_HRV(i__),n_HRV(i__),Sep(i__),dV(i__),
     +  u_dV(i__),vsini(i__),e_vsini(i__),r_vsini(i__),Plx(i__),
     +  e_Plx(i__),r_Plx(i__),Vmag(i__),f_Vmag(i__),B_V(i__),V_I(i__),
     +  r_Vmag(i__),n_Vmag(i__),SpType(i__),EWLi(i__),f_EWLi(i__),
     +  EWHa(i__),r_EWLi(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 'table4.dat'	! Astrometric, spectroscopic and photometric data
*                              of other observed stars

C  Format for file interpretation

    2 format(
     +  I4,2X,A16,2X,A16,1X,A13,1X,I2,1X,I2,1X,F4.1,2X,A1,I2,1X,I2,1X,
     +  I2,2X,F7.1,1X,F7.1,2X,A3,1X,F6.1,A1,1X,F4.1,1X,I3,A1,1X,A3,2X,
     +  A5,1X,F4.1,2X,F4.1,A1,2X,F5.1,1X,F4.1,2X,A3,2X,F5.1,2X,F4.1,
     +  2X,A3,2X,F5.2,A1,F5.2,1X,F5.2,2X,A3,A1,3X,A8,1X,I3,A1,1X,F5.1,
     +  2X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,165
        read(1,'(A228)')ar__1
        read(ar__1,2)
     +  Seq_1(i__),v1RXS_1(i__),Name_1(i__),OName_1(i__),RAh_1(i__),
     +  RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),pmRA_1(i__),pmDE_1(i__),r_pmRA_1(i__),HRV_1(i__),
     +  f_HRV_1(i__),e_HRV_1(i__),Nobs_1(i__),f_Nobs_1(i__),
     +  r_HRV_1(i__),n_HRV_1(i__),Sep_1(i__),dV_1(i__),u_dV_1(i__),
     +  vsini_1(i__),e_vsini_1(i__),r_vsini_1(i__),Plx_1(i__),
     +  e_Plx_1(i__),r_Plx_1(i__),Vmag_1(i__),f_Vmag_1(i__),
     +  B_V_1(i__),V_I_1(i__),r_Vmag_1(i__),n_Vmag_1(i__),
     +  SpType_1(i__),EWLi_1(i__),f_EWLi_1(i__),EWHa_1(i__),
     +  r_EWLi_1(i__)
        if(ar__1(108:111) .EQ. '') e_HRV_1(i__) = rNULL__
        if(ar__1(113:115) .EQ. '') Nobs_1(i__) = iNULL__
        if(ar__1(129:132) .EQ. '') Sep_1(i__) = rNULL__
        if(ar__1(135:138) .EQ. '') dV_1(i__) = rNULL__
        if(ar__1(142:146) .EQ. '') vsini_1(i__) = rNULL__
        if(ar__1(148:151) .EQ. '') e_vsini_1(i__) = rNULL__
        if(ar__1(159:163) .EQ. '') Plx_1(i__) = rNULL__
        if(ar__1(166:169) .EQ. '') e_Plx_1(i__) = rNULL__
        if(ar__1(183:187) .EQ. '') B_V_1(i__) = rNULL__
        if(ar__1(189:193) .EQ. '') V_I_1(i__) = rNULL__
        if(ar__1(212:214) .EQ. '') EWLi_1(i__) = iNULL__
        if(ar__1(217:221) .EQ. '') EWHa_1(i__) = rNULL__
        RAdeg_1(i__) = rNULL__
        DEdeg_1(i__) = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_1 from input data
c  (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
        if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
        if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
        if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
        if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
        if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
        if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
        if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Seq_1(i__),v1RXS_1(i__),Name_1(i__),OName_1(i__),RAh_1(i__),
     +  RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
     +  DEs_1(i__),pmRA_1(i__),pmDE_1(i__),r_pmRA_1(i__),HRV_1(i__),
     +  f_HRV_1(i__),e_HRV_1(i__),Nobs_1(i__),f_Nobs_1(i__),
     +  r_HRV_1(i__),n_HRV_1(i__),Sep_1(i__),dV_1(i__),u_dV_1(i__),
     +  vsini_1(i__),e_vsini_1(i__),r_vsini_1(i__),Plx_1(i__),
     +  e_Plx_1(i__),r_Plx_1(i__),Vmag_1(i__),f_Vmag_1(i__),
     +  B_V_1(i__),V_I_1(i__),r_Vmag_1(i__),n_Vmag_1(i__),
     +  SpType_1(i__),EWLi_1(i__),f_EWLi_1(i__),EWHa_1(i__),
     +  r_EWLi_1(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References (table5 of the paper)

C  Format for file interpretation

    3 format(A3,1X,A19,1X,A24,A48)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,107
        read(1,'(A96)')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