FORTRAN Generation
(/./ftp/cats/I/239)

Conversion of standardized ReadMe file for file /./ftp/cats/I/239 into FORTRAN code for reading data files line by line.

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 2018-Jul-21
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. I/239           The Hipparcos and Tycho Catalogues                    (ESA 1997)
*================================================================================
*The Hipparcos and Tycho Catalogues
*    ESA 1997
*   <ESA, 1997, The Hipparcos Catalogue, ESA SP-1200>
*   <ESA, 1997, The Tycho Catalogue, ESA SP-1200>
*   =1997HIP...C......0E
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'hip_main.dat'	! The Hipparcos Main Catalogue

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

C  Position composed of: RAdeg DEdeg (Epoch=J1991.25)
      character*1   Catalog     ! [H] Catalogue (H=Hipparcos)               (H0)
      integer*4     HIP         ! Identifier (HIP number)                   (H1)
      character*1   Proxy       ! *[HT] Proximity flag                       (H2)
      character*11  RAhms       ! Right ascension in h m s, ICRS (J1991.25) (H3)
      character*11  DEdms       ! Declination in deg ' ", ICRS (J1991.25)   (H4)
      real*4        Vmag        ! (mag) ? Magnitude in Johnson V                  (H5)
      integer*4     VarFlag     ! *[1,3]? Coarse variability flag            (H6)
      character*1   r_Vmag      ! *[GHT] Source of magnitude                 (H7)
      real*8        RAdeg       ! (deg) *? alpha, degrees (ICRS, Epoch=J1991.25)   (H8)
      real*8        DEdeg       ! (deg) *? delta, degrees (ICRS, Epoch=J1991.25)   (H9)
      character*1   AstroRef    ! *[*+A-Z] Reference flag for astrometry    (H10)
      real*8        Plx         ! (mas) ? Trigonometric parallax                 (H11)
      real*8        pmRA        ! (mas/yr) *? Proper motion mu_alpha.cos(delta), ICRS(H12)
      real*8        pmDE        ! (mas/yr) *? Proper motion mu_delta, ICRS           (H13)
      real*4        e_RAdeg     ! (mas) *? Standard error in RA*cos(DEdeg)        (H14)
      real*4        e_DEdeg     ! (mas) *? Standard error in DE                   (H15)
      real*4        e_Plx       ! (mas) ? Standard error in Plx                  (H16)
      real*4        e_pmRA      ! (mas/yr) ? Standard error in pmRA                 (H17)
      real*4        e_pmDE      ! (mas/yr) ? Standard error in pmDE                 (H18)
      real*4        DE_RA       ! [-1/1]? Correlation, DE/RA*cos(delta)    (H19)
      real*4        Plx_RA      ! [-1/1]? Correlation, Plx/RA*cos(delta)   (H20)
      real*4        Plx_DE      ! [-1/1]? Correlation, Plx/DE              (H21)
      real*4        pmRA_RA     ! [-1/1]? Correlation, pmRA/RA*cos(delta)  (H22)
      real*4        pmRA_DE     ! [-1/1]? Correlation, pmRA/DE             (H23)
      real*4        pmRA_Plx    ! [-1/1]? Correlation, pmRA/Plx            (H24)
      real*4        pmDE_RA     ! [-1/1]? Correlation, pmDE/RA*cos(delta)  (H25)
      real*4        pmDE_DE     ! [-1/1]? Correlation, pmDE/DE             (H26)
      real*4        pmDE_Plx    ! [-1/1]? Correlation, pmDE/Plx            (H27)
      real*4        pmDE_pmRA   ! [-1/1]? Correlation, pmDE/pmRA           (H28)
      integer*4     F1          ! (%) ? Percentage of rejected data            (H29)
      real*4        F2          ! *? Goodness-of-fit parameter              (H30)
      real*4        BTmag       ! (mag) ? Mean BT magnitude                      (H32)
      real*4        e_BTmag     ! (mag) ? Standard error on BTmag                (H33)
      real*4        VTmag       ! (mag) ? Mean VT magnitude                      (H34)
      real*4        e_VTmag     ! (mag) ? Standard error on VTmag                (H35)
      character*1   m_BTmag     ! *[A-Z*-] Reference flag for BT and VTmag  (H36)
      real*4        B_V         ! (mag) ? Johnson B-V colour                     (H37)
      real*4        e_B_V       ! (mag) ? Standard error on B-V                  (H38)
      character*1   r_B_V       ! [GT] Source of B-V from Ground or Tycho  (H39)
      real*4        V_I         ! (mag) ? Colour index in Cousins' system        (H40)
      real*4        e_V_I       ! (mag) ? Standard error on V-I                  (H41)
      character*1   r_V_I       ! *[A-T] Source of V-I                      (H42)
      character*1   CombMag     ! [*] Flag for combined Vmag, B-V, V-I     (H43)
      real*8        Hpmag       ! (mag) *? Median magnitude in Hipparcos system   (H44)
      real*4        e_Hpmag     ! (mag) *? Standard error on Hpmag                (H45)
      real*4        Hpscat      ! (mag) ? Scatter on Hpmag                       (H46)
      integer*4     o_Hpmag     ! ? Number of observations for Hpmag       (H47)
      character*1   m_Hpmag     ! *[A-Z*-] Reference flag for Hpmag         (H48)
      real*4        Hpmax       ! (mag) ? Hpmag at maximum (5th percentile)      (H49)
      real*4        HPmin       ! (mag) ? Hpmag at minimum (95th percentile)     (H50)
      real*8        Period      ! (d) ? Variability period (days)              (H51)
      character*1   HvarType    ! *[CDMPRU]? variability type               (H52)
      character*1   moreVar     ! *[12] Additional data about variability   (H53)
      character*1   morePhoto   ! [ABC] Light curve Annex                  (H54)
      character*10  CCDM        ! CCDM identifier                          (H55)
      character*1   n_CCDM      ! *[HIM] Historical status flag             (H56)
      integer*4     Nsys        ! ? Number of entries with same CCDM       (H57)
      integer*4     Ncomp       ! ? Number of components in this entry     (H58)
      character*1   MultFlag    ! *[CGOVX] Double/Multiple Systems flag     (H59)
      character*1   Source      ! *[PFILS] Astrometric source flag          (H60)
      character*1   Qual        ! *[ABCDS] Solution quality                 (H61)
      character*2   m_HIP       ! Component identifiers                    (H62)
      integer*4     theta       ! (deg) ? Position angle between components      (H63)
      real*8        rho         ! (arcsec) ? Angular separation between components  (H64)
      real*4        e_rho       ! (arcsec) ? Standard error on rho                  (H65)
      real*4        dHp         ! (mag) ? Magnitude difference of components     (H66)
      real*4        e_dHp       ! (mag) ? Standard error on dHp                  (H67)
      character*1   Survey      ! [S] Flag indicating a Survey Star        (H68)
      character*1   Chart       ! *[DG] Identification Chart                (H69)
      character*1   Notes       ! *[DGPWXYZ] Existence of notes             (H70)
      integer*4     HD          ! [1/359083]? HD number <III/135>          (H71)
      character*10  BD          ! Bonner DM <I/119>, <I/122>               (H72)
      character*10  CoD         ! Cordoba Durchmusterung (DM) <I/114>      (H73)
      character*10  CPD         ! Cape Photographic DM <I/108>             (H74)
      real*4        v_V_I_red   ! (mag) V-I used for reductions                  (H75)
      character*12  SpType      ! Spectral type                            (H76)
      character*1   r_SpType    ! *[1234GKSX]? Source of spectral type      (H77)
*Note on Proxy: this flag provides a coarse indication of the presence
*     of nearby objects within 10arcsec of the given entry.
*     If non-blank, it indicates that 
*     'H' there is one or more distinct Hipparcos Catalogue entries, 
*         or distinct components of system from h_dm_com.dat
*     'T' there is one or more distinct Tycho entries
*     If 'H' and 'T' apply, 'H' is adopted.
*     The 'T' flag implies either an inconsistency between the Hipparcos
*     and Tycho catalogues, or a deficiency in one or both of the 
*     catalogues.
*Note on RAdeg, DEdeg: right ascension and declination are
*     expressed in degrees for epoch J1991.25 (JD2448349.0625 (TT)) in the
*     ICRS (International Celestial Reference System, consistent with
*     J2000) reference system.
*     There are 263 cases where these fields are missing (no astrometric
*     solution could be found)
*Note on pmRA, pmDE:
*     The proper motions refer to the ICRS and to the epoch J1991.25.
*Note on e_RAdeg, e_DEdeg:
*     The standard errors refer to the epoch J1991.25, and represent a
*     minimum of the error on the position. The actual standard error
*     on the positions is increasing for epochs increasingly differing
*     from the nominal J1991.25 epoch.
*Note on VarFlag: the values are
*     1: < 0.06mag ; 2: 0.06-0.6mag ; 3: >0.6mag
*Note on r_Vmag: the source is
*     G = ground-based, H=HIP, T=Tycho
*Note on AstroRef: this flag indicates that the astrometric parameters in H3-4
*     and H8-30 refer to:
*     A to Z: the letter indicates the component of a double or multiple system
*     *: the photocentre of a double or multiple system
*     +: the centre of mass
*Note on F2: values exceeding +3 indicate a bad fit to the data.
*Note on m_BTmag: this flag indicates the component or combined photometry:
*     A to Z : the letter indicates the component measured in Tycho
*              (non-single star)
*     * : the photometry refers to all components of the Hipparcos entry
*     - : single-pointing triple or quadruple system
*Note on r_V-I: the origin of the V-I colour, in summary:
*     'A'        for an observation of V-I in Cousins' system;
*     'B' to 'K' when V-I derived from measurements in other
*                bands/photoelectric systems
*     'L' to 'P' when V-I derived from Hipparcos and Star Mapper photometry
*     'Q'        for long-period variables
*     'R' to 'T' when colours are unknown
*Note on Hpmag, e_Hpmag:
*     the Hipparcos magnitude could not be determined for 14 stars.
*Note on m_Hpmag: this flag indicates for double or multiple entries:
*     A to Z : the letter indicates the specified component measured
*     * : combined Hpmag of a double system, corrected for attenuation
*     - : combined Hpmag of a multiple system, not corrected for attenuation
*Note on HvarType: Hipparcos-defined type of variability (a blank entry
*     signifies that the entry could not be classified as variable or constant):
*     C : no variability detected ("constant")
*     D : duplicity-induced variability
*     M : possibly micro-variable (amplitude < 0.03mag)
*     P : periodic variable
*     R : V-I colour index was revised due to variability analysis
*     U : unsolved variable which does not fall in the other categories
*Note on moreVar: more data about periodic variability are provided
*Note on n_CCDM: the flag takes the following values:
*     H : determined multiple by Hipparcos, previously unknown
*     I : system previously identified as multiple in HIC <I/196> (annex1)
*     M : miscellaneous (system identified after publication of HIC)
*Note on MultFlag: indicates that further details are given in the Double
*     and Multiple Systems Annex:
*     C : solutions for the components
*     G : acceleration or higher order terms
*     O : orbital solutions
*     V : variability-induced movers (apparent motion arises from variability)
*     X : stochastic solution (probably astrometric binaries with short period)
*Note on Source: qualifies the source of the astrometric parameters H8-30
*        with a 'C' in MultFlag:
*     P : primary target of a 2- or 3-pointing system
*     F : secondary or tertiary of a 2- or 3-pointing 'fixed' system
*         (common parallax and proper motions)
*     I : secondary or tertiary of a 2- or 3-pointing 'independent' system
*         (no constraints on parallax or proper motions)
*     L : secondary or tertiary of a 2- or 3-pointing 'linear' system
*         (common parallax)
*     S : astrometric parameters from 'single-star merging' process.
*Note on Qual: Reliability of the double or multiple star solution:
*        A=good, B=fair, C=poor, D=uncertain, S=suspected non-single
*Note on Chart: the chart was produced:
*     D : from the STScI Digitized Sky Survey
*     G : from the Guide Star Catalog
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on r_SpType: the flag indicates the source, as:
*     1 : Michigan catalogue for the HD stars, vol. 1 (Houk+, 1975) <III/31>
*     2 : Michigan catalogue for the HD stars, vol. 2 (Houk, 1978)  <III/51>
*     3 : Michigan Catalogue for the HD stars, vol. 3 (Houk, 1982)  <III/80>
*     4 : Michigan Catalogue for the HD stars, vol. 4 (Houk+, 1988) <III/133>
*     G : updated after publication of the HIC <I/196>
*     K : General Catalog of Variable Stars, 4th Ed. (Kholopov+ 1988) <II/214>
*     S : SIMBAD data-base http://cdsweb.u-strasbg.fr/Simbad.html
*     X : Miscellaneous
*     A blank entry has no corresponding information.

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

C  Declarations for 'h_dm_com.dat'	! Double and Multiples: Component solutions -COMP

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

C  Position composed of: RAdeg DEdeg (Epoch=J1991.25)
      character*10  CCDM_1      ! CCDM number                              (DC1)
      integer*4     S           ! *Solution identifier                      (DC2)
      character*1   Type        ! *[FIL]Type of solution                    (DC3)
      character*1   Source_1    ! *[CFN] Solution source                    (DC4)
      character*1   Qual_1      ! *[ABCD] Solution quality                  (DC5)
      character*1   Notes_1     ! *[DGPWXYZ] Existence of notes             (DC6)
      integer*4     Nsys_1      ! Number of solutions for the system      (DCM1)
      integer*4     Ncomp_1     ! Number of components in this solution   (DCM2)
      integer*4     Nparm       ! Number of free parameters in solution   (DCM3)
      integer*4     Ncorr       ! Number of correlation records           (DCM4)
      integer*4     seq         ! Sequential component number             (DCM6)
      character*1   comp_id     ! Component identifier                     (DC7)
      integer*4     HIP_1       ! HIP number                               (DC8)
      real*4        Hp          ! (mag) Magnitude of component                   (DC9)
      real*4        e_Hp        ! (mag) Standard error of Hp magnitude          (DC10)
      real*4        BT          ! (mag) ? Magnitude of component, BT            (DC11)
      real*4        e_BT        ! (mag) ? Standard error of BT                  (DC12)
      real*4        VT          ! (mag) ? Magnitude of component, VT            (DC13)
      real*4        e_VT        ! (mag) ? Standard error of VT                  (DC14)
      real*8        RAdeg_1     ! (deg) alpha, degrees (ICRS, Epoch=J1991.25)   (DC15)
      real*8        DEdeg_1     ! (deg) delta, degrees (ICRS, Epoch=J1991.25)   (DC16)
      real*8        Plx_1       ! (mas) Trigonometric parallax                  (DC17)
      real*8        pmRA_1      ! (mas/yr) Proper motion in mu_alpha.cos(delta)ICRS(DC18)
      real*8        pmDE_1      ! (mas/yr) Proper motion in mu_delta in ICRS       (DC19)
      real*4        e_RAdeg_1   ! (mas) Standard error in RA*cos(DEdeg)         (DC20)
      real*4        e_DEdeg_1   ! (mas) Standard error in DE                    (DC21)
      real*4        e_Plx_1     ! (mas) Standard error in Plx                   (DC22)
      real*4        e_pmRA_1    ! (mas/yr) Standard error in pmRA                  (DC23)
      real*4        e_pmDE_1    ! (mas/yr) Standard error in pmDE                  (DC24)
      character*1   ref         ! Reference component for following data  (DC25)
      real*8        theta_1     ! (deg) *? Position angle                        (DC26)
      real*8        rho_1       ! (arcsec) ? Separation from reference component   (DC27)
      real*8        d_theta     ! (deg/yr) ? Rate of change of theta               (DC28)
      real*4        d_rho       ! (arcsec/yr) ? Rate of change of separation          (DC29)
      integer*4     seq_ref     ! *Sequential record number                (DCM7)
      character*25  flag        ! *Status flags for parameters             (DCM8)
*Note on S: a digit identifies different solutions pertaining to the
*     same CCDM number.
*Note on Type: Summary of double or multiple star solution:
*     F: fixed double or multiple system
*        (identical proper motions and parallaxes)
*     I: individual parallaxes and linear (relative) motion
*        (possible optical double star)
*     L: linear double or multiple system
*        (may have different proper motions but assumed to have same parallax)
*Note on Source: the source of the solution is given by this flag:
*     C: combined FAST and NDAC solution
*     F: solution taken from the FAST Consortium only
*     N: solution taken from the NDAC Consortium only
*Note on Qual: Reliability of the double or multiple star solution:
*        A=good  B=fair  C=poor   D=uncertain
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on theta: position angle relative to reference component
*Note on seq_ref: Sequential record number for the reference component, this
*     field is set to zero if DC25 is blank.
*Note on flag: Status flags for Hp, RAdeg, DEdeg, Plx, pmRA, pmDE.
*     1 = estimated,
*     0 = constrained to the value of the first component.

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

C  Declarations for 'h_dm_cor.dat'	! Double and Multiples: Component solutions -CORR

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

      character*10  CCDM_2      ! CCDM number                               (DC1)
      integer*4     S_1         ! *Solution identifier                       (DC2)
      character*1   Type_1      ! *[FIL]Type of solution                     (DC3)
      character*1   Source_2    ! *[CFN] Solution source                     (DC4)
      character*1   Qual_2      ! *[ABCD] Solution quality                   (DC5)
      character*1   Notes_2     ! *[DGPWXYZ] Existence of notes              (DC6)
      integer*4     Nsys_2      ! Number of solutions for the system       (DCM1)
      integer*4     Ncomp_2     ! Number of components in this solution    (DCM2)
      integer*4     Nparm_1     ! Number of free parameters in solution    (DCM3)
      integer*4     Ncorr_1     ! Number of correlation records            (DCM4)
      integer*4     seq_1       ! Sequential component number              (DCM6)
      integer*4     corr       (66) ! *[-99/999]? Correlation coefficients      (DCM7)
*Note on S: a digit identifies different solutions pertaining to the
*     same CCDM number.
*Note on Type: Summary of double or multiple star solution:
*     F: fixed double or multiple system
*        (identical proper motions and parallaxes)
*     I: individual parallaxes and linear (relative) motion
*        (possible optical double star)
*     L: linear double or multiple system
*        (may have different proper motions but assumed to have same parallax)
*Note on Source: the source of the solution is given by this flag:
*     C: combined FAST and NDAC solution
*     F: solution taken from the FAST Consortium only
*     N: solution taken from the NDAC Consortium only
*Note on Qual: Reliability of the double or multiple star solution:
*        A=good  B=fair  C=poor   D=uncertain
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on corr: the correlation records contain the correlation coefficients among
*     the maximum set of 6N_C parameters. The number of correlation coefficients
*     is 6N_C(6N_C-1)/2=66, 153 and 276 for N_C=2,3,4. A correlation record
*     contains up to 66 coefficients; thus respectively 1,3 and 5 records are
*     needed for a double, triple or quadruple star. The correlation coefficients
*     are coded as integers between -99 and 999 using the arcsine
*     transformation. The order of the correlation coefficients is indicated by
*     the following table:
*            Hp1 RA1 DE1 Plx1 pmRA1 pmDE1  Hp2 RA2 DE2 Plx2 pmRA2 pmDE2  Hp3 ...
*     Hp1    1   r1  r2   r4   r7    r11   r16 r22 r29  r37  r46   r56   r67 ...
*     RA1    r1  1   r3   r5   r8    r12   r17 r23 r30  r38  r47   r57   r68 ...
*     DE1    r2  r3  1    r6   r9    r13   r18 r24 r31  r39  r48   r58   r69 ...
*     Plx1   r4  r5  r6   1    r10   r14   r19 r25 r32  r40  r49   r59   r70 ...
*     pmRA1  r7  r8  r9   r10  1     r15   r20 r26 r33  r41  r50   r60   r71 ...
*     pmDE1  r11 r12 r13  r14  r15   1     r21 r27 r34  r42  r51   r61   r72 ...
*     Hp2    r16 r17 r18  r19  r20   r21   1   r28 r35  r43  r52   r62   r73 ...
*     RA2    r22 r23 r24  r25  r26   r27   r28 1   r36  r44  r53   r63   r74 ...
*     DE2    r29 r30 r31  r32  r33   r34   r35 r36 1    r45  r54   r64   r75 ...
*     Plx2   r37 r38 r39  r40  r41   r42   r43 r44 r45  1    r55   r65   r76 ...
*     pmRA2  r46 r47 r48  r49  r50   r51   r52 r53 r54  r55  1     r66   r77 ...
*     pmDE2  r56 r57 r58  r59  r60   r61   r62 r63 r64  r65  r66   1     r78 ...
*     Hp3    r67 r68 r69  r70  r71   r72   r73 r74 r75  r76  r77   r78   1  ...
*     ...    ... ... ...  ...  ...   ...   ... ... ...  ...  ...  ...   ... ...
*     These correlation coefficients are written as 66I3

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

C  Declarations for 'hip_dm_g.dat'	! Double and Multiples: Acceleration solutions

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

      integer*4     HIP_2       ! Identifier (HIP number)                   (DG1)
      real*8        gRA         ! (mas/yr2) *Acceleration gRA = d(pmRA)/dt             (DG2)
      real*8        gDE         ! (mas/yr2) *Acceleration gDE = d(pmDE)/dt             (DG3)
      real*8        e_gRA       ! (mas/yr2) Standard error of gRA                     (DG4)
      real*8        e_gDE       ! (mas/yr2) Standard error of gDE                     (DG5)
      real*4        Fg          ! *Significance of the g terms               (DG6)
      real*8        dgRA        ! (mas/yr3) *? dgRA = d2(pmRA)/dt2                     (DG7)
      real*8        dgDE        ! (mas/yr3) *? dgDE = d2(pmDE)/dt2                     (DG8)
      real*8        e_dgRA      ! (mas/yr3) ? Standard error of dgRA                  (DG9)
      real*8        e_dgDE      ! (mas/yr3) ? Standard error of dgDE                 (DG10)
      real*4        Fdg         ! *? Signifance of the dg terms             (DG11)
      character*1   Notes_3     ! *[DGPWXYZ] Existence of notes             (DG12)
      integer*4     num         ! *[7,9] Number of astrometric parameters   (DGM1)
      integer*4     corr_1     (36) ! *[-99/999]? Correlation coefficients      (DGM2)
*Note on gRA: component in right ascension of the apparent acceleration
*     of the photocentre at epoch J1991.25.
*Note on gDE: component in declination of the apparent acceleration
*     of the photocentre at epoch J1991.25.
*Note on Fg: the quadratic model is only adopted if the g terms
*     are significant (Fg>3.44)
*Note on dgRA: component in right ascension of the rate of change of the
*     apparent acceleration of the photocentre at epoch J1991.25.
*Note on dgDE: component in declination of the rate of change of the
*     apparent acceleration of the photocentre at epoch J1991.25.
*Note on Fdg: The cubic model is only adopted if the dg terms are
*     significant (Fdg>3.44)
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on num: for a quadratic model of the photocentric motion this
*     parameter is 7, for a cubic model it is 9.
*Note on corr: the complete set of n(n-1)*0.5 correlation coefficients
*     (where n=7 for a quadratic  and n=9 for a cubic model of the
*     photocentric motion) is given in the order indicated by the
*     following table:
*              RA   Dec  Plx  pmRA pmDE gRA  gDE  dgRA dgDE
*     RA       1    r1   r2   r4   r7   r11  r16  r22  r29
*     Dec      r1   1    r3   r5   r8   r12  r17  r23  r30
*     Plx      r2   r3   1    r6   r9   r13  r18  r24  r31
*     pmRA     r4   r5   r6   1    r10  r14  r19  r25  r32
*     pmDE     r7   r8   r9   r10  1    r15  r20  r26  r33
*     gRA      r11  r12  r13  r14  r15  1    r21  r27  r34
*     gDE      r16  r17  r18  r19  r20  r21  1    r28  r35
*     dgRA     r22  r23  r24  r25  r26  r27  r28  1    r36
*     dgDE     r29  r30  r31  r32  r33  r34  r35  r36  1

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

C  Declarations for 'hip_dm_o.dat'	! Double and Multiples: Orbital solutions

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

      integer*4     HIP_3       ! Identifier (HIP)                         (D01)
      real*8        P           ! (d) Orbital period                           (DO2)
      real*8        T           ! (d) *Time of periastron passage               (DO3)
      real*8        a0          ! (mas) Semi-major axis of photocentric orbit    (DO4)
      real*4        ecc         ! [0,1] Eccentricity                       (DO5)
      real*4        w           ! (deg) *[0,360] Argument of periastron           (DO6)
      real*4        i           ! (deg) *[0,180] Inclination                      (DO7)
      real*4        Omega       ! (deg) *[0,360] Position angle of the node       (DO8)
      real*8        e_P         ! (d) ? Standard error of P                    (DO9)
      real*8        e_T         ! (d) ? Standard error of T                   (DO10)
      real*4        e_a0        ! (mas) ? Standard error of a0                  (DO11)
      real*4        e_ecc       ! ? Standard error of ecc                 (DO12)
      real*4        e_w         ! (deg) ? Standard error of w                   (DO13)
      real*4        e_i         ! (deg) ? Standard error of i                   (DO14)
      real*4        e_Omega     ! (deg) ? Standard error of Omega               (DO15)
      integer*4     dmRef       ! *? reference to the literature           (DO16)
      character*1   Notes_4     ! *[DGPWXYZ] Existence of notes            (DO17)
      character*12  flag_1      ! *Status flags for the parameters         (DOM1)
      integer*4     corr_2     (66) ! *[-99/999]?=450 Correlation coefficients (DOM2)
*Note on T: this is the date when the photocentre is closest to the
*     centre of mass in the orbital plane. It is equivalent to the
*     closest approach of the stellar components.
*Note on w: this is the angle in the orbital plane from the line of
*     nodes to the major axis, measured from the nodal point (DO8) to
*     the periastron in the direction of motion.
*Note on i: the inclination of the orbital plane to the tangent plane
*     of the sky.  Taken to be in the first quadrant if the apparent
*     motion is direct (counter-clockwise) and in the second quadrant
*     for retrograde (clockwise) apparent motion.
*Note on Omega: This is the position angle (measured counter-clockwise
*     as seen on the sky from the +Dec direction) of the line of nodes,
*     or the intersection of the orbital and tangent planes. If the
*     radial motion of the components is known from spectroscopic studies,
*     then Omega should give the position angle of the ascending node,
*     at which the primary star crosses the tangent plane while receding
*     from the observer. In the absence of spectroscopic information
*     Omega refers to the mode with the smallest positive position angle.
*Note on dmRef: A reference number in this field points to references in the
*     printed catalogue, explained in dmsa_o.doc file.
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on flag: status flags for the 12 astrometric and orbital
*     parameters taken in the order indicated below.
*     1= estimated; 0= not estimated.
*Note on corr: the correlation coefficients in bytes 140-337 are given
*     in the table below. Correlation coefficients which are undefined
*     (corresponding to blanks in DO9-15, and a status flag=0 in
*     DOB127) are set to zero. All correlation coefficients are coded
*     as integers between -99 and 999 using the arcsine transformation.
*     Undefined values are coded as 450.
*           RA   Dec  Plx  pmRA pmDE  P    T   a0    e    w    i Omega
*     RA    1    r1   r2   r4   r7   r11  r16  r22  r29  r37  r46 r56
*     Dec   r1   1    r3   r5   r8   r12  r17  r23  r30  r36  r47 r57
*     Plx   r2   r3   1    r6   r9   r13  r18  r24  r31  r37  r46 r58
*     pmRA  r4   r5   r6   1    r10  r14  r19  r25  r32  r38  r47 r59
*     pmDE  r7   r8   r9   r10  1    r15  r20  r26  r33  r39  r48 r60
*     P     r11  r12  r13  r14  r15  1    r21  r27  r34  r40  r49 r61
*     T     r16  r17  r18  r19  r20  r21  1    r28  r35  r41  r50 r62
*     a0    r22  r23  r24  r25  r26  r27  r28  1    r36  r42  r51 r63
*     e     r29  r30  r31  r32  r33  r34  r35  r36  1    r43  r52 r64
*     w     r37  r38  r39  r40  r41  r42  r43  r44  r45  1    r53 r65
*     i     r46  r47  r48  r49  r50  r51  r52  r53  r54  r55  1   r66
*     Omega r56  r57  r59  r59  r60  r61  r62  r63  r64  r65  r66 1
*     These correlation coefficients are written as 66I3

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

C  Declarations for 'hip_dm_v.dat'	! Double and Multiples: VIM solutions

      integer*4 nr__5
      parameter (nr__5=288)	! Number of records
      character*144 ar__5  	! Full-size record

      integer*4     HIP_4       ! Identifier (HIP)                         (DV1)
      real*4        Hp_ref      ! (mag) *Reference magnitude                      (DV2)
      real*8        DRA         ! (mas) VIM element in right ascension           (DV3)
      real*8        DDE         ! (mas) VIM element in declination               (DV4)
      real*8        e_DRA       ! (mas) Standard error of DRA                    (DV5)
      real*8        e_DDE       ! (mas) Standard error of DDE                    (DV6)
      real*4        FD          ! *Significance of the VIM elements         (DV7)
      real*4        theta_C     ! (deg) *Position angle of the constant component (DV8)
      real*4        e_theta_C   ! (deg) Standard error of theta_C                (DV9)
      real*4        minSep      ! (mas) Lower limit for separation of binary    (DV10)
      real*4        dvar        ! (mas) *Displacement of photocentre             (DV11)
      character*1   Notes_5     ! *[DGPWXYZ] Existence of notes            (DV12)
      integer*4     corr_3     (21) ! *[-99/999]Correlation coefficients       (DVM1)
*Note on Hp_ref: the reference magnitude is freely chosen and defines the
*     reference point for the object. The positional parameters RA and Dec (H8,
*     H9 in hip_main.dat) and the VIM elements DV3,DV4 depend on the chosen
*     Hp_ref. DV7-11 do not depend on Hp_ref.
*Note on FD: The VIM solution is only accepted if FD>2.15
*Note on theta_C: position angle of the constant component of the binary with
*     respect to the variable component, measured counterclockwise as seen on
*     the sky from the +Dec direction.
*Note on dvar: displacement of photocentre between minimum and maximum
*     luminosity of the system. Indicates the size of the VIM effect.
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on corr: the correlation coefficients in bytes 82-144 are given in the
*     following sequence:
*           RA   Dec  Plx  pmRA pmDE DRA  DDE
*     RA    1    r1   r2   r4   r7   r11  r16
*     Dec   r1   1    r3   r5   r8   r12  r17
*     Plx   r2   r3   1    r6   r9   r13  r18
*     pmRA  r4   r5   r6   1    r10  r14  r19
*     pmDE  r7   r8   r9   r10  1    r15  r20
*     DRA   r11  r12  r13  r14  r15  1    r21
*     DDE   r16  r17  r18  r19  r20  r21  1
*     These correlation coefficients are written as 21I3

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

C  Declarations for 'hip_dm_x.dat'	! Double and Multiples: Stochastic solutions

      integer*4 nr__6
      parameter (nr__6=1561)	! Number of records
      character*22 ar__6  	! Full-size record

      integer*4     HIP_5       ! Identifier (HIP)                         (DX1)
      real*4        epsilon     ! (mas) Cosmic error, epsilon                    (DX2)
      real*4        e_epsilon   ! (mas) Standard error of epsilon                (DX3)
      character*1   Notes_6     ! *[DGPWXYZ] Existence of notes             (DX4)
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P

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

C  Declarations for 'hip_va_1.dat'	! Variability Annex: Periodic variables

      integer*4 nr__7
      parameter (nr__7=2712)	! Number of records
      character*142 ar__7  	! Full-size record

      integer*4     HIP_6       ! Identifier (HIP)                          (P1)
      character*1   NewVar      ! [*] Flag if new variable                  (P2)
      character*12  SpType_1    ! Spectral type                             (P3)
      character*1   HvarType_1  ! *[CDMPRU] Variability type (1-letter)      (P4)
      character*6   VarType     ! *Variability type as in GCVS/NSV           (P5)
      character*1   n_VarType   ! [*] Flag if newly classified by Hipparcos (P6)
      real*4        maxMag      ! (mag) Magnitude at max from curve fitting       (P7)
      character*1   l_minMag    ! *[>] Limit flag (>)                        (P8)
      real*4        minMag      ! (mag) Magnitude at min from curve fitting       (P9)
      real*4        log_sA_A    ! ? log_10(sigma_A / A)                    (P10)
      real*8        Period_1    ! (d) ? Mean period in days                    (P11)
      real*4        log_sP      ! ([d]) ? log_10(sigma_P)                        (P12)
      real*8        Ep_2440000  ! (d) ? Epoch (JD-2440000) of zero phase       (P13)
      integer*4     q_Ep_2440000 ! *[0/5]? Precision flag                    (P14)
      character*1   morePhoto_1 ! [ABC] Light curve Annex                  (P15)
      character*1   Notes_7     ! *[DGPWXYZ] Existence of notes             (P16)
      character*12  VarName     ! Variable star name                       (P17)
      real*8        period_2    ! (d) ? Period from literature                 (P18)
      real*8        epoch       ! (d) ? Epoch  from literature                 (P19)
      real*4        max         ! (mag) ? Magnitude at max from literature       (P20)
      real*4        min         ! (mag) ? Magnitude at min from literature       (P21)
      character*1   Band        ! *[UBVKIRPYb] Photometric band             (P22)
      character*1   refFlag     ! [R] Reference in printed catalogue       (P23)
*Note on HvarType: Hipparcos-defined type of variability:
*     C : no variability detected ("constant")
*     D : duplicity-induced variability
*     M : possibly micro-variable (amplitude < 0.03mag)
*     P : periodic variable
*     R : V-I colour index was revised due to variability analysis
*     U : unsolved variable which does not fall in the other categories
*Note on VarType: this is the 6-letter type defined
*        in the General Catalog of Variable Stars <II/214>
*        (see also the "Types of Variability" section below)
*Note on l_minMag: the flag (>) indicates that the true magnitude at minimum
*     luminosity is likely to be larger than the value of minMag given in P9.
*Note on q_Ep-2440000: the value is 1-log10(sigma_epoch), i.e.
*     0 for an accuracy of about 10 days,
*     1 for an accuracy of about 1 day,
*     2 for an accuracy of about 0.1day,
*     3 for an accuracy of about 0.01day,
*     4 for an accuracy of about 0.001day,
*     5 for an accuracy of about 0.0001day.
*Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on Band: U,B,V,K,I,R refer to Johnson system, or closely related
*     systems; P refers to photographic magnitudes, Y and b for the
*     Stroemgren y and b bands.

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

C  Declarations for 'hip_va_2.dat'	! Variability Annex: Unsolved variables

      integer*4 nr__8
      parameter (nr__8=5542)	! Number of records
      character*142 ar__8  	! Full-size record

      integer*4     HIP_7       ! Identifier (HIP number)                   (U1)
      character*1   NewVar_1    ! [*] Flag if new variable                  (U2)
      character*12  SpType_2    ! Spectral type                             (U3)
      character*1   HvarType_2  ! *[CDMPRU] Variability type (1-letter)      (U4)
      character*6   VarType_1   ! *Variability type (6-letter) as in GCVS/NSV(U5)
      character*1   n_VarType_1 ! [*] Flag if newly classified by Hipparcos (U6)
      real*4        maxMag_1    ! (mag) *Magnitude at max from curve fitting       (U7)
      character*1   l_minMag_1  ! [>] Limit flag (>)                        (U8)
      real*4        minMag_1    ! (mag) *Magnitude at min from curve fitting       (U9)
      real*4        med_Hp      ! (mag) Median Hp                                (U10)
      real*8        A           ! (mag) Intrinsic variability amplitude          (U11)
      real*4        e_A         ! (mag) Standard error of A                      (U12)
      character*1   morePhoto_2 ! [ABC] Light curve Annex                  (U15)
      character*1   Notes_8     ! *[DGPWXYZ] Notes                          (U16)
      character*12  VarName_1   ! Variable star name                       (U17)
      real*8        period_3    ! (d) ? Period from literature                 (U18)
      real*8        Ep_2440000_1 ! (d) ? Epoch (JD-2440000) from literature     (U19)
      real*4        max_1       ! (mag) ? Magnitude at max from literature       (U20)
      real*4        min_1       ! (mag) ? Magnitude at min from literature       (U21)
      character*1   Band_1      ! *[UBVKIRPYb] Photometric band             (U22)
      character*1   refFlag_1   ! [R] Reference in printed catalogue       (U23)
*Note on HvarType: Hipparcos-defined type of variability:
*     C : no variability detected ("constant")
*     D : duplicity-induced variability
*     M : possibly micro-variable (amplitude < 0.03mag)
*     P : periodic variable
*     R : V-I colour index was revised due to variability analysis
*     U : unsolved variable which does not fall in the other categories
*Note on VarType: this is the 6-letter type defined
*        in the General Catalog of Variable Stars <II/214>
*        (see also the "Types of Variability" section below)
*Note on maxMag: this is given if significant
*Note on minMag: this is given if significant
*Note on Notes: Note on Notes: the flag has the following meaning:
*     D : double and multiple systems note only (note in hd_notes.doc file)
*     G : general note only (note in hg_notes.doc file)
*     P : photometric notes only (note in hp_notes.doc file)
*     W : D + P
*     X : D + G
*     Y : G + P
*     Z : D + G + P
*Note on Band: U,B,V,K,I,R refer to Johnson system, or closely related
*     systems; P refers to photographic magnitudes, Y and b for the
*     Stroemgren y and b bands.

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

C  Declarations for 'solar_ha.dat'	! Solar System Annex: Astrometric catalogue

      integer*4 nr__9
      parameter (nr__9=5609)	! Number of records
      character*64 ar__9  	! Full-size record

C  Position composed of: RAdeg DEdeg
      integer*4     ID          ! Object number                         (SHA1)
      real*8        RAdeg_2     ! (deg) Reference point RA  (ICRS system)     (SHA2)
      real*8        DEdeg_2     ! (deg) Reference point Dec (ICRS system)     (SHA3)
      real*8        Ep_2440000_2 ! (d) *Measurement epoch                     (SHA4)
      real*4        delay       ! (s) *Light delay time                      (SHA5)
      real*8        theta_2     ! (deg) *Position angle, theta                 (SHA6)
      real*4        e_lambda    ! (mas) Estimated standard error of abscissa  (SHA7)
      integer*4     flag_2      ! *[1,2]FAST or NDAC flag                (SHA8)
*Note on Ep-2440000: the measurement epoch is specified in JD with respect to
*     JD(TT)2440000.0 and is corrected to the geocentre.
*Note on delay: this gives the applied light time delay in the geocentric
*     direction of the observed object between the satellite and the Earth.
*Note on theta: The position angle of the slit coordinate direction w is
*     reckoned positive from North through East.
*Note on flag: If the transit corresponds to an NDAC record the flag is 1.
*     For a FAST record the flag is 2.

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

C  Declarations for 'solar_hp.dat'	! Solar System Annex: Photometric catalogue

      integer*4 nr__10
      parameter (nr__10=2639)	! Number of records
      character*63 ar__10 	! Full-size record

      integer*4     ID_1        ! Object number                         (SHP1)
      real*8        Ep_2440000_3 ! (d) *Measurement epoch                     (SHP2)
      real*8        Hp_dc       ! (mag) Magnitude from unmodulated signal     (SHP3)
      real*4        e_Hp_dc     ! (mag) Standard error on Hp_dc               (SHP4)
      real*8        Hp_ac       ! (mag) Magnitude from modulated signal       (SHP5)
      real*4        e_Hp_ac     ! (mag) Standard error on Hp_ac               (SHP6)
      real*4        r           ! (AU) Distance: Sun-asteroid                (SHP7)
      real*4        Delta       ! (AU) Distance: satellite-asteroid          (SHP8)
      real*4        alpha       ! (deg) *Solar phase angle                     (SHP9)
*Note on Ep-2440000: the measurement epoch is specified in JD with respect
*     to JD(TT)2440000.0 and is not corrected to the geocentre  (i.e. the
*     light-time delay between the satellite and the earth is also neglected.
*Note on alpha: the solar phase angle is the angle between the Sun and the
*     satellite as viewed from the asteroid.

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

C  Declarations for 'solar_t.dat'	! Solar System Annex: Tycho astrometry/photometry

      integer*4 nr__11
      parameter (nr__11=291)	! Number of records
      character*95 ar__11 	! Full-size record

C  Position composed of: RAdeg DEdeg
      integer*4     ID_2        ! Object number                           (ST1)
      real*8        Ep_2440000_4 ! (d) *Measurement epoch                       (ST2)
      real*8        RAdeg_3     ! (deg) Right ascension (ICRS system)           (ST3)
      real*8        DEdeg_3     ! (deg) Declination (ICRS system)               (ST4)
      real*4        BTmag_1     ! (mag) ? Mean BT magnitude                     (ST5)
      real*4        VTmag_1     ! (mag) ? Mean VT magnitude                     (ST6)
      integer*4     Flag_3      ! *[1,2] Transit flag                      (ST7)
      real*4        e_RAdeg_2   ! (mas) Standard error on RA                    (ST8)
      real*4        e_DEdeg_2   ! (mas) Standard error on Dec                   (ST9)
      real*4        DE_RA_1     ! Correlation, Dec/RA*cos(delta)         (ST10)
      real*4        theta_3     ! (deg) Position angle of slit (direction w)   (ST11)
      integer*4     sign_z      ! *[-1/1] Inclined slit flag              (ST12)
      real*4        e_incl      ! (mas) Standard error on slit pos. (inclined) (ST13)
      real*4        e_vert      ! (mas) Standard error on slit pos. (vertical) (ST14)
*Note on Ep-2440000: the measurement epoch is specified in JD with respect
*     to JD(TT)2440000.0 and is not corrected to the geocentre (i.e. the
*     light-time delay between the satellite and the earth is also neglected.
*Note on Flag: if flag =1 then only one crossing of the field of view has
*     been detected or retained. If Flag n>1, then n candidate observed
*     transits have been detected or retained for that predicted
*     observation. ST1, ST2, ST7, ST11, ST12 are identical in those cases.
*Note on sign_z: if the transit occurred in the upper branch of the
*     inclined slits, the flag is +1, if the transit occurred in the
*     lower branch of the inclined slits, the flag is -1.

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

C  Declarations for 'hd_notes.doc'	! Hipparcos notes: Double and multiple systems

      integer*4 nr__12
      parameter (nr__12=2622)	! Number of records
      character*97 ar__12 	! Full-size record

      integer*4     HIP_8       ! HIP number
      character*1   Note1       ! [DG] Double and multiple or General note
      character*1   Note2       ! [DP] Double and multiple or Photometric note
      integer*4     Ntot        ! Total number of lines for the HIP object
      integer*4     Nline       ! Running line number in range [1,Ntot]
      character*80  Text        ! Text of note.

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

C  Declarations for 'hg_notes.doc'	! Hipparcos notes: General notes

      integer*4 nr__13
      parameter (nr__13=3898)	! Number of records
      character*97 ar__13 	! Full-size record

      integer*4     HIP_9       ! HIP number
      character*1   Note1_1     ! [DG] Double and multiple or General note
      character*1   Note2_1     ! [DP] Double and multiple or Photometric note
      integer*4     Ntot_1      ! Total number of lines for the HIP object
      integer*4     Nline_1     ! Running line number in range [1,Ntot]
      character*80  Text_1      ! Text of note.

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

C  Declarations for 'hp_notes.doc'	! Hipparcos notes: Photometric notes

      integer*4 nr__14
      parameter (nr__14=2444)	! Number of records
      character*97 ar__14 	! Full-size record

      integer*4     HIP_10      ! HIP number
      character*1   Note1_2     ! [DG] Double and multiple or General note
      character*1   Note2_2     ! [DP] Double and multiple or Photometric note
      integer*4     Ntot_2      ! Total number of lines for the HIP object
      integer*4     Nline_2     ! Running line number in range [1,Ntot]
      character*80  Text_2      ! Text of note.

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

C  Declarations for 'hp_refs.doc'	! References Hipparcos stars

      integer*4 nr__15
      parameter (nr__15=33769)	! Number of records
      character*19 ar__15 	! Full-size record

      integer*4     HIP_11      ! HIP number
      integer*4     Ntot_3      ! Total number of lines for the HIP object
      integer*4     Nline_3     ! Running line number in range [1,Ntot]
      real*4        nRef        ! Reference number (explained in hp_auth.doc)

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

C  Declarations for 'hp_auth.doc'	! References of hp_notes.doc

      integer*4 nr__16
      parameter (nr__16=4335)	! Number of records
      character*80 ar__16 	! Full-size record

      real*4        nRef_1      ! Reference number from hp_refs.doc
      character*70  Text_3      ! Text of reference

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

C  Declarations for 'dmsa_o.doc'	! References of hip_dm_o.dat

      integer*4 nr__17
      parameter (nr__17=118)	! Number of records
      character*80 ar__17 	! Full-size record

      integer*4     dmRef_1     ! Reference number (hip_dm_o.dat file)
      integer*4     Ntot_4      ! Total number of lines for the reference
      integer*4     Nline_4     ! Running line number in range [1,Ntot]
      character*72  Text_4      ! Text of reference

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

C  Declarations for 'tyc_main.dat'	! The main part of Tycho Catalogue

      integer*4 nr__18
      parameter (nr__18=1058332)	! Number of records
      character*350 ar__18 	! Full-size record

C  Position composed of: RAdeg DEdeg (Epoch=J1991.25)
      character*1   Catalog_1   ! [T] Catalogue (T = Tycho)                 (T0)
      character*12  TYC         ! *TYC1-3 (TYC number)                       (T1)
      character*1   Proxy_1     ! [HT]? Proximity flag                      (T2)
      character*11  RAhms_1     ! Right ascension in h m s, ICRS (J1991.25) (T3)
      character*11  DEdms_1     ! Declination in deg ' ", ICRS (J1991.25)   (T4)
      real*4        Vmag_1      ! (mag) ? Magnitude in Johnson V                  (T5)
      character*1   r_Vmag_1    ! *[BDTV] Source of magnitude                (T7)
      real*8        RAdeg_4     ! (deg) *alpha, degrees (ICRS, Epoch=J1991.25)     (T8)
      real*8        DEdeg_4     ! (deg) *delta, degrees (ICRS, Epoch=J1991.25)     (T9)
      character*1   AstroRef_1  ! *[X]? Reference flag for astrometry       (T10)
      real*8        Plx_2       ! (mas) *? Trigonometric parallax                 (T11)
      real*8        pmRA_2      ! (mas/yr) *? Proper motion mu_alpha.cos(delta), ICRS(T12)
      real*8        pmDE_2      ! (mas/yr) *? Proper motion mu_delta, ICRS           (T13)
      real*4        e_RAdeg_3   ! (mas) *? Standard error in RA*cos(delta)        (T14)
      real*4        e_DEdeg_3   ! (mas) *? Standard error in DE                   (T15)
      real*4        e_Plx_2     ! (mas) *? Standard error in Plx                  (T16)
      real*4        e_pmRA_2    ! (mas/yr) *? Standard error in pmRA                 (T17)
      real*4        e_pmDE_2    ! (mas/yr) *? Standard error in pmDE                 (T18)
      real*4        DE_RA_2     ! [-1/1]? Correlation, DE/RA*cos(delta)    (T19)
      real*4        Plx_RA_1    ! [-1/1]? Correlation, Plx/RA*cos(delta)   (T20)
      real*4        Plx_DE_1    ! [-1/1]? Correlation, Plx/DE              (T21)
      real*4        pmRA_RA_1   ! [-1/1]? Correlation, pmRA/RA*cos(delta)  (T22)
      real*4        pmRA_DE_1   ! [-1/1]? Correlation, pmRA/DE             (T23)
      real*4        pmRA_Plx_1  ! [-1/1]? Correlation, pmRA/Plx            (T24)
      real*4        pmDE_RA_1   ! [-1/1]? Correlation, pmDE/RA*cos(delta)  (T25)
      real*4        pmDE_DE_1   ! [-1/1]? Correlation, pmDE/DE             (T26)
      real*4        pmDE_Plx_1  ! [-1/1]? Correlation, pmDE/Plx            (T27)
      real*4        pmDE_pmRA_1 ! [-1/1]? Correlation, pmDE/pmRA           (T28)
      integer*4     Nastro      ! ? Number of transits for astrometry      (T29)
      real*4        F2_1        ! *? Goodness-of-fit parameter              (T30)
      integer*4     HIP_12      ! ? Hipparcos HIP number                   (T31)
      real*4        BTmag_2     ! (mag) ? Mean BT magnitude                      (T32)
      real*4        e_BTmag_1   ! (mag) ? Standard error in BTmag                (T33)
      real*4        VTmag_2     ! (mag) ? Mean VT magnitude                      (T34)
      real*4        e_VTmag_1   ! (mag) ? Standard error in VTmag                (T35)
      character*1   r_BTmag     ! *[DMNT] Source of photometry              (T36)
      real*4        B_V_1       ! (mag) ? Johnson B-V colour                     (T37)
      real*4        e_B_V_1     ! (mag) ? Standard error on B-V                  (T38)
      integer*4     Q           ! *? Astrometric quality flag, Q            (T40)
      real*4        Fs          ! ? Signal-to-noise ratio of the star image(T41)
      character*1   Source_3    ! *[HPR] Source of astrometric data         (T42)
      integer*4     Nphoto      ! ? Number of transits for photometry      (T43)
      real*4        VTscat      ! (mag) ? Estimate of VTmag scatter              (T44)
      real*4        VTmax       ! (mag) ? VTmag at maximum (15th percentile)     (T45)
      real*4        VTmin       ! (mag) ? VTmag at minimum (85th percentile)     (T46)
      character*1   Var         ! *[GN]? Known variability from GCVS/NSV    (T47)
      character*1   VarFlag_1   ! *[UVW]? Variability from Tycho            (T48)
      character*1   MultFlag_1  ! *[DRSYZ]? Duplicity from Tycho            (T49)
      character*1   morePhoto_3 ! [AB]  Epoch photometry in Annex A or B   (T50)
      character*2   m_HIP_1     ! CCDM component identifier                (T51)
      integer*4     PPM         ! *[1/789676]? PPM and Supplement           (T52)
      integer*4     HD_1        ! [1/359083]? HD cat. <III/135>            (T53)
      character*10  BD_1        ! Bonner DM <I/119>, <I/122>               (T54)
      character*10  CoD_1       ! Cordoba DM <I/114>                       (T55)
      character*10  CPD_1       ! Cape Photographic DM <I/108>             (T56)
      character*1   Remark      ! *[JKLM] Notes                             (T57)
*Note on TYC: the designation of an object in the Tycho Catalogue uses the
*     Guide Star Catalog numbering system (a region number (TYC1) and a number
*     within the region (TYC2)) followed by a Tycho specific component number
*     (TYC3).
*Note on r_Vmag: if non blank, the field has the following meaning:
*     B : no VTmag available, therefore BTmag was adopted
*     D : derived from approximate BTmag and VTmag (r_BTmag field T36 is 'D')
*     T : derived from approximate VTmag (r_BTmag field T36 is 'T')
*     V : no BTmag available, therefore VTmag was adopted in Vmag
*Note on RAdeg, DEdeg: right ascension and declination are
*     expressed in degrees for epoch J1991.25 (JD2448349.0625 (TT)) in the
*      ICRS (International Celestial Reference System, close to
*     J2000) reference system.
*Note on AstroRef: 'X' indicates a dubious astrometric reference star
*     in the context of the Tycho catalogue.
*Note on Plx, pmRA, pmDE, e_RAdeg, e_DEdeg, e_Plx, e_pmRA, e_pmDE:
*     For these fields the second decimal digit is always blank
*Note on F2: values exceeding +2.5 to +3 indicate a bad fit to the data.
*Note on r_BTmag: the source flag has the following meaning:
*     D : approximate BTmag, VTmag obtained for resolved double stars
*     M : BTmag and VTmag are median values rather than de-censored means
*     N : BTmag and VTmag are de-censored means
*     T : BTmag is not given, and VTmag is an estimate; these magnitudes
*         are systematically too bright by up to 1mag.
*     A blank indicates an Hipparcos star not observed by Tycho (T42 = H).
*Note on Q: the astrometric quality flag is defined by the following table,
*     where sigma(max) is the largest of the 5 astrometric standard errors:
*   ---------------------------------------------------------
*     Q  sigma(max)   Astrometric quality
*   ---------------------------------------------------------
*     1     <   5     very high
*     2     5- 10     very high
*     3    10- 25     high
*     4    25- 50     high
*     5    50-150     medium
*     6     < 150     perhaps non-single
*     7     < 150     low
*     8     < 150     perhaps non-stellar
*     9     ~ 200     low, position derived from TICR
*                     ('R' in Source field T42)
*   blank     ---     unassigned, 'H' in Source field T42
*                     (Hipparcos star not observed by Tycho)
*   ---------------------------------------------------------
*Note on Source:
*     H : Hipparcos star not observed by Tycho
*     P : only the position was determined (no proper motion, no parallax)
*     R : the position is derived from TICR catalog <I/197>
*Note on Var: this flag is set when the variability is known:
*     G : variable known in General Catalog of Variable Stars <II/214>
*     N : variable known in New Suspected Variables catalog <II/140>
*Note on VarFlag: this flag has the following meaning:
*     U : apparent variability in the Tycho data; may be due to duplicity
*     V : strong evidence of intrinsic variability
*     W : suspected intrinsic variability
*Note on MultFlag: unresolved duplicity status from Tycho data analysis:
*     D : duplicity clearly indicated (BTmag and VTmag refer to combined light)
*     R : duplicity weakly indicated, combined with indication of variability
*     S : duplicity suspected
*     Y : investigation for duplicity carried out on Tycho data,
*         no indication of duplicity was found
*     Z : investigation for duplicity not carried out
* blank : Hipparcos star not observed by Tycho
*Note on PPM: from the 3 parts of the PPM catalogue
*    North <I/146>, South <I/193>, and 90000Supplement <I/208>
*Note on Remark: if non-blank, the remark has the following meaning:
*     J : disagreement with position of magnitude in GSC1.1 catalog <I/220>
*     K : dubious Tycho parallax (Plx)
*     L : dubious Tycho proper motion (disagrees with PPM catalogue)
*     M : very uncertain Tycho magnitude (standard error larger than 0.3mag)

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

C  Loading file 'hip_main.dat'	! The Hipparcos Main Catalogue

C  Format for file interpretation

    1 format(
     +  A1,7X,I6,1X,A1,1X,A11,1X,A11,1X,F5.2,1X,I1,1X,A1,1X,F12.8,1X,
     +  F12.8,1X,A1,1X,F7.2,1X,F8.2,1X,F8.2,1X,F6.2,1X,F6.2,1X,F6.2,
     +  1X,F6.2,1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
     +  F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F5.2,8X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F5.3,1X,A1,1X,F4.2,1X,
     +  F4.2,1X,A1,1X,A1,1X,F7.4,1X,F6.4,1X,F5.3,1X,I3,1X,A1,1X,F5.2,
     +  1X,F5.2,1X,F7.2,1X,A1,1X,A1,1X,A1,1X,A10,1X,A1,1X,I2,1X,I2,1X,
     +  A1,1X,A1,1X,A1,1X,A2,1X,I3,1X,F7.3,1X,F5.3,1X,F5.2,1X,F4.2,1X,
     +  A1,1X,A1,1X,A1,1X,I6,1X,A10,1X,A10,1X,A10,1X,F4.2,1X,A12,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hip_main.dat')
      write(6,*) '....Loading file: hip_main.dat'
      do i__=1,118218
        read(1,'(A450)')ar__
        read(ar__,1)
     +  Catalog,HIP,Proxy,RAhms,DEdms,Vmag,VarFlag,r_Vmag,RAdeg,DEdeg,
     +  AstroRef,Plx,pmRA,pmDE,e_RAdeg,e_DEdeg,e_Plx,e_pmRA,e_pmDE,
     +  DE_RA,Plx_RA,Plx_DE,pmRA_RA,pmRA_DE,pmRA_Plx,pmDE_RA,pmDE_DE,
     +  pmDE_Plx,pmDE_pmRA,F1,F2,BTmag,e_BTmag,VTmag,e_VTmag,m_BTmag,
     +  B_V,e_B_V,r_B_V,V_I,e_V_I,r_V_I,CombMag,Hpmag,e_Hpmag,Hpscat,
     +  o_Hpmag,m_Hpmag,Hpmax,HPmin,Period,HvarType,moreVar,morePhoto,
     +  CCDM,n_CCDM,Nsys,Ncomp,MultFlag,Source,Qual,m_HIP,theta,rho,
     +  e_rho,dHp,e_dHp,Survey,Chart,Notes,HD,BD,CoD,CPD,v_V_I_red,
     +  SpType,r_SpType
        if(ar__(42:46) .EQ. '') Vmag = rNULL__
        if(ar__(48:48) .EQ. '') VarFlag = iNULL__
        if(ar__(52:63) .EQ. '') RAdeg = rNULL__
        if(ar__(65:76) .EQ. '') DEdeg = rNULL__
        if(ar__(80:86) .EQ. '') Plx = rNULL__
        if(ar__(88:95) .EQ. '') pmRA = rNULL__
        if(ar__(97:104) .EQ. '') pmDE = rNULL__
        if(ar__(106:111) .EQ. '') e_RAdeg = rNULL__
        if(ar__(113:118) .EQ. '') e_DEdeg = rNULL__
        if(ar__(120:125) .EQ. '') e_Plx = rNULL__
        if(ar__(127:132) .EQ. '') e_pmRA = rNULL__
        if(ar__(134:139) .EQ. '') e_pmDE = rNULL__
        if(ar__(141:145) .EQ. '') DE_RA = rNULL__
        if(ar__(147:151) .EQ. '') Plx_RA = rNULL__
        if(ar__(153:157) .EQ. '') Plx_DE = rNULL__
        if(ar__(159:163) .EQ. '') pmRA_RA = rNULL__
        if(ar__(165:169) .EQ. '') pmRA_DE = rNULL__
        if(ar__(171:175) .EQ. '') pmRA_Plx = rNULL__
        if(ar__(177:181) .EQ. '') pmDE_RA = rNULL__
        if(ar__(183:187) .EQ. '') pmDE_DE = rNULL__
        if(ar__(189:193) .EQ. '') pmDE_Plx = rNULL__
        if(ar__(195:199) .EQ. '') pmDE_pmRA = rNULL__
        if(ar__(201:203) .EQ. '') F1 = iNULL__
        if(ar__(205:209) .EQ. '') F2 = rNULL__
        if(ar__(218:223) .EQ. '') BTmag = rNULL__
        if(ar__(225:229) .EQ. '') e_BTmag = rNULL__
        if(ar__(231:236) .EQ. '') VTmag = rNULL__
        if(ar__(238:242) .EQ. '') e_VTmag = rNULL__
        if(ar__(246:251) .EQ. '') B_V = rNULL__
        if(ar__(253:257) .EQ. '') e_B_V = rNULL__
        if(ar__(261:264) .EQ. '') V_I = rNULL__
        if(ar__(266:269) .EQ. '') e_V_I = rNULL__
        if(ar__(275:281) .EQ. '') Hpmag = rNULL__
        if(ar__(283:288) .EQ. '') e_Hpmag = rNULL__
        if(ar__(290:294) .EQ. '') Hpscat = rNULL__
        if(ar__(296:298) .EQ. '') o_Hpmag = iNULL__
        if(ar__(302:306) .EQ. '') Hpmax = rNULL__
        if(ar__(308:312) .EQ. '') HPmin = rNULL__
        if(ar__(314:320) .EQ. '') Period = rNULL__
        if(ar__(341:342) .EQ. '') Nsys = iNULL__
        if(ar__(344:345) .EQ. '') Ncomp = iNULL__
        if(ar__(356:358) .EQ. '') theta = iNULL__
        if(ar__(360:366) .EQ. '') rho = rNULL__
        if(ar__(368:372) .EQ. '') e_rho = rNULL__
        if(ar__(374:378) .EQ. '') dHp = rNULL__
        if(ar__(380:383) .EQ. '') e_dHp = rNULL__
        if(ar__(391:396) .EQ. '') HD = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Catalog,HIP,Proxy,RAhms,DEdms,Vmag,VarFlag,r_Vmag,RAdeg,DEdeg,
     +  AstroRef,Plx,pmRA,pmDE,e_RAdeg,e_DEdeg,e_Plx,e_pmRA,e_pmDE,
     +  DE_RA,Plx_RA,Plx_DE,pmRA_RA,pmRA_DE,pmRA_Plx,pmDE_RA,pmDE_DE,
     +  pmDE_Plx,pmDE_pmRA,F1,F2,BTmag,e_BTmag,VTmag,e_VTmag,m_BTmag,
     +  B_V,e_B_V,r_B_V,V_I,e_V_I,r_V_I,CombMag,Hpmag,e_Hpmag,Hpscat,
     +  o_Hpmag,m_Hpmag,Hpmax,HPmin,Period,HvarType,moreVar,morePhoto,
     +  CCDM,n_CCDM,Nsys,Ncomp,MultFlag,Source,Qual,m_HIP,theta,rho,
     +  e_rho,dHp,e_dHp,Survey,Chart,Notes,HD,BD,CoD,CPD,v_V_I_red,
     +  SpType,r_SpType
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'h_dm_com.dat'	! Double and Multiples: Component solutions -COMP

C  Format for file interpretation

    2 format(
     +  A10,1X,I1,1X,A1,1X,A1,1X,A1,1X,A1,1X,I1,1X,I2,1X,I2,1X,I2,6X,
     +  I2,1X,A1,1X,I6,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
     +  F5.3,1X,F12.8,1X,F12.8,1X,F7.2,1X,F8.2,1X,F8.2,1X,F6.2,1X,
     +  F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,A1,1X,F7.3,1X,F8.3,1X,F8.3,1X,
     +  F6.3,1X,I2,1X,A25)

C  Effective file loading

      open(unit=1,status='old',file=
     +'h_dm_com.dat')
      write(6,*) '....Loading file: h_dm_com.dat'
      do i__=1,24588
        read(1,'(A238)')ar__1
        read(ar__1,2)
     +  CCDM_1,S,Type,Source_1,Qual_1,Notes_1,Nsys_1,Ncomp_1,Nparm,
     +  Ncorr,seq,comp_id,HIP_1,Hp,e_Hp,BT,e_BT,VT,e_VT,RAdeg_1,
     +  DEdeg_1,Plx_1,pmRA_1,pmDE_1,e_RAdeg_1,e_DEdeg_1,e_Plx_1,
     +  e_pmRA_1,e_pmDE_1,ref,theta_1,rho_1,d_theta,d_rho,seq_ref,flag
        if(ar__1(63:68) .EQ. '') BT = rNULL__
        if(ar__1(70:74) .EQ. '') e_BT = rNULL__
        if(ar__1(76:81) .EQ. '') VT = rNULL__
        if(ar__1(83:87) .EQ. '') e_VT = rNULL__
        if(ar__1(178:184) .EQ. '') theta_1 = rNULL__
        if(ar__1(186:193) .EQ. '') rho_1 = rNULL__
        if(ar__1(195:202) .EQ. '') d_theta = rNULL__
        if(ar__1(204:209) .EQ. '') d_rho = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  CCDM_1,S,Type,Source_1,Qual_1,Notes_1,Nsys_1,Ncomp_1,Nparm,
     +  Ncorr,seq,comp_id,HIP_1,Hp,e_Hp,BT,e_BT,VT,e_VT,RAdeg_1,
     +  DEdeg_1,Plx_1,pmRA_1,pmDE_1,e_RAdeg_1,e_DEdeg_1,e_Plx_1,
     +  e_pmRA_1,e_pmDE_1,ref,theta_1,rho_1,d_theta,d_rho,seq_ref,flag
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'h_dm_cor.dat'	! Double and Multiples: Component solutions -CORR

C  Format for file interpretation

    3 format(
     +  A10,1X,I1,1X,A1,1X,A1,1X,A1,1X,A1,1X,I1,1X,I2,1X,I2,1X,I2,6X,
     +  I2,1X,66I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'h_dm_cor.dat')
      write(6,*) '....Loading file: h_dm_cor.dat'
      do i__=1,12591
        read(1,'(A238)')ar__2
        read(ar__2,3)
     +  CCDM_2,S_1,Type_1,Source_2,Qual_2,Notes_2,Nsys_2,Ncomp_2,
     +  Nparm_1,Ncorr_1,seq_1,(corr(j__,),j__=1,66)
        if(ar__2(41:238) .EQ. '') corr = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  CCDM_2,S_1,Type_1,Source_2,Qual_2,Notes_2,Nsys_2,Ncomp_2,
     +  Nparm_1,Ncorr_1,seq_1,(corr(j__,),j__=1,66)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hip_dm_g.dat'	! Double and Multiples: Acceleration solutions

C  Format for file interpretation

    4 format(
     +  I6,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F7.2,1X,F7.2,1X,
     +  F7.2,1X,F7.2,1X,F5.2,1X,A1,1X,I1,1X,36I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hip_dm_g.dat')
      write(6,*) '....Loading file: hip_dm_g.dat'
      do i__=1,2622
        read(1,'(A195)')ar__3
        read(ar__3,4)
     +  HIP_2,gRA,gDE,e_gRA,e_gDE,Fg,dgRA,dgDE,e_dgRA,e_dgDE,Fdg,
     +  Notes_3,num,(corr_1(j__,),j__=1,36)
        if(ar__3(46:52) .EQ. '') dgRA = rNULL__
        if(ar__3(54:60) .EQ. '') dgDE = rNULL__
        if(ar__3(62:68) .EQ. '') e_dgRA = rNULL__
        if(ar__3(70:76) .EQ. '') e_dgDE = rNULL__
        if(ar__3(78:82) .EQ. '') Fdg = rNULL__
        if(ar__3(88:195) .EQ. '') corr_1 = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  HIP_2,gRA,gDE,e_gRA,e_gDE,Fg,dgRA,dgDE,e_dgRA,e_dgDE,Fdg,
     +  Notes_3,num,(corr_1(j__,),j__=1,36)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hip_dm_o.dat'	! Double and Multiples: Orbital solutions

C  Format for file interpretation

    5 format(
     +  I6,1X,F10.4,1X,F11.4,1X,F8.2,1X,F6.4,1X,F6.2,1X,F6.2,1X,F6.2,
     +  1X,F8.4,1X,F9.4,1X,F5.2,1X,F6.4,1X,F6.2,1X,F6.2,1X,F6.2,1X,I3,
     +  1X,A1,1X,A12,1X,66I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hip_dm_o.dat')
      write(6,*) '....Loading file: hip_dm_o.dat'
      do i__=1,235
        read(1,'(A337)')ar__4
        read(ar__4,5)
     +  HIP_3,P,T,a0,ecc,w,i,Omega,e_P,e_T,e_a0,e_ecc,e_w,e_i,e_Omega,
     +  dmRef,Notes_4,flag_1,(corr_2(j__,),j__=1,66)
        if(ar__4(68:75) .EQ. '') e_P = rNULL__
        if(ar__4(77:85) .EQ. '') e_T = rNULL__
        if(ar__4(87:91) .EQ. '') e_a0 = rNULL__
        if(ar__4(93:98) .EQ. '') e_ecc = rNULL__
        if(ar__4(100:105) .EQ. '') e_w = rNULL__
        if(ar__4(107:112) .EQ. '') e_i = rNULL__
        if(ar__4(114:119) .EQ. '') e_Omega = rNULL__
        if(ar__4(121:123) .EQ. '') dmRef = iNULL__
c    ..............Just test output...........
        write(6,5)
     +  HIP_3,P,T,a0,ecc,w,i,Omega,e_P,e_T,e_a0,e_ecc,e_w,e_i,e_Omega,
     +  dmRef,Notes_4,flag_1,(corr_2(j__,),j__=1,66)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hip_dm_v.dat'	! Double and Multiples: VIM solutions

C  Format for file interpretation

    6 format(
     +  I6,1X,F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F6.2,1X,
     +  F6.2,1X,F6.1,1X,F6.1,1X,A1,1X,21I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hip_dm_v.dat')
      write(6,*) '....Loading file: hip_dm_v.dat'
      do i__=1,288
        read(1,'(A144)')ar__5
        read(ar__5,6)
     +  HIP_4,Hp_ref,DRA,DDE,e_DRA,e_DDE,FD,theta_C,e_theta_C,minSep,
     +  dvar,Notes_5,(corr_3(j__,),j__=1,21)
c    ..............Just test output...........
        write(6,6)
     +  HIP_4,Hp_ref,DRA,DDE,e_DRA,e_DDE,FD,theta_C,e_theta_C,minSep,
     +  dvar,Notes_5,(corr_3(j__,),j__=1,21)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hip_dm_x.dat'	! Double and Multiples: Stochastic solutions

C  Format for file interpretation

    7 format(I6,1X,F6.2,1X,F6.2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hip_dm_x.dat')
      write(6,*) '....Loading file: hip_dm_x.dat'
      do i__=1,1561
        read(1,'(A22)')ar__6
        read(ar__6,7)HIP_5,epsilon,e_epsilon,Notes_6
c    ..............Just test output...........
        write(6,7)HIP_5,epsilon,e_epsilon,Notes_6
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hip_va_1.dat'	! Variability Annex: Periodic variables

C  Format for file interpretation

    8 format(
     +  I6,1X,A1,1X,A12,1X,A1,1X,A6,1X,A1,1X,F6.3,1X,A1,1X,F6.3,1X,
     +  F6.1,1X,F12.7,1X,F6.1,1X,F9.4,1X,I1,1X,A1,1X,A1,1X,A12,1X,
     +  F10.5,1X,F10.2,1X,F5.2,1X,F5.2,1X,A1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hip_va_1.dat')
      write(6,*) '....Loading file: hip_va_1.dat'
      do i__=1,2712
        read(1,'(A142)')ar__7
        read(ar__7,8)
     +  HIP_6,NewVar,SpType_1,HvarType_1,VarType,n_VarType,maxMag,
     +  l_minMag,minMag,log_sA_A,Period_1,log_sP,Ep_2440000,
     +  q_Ep_2440000,morePhoto_1,Notes_7,VarName,period_2,epoch,max,
     +  min,Band,refFlag
        if(ar__7(50:55) .EQ. '') log_sA_A = rNULL__
        if(ar__7(57:68) .EQ. '') Period_1 = rNULL__
        if(ar__7(70:75) .EQ. '') log_sP = rNULL__
        if(ar__7(77:85) .EQ. '') Ep_2440000 = rNULL__
        if(ar__7(87:87) .EQ. '') q_Ep_2440000 = iNULL__
        if(ar__7(106:115) .EQ. '') period_2 = rNULL__
        if(ar__7(117:126) .EQ. '') epoch = rNULL__
        if(ar__7(128:132) .EQ. '') max = rNULL__
        if(ar__7(134:138) .EQ. '') min = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  HIP_6,NewVar,SpType_1,HvarType_1,VarType,n_VarType,maxMag,
     +  l_minMag,minMag,log_sA_A,Period_1,log_sP,Ep_2440000,
     +  q_Ep_2440000,morePhoto_1,Notes_7,VarName,period_2,epoch,max,
     +  min,Band,refFlag
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hip_va_2.dat'	! Variability Annex: Unsolved variables

C  Format for file interpretation

    9 format(
     +  I6,1X,A1,1X,A12,1X,A1,1X,A6,1X,A1,1X,F6.3,1X,A1,1X,F6.3,1X,
     +  F6.3,1X,F12.3,1X,F6.3,13X,A1,1X,A1,1X,A12,1X,F10.5,1X,F10.2,
     +  1X,F5.2,1X,F5.2,1X,A1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hip_va_2.dat')
      write(6,*) '....Loading file: hip_va_2.dat'
      do i__=1,5542
        read(1,'(A142)')ar__8
        read(ar__8,9)
     +  HIP_7,NewVar_1,SpType_2,HvarType_2,VarType_1,n_VarType_1,
     +  maxMag_1,l_minMag_1,minMag_1,med_Hp,A,e_A,morePhoto_2,Notes_8,
     +  VarName_1,period_3,Ep_2440000_1,max_1,min_1,Band_1,refFlag_1
        if(ar__8(106:115) .EQ. '') period_3 = rNULL__
        if(ar__8(117:126) .EQ. '') Ep_2440000_1 = rNULL__
        if(ar__8(128:132) .EQ. '') max_1 = rNULL__
        if(ar__8(134:138) .EQ. '') min_1 = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  HIP_7,NewVar_1,SpType_2,HvarType_2,VarType_1,n_VarType_1,
     +  maxMag_1,l_minMag_1,minMag_1,med_Hp,A,e_A,morePhoto_2,Notes_8,
     +  VarName_1,period_3,Ep_2440000_1,max_1,min_1,Band_1,refFlag_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'solar_ha.dat'	! Solar System Annex: Astrometric catalogue

C  Format for file interpretation

   10 format(
     +  I3,1X,F11.7,1X,F11.7,1X,F13.7,1X,F5.2,1X,F7.3,1X,F6.2,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'solar_ha.dat')
      write(6,*) '....Loading file: solar_ha.dat'
      do i__=1,5609
        read(1,'(A64)')ar__9
        read(ar__9,10)
     +  ID,RAdeg_2,DEdeg_2,Ep_2440000_2,delay,theta_2,e_lambda,flag_2
c    ..............Just test output...........
        write(6,10)
     +  ID,RAdeg_2,DEdeg_2,Ep_2440000_2,delay,theta_2,e_lambda,flag_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'solar_hp.dat'	! Solar System Annex: Photometric catalogue

C  Format for file interpretation

   11 format(
     +  I3,1X,F11.5,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F5.3,1X,F5.3,
     +  1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'solar_hp.dat')
      write(6,*) '....Loading file: solar_hp.dat'
      do i__=1,2639
        read(1,'(A63)')ar__10
        read(ar__10,11)
     +  ID_1,Ep_2440000_3,Hp_dc,e_Hp_dc,Hp_ac,e_Hp_ac,r,Delta,alpha
c    ..............Just test output...........
        write(6,11)
     +  ID_1,Ep_2440000_3,Hp_dc,e_Hp_dc,Hp_ac,e_Hp_ac,r,Delta,alpha
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'solar_t.dat'	! Solar System Annex: Tycho astrometry/photometry

C  Format for file interpretation

   12 format(
     +  I3,1X,F13.7,1X,F11.7,1X,F11.7,1X,F5.2,1X,F5.2,1X,I1,1X,F5.1,
     +  1X,F5.1,1X,F5.2,1X,F6.2,1X,I2,1X,F5.1,1X,F5.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'solar_t.dat')
      write(6,*) '....Loading file: solar_t.dat'
      do i__=1,291
        read(1,'(A95)')ar__11
        read(ar__11,12)
     +  ID_2,Ep_2440000_4,RAdeg_3,DEdeg_3,BTmag_1,VTmag_1,Flag_3,
     +  e_RAdeg_2,e_DEdeg_2,DE_RA_1,theta_3,sign_z,e_incl,e_vert
        if(ar__11(43:47) .EQ. '') BTmag_1 = rNULL__
        if(ar__11(49:53) .EQ. '') VTmag_1 = rNULL__
c    ..............Just test output...........
        write(6,12)
     +  ID_2,Ep_2440000_4,RAdeg_3,DEdeg_3,BTmag_1,VTmag_1,Flag_3,
     +  e_RAdeg_2,e_DEdeg_2,DE_RA_1,theta_3,sign_z,e_incl,e_vert
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hd_notes.doc'	! Hipparcos notes: Double and multiple systems

C  Format for file interpretation

   13 format(I6,1X,A1,1X,A1,1X,I2,1X,I2,1X,A80)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hd_notes.doc')
      write(6,*) '....Loading file: hd_notes.doc'
      do i__=1,2622
        read(1,'(A97)')ar__12
        read(ar__12,13)HIP_8,Note1,Note2,Ntot,Nline,Text
c    ..............Just test output...........
        write(6,13)HIP_8,Note1,Note2,Ntot,Nline,Text
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hg_notes.doc'	! Hipparcos notes: General notes

C  Format for file interpretation

   14 format(I6,1X,A1,1X,A1,1X,I2,1X,I2,1X,A80)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hg_notes.doc')
      write(6,*) '....Loading file: hg_notes.doc'
      do i__=1,3898
        read(1,'(A97)')ar__13
        read(ar__13,14)HIP_9,Note1_1,Note2_1,Ntot_1,Nline_1,Text_1
c    ..............Just test output...........
        write(6,14)HIP_9,Note1_1,Note2_1,Ntot_1,Nline_1,Text_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hp_notes.doc'	! Hipparcos notes: Photometric notes

C  Format for file interpretation

   15 format(I6,1X,A1,1X,A1,1X,I2,1X,I2,1X,A80)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hp_notes.doc')
      write(6,*) '....Loading file: hp_notes.doc'
      do i__=1,2444
        read(1,'(A97)')ar__14
        read(ar__14,15)HIP_10,Note1_2,Note2_2,Ntot_2,Nline_2,Text_2
c    ..............Just test output...........
        write(6,15)HIP_10,Note1_2,Note2_2,Ntot_2,Nline_2,Text_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hp_refs.doc'	! References Hipparcos stars

C  Format for file interpretation

   16 format(I6,1X,I2,1X,I2,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hp_refs.doc')
      write(6,*) '....Loading file: hp_refs.doc'
      do i__=1,33769
        read(1,'(A19)')ar__15
        read(ar__15,16)HIP_11,Ntot_3,Nline_3,nRef
c    ..............Just test output...........
        write(6,16)HIP_11,Ntot_3,Nline_3,nRef
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'hp_auth.doc'	! References of hp_notes.doc

C  Format for file interpretation

   17 format(F6.3,1X,A70)

C  Effective file loading

      open(unit=1,status='old',file=
     +'hp_auth.doc')
      write(6,*) '....Loading file: hp_auth.doc'
      do i__=1,4335
        read(1,'(A80)')ar__16
        read(ar__16,17)nRef_1,Text_3
c    ..............Just test output...........
        write(6,17)nRef_1,Text_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'dmsa_o.doc'	! References of hip_dm_o.dat

C  Format for file interpretation

   18 format(I3,1X,I1,1X,I1,1X,A72)

C  Effective file loading

      open(unit=1,status='old',file=
     +'dmsa_o.doc')
      write(6,*) '....Loading file: dmsa_o.doc'
      do i__=1,118
        read(1,'(A80)')ar__17
        read(ar__17,18)dmRef_1,Ntot_4,Nline_4,Text_4
c    ..............Just test output...........
        write(6,18)dmRef_1,Ntot_4,Nline_4,Text_4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tyc_main.dat'	! The main part of Tycho Catalogue

C  Format for file interpretation

   19 format(
     +  A1,1X,A12,1X,A1,1X,A11,1X,A11,1X,F5.2,3X,A1,1X,F12.8,1X,F12.8,
     +  1X,A1,1X,F7.2,1X,F8.2,1X,F8.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
     +  1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
     +  F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F5.2,1X,I6,1X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F5.3,3X,I1,1X,F4.1,1X,
     +  A1,1X,I3,1X,F5.3,1X,F5.2,1X,F5.2,1X,A1,1X,A1,1X,A1,1X,A1,1X,
     +  A2,1X,I6,1X,I6,1X,A10,1X,A10,1X,A10,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tyc_main.dat')
      write(6,*) '....Loading file: tyc_main.dat'
      do i__=1,1058332
        read(1,'(A350)')ar__18
        read(ar__18,19)
     +  Catalog_1,TYC,Proxy_1,RAhms_1,DEdms_1,Vmag_1,r_Vmag_1,RAdeg_4,
     +  DEdeg_4,AstroRef_1,Plx_2,pmRA_2,pmDE_2,e_RAdeg_3,e_DEdeg_3,
     +  e_Plx_2,e_pmRA_2,e_pmDE_2,DE_RA_2,Plx_RA_1,Plx_DE_1,pmRA_RA_1,
     +  pmRA_DE_1,pmRA_Plx_1,pmDE_RA_1,pmDE_DE_1,pmDE_Plx_1,
     +  pmDE_pmRA_1,Nastro,F2_1,HIP_12,BTmag_2,e_BTmag_1,VTmag_2,
     +  e_VTmag_1,r_BTmag,B_V_1,e_B_V_1,Q,Fs,Source_3,Nphoto,VTscat,
     +  VTmax,VTmin,Var,VarFlag_1,MultFlag_1,morePhoto_3,m_HIP_1,PPM,
     +  HD_1,BD_1,CoD_1,CPD_1,Remark
        if(ar__18(42:46) .EQ. '') Vmag_1 = rNULL__
        if(ar__18(80:86) .EQ. '') Plx_2 = rNULL__
        if(ar__18(88:95) .EQ. '') pmRA_2 = rNULL__
        if(ar__18(97:104) .EQ. '') pmDE_2 = rNULL__
        if(ar__18(106:111) .EQ. '') e_RAdeg_3 = rNULL__
        if(ar__18(113:118) .EQ. '') e_DEdeg_3 = rNULL__
        if(ar__18(120:125) .EQ. '') e_Plx_2 = rNULL__
        if(ar__18(127:132) .EQ. '') e_pmRA_2 = rNULL__
        if(ar__18(134:139) .EQ. '') e_pmDE_2 = rNULL__
        if(ar__18(141:145) .EQ. '') DE_RA_2 = rNULL__
        if(ar__18(147:151) .EQ. '') Plx_RA_1 = rNULL__
        if(ar__18(153:157) .EQ. '') Plx_DE_1 = rNULL__
        if(ar__18(159:163) .EQ. '') pmRA_RA_1 = rNULL__
        if(ar__18(165:169) .EQ. '') pmRA_DE_1 = rNULL__
        if(ar__18(171:175) .EQ. '') pmRA_Plx_1 = rNULL__
        if(ar__18(177:181) .EQ. '') pmDE_RA_1 = rNULL__
        if(ar__18(183:187) .EQ. '') pmDE_DE_1 = rNULL__
        if(ar__18(189:193) .EQ. '') pmDE_Plx_1 = rNULL__
        if(ar__18(195:199) .EQ. '') pmDE_pmRA_1 = rNULL__
        if(ar__18(201:203) .EQ. '') Nastro = iNULL__
        if(ar__18(205:209) .EQ. '') F2_1 = rNULL__
        if(ar__18(211:216) .EQ. '') HIP_12 = iNULL__
        if(ar__18(218:223) .EQ. '') BTmag_2 = rNULL__
        if(ar__18(225:229) .EQ. '') e_BTmag_1 = rNULL__
        if(ar__18(231:236) .EQ. '') VTmag_2 = rNULL__
        if(ar__18(238:242) .EQ. '') e_VTmag_1 = rNULL__
        if(ar__18(246:251) .EQ. '') B_V_1 = rNULL__
        if(ar__18(253:257) .EQ. '') e_B_V_1 = rNULL__
        if(ar__18(261:261) .EQ. '') Q = iNULL__
        if(ar__18(263:266) .EQ. '') Fs = rNULL__
        if(ar__18(270:272) .EQ. '') Nphoto = iNULL__
        if(ar__18(274:278) .EQ. '') VTscat = rNULL__
        if(ar__18(280:284) .EQ. '') VTmax = rNULL__
        if(ar__18(286:290) .EQ. '') VTmin = rNULL__
        if(ar__18(303:308) .EQ. '') PPM = iNULL__
        if(ar__18(310:315) .EQ. '') HD_1 = iNULL__
c    ..............Just test output...........
        write(6,19)
     +  Catalog_1,TYC,Proxy_1,RAhms_1,DEdms_1,Vmag_1,r_Vmag_1,RAdeg_4,
     +  DEdeg_4,AstroRef_1,Plx_2,pmRA_2,pmDE_2,e_RAdeg_3,e_DEdeg_3,
     +  e_Plx_2,e_pmRA_2,e_pmDE_2,DE_RA_2,Plx_RA_1,Plx_DE_1,pmRA_RA_1,
     +  pmRA_DE_1,pmRA_Plx_1,pmDE_RA_1,pmDE_DE_1,pmDE_Plx_1,
     +  pmDE_pmRA_1,Nastro,F2_1,HIP_12,BTmag_2,e_BTmag_1,VTmag_2,
     +  e_VTmag_1,r_BTmag,B_V_1,e_B_V_1,Q,Fs,Source_3,Nphoto,VTscat,
     +  VTmax,VTmin,Var,VarFlag_1,MultFlag_1,morePhoto_3,m_HIP_1,PPM,
     +  HD_1,BD_1,CoD_1,CPD_1,Remark
c    .......End.of.Just test output...........
      end do
      close(1)

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