FORTRAN Generation
(/./ftp/cats/J/AJ/153/14)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/153/14 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/AJ/153/14  Solar neighborhood XXXVIII. Nearby M dwarf systems (Winters+, 2017)
*================================================================================
*The solar neighborhood XXXVIII. Results from the CTIO/SMARTS 0.9M: trigonometric
*parallaxes for 151 nearby M dwarf systems.
*    Winters J.G., Sevrinsky R.A., Jao W.-C., Henry T.J., Riedel A.R.,
*    Subasavage J.P., Lurie J.C., Ianna P.A., Finch C.T.
*   <Astron. J., 153, 14-14 (2017)>
*   =2017AJ....153...14W    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Photometric data

      integer*4 nr__
      parameter (nr__=160)	! Number of records
      character*99 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*13  Name       (nr__) ! Name of the M dwarf (1)
      character*3   m_Name     (nr__) ! Components of multiple system (noted with
*                                   capital letters A and B after the name) (G1)
      integer*4     RAh        (nr__) ! (h) Hour of Right Ascension (J2000)
      integer*4     RAm        (nr__) ! (min) Minute of Right Ascension (J2000)
      real*4        RAs        (nr__) ! (s) Second of Right Ascension (J2000)
      character*1   DE_        (nr__) ! Sign of the Declination (J2000)
      integer*4     DEd        (nr__) ! (deg) Degree of Declination (J2000)
      integer*4     DEm        (nr__) ! (arcmin) Arcminute of Declination (J2000)
      integer*4     DEs        (nr__) ! (arcsec) Arcsecond of Declination (J2000)
      real*4        Vmag       (nr__) ! (mag) [9.8/20] Johnson V-band magnitude (2)
      character*1   f_Vmag     (nr__) ! [Jd] Flag on Vmag (J or d) (3)
      real*4        Rmag       (nr__) ! (mag) [8.9/19.2] Kron-Cousins R-band magnitude (2)
      character*1   f_Rmag     (nr__) ! [J] Flag 'J' on Rmag (3)
      real*4        Imag       (nr__) ! (mag) [7.93/18.26] Kron-Cousins I-band magnitude (2)
      character*1   f_Imag     (nr__) ! [J] Flag 'J' on Imag (3)
      integer*4     Nobs       (nr__) ! [1/5]? Number of observations (1-5)
      integer*4     Ref        (nr__) ! [1/3] Reference (1-3) (4)
      character*1   f_Ref      (nr__) ! [u] Flag 'u' on Ref for previously reported,
*                                   but updated, photometry (5)
      real*4        Jmag       (nr__) ! (mag) [6.84/12.02]? 2MASS J band magnitude
      character*1   f_Jmag     (nr__) ! [Jc] Flag on Jmag (J or c) (3)
      real*4        Hmag       (nr__) ! (mag) [6.21/11.45]? 2MASS H band magnitude
      character*3   f_Hmag     (nr__) ! Flag on Hmag (J and/or c) (3)
      real*4        Ksmag      (nr__) ! (mag) [5.96/11.24]? 2MASS Ks band magnitude
      character*1   f_Ksmag    (nr__) ! [Jc] Flag on Ksmag (J or c) (3)
      real*4        Dist       (nr__) ! (pc) [6.4/67.9]? Photometric distance estimate
*                                   (d_CCD_) (6)
      character*3   f_Dist     (nr__) ! Flag on Dist (a and/or b) (7)
      real*4        e_Dist     (nr__) ! (pc) [1/14]? Distance uncertainty based on the
*                                   photometry ({sigma}_CCD_)
*Note (1): This table provides the VRIJHK photometry data for the 160 stars in
*     151 M dwarf systems.
*Note (2): As outlined in Winters et al. 2011 (Cat. J/AJ/141/21), photometric
*     uncertainties are typically 30mmag for the V-band and 20mmag for the R-
*     and I-bands.
*Note (3): The flags on photometry are defined as follows:
*     c = 2MASS magnitude uncertainty greater than 50mmag.
*     d = VRI photometry for individual components has been calculated using
*         Point-Spread Function (PSF) reduction;
*     J = The magnitude is blended due to one or more close companions.
*Note (4): The reference codes are defined as follows:
*     1 = This work;
*     2 = Winters et al. 2011 (Cat. J/AJ/141/21);
*     3 = Winters et al. 2015 (Cat. J/AJ/149/5).
*Note (5): A "u" next to the reference indicates that REsearch Consortium On
*     Nearby Stars (RECONS) photometry was reported in a previous publication,
*     but that the value listed here is updated due to additional data and
*     supersedes previous results.
*Note (6): Photometric distance estimates are calculated using the relations
*     described in Henry et al. (2004AJ....128.2460H). Please refer to Section
*     4.1 in the paper for further details.
*Note (7): The flags on the distance are defined as follows:
*     a = The {pi}_trig_ for this object places it beyond 25pc. Nearly all of
*         these stars were originally thought to be within 25pc, but our new
*         astrometry indicates that they are not. This in itself reveals
*         important information about the stars, as discussed in Section 4.1;
*     b = The photometric distance estimate is based upon the joint photometry
*         for the system and is thus likely an underestimate.

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

C  Declarations for 'table2.dat'	! Astrometry data

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

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

      character*13  Name_1     (nr__1) ! Name of the M dwarf (1)
      character*3   m_Name_1   (nr__1) ! Components of multiple system (noted with
*                                 capital letters A and B after the name) (G1)
      integer*4     RAh_1      (nr__1) ! (h) Hour of Right Ascension (J2000)
      integer*4     RAm_1      (nr__1) ! (min) Minute of Right Ascension (J2000)
      real*4        RAs_1      (nr__1) ! (s) Second of Right Ascension (J2000)
      character*1   DE__1      (nr__1) ! Sign of the Declination (J2000)
      integer*4     DEd_1      (nr__1) ! (deg) Degree of Declination (J2000)
      integer*4     DEm_1      (nr__1) ! (arcmin) Arcminute of Declination (J2000)
      real*4        DEs_1      (nr__1) ! (arcsec) Arcsecond of Declination (J2000)
      character*1   Flt        (nr__1) ! [VRI] Filter in which {pi}_trig_ was measured
*                                 (V, R, or I)
      integer*4     Nsea       (nr__1) ! [3/15] Number of seasons covered (N_sea_) (2)
      character*1   f_Nsea     (nr__1) ! [cs] Flag on Nsea (c or s) (3)
      integer*4     Nfr        (nr__1) ! [33/200] Number of frames acquired for each star
*                                 (N_fr_)
      real*8        Date0      (nr__1) ! (yr) [1999.9/2013.4] Starting date of the time
*                                 interval covered by the data sequence
      real*8        Date1      (nr__1) ! (yr) [2002.46/2016.21] Ending date of the time
*                                 interval covered by the data sequence
      real*4        Time       (nr__1) ! (yr) [2.32/15.61] Number of years of coverage (T)
      integer*4     Nref       (nr__1) ! [5/16] Number of reference stars used in the
*                                 field
      real*4        Var        (nr__1) ! (mmag) [6.4/34]? Result for photometric variability in
*                                 the filter used for parallax observation
      real*4        plxrel     (nr__1) ! (mas) [15.74/97.82] Relative parallax ({pi}_rel_)
      real*4        e_plxrel   (nr__1) ! (mas) [0.55/3.5] Uncertainty in relative parallax
      real*4        pcorr      (nr__1) ! (mas) [0.14/3.79] Correction to absolute parallax
*                                 ({pi}_cor_) (4)
      real*4        e_pcorr    (nr__1) ! (mas) [0.02/0.79] Uncertainty in pcorr
      character*1   f_pcorr    (nr__1) ! [b] Flag 'b' on pcorr when a generic
*                                 correction of 1.50+/-0.50mas was adopted (5)
      real*4        plxabs     (nr__1) ! (mas) [16.43/99.16] Final absolute parallax {pi}_abs_
      real*4        e_plxabs   (nr__1) ! (mas) [0.59/3.51] Uncertainty in plxabs
      character*1   f_plxabs   (nr__1) ! [a] Flag 'a' on plxabs when the parallax value
*                                 have had the orbital motion of the system
*                                 subtracted
      real*4        pm         (nr__1) ! (mas/yr) [118.5/827.7] Magnitude of the proper motion
      real*4        e_pm       (nr__1) ! (mas/yr) [0.1/2.9] Uncertainty in pm
      real*4        pmPA       (nr__1) ! (deg) [5.8/354] Position angle (measured east from
*                                 north) of the proper motion
      real*4        e_pmPA     (nr__1) ! (deg) [0.02/1.72] Uncertainty in pmPA
      real*4        Vtan       (nr__1) ! (km/s) [11.2/111.4] Derived tangential velocity V_tan_
      character*1   Note       (nr__1) ! [N] Note (indicated with a 'N') from Section 4.5
*Note (1): The object name "EHPM 2-801" was probably a mistake and was modified
*     into "SCR 2105-5503" at CDS.
*Note (2): Where the six months when a star is observable counts as one season.
*Note (3): The observational flags are defined as follows:
*     c = Continuous set of observations where multiple nights of data were taken
*         in each season;
*     s = Scattered observations when some seasons have only a single night of
*         observations.
*Note (4): Needed because the reference stars are not at an infinite distance.
*Note (5): Because this field is in the direction of the Galactic center, the
*     initial correction from the relative parallax to the absolute parallax for
*     this object was erroneously large due to reddened reference stars. We have
*     adopted a generic correction of 1.50+/-0.50mas.

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

C  Declarations for 'table4.dat'	! Information for confirmed multiples

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

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

      character*13  Name_2     (nr__2) ! Name of the M dwarf
      character*2   m_Name_2   (nr__2) ! [AB] Components of multiple system (noted with
*                                   capital letters A and B after the name) (G1)
      integer*4     RAh_2      (nr__2) ! (h) Hour of Right Ascension (J2000)
      integer*4     RAm_2      (nr__2) ! (min) Minute of Right Ascension (J2000)
      real*4        RAs_2      (nr__2) ! (s) Second of Right Ascension (J2000)
      character*1   DE__2      (nr__2) ! Sign of the Declination (J2000)
      integer*4     DEd_2      (nr__2) ! (deg) Degree of Declination (J2000)
      integer*4     DEm_2      (nr__2) ! (arcmin) Arcminute of Declination (J2000)
      integer*4     DEs_2      (nr__2) ! (arcsec) Arcsecond of Declination (J2000)
      character*1   l_Sep      (nr__2) ! [<] Upper limit flag on Sep
      real*4        Sep        (nr__2) ! (arcsec) [0.02/49.58] Separation {rho} (1)
      integer*4     PA         (nr__2) ! (deg) [32/354]? Position angle {theta} (1)
      integer*4     Year       (nr__2) ! (yr) [2001/2016]? Date of the measurement
      character*6   Tech       (nr__2) ! Technique used to identify the binary (2)
      integer*4     r_Tech     (nr__2) ! [1/5] Reference for the separation measurement
*                                   and the detection technique (3)
      real*4        Dmag       (nr__2) ! (mag) [0.1/6]? Magnitude difference between the
*                                   components, where available ({Delta}mag)
      character*5   Flt_1      (nr__2) ! Filter in which the measurement was made
*                                   (either 754nm, H, I_KC_, or V_J_) (4)
      integer*4     r_Dmag     (nr__2) ! [1/5]? Reference for the {Delta}mag (3)
*Note (1): Of the companion from the primary.
*Note (2): Codes for the techniques used to identify companions are as follows:
*     astdet = Astrometric detection of a perturbation of the photocenter, as
*              discussed in Section 5.1;
*     astorb = Astrometric orbit;
*        SB2 = Double-lined spectroscopic binary;
*     spkdet = Detection using speckle interferometry;
*     visdet = Visual detection.
*Note (3): The reference codes are defined as follows:
*     1 = This work;
*     2 = Bonfils et al. 2013 (Cat. J/A+A/549/A109);
*     3 = Forveille et al. 2005A&A...435L...5F;
*     4 = Horch et al. 2010 (Cat. J/AJ/139/205);
*     5 = Jao et al. 2003AJ....125..332J;
*     6 = Winters et al. 2011 (Cat. J/AJ/141/21).
*Note (4): The filters in which the {Delta}mags were measured are either V or I
*     on the Johnson-Kron-Cousins system, the filter at 754nm used on the DSSI
*     speckle instrument, or the 2MASS H-band.

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

C  Declarations for 'table5.dat'	! Candidate multiples and young stars

      integer*4 nr__3
      parameter (nr__3=12)	! Number of records
      character*53 ar__3  	! Full-size record

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

      character*13  Name_3     (nr__3) ! Name of the M dwarf
      integer*4     RAh_3      (nr__3) ! (h) Hour of Right Ascension (J2000)
      integer*4     RAm_3      (nr__3) ! (min) Minute of Right Ascension (J2000)
      real*4        RAs_3      (nr__3) ! (s) Second of Right Ascension (J2000)
      character*1   DE__3      (nr__3) ! Sign of the Declination (J2000)
      integer*4     DEd_3      (nr__3) ! (deg) Degree of Declination (J2000)
      integer*4     DEm_3      (nr__3) ! (arcmin) Arcminute of Declination (J2000)
      integer*4     DEs_3      (nr__3) ! (arcsec) Arcsecond of Declination (J2000)
      real*4        Vtan_1     (nr__3) ! (km/s) [22.6/64.8] Tangential velocity
      character*12  Youth      (nr__3) ! Youth indicator (1)
*Note (1): The youth indicators "GALEX" (Galaxy Evolution Explorer DR5 UV
*     catalog; Bianchi et al. 2011, Cat. II/312) and "ROSAT" (ROentgen SATellite
*     All Sky Survey X-ray catalog; Voges et al. 1999, Cat. IX/10; Voges et al.
*     2000, Cat. IX/29) indicate detections in each respective catalog, while
*     "Tuc-Hor," "AB Dor," "Carina-Near," and "Argus" indicate the young nearby
*     moving groups to which the star in question has a probability of belonging,
*     based on its kinematic information.

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

C  Declarations for 'notes.dat'	! Notes

      integer*4 nr__4
      parameter (nr__4=29)	! Number of records
      character*605 ar__4  	! Full-size record

      character*13  Name_4     (nr__4) ! Name of the M dwarf
      character*2   m_Name_3   (nr__4) ! Components of multiple system
      character*9   OName      (nr__4) ! Other name (first four digits of both the R.A.
*                                   and decl.)
      character*579 Note_1     (nr__4) ! Notes (figure numbers are listed if plots of
*                                   the astrometry series are included in the
*                                   paper)

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

C  Loading file 'table1.dat'	! Photometric data

C  Format for file interpretation

    1 format(
     +  A13,A3,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F5.2,A1,1X,
     +  F5.2,A1,1X,F5.2,A1,1X,I1,1X,I1,A1,1X,F5.2,A1,1X,F5.2,A3,1X,
     +  F5.2,A1,1X,F4.1,A3,1X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,160
        read(1,'(A99)')ar__
        read(ar__,1)
     +  Name(i__),m_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Vmag(i__),f_Vmag(i__),Rmag(i__),
     +  f_Rmag(i__),Imag(i__),f_Imag(i__),Nobs(i__),Ref(i__),
     +  f_Ref(i__),Jmag(i__),f_Jmag(i__),Hmag(i__),f_Hmag(i__),
     +  Ksmag(i__),f_Ksmag(i__),Dist(i__),f_Dist(i__),e_Dist(i__)
        if(ar__(60:60) .EQ. '') Nobs(i__) = iNULL__
        if(ar__(65:69) .EQ. '') Jmag(i__) = rNULL__
        if(ar__(72:76) .EQ. '') Hmag(i__) = rNULL__
        if(ar__(81:85) .EQ. '') Ksmag(i__) = rNULL__
        if(ar__(88:91) .EQ. '') Dist(i__) = rNULL__
        if(ar__(96:99) .EQ. '') e_Dist(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)
     +  Name(i__),m_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
     +  DEd(i__),DEm(i__),DEs(i__),Vmag(i__),f_Vmag(i__),Rmag(i__),
     +  f_Rmag(i__),Imag(i__),f_Imag(i__),Nobs(i__),Ref(i__),
     +  f_Ref(i__),Jmag(i__),f_Jmag(i__),Hmag(i__),f_Hmag(i__),
     +  Ksmag(i__),f_Ksmag(i__),Dist(i__),f_Dist(i__),e_Dist(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 'table2.dat'	! Astrometry data

C  Format for file interpretation

    2 format(
     +  A13,A3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A1,1X,I2,
     +  A1,1X,I3,1X,F7.2,1X,F7.2,1X,F5.2,1X,I2,1X,F4.1,1X,F5.2,1X,
     +  F4.2,1X,F4.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F5.1,1X,
     +  F3.1,1X,F5.1,1X,F4.2,1X,F5.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,160
        read(1,'(A145)')ar__1
        read(ar__1,2)
     +  Name_1(i__),m_Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Flt(i__),
     +  Nsea(i__),f_Nsea(i__),Nfr(i__),Date0(i__),Date1(i__),
     +  Time(i__),Nref(i__),Var(i__),plxrel(i__),e_plxrel(i__),
     +  pcorr(i__),e_pcorr(i__),f_pcorr(i__),plxabs(i__),
     +  e_plxabs(i__),f_plxabs(i__),pm(i__),e_pm(i__),pmPA(i__),
     +  e_pmPA(i__),Vtan(i__),Note(i__)
        if(ar__1(77:80) .EQ. '') Var(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)
     +  Name_1(i__),m_Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
     +  DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Flt(i__),
     +  Nsea(i__),f_Nsea(i__),Nfr(i__),Date0(i__),Date1(i__),
     +  Time(i__),Nref(i__),Var(i__),plxrel(i__),e_plxrel(i__),
     +  pcorr(i__),e_pcorr(i__),f_pcorr(i__),plxabs(i__),
     +  e_plxabs(i__),f_plxabs(i__),pm(i__),e_pm(i__),pmPA(i__),
     +  e_pmPA(i__),Vtan(i__),Note(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 'table4.dat'	! Information for confirmed multiples

C  Format for file interpretation

    3 format(
     +  A13,A2,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A1,F6.3,1X,
     +  I3,1X,I4,1X,A6,1X,I1,1X,F4.2,1X,A5,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,32
        read(1,'(A75)')ar__2
        read(ar__2,3)
     +  Name_2(i__),m_Name_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),l_Sep(i__),
     +  Sep(i__),PA(i__),Year(i__),Tech(i__),r_Tech(i__),Dmag(i__),
     +  Flt_1(i__),r_Dmag(i__)
        if(ar__2(46:48) .EQ. '') PA(i__) = iNULL__
        if(ar__2(50:53) .EQ. '') Year(i__) = iNULL__
        if(ar__2(64:67) .EQ. '') Dmag(i__) = rNULL__
        if(ar__2(75:75) .EQ. '') r_Dmag(i__) = iNULL__
        RAdeg_2(i__) = rNULL__
        DEdeg_2(i__) = rNULL__
c  Derive coordinates RAdeg_2 and DEdeg_2 from input data
c  (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
        if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
        if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
        if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
        if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
        if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
        if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
        if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  Name_2(i__),m_Name_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),
     +  DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),l_Sep(i__),
     +  Sep(i__),PA(i__),Year(i__),Tech(i__),r_Tech(i__),Dmag(i__),
     +  Flt_1(i__),r_Dmag(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Candidate multiples and young stars

C  Format for file interpretation

    4 format(
     +  A13,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F4.1,1X,A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,12
        read(1,'(A53)')ar__3
        read(ar__3,4)
     +  Name_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),
     +  DEd_3(i__),DEm_3(i__),DEs_3(i__),Vtan_1(i__),Youth(i__)
        RAdeg_3(i__) = rNULL__
        DEdeg_3(i__) = rNULL__
c  Derive coordinates RAdeg_3 and DEdeg_3 from input data
c  (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
        if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
        if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
        if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
        if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
        if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
        if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
        if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  Name_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),DE__3(i__),
     +  DEd_3(i__),DEm_3(i__),DEs_3(i__),Vtan_1(i__),Youth(i__)
        write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'notes.dat'	! Notes

C  Format for file interpretation

    5 format(A13,A2,1X,A9,1X,A579)

C  Effective file loading

      open(unit=1,status='old',file=
     +'notes.dat')
      write(6,*) '....Loading file: notes.dat'
      do i__=1,29
        read(1,'(A605)')ar__4
        read(ar__4,5)Name_4(i__),m_Name_3(i__),OName(i__),Note_1(i__)
c    ..............Just test output...........
        write(6,5)Name_4(i__),m_Name_3(i__),OName(i__),Note_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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