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

Conversion of standardized ReadMe file for file /./ftp/cats/I/239 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 2018-Apr-19
C=============================================================================
*  This code was generated from the ReadMe file documenting a catalogue
*  according to the "Standard for Documentation of Astronomical Catalogues"
*  currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
*  (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
*  Please report problems or questions to   
C=============================================================================

      implicit none
*  Unspecified or NULL values, generally corresponding to blank columns,
*  are assigned one of the following special values:
*     rNULL__    for unknown or NULL floating-point values
*     iNULL__    for unknown or NULL   integer      values
      real*4     rNULL__
      integer*4  iNULL__
      parameter  (rNULL__=--2147483648.)  	! NULL real number
      parameter  (iNULL__=(-2147483647-1))	! NULL int  number
      integer    idig			! testing NULL number

C=============================================================================
Cat. 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    (nr__) ! [H] Catalogue (H=Hipparcos)               (H0)
      integer*4     HIP        (nr__) ! Identifier (HIP number)                   (H1)
      character*1   Proxy      (nr__) ! *[HT] Proximity flag                       (H2)
      character*11  RAhms      (nr__) ! Right ascension in h m s, ICRS (J1991.25) (H3)
      character*11  DEdms      (nr__) ! Declination in deg ' ", ICRS (J1991.25)   (H4)
      real*4        Vmag       (nr__) ! (mag) ? Magnitude in Johnson V                  (H5)
      integer*4     VarFlag    (nr__) ! *[1,3]? Coarse variability flag            (H6)
      character*1   r_Vmag     (nr__) ! *[GHT] Source of magnitude                 (H7)
      real*8        RAdeg      (nr__) ! (deg) *? alpha, degrees (ICRS, Epoch=J1991.25)   (H8)
      real*8        DEdeg      (nr__) ! (deg) *? delta, degrees (ICRS, Epoch=J1991.25)   (H9)
      character*1   AstroRef   (nr__) ! *[*+A-Z] Reference flag for astrometry    (H10)
      real*8        Plx        (nr__) ! (mas) ? Trigonometric parallax                 (H11)
      real*8        pmRA       (nr__) ! (mas/yr) *? Proper motion mu_alpha.cos(delta), ICRS(H12)
      real*8        pmDE       (nr__) ! (mas/yr) *? Proper motion mu_delta, ICRS           (H13)
      real*4        e_RAdeg    (nr__) ! (mas) *? Standard error in RA*cos(DEdeg)        (H14)
      real*4        e_DEdeg    (nr__) ! (mas) *? Standard error in DE                   (H15)
      real*4        e_Plx      (nr__) ! (mas) ? Standard error in Plx                  (H16)
      real*4        e_pmRA     (nr__) ! (mas/yr) ? Standard error in pmRA                 (H17)
      real*4        e_pmDE     (nr__) ! (mas/yr) ? Standard error in pmDE                 (H18)
      real*4        DE_RA      (nr__) ! [-1/1]? Correlation, DE/RA*cos(delta)    (H19)
      real*4        Plx_RA     (nr__) ! [-1/1]? Correlation, Plx/RA*cos(delta)   (H20)
      real*4        Plx_DE     (nr__) ! [-1/1]? Correlation, Plx/DE              (H21)
      real*4        pmRA_RA    (nr__) ! [-1/1]? Correlation, pmRA/RA*cos(delta)  (H22)
      real*4        pmRA_DE    (nr__) ! [-1/1]? Correlation, pmRA/DE             (H23)
      real*4        pmRA_Plx   (nr__) ! [-1/1]? Correlation, pmRA/Plx            (H24)
      real*4        pmDE_RA    (nr__) ! [-1/1]? Correlation, pmDE/RA*cos(delta)  (H25)
      real*4        pmDE_DE    (nr__) ! [-1/1]? Correlation, pmDE/DE             (H26)
      real*4        pmDE_Plx   (nr__) ! [-1/1]? Correlation, pmDE/Plx            (H27)
      real*4        pmDE_pmRA  (nr__) ! [-1/1]? Correlation, pmDE/pmRA           (H28)
      integer*4     F1         (nr__) ! (%) ? Percentage of rejected data            (H29)
      real*4        F2         (nr__) ! *? Goodness-of-fit parameter              (H30)
      real*4        BTmag      (nr__) ! (mag) ? Mean BT magnitude                      (H32)
      real*4        e_BTmag    (nr__) ! (mag) ? Standard error on BTmag                (H33)
      real*4        VTmag      (nr__) ! (mag) ? Mean VT magnitude                      (H34)
      real*4        e_VTmag    (nr__) ! (mag) ? Standard error on VTmag                (H35)
      character*1   m_BTmag    (nr__) ! *[A-Z*-] Reference flag for BT and VTmag  (H36)
      real*4        B_V        (nr__) ! (mag) ? Johnson B-V colour                     (H37)
      real*4        e_B_V      (nr__) ! (mag) ? Standard error on B-V                  (H38)
      character*1   r_B_V      (nr__) ! [GT] Source of B-V from Ground or Tycho  (H39)
      real*4        V_I        (nr__) ! (mag) ? Colour index in Cousins' system        (H40)
      real*4        e_V_I      (nr__) ! (mag) ? Standard error on V-I                  (H41)
      character*1   r_V_I      (nr__) ! *[A-T] Source of V-I                      (H42)
      character*1   CombMag    (nr__) ! [*] Flag for combined Vmag, B-V, V-I     (H43)
      real*8        Hpmag      (nr__) ! (mag) *? Median magnitude in Hipparcos system   (H44)
      real*4        e_Hpmag    (nr__) ! (mag) *? Standard error on Hpmag                (H45)
      real*4        Hpscat     (nr__) ! (mag) ? Scatter on Hpmag                       (H46)
      integer*4     o_Hpmag    (nr__) ! ? Number of observations for Hpmag       (H47)
      character*1   m_Hpmag    (nr__) ! *[A-Z*-] Reference flag for Hpmag         (H48)
      real*4        Hpmax      (nr__) ! (mag) ? Hpmag at maximum (5th percentile)      (H49)
      real*4        HPmin      (nr__) ! (mag) ? Hpmag at minimum (95th percentile)     (H50)
      real*8        Period     (nr__) ! (d) ? Variability period (days)              (H51)
      character*1   HvarType   (nr__) ! *[CDMPRU]? variability type               (H52)
      character*1   moreVar    (nr__) ! *[12] Additional data about variability   (H53)
      character*1   morePhoto  (nr__) ! [ABC] Light curve Annex                  (H54)
      character*10  CCDM       (nr__) ! CCDM identifier                          (H55)
      character*1   n_CCDM     (nr__) ! *[HIM] Historical status flag             (H56)
      integer*4     Nsys       (nr__) ! ? Number of entries with same CCDM       (H57)
      integer*4     Ncomp      (nr__) ! ? Number of components in this entry     (H58)
      character*1   MultFlag   (nr__) ! *[CGOVX] Double/Multiple Systems flag     (H59)
      character*1   Source     (nr__) ! *[PFILS] Astrometric source flag          (H60)
      character*1   Qual       (nr__) ! *[ABCDS] Solution quality                 (H61)
      character*2   m_HIP      (nr__) ! Component identifiers                    (H62)
      integer*4     theta      (nr__) ! (deg) ? Position angle between components      (H63)
      real*8        rho        (nr__) ! (arcsec) ? Angular separation between components  (H64)
      real*4        e_rho      (nr__) ! (arcsec) ? Standard error on rho                  (H65)
      real*4        dHp        (nr__) ! (mag) ? Magnitude difference of components     (H66)
      real*4        e_dHp      (nr__) ! (mag) ? Standard error on dHp                  (H67)
      character*1   Survey     (nr__) ! [S] Flag indicating a Survey Star        (H68)
      character*1   Chart      (nr__) ! *[DG] Identification Chart                (H69)
      character*1   Notes      (nr__) ! *[DGPWXYZ] Existence of notes             (H70)
      integer*4     HD         (nr__) ! [1/359083]? HD number <III/135>          (H71)
      character*10  BD         (nr__) ! Bonner DM <I/119>, <I/122>               (H72)
      character*10  CoD        (nr__) ! Cordoba Durchmusterung (DM) <I/114>      (H73)
      character*10  CPD        (nr__) ! Cape Photographic DM <I/108>             (H74)
      real*4        v_V_I_red  (nr__) ! (mag) V-I used for reductions                  (H75)
      character*12  SpType     (nr__) ! Spectral type                            (H76)
      character*1   r_SpType   (nr__) ! *[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     (nr__1) ! CCDM number                              (DC1)
      integer*4     S          (nr__1) ! *Solution identifier                      (DC2)
      character*1   Type       (nr__1) ! *[FIL]Type of solution                    (DC3)
      character*1   Source_1   (nr__1) ! *[CFN] Solution source                    (DC4)
      character*1   Qual_1     (nr__1) ! *[ABCD] Solution quality                  (DC5)
      character*1   Notes_1    (nr__1) ! *[DGPWXYZ] Existence of notes             (DC6)
      integer*4     Nsys_1     (nr__1) ! Number of solutions for the system      (DCM1)
      integer*4     Ncomp_1    (nr__1) ! Number of components in this solution   (DCM2)
      integer*4     Nparm      (nr__1) ! Number of free parameters in solution   (DCM3)
      integer*4     Ncorr      (nr__1) ! Number of correlation records           (DCM4)
      integer*4     seq        (nr__1) ! Sequential component number             (DCM6)
      character*1   comp_id    (nr__1) ! Component identifier                     (DC7)
      integer*4     HIP_1      (nr__1) ! HIP number                               (DC8)
      real*4        Hp         (nr__1) ! (mag) Magnitude of component                   (DC9)
      real*4        e_Hp       (nr__1) ! (mag) Standard error of Hp magnitude          (DC10)
      real*4        BT         (nr__1) ! (mag) ? Magnitude of component, BT            (DC11)
      real*4        e_BT       (nr__1) ! (mag) ? Standard error of BT                  (DC12)
      real*4        VT         (nr__1) ! (mag) ? Magnitude of component, VT            (DC13)
      real*4        e_VT       (nr__1) ! (mag) ? Standard error of VT                  (DC14)
      real*8        RAdeg_1    (nr__1) ! (deg) alpha, degrees (ICRS, Epoch=J1991.25)   (DC15)
      real*8        DEdeg_1    (nr__1) ! (deg) delta, degrees (ICRS, Epoch=J1991.25)   (DC16)
      real*8        Plx_1      (nr__1) ! (mas) Trigonometric parallax                  (DC17)
      real*8        pmRA_1     (nr__1) ! (mas/yr) Proper motion in mu_alpha.cos(delta)ICRS(DC18)
      real*8        pmDE_1     (nr__1) ! (mas/yr) Proper motion in mu_delta in ICRS       (DC19)
      real*4        e_RAdeg_1  (nr__1) ! (mas) Standard error in RA*cos(DEdeg)         (DC20)
      real*4        e_DEdeg_1  (nr__1) ! (mas) Standard error in DE                    (DC21)
      real*4        e_Plx_1    (nr__1) ! (mas) Standard error in Plx                   (DC22)
      real*4        e_pmRA_1   (nr__1) ! (mas/yr) Standard error in pmRA                  (DC23)
      real*4        e_pmDE_1   (nr__1) ! (mas/yr) Standard error in pmDE                  (DC24)
      character*1   ref        (nr__1) ! Reference component for following data  (DC25)
      real*8        theta_1    (nr__1) ! (deg) *? Position angle                        (DC26)
      real*8        rho_1      (nr__1) ! (arcsec) ? Separation from reference component   (DC27)
      real*8        d_theta    (nr__1) ! (deg/yr) ? Rate of change of theta               (DC28)
      real*4        d_rho      (nr__1) ! (arcsec/yr) ? Rate of change of separation          (DC29)
      integer*4     seq_ref    (nr__1) ! *Sequential record number                (DCM7)
      character*25  flag       (nr__1) ! *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     (nr__2) ! CCDM number                               (DC1)
      integer*4     S_1        (nr__2) ! *Solution identifier                       (DC2)
      character*1   Type_1     (nr__2) ! *[FIL]Type of solution                     (DC3)
      character*1   Source_2   (nr__2) ! *[CFN] Solution source                     (DC4)
      character*1   Qual_2     (nr__2) ! *[ABCD] Solution quality                   (DC5)
      character*1   Notes_2    (nr__2) ! *[DGPWXYZ] Existence of notes              (DC6)
      integer*4     Nsys_2     (nr__2) ! Number of solutions for the system       (DCM1)
      integer*4     Ncomp_2    (nr__2) ! Number of components in this solution    (DCM2)
      integer*4     Nparm_1    (nr__2) ! Number of free parameters in solution    (DCM3)
      integer*4     Ncorr_1    (nr__2) ! Number of correlation records            (DCM4)
      integer*4     seq_1      (nr__2) ! Sequential component number              (DCM6)
      integer*4     corr       (66,nr__2) ! *[-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      (nr__3) ! Identifier (HIP number)                   (DG1)
      real*8        gRA        (nr__3) ! (mas/yr2) *Acceleration gRA = d(pmRA)/dt             (DG2)
      real*8        gDE        (nr__3) ! (mas/yr2) *Acceleration gDE = d(pmDE)/dt             (DG3)
      real*8        e_gRA      (nr__3) ! (mas/yr2) Standard error of gRA                     (DG4)
      real*8        e_gDE      (nr__3) ! (mas/yr2) Standard error of gDE                     (DG5)
      real*4        Fg         (nr__3) ! *Significance of the g terms               (DG6)
      real*8        dgRA       (nr__3) ! (mas/yr3) *? dgRA = d2(pmRA)/dt2                     (DG7)
      real*8        dgDE       (nr__3) ! (mas/yr3) *? dgDE = d2(pmDE)/dt2                     (DG8)
      real*8        e_dgRA     (nr__3) ! (mas/yr3) ? Standard error of dgRA                  (DG9)
      real*8        e_dgDE     (nr__3) ! (mas/yr3) ? Standard error of dgDE                 (DG10)
      real*4        Fdg        (nr__3) ! *? Signifance of the dg terms             (DG11)
      character*1   Notes_3    (nr__3) ! *[DGPWXYZ] Existence of notes             (DG12)
      integer*4     num        (nr__3) ! *[7,9] Number of astrometric parameters   (DGM1)
      integer*4     corr_1     (36,nr__3) ! *[-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      (nr__4) ! Identifier (HIP)                         (D01)
      real*8        P          (nr__4) ! (d) Orbital period                           (DO2)
      real*8        T          (nr__4) ! (d) *Time of periastron passage               (DO3)
      real*8        a0         (nr__4) ! (mas) Semi-major axis of photocentric orbit    (DO4)
      real*4        ecc        (nr__4) ! [0,1] Eccentricity                       (DO5)
      real*4        w          (nr__4) ! (deg) *[0,360] Argument of periastron           (DO6)
      real*4        i          (nr__4) ! (deg) *[0,180] Inclination                      (DO7)
      real*4        Omega      (nr__4) ! (deg) *[0,360] Position angle of the node       (DO8)
      real*8        e_P        (nr__4) ! (d) ? Standard error of P                    (DO9)
      real*8        e_T        (nr__4) ! (d) ? Standard error of T                   (DO10)
      real*4        e_a0       (nr__4) ! (mas) ? Standard error of a0                  (DO11)
      real*4        e_ecc      (nr__4) ! ? Standard error of ecc                 (DO12)
      real*4        e_w        (nr__4) ! (deg) ? Standard error of w                   (DO13)
      real*4        e_i        (nr__4) ! (deg) ? Standard error of i                   (DO14)
      real*4        e_Omega    (nr__4) ! (deg) ? Standard error of Omega               (DO15)
      integer*4     dmRef      (nr__4) ! *? reference to the literature           (DO16)
      character*1   Notes_4    (nr__4) ! *[DGPWXYZ] Existence of notes            (DO17)
      character*12  flag_1     (nr__4) ! *Status flags for the parameters         (DOM1)
      integer*4     corr_2     (66,nr__4) ! *[-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      (nr__5) ! Identifier (HIP)                         (DV1)
      real*4        Hp_ref     (nr__5) ! (mag) *Reference magnitude                      (DV2)
      real*8        DRA        (nr__5) ! (mas) VIM element in right ascension           (DV3)
      real*8        DDE        (nr__5) ! (mas) VIM element in declination               (DV4)
      real*8        e_DRA      (nr__5) ! (mas) Standard error of DRA                    (DV5)
      real*8        e_DDE      (nr__5) ! (mas) Standard error of DDE                    (DV6)
      real*4        FD         (nr__5) ! *Significance of the VIM elements         (DV7)
      real*4        theta_C    (nr__5) ! (deg) *Position angle of the constant component (DV8)
      real*4        e_theta_C  (nr__5) ! (deg) Standard error of theta_C                (DV9)
      real*4        minSep     (nr__5) ! (mas) Lower limit for separation of binary    (DV10)
      real*4        dvar       (nr__5) ! (mas) *Displacement of photocentre             (DV11)
      character*1   Notes_5    (nr__5) ! *[DGPWXYZ] Existence of notes            (DV12)
      integer*4     corr_3     (21,nr__5) ! *[-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      (nr__6) ! Identifier (HIP)                         (DX1)
      real*4        epsilon    (nr__6) ! (mas) Cosmic error, epsilon                    (DX2)
      real*4        e_epsilon  (nr__6) ! (mas) Standard error of epsilon                (DX3)
      character*1   Notes_6    (nr__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      (nr__7) ! Identifier (HIP)                          (P1)
      character*1   NewVar     (nr__7) ! [*] Flag if new variable                  (P2)
      character*12  SpType_1   (nr__7) ! Spectral type                             (P3)
      character*1   HvarType_1 (nr__7) ! *[CDMPRU] Variability type (1-letter)      (P4)
      character*6   VarType    (nr__7) ! *Variability type as in GCVS/NSV           (P5)
      character*1   n_VarType  (nr__7) ! [*] Flag if newly classified by Hipparcos (P6)
      real*4        maxMag     (nr__7) ! (mag) Magnitude at max from curve fitting       (P7)
      character*1   l_minMag   (nr__7) ! *[>] Limit flag (>)                        (P8)
      real*4        minMag     (nr__7) ! (mag) Magnitude at min from curve fitting       (P9)
      real*4        log_sA_A   (nr__7) ! ? log_10(sigma_A / A)                    (P10)
      real*8        Period_1   (nr__7) ! (d) ? Mean period in days                    (P11)
      real*4        log_sP     (nr__7) ! ([d]) ? log_10(sigma_P)                        (P12)
      real*8        Ep_2440000 (nr__7) ! (d) ? Epoch (JD-2440000) of zero phase       (P13)
      integer*4     q_Ep_2440000(nr__7) ! *[0/5]? Precision flag                    (P14)
      character*1   morePhoto_1(nr__7) ! [ABC] Light curve Annex                  (P15)
      character*1   Notes_7    (nr__7) ! *[DGPWXYZ] Existence of notes             (P16)
      character*12  VarName    (nr__7) ! Variable star name                       (P17)
      real*8        period_2   (nr__7) ! (d) ? Period from literature                 (P18)
      real*8        epoch      (nr__7) ! (d) ? Epoch  from literature                 (P19)
      real*4        max        (nr__7) ! (mag) ? Magnitude at max from literature       (P20)
      real*4        min        (nr__7) ! (mag) ? Magnitude at min from literature       (P21)
      character*1   Band       (nr__7) ! *[UBVKIRPYb] Photometric band             (P22)
      character*1   refFlag    (nr__7) ! [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      (nr__8) ! Identifier (HIP number)                   (U1)
      character*1   NewVar_1   (nr__8) ! [*] Flag if new variable                  (U2)
      character*12  SpType_2   (nr__8) ! Spectral type                             (U3)
      character*1   HvarType_2 (nr__8) ! *[CDMPRU] Variability type (1-letter)      (U4)
      character*6   VarType_1  (nr__8) ! *Variability type (6-letter) as in GCVS/NSV(U5)
      character*1   n_VarType_1(nr__8) ! [*] Flag if newly classified by Hipparcos (U6)
      real*4        maxMag_1   (nr__8) ! (mag) *Magnitude at max from curve fitting       (U7)
      character*1   l_minMag_1 (nr__8) ! [>] Limit flag (>)                        (U8)
      real*4        minMag_1   (nr__8) ! (mag) *Magnitude at min from curve fitting       (U9)
      real*4        med_Hp     (nr__8) ! (mag) Median Hp                                (U10)
      real*8        A          (nr__8) ! (mag) Intrinsic variability amplitude          (U11)
      real*4        e_A        (nr__8) ! (mag) Standard error of A                      (U12)
      character*1   morePhoto_2(nr__8) ! [ABC] Light curve Annex                  (U15)
      character*1   Notes_8    (nr__8) ! *[DGPWXYZ] Notes                          (U16)
      character*12  VarName_1  (nr__8) ! Variable star name                       (U17)
      real*8        period_3   (nr__8) ! (d) ? Period from literature                 (U18)
      real*8        Ep_2440000_1(nr__8) ! (d) ? Epoch (JD-2440000) from literature     (U19)
      real*4        max_1      (nr__8) ! (mag) ? Magnitude at max from literature       (U20)
      real*4        min_1      (nr__8) ! (mag) ? Magnitude at min from literature       (U21)
      character*1   Band_1     (nr__8) ! *[UBVKIRPYb] Photometric band             (U22)
      character*1   refFlag_1  (nr__8) ! [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         (nr__9) ! Object number                         (SHA1)
      real*8        RAdeg_2    (nr__9) ! (deg) Reference point RA  (ICRS system)     (SHA2)
      real*8        DEdeg_2    (nr__9) ! (deg) Reference point Dec (ICRS system)     (SHA3)
      real*8        Ep_2440000_2(nr__9) ! (d) *Measurement epoch                     (SHA4)
      real*4        delay      (nr__9) ! (s) *Light delay time                      (SHA5)
      real*8        theta_2    (nr__9) ! (deg) *Position angle, theta                 (SHA6)
      real*4        e_lambda   (nr__9) ! (mas) Estimated standard error of abscissa  (SHA7)
      integer*4     flag_2     (nr__9) ! *[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       (nr__10) ! Object number                         (SHP1)
      real*8        Ep_2440000_3(nr__10) ! (d) *Measurement epoch                     (SHP2)
      real*8        Hp_dc      (nr__10) ! (mag) Magnitude from unmodulated signal     (SHP3)
      real*4        e_Hp_dc    (nr__10) ! (mag) Standard error on Hp_dc               (SHP4)
      real*8        Hp_ac      (nr__10) ! (mag) Magnitude from modulated signal       (SHP5)
      real*4        e_Hp_ac    (nr__10) ! (mag) Standard error on Hp_ac               (SHP6)
      real*4        r          (nr__10) ! (AU) Distance: Sun-asteroid                (SHP7)
      real*4        Delta      (nr__10) ! (AU) Distance: satellite-asteroid          (SHP8)
      real*4        alpha      (nr__10) ! (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       (nr__11) ! Object number                           (ST1)
      real*8        Ep_2440000_4(nr__11) ! (d) *Measurement epoch                       (ST2)
      real*8        RAdeg_3    (nr__11) ! (deg) Right ascension (ICRS system)           (ST3)
      real*8        DEdeg_3    (nr__11) ! (deg) Declination (ICRS system)               (ST4)
      real*4        BTmag_1    (nr__11) ! (mag) ? Mean BT magnitude                     (ST5)
      real*4        VTmag_1    (nr__11) ! (mag) ? Mean VT magnitude                     (ST6)
      integer*4     Flag_3     (nr__11) ! *[1,2] Transit flag                      (ST7)
      real*4        e_RAdeg_2  (nr__11) ! (mas) Standard error on RA                    (ST8)
      real*4        e_DEdeg_2  (nr__11) ! (mas) Standard error on Dec                   (ST9)
      real*4        DE_RA_1    (nr__11) ! Correlation, Dec/RA*cos(delta)         (ST10)
      real*4        theta_3    (nr__11) ! (deg) Position angle of slit (direction w)   (ST11)
      integer*4     sign_z     (nr__11) ! *[-1/1] Inclined slit flag              (ST12)
      real*4        e_incl     (nr__11) ! (mas) Standard error on slit pos. (inclined) (ST13)
      real*4        e_vert     (nr__11) ! (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      (nr__12) ! HIP number
      character*1   Note1      (nr__12) ! [DG] Double and multiple or General note
      character*1   Note2      (nr__12) ! [DP] Double and multiple or Photometric note
      integer*4     Ntot       (nr__12) ! Total number of lines for the HIP object
      integer*4     Nline      (nr__12) ! Running line number in range [1,Ntot]
      character*80  Text       (nr__12) ! 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      (nr__13) ! HIP number
      character*1   Note1_1    (nr__13) ! [DG] Double and multiple or General note
      character*1   Note2_1    (nr__13) ! [DP] Double and multiple or Photometric note
      integer*4     Ntot_1     (nr__13) ! Total number of lines for the HIP object
      integer*4     Nline_1    (nr__13) ! Running line number in range [1,Ntot]
      character*80  Text_1     (nr__13) ! 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     (nr__14) ! HIP number
      character*1   Note1_2    (nr__14) ! [DG] Double and multiple or General note
      character*1   Note2_2    (nr__14) ! [DP] Double and multiple or Photometric note
      integer*4     Ntot_2     (nr__14) ! Total number of lines for the HIP object
      integer*4     Nline_2    (nr__14) ! Running line number in range [1,Ntot]
      character*80  Text_2     (nr__14) ! 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     (nr__15) ! HIP number
      integer*4     Ntot_3     (nr__15) ! Total number of lines for the HIP object
      integer*4     Nline_3    (nr__15) ! Running line number in range [1,Ntot]
      real*4        nRef       (nr__15) ! 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     (nr__16) ! Reference number from hp_refs.doc
      character*70  Text_3     (nr__16) ! 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    (nr__17) ! Reference number (hip_dm_o.dat file)
      integer*4     Ntot_4     (nr__17) ! Total number of lines for the reference
      integer*4     Nline_4    (nr__17) ! Running line number in range [1,Ntot]
      character*72  Text_4     (nr__17) ! 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  (nr__18) ! [T] Catalogue (T = Tycho)                 (T0)
      character*12  TYC        (nr__18) ! *TYC1-3 (TYC number)                       (T1)
      character*1   Proxy_1    (nr__18) ! [HT]? Proximity flag                      (T2)
      character*11  RAhms_1    (nr__18) ! Right ascension in h m s, ICRS (J1991.25) (T3)
      character*11  DEdms_1    (nr__18) ! Declination in deg ' ", ICRS (J1991.25)   (T4)
      real*4        Vmag_1     (nr__18) ! (mag) ? Magnitude in Johnson V                  (T5)
      character*1   r_Vmag_1   (nr__18) ! *[BDTV] Source of magnitude                (T7)
      real*8        RAdeg_4    (nr__18) ! (deg) *alpha, degrees (ICRS, Epoch=J1991.25)     (T8)
      real*8        DEdeg_4    (nr__18) ! (deg) *delta, degrees (ICRS, Epoch=J1991.25)     (T9)
      character*1   AstroRef_1 (nr__18) ! *[X]? Reference flag for astrometry       (T10)
      real*8        Plx_2      (nr__18) ! (mas) *? Trigonometric parallax                 (T11)
      real*8        pmRA_2     (nr__18) ! (mas/yr) *? Proper motion mu_alpha.cos(delta), ICRS(T12)
      real*8        pmDE_2     (nr__18) ! (mas/yr) *? Proper motion mu_delta, ICRS           (T13)
      real*4        e_RAdeg_3  (nr__18) ! (mas) *? Standard error in RA*cos(delta)        (T14)
      real*4        e_DEdeg_3  (nr__18) ! (mas) *? Standard error in DE                   (T15)
      real*4        e_Plx_2    (nr__18) ! (mas) *? Standard error in Plx                  (T16)
      real*4        e_pmRA_2   (nr__18) ! (mas/yr) *? Standard error in pmRA                 (T17)
      real*4        e_pmDE_2   (nr__18) ! (mas/yr) *? Standard error in pmDE                 (T18)
      real*4        DE_RA_2    (nr__18) ! [-1/1]? Correlation, DE/RA*cos(delta)    (T19)
      real*4        Plx_RA_1   (nr__18) ! [-1/1]? Correlation, Plx/RA*cos(delta)   (T20)
      real*4        Plx_DE_1   (nr__18) ! [-1/1]? Correlation, Plx/DE              (T21)
      real*4        pmRA_RA_1  (nr__18) ! [-1/1]? Correlation, pmRA/RA*cos(delta)  (T22)
      real*4        pmRA_DE_1  (nr__18) ! [-1/1]? Correlation, pmRA/DE             (T23)
      real*4        pmRA_Plx_1 (nr__18) ! [-1/1]? Correlation, pmRA/Plx            (T24)
      real*4        pmDE_RA_1  (nr__18) ! [-1/1]? Correlation, pmDE/RA*cos(delta)  (T25)
      real*4        pmDE_DE_1  (nr__18) ! [-1/1]? Correlation, pmDE/DE             (T26)
      real*4        pmDE_Plx_1 (nr__18) ! [-1/1]? Correlation, pmDE/Plx            (T27)
      real*4        pmDE_pmRA_1(nr__18) ! [-1/1]? Correlation, pmDE/pmRA           (T28)
      integer*4     Nastro     (nr__18) ! ? Number of transits for astrometry      (T29)
      real*4        F2_1       (nr__18) ! *? Goodness-of-fit parameter              (T30)
      integer*4     HIP_12     (nr__18) ! ? Hipparcos HIP number                   (T31)
      real*4        BTmag_2    (nr__18) ! (mag) ? Mean BT magnitude                      (T32)
      real*4        e_BTmag_1  (nr__18) ! (mag) ? Standard error in BTmag                (T33)
      real*4        VTmag_2    (nr__18) ! (mag) ? Mean VT magnitude                      (T34)
      real*4        e_VTmag_1  (nr__18) ! (mag) ? Standard error in VTmag                (T35)
      character*1   r_BTmag    (nr__18) ! *[DMNT] Source of photometry              (T36)
      real*4        B_V_1      (nr__18) ! (mag) ? Johnson B-V colour                     (T37)
      real*4        e_B_V_1    (nr__18) ! (mag) ? Standard error on B-V                  (T38)
      integer*4     Q          (nr__18) ! *? Astrometric quality flag, Q            (T40)
      real*4        Fs         (nr__18) ! ? Signal-to-noise ratio of the star image(T41)
      character*1   Source_3   (nr__18) ! *[HPR] Source of astrometric data         (T42)
      integer*4     Nphoto     (nr__18) ! ? Number of transits for photometry      (T43)
      real*4        VTscat     (nr__18) ! (mag) ? Estimate of VTmag scatter              (T44)
      real*4        VTmax      (nr__18) ! (mag) ? VTmag at maximum (15th percentile)     (T45)
      real*4        VTmin      (nr__18) ! (mag) ? VTmag at minimum (85th percentile)     (T46)
      character*1   Var        (nr__18) ! *[GN]? Known variability from GCVS/NSV    (T47)
      character*1   VarFlag_1  (nr__18) ! *[UVW]? Variability from Tycho            (T48)
      character*1   MultFlag_1 (nr__18) ! *[DRSYZ]? Duplicity from Tycho            (T49)
      character*1   morePhoto_3(nr__18) ! [AB]  Epoch photometry in Annex A or B   (T50)
      character*2   m_HIP_1    (nr__18) ! CCDM component identifier                (T51)
      integer*4     PPM        (nr__18) ! *[1/789676]? PPM and Supplement           (T52)
      integer*4     HD_1       (nr__18) ! [1/359083]? HD cat. <III/135>            (T53)
      character*10  BD_1       (nr__18) ! Bonner DM <I/119>, <I/122>               (T54)
      character*10  CoD_1      (nr__18) ! Cordoba DM <I/114>                       (T55)
      character*10  CPD_1      (nr__18) ! Cape Photographic DM <I/108>             (T56)
      character*1   Remark     (nr__18) ! *[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(i__),HIP(i__),Proxy(i__),RAhms(i__),DEdms(i__),
     +  Vmag(i__),VarFlag(i__),r_Vmag(i__),RAdeg(i__),DEdeg(i__),
     +  AstroRef(i__),Plx(i__),pmRA(i__),pmDE(i__),e_RAdeg(i__),
     +  e_DEdeg(i__),e_Plx(i__),e_pmRA(i__),e_pmDE(i__),DE_RA(i__),
     +  Plx_RA(i__),Plx_DE(i__),pmRA_RA(i__),pmRA_DE(i__),
     +  pmRA_Plx(i__),pmDE_RA(i__),pmDE_DE(i__),pmDE_Plx(i__),
     +  pmDE_pmRA(i__),F1(i__),F2(i__),BTmag(i__),e_BTmag(i__),
     +  VTmag(i__),e_VTmag(i__),m_BTmag(i__),B_V(i__),e_B_V(i__),
     +  r_B_V(i__),V_I(i__),e_V_I(i__),r_V_I(i__),CombMag(i__),
     +  Hpmag(i__),e_Hpmag(i__),Hpscat(i__),o_Hpmag(i__),m_Hpmag(i__),
     +  Hpmax(i__),HPmin(i__),Period(i__),HvarType(i__),moreVar(i__),
     +  morePhoto(i__),CCDM(i__),n_CCDM(i__),Nsys(i__),Ncomp(i__),
     +  MultFlag(i__),Source(i__),Qual(i__),m_HIP(i__),theta(i__),
     +  rho(i__),e_rho(i__),dHp(i__),e_dHp(i__),Survey(i__),
     +  Chart(i__),Notes(i__),HD(i__),BD(i__),CoD(i__),CPD(i__),
     +  v_V_I_red(i__),SpType(i__),r_SpType(i__)
        if(ar__(42:46) .EQ. '') Vmag(i__) = rNULL__
        if(ar__(48:48) .EQ. '') VarFlag(i__) = iNULL__
        if(ar__(52:63) .EQ. '') RAdeg(i__) = rNULL__
        if(ar__(65:76) .EQ. '') DEdeg(i__) = rNULL__
        if(ar__(80:86) .EQ. '') Plx(i__) = rNULL__
        if(ar__(88:95) .EQ. '') pmRA(i__) = rNULL__
        if(ar__(97:104) .EQ. '') pmDE(i__) = rNULL__
        if(ar__(106:111) .EQ. '') e_RAdeg(i__) = rNULL__
        if(ar__(113:118) .EQ. '') e_DEdeg(i__) = rNULL__
        if(ar__(120:125) .EQ. '') e_Plx(i__) = rNULL__
        if(ar__(127:132) .EQ. '') e_pmRA(i__) = rNULL__
        if(ar__(134:139) .EQ. '') e_pmDE(i__) = rNULL__
        if(ar__(141:145) .EQ. '') DE_RA(i__) = rNULL__
        if(ar__(147:151) .EQ. '') Plx_RA(i__) = rNULL__
        if(ar__(153:157) .EQ. '') Plx_DE(i__) = rNULL__
        if(ar__(159:163) .EQ. '') pmRA_RA(i__) = rNULL__
        if(ar__(165:169) .EQ. '') pmRA_DE(i__) = rNULL__
        if(ar__(171:175) .EQ. '') pmRA_Plx(i__) = rNULL__
        if(ar__(177:181) .EQ. '') pmDE_RA(i__) = rNULL__
        if(ar__(183:187) .EQ. '') pmDE_DE(i__) = rNULL__
        if(ar__(189:193) .EQ. '') pmDE_Plx(i__) = rNULL__
        if(ar__(195:199) .EQ. '') pmDE_pmRA(i__) = rNULL__
        if(ar__(201:203) .EQ. '') F1(i__) = iNULL__
        if(ar__(205:209) .EQ. '') F2(i__) = rNULL__
        if(ar__(218:223) .EQ. '') BTmag(i__) = rNULL__
        if(ar__(225:229) .EQ. '') e_BTmag(i__) = rNULL__
        if(ar__(231:236) .EQ. '') VTmag(i__) = rNULL__
        if(ar__(238:242) .EQ. '') e_VTmag(i__) = rNULL__
        if(ar__(246:251) .EQ. '') B_V(i__) = rNULL__
        if(ar__(253:257) .EQ. '') e_B_V(i__) = rNULL__
        if(ar__(261:264) .EQ. '') V_I(i__) = rNULL__
        if(ar__(266:269) .EQ. '') e_V_I(i__) = rNULL__
        if(ar__(275:281) .EQ. '') Hpmag(i__) = rNULL__
        if(ar__(283:288) .EQ. '') e_Hpmag(i__) = rNULL__
        if(ar__(290:294) .EQ. '') Hpscat(i__) = rNULL__
        if(ar__(296:298) .EQ. '') o_Hpmag(i__) = iNULL__
        if(ar__(302:306) .EQ. '') Hpmax(i__) = rNULL__
        if(ar__(308:312) .EQ. '') HPmin(i__) = rNULL__
        if(ar__(314:320) .EQ. '') Period(i__) = rNULL__
        if(ar__(341:342) .EQ. '') Nsys(i__) = iNULL__
        if(ar__(344:345) .EQ. '') Ncomp(i__) = iNULL__
        if(ar__(356:358) .EQ. '') theta(i__) = iNULL__
        if(ar__(360:366) .EQ. '') rho(i__) = rNULL__
        if(ar__(368:372) .EQ. '') e_rho(i__) = rNULL__
        if(ar__(374:378) .EQ. '') dHp(i__) = rNULL__
        if(ar__(380:383) .EQ. '') e_dHp(i__) = rNULL__
        if(ar__(391:396) .EQ. '') HD(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Catalog(i__),HIP(i__),Proxy(i__),RAhms(i__),DEdms(i__),
     +  Vmag(i__),VarFlag(i__),r_Vmag(i__),RAdeg(i__),DEdeg(i__),
     +  AstroRef(i__),Plx(i__),pmRA(i__),pmDE(i__),e_RAdeg(i__),
     +  e_DEdeg(i__),e_Plx(i__),e_pmRA(i__),e_pmDE(i__),DE_RA(i__),
     +  Plx_RA(i__),Plx_DE(i__),pmRA_RA(i__),pmRA_DE(i__),
     +  pmRA_Plx(i__),pmDE_RA(i__),pmDE_DE(i__),pmDE_Plx(i__),
     +  pmDE_pmRA(i__),F1(i__),F2(i__),BTmag(i__),e_BTmag(i__),
     +  VTmag(i__),e_VTmag(i__),m_BTmag(i__),B_V(i__),e_B_V(i__),
     +  r_B_V(i__),V_I(i__),e_V_I(i__),r_V_I(i__),CombMag(i__),
     +  Hpmag(i__),e_Hpmag(i__),Hpscat(i__),o_Hpmag(i__),m_Hpmag(i__),
     +  Hpmax(i__),HPmin(i__),Period(i__),HvarType(i__),moreVar(i__),
     +  morePhoto(i__),CCDM(i__),n_CCDM(i__),Nsys(i__),Ncomp(i__),
     +  MultFlag(i__),Source(i__),Qual(i__),m_HIP(i__),theta(i__),
     +  rho(i__),e_rho(i__),dHp(i__),e_dHp(i__),Survey(i__),
     +  Chart(i__),Notes(i__),HD(i__),BD(i__),CoD(i__),CPD(i__),
     +  v_V_I_red(i__),SpType(i__),r_SpType(i__)
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(i__),S(i__),Type(i__),Source_1(i__),Qual_1(i__),
     +  Notes_1(i__),Nsys_1(i__),Ncomp_1(i__),Nparm(i__),Ncorr(i__),
     +  seq(i__),comp_id(i__),HIP_1(i__),Hp(i__),e_Hp(i__),BT(i__),
     +  e_BT(i__),VT(i__),e_VT(i__),RAdeg_1(i__),DEdeg_1(i__),
     +  Plx_1(i__),pmRA_1(i__),pmDE_1(i__),e_RAdeg_1(i__),
     +  e_DEdeg_1(i__),e_Plx_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),
     +  ref(i__),theta_1(i__),rho_1(i__),d_theta(i__),d_rho(i__),
     +  seq_ref(i__),flag(i__)
        if(ar__1(63:68) .EQ. '') BT(i__) = rNULL__
        if(ar__1(70:74) .EQ. '') e_BT(i__) = rNULL__
        if(ar__1(76:81) .EQ. '') VT(i__) = rNULL__
        if(ar__1(83:87) .EQ. '') e_VT(i__) = rNULL__
        if(ar__1(178:184) .EQ. '') theta_1(i__) = rNULL__
        if(ar__1(186:193) .EQ. '') rho_1(i__) = rNULL__
        if(ar__1(195:202) .EQ. '') d_theta(i__) = rNULL__
        if(ar__1(204:209) .EQ. '') d_rho(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  CCDM_1(i__),S(i__),Type(i__),Source_1(i__),Qual_1(i__),
     +  Notes_1(i__),Nsys_1(i__),Ncomp_1(i__),Nparm(i__),Ncorr(i__),
     +  seq(i__),comp_id(i__),HIP_1(i__),Hp(i__),e_Hp(i__),BT(i__),
     +  e_BT(i__),VT(i__),e_VT(i__),RAdeg_1(i__),DEdeg_1(i__),
     +  Plx_1(i__),pmRA_1(i__),pmDE_1(i__),e_RAdeg_1(i__),
     +  e_DEdeg_1(i__),e_Plx_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),
     +  ref(i__),theta_1(i__),rho_1(i__),d_theta(i__),d_rho(i__),
     +  seq_ref(i__),flag(i__)
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(i__),S_1(i__),Type_1(i__),Source_2(i__),Qual_2(i__),
     +  Notes_2(i__),Nsys_2(i__),Ncomp_2(i__),Nparm_1(i__),
     +  Ncorr_1(i__),seq_1(i__),(corr(j__,i__),j__=1,66)
        if(ar__2(41:238) .EQ. '') corr(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  CCDM_2(i__),S_1(i__),Type_1(i__),Source_2(i__),Qual_2(i__),
     +  Notes_2(i__),Nsys_2(i__),Ncomp_2(i__),Nparm_1(i__),
     +  Ncorr_1(i__),seq_1(i__),(corr(j__,i__),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(i__),gRA(i__),gDE(i__),e_gRA(i__),e_gDE(i__),Fg(i__),
     +  dgRA(i__),dgDE(i__),e_dgRA(i__),e_dgDE(i__),Fdg(i__),
     +  Notes_3(i__),num(i__),(corr_1(j__,i__),j__=1,36)
        if(ar__3(46:52) .EQ. '') dgRA(i__) = rNULL__
        if(ar__3(54:60) .EQ. '') dgDE(i__) = rNULL__
        if(ar__3(62:68) .EQ. '') e_dgRA(i__) = rNULL__
        if(ar__3(70:76) .EQ. '') e_dgDE(i__) = rNULL__
        if(ar__3(78:82) .EQ. '') Fdg(i__) = rNULL__
        if(ar__3(88:195) .EQ. '') corr_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  HIP_2(i__),gRA(i__),gDE(i__),e_gRA(i__),e_gDE(i__),Fg(i__),
     +  dgRA(i__),dgDE(i__),e_dgRA(i__),e_dgDE(i__),Fdg(i__),
     +  Notes_3(i__),num(i__),(corr_1(j__,i__),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(i__),P(i__),T(i__),a0(i__),ecc(i__),w(i__),i(i__),
     +  Omega(i__),e_P(i__),e_T(i__),e_a0(i__),e_ecc(i__),e_w(i__),
     +  e_i(i__),e_Omega(i__),dmRef(i__),Notes_4(i__),flag_1(i__),
     +  (corr_2(j__,i__),j__=1,66)
        if(ar__4(68:75) .EQ. '') e_P(i__) = rNULL__
        if(ar__4(77:85) .EQ. '') e_T(i__) = rNULL__
        if(ar__4(87:91) .EQ. '') e_a0(i__) = rNULL__
        if(ar__4(93:98) .EQ. '') e_ecc(i__) = rNULL__
        if(ar__4(100:105) .EQ. '') e_w(i__) = rNULL__
        if(ar__4(107:112) .EQ. '') e_i(i__) = rNULL__
        if(ar__4(114:119) .EQ. '') e_Omega(i__) = rNULL__
        if(ar__4(121:123) .EQ. '') dmRef(i__) = iNULL__
c    ..............Just test output...........
        write(6,5)
     +  HIP_3(i__),P(i__),T(i__),a0(i__),ecc(i__),w(i__),i(i__),
     +  Omega(i__),e_P(i__),e_T(i__),e_a0(i__),e_ecc(i__),e_w(i__),
     +  e_i(i__),e_Omega(i__),dmRef(i__),Notes_4(i__),flag_1(i__),
     +  (corr_2(j__,i__),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(i__),Hp_ref(i__),DRA(i__),DDE(i__),e_DRA(i__),
     +  e_DDE(i__),FD(i__),theta_C(i__),e_theta_C(i__),minSep(i__),
     +  dvar(i__),Notes_5(i__),(corr_3(j__,i__),j__=1,21)
c    ..............Just test output...........
        write(6,6)
     +  HIP_4(i__),Hp_ref(i__),DRA(i__),DDE(i__),e_DRA(i__),
     +  e_DDE(i__),FD(i__),theta_C(i__),e_theta_C(i__),minSep(i__),
     +  dvar(i__),Notes_5(i__),(corr_3(j__,i__),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(i__),epsilon(i__),e_epsilon(i__),Notes_6(i__)
c    ..............Just test output...........
        write(6,7)
     +  HIP_5(i__),epsilon(i__),e_epsilon(i__),Notes_6(i__)
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(i__),NewVar(i__),SpType_1(i__),HvarType_1(i__),
     +  VarType(i__),n_VarType(i__),maxMag(i__),l_minMag(i__),
     +  minMag(i__),log_sA_A(i__),Period_1(i__),log_sP(i__),
     +  Ep_2440000(i__),q_Ep_2440000(i__),morePhoto_1(i__),
     +  Notes_7(i__),VarName(i__),period_2(i__),epoch(i__),max(i__),
     +  min(i__),Band(i__),refFlag(i__)
        if(ar__7(50:55) .EQ. '') log_sA_A(i__) = rNULL__
        if(ar__7(57:68) .EQ. '') Period_1(i__) = rNULL__
        if(ar__7(70:75) .EQ. '') log_sP(i__) = rNULL__
        if(ar__7(77:85) .EQ. '') Ep_2440000(i__) = rNULL__
        if(ar__7(87:87) .EQ. '') q_Ep_2440000(i__) = iNULL__
        if(ar__7(106:115) .EQ. '') period_2(i__) = rNULL__
        if(ar__7(117:126) .EQ. '') epoch(i__) = rNULL__
        if(ar__7(128:132) .EQ. '') max(i__) = rNULL__
        if(ar__7(134:138) .EQ. '') min(i__) = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  HIP_6(i__),NewVar(i__),SpType_1(i__),HvarType_1(i__),
     +  VarType(i__),n_VarType(i__),maxMag(i__),l_minMag(i__),
     +  minMag(i__),log_sA_A(i__),Period_1(i__),log_sP(i__),
     +  Ep_2440000(i__),q_Ep_2440000(i__),morePhoto_1(i__),
     +  Notes_7(i__),VarName(i__),period_2(i__),epoch(i__),max(i__),
     +  min(i__),Band(i__),refFlag(i__)
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(i__),NewVar_1(i__),SpType_2(i__),HvarType_2(i__),
     +  VarType_1(i__),n_VarType_1(i__),maxMag_1(i__),l_minMag_1(i__),
     +  minMag_1(i__),med_Hp(i__),A(i__),e_A(i__),morePhoto_2(i__),
     +  Notes_8(i__),VarName_1(i__),period_3(i__),Ep_2440000_1(i__),
     +  max_1(i__),min_1(i__),Band_1(i__),refFlag_1(i__)
        if(ar__8(106:115) .EQ. '') period_3(i__) = rNULL__
        if(ar__8(117:126) .EQ. '') Ep_2440000_1(i__) = rNULL__
        if(ar__8(128:132) .EQ. '') max_1(i__) = rNULL__
        if(ar__8(134:138) .EQ. '') min_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,9)
     +  HIP_7(i__),NewVar_1(i__),SpType_2(i__),HvarType_2(i__),
     +  VarType_1(i__),n_VarType_1(i__),maxMag_1(i__),l_minMag_1(i__),
     +  minMag_1(i__),med_Hp(i__),A(i__),e_A(i__),morePhoto_2(i__),
     +  Notes_8(i__),VarName_1(i__),period_3(i__),Ep_2440000_1(i__),
     +  max_1(i__),min_1(i__),Band_1(i__),refFlag_1(i__)
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(i__),RAdeg_2(i__),DEdeg_2(i__),Ep_2440000_2(i__),
     +  delay(i__),theta_2(i__),e_lambda(i__),flag_2(i__)
c    ..............Just test output...........
        write(6,10)
     +  ID(i__),RAdeg_2(i__),DEdeg_2(i__),Ep_2440000_2(i__),
     +  delay(i__),theta_2(i__),e_lambda(i__),flag_2(i__)
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(i__),Ep_2440000_3(i__),Hp_dc(i__),e_Hp_dc(i__),
     +  Hp_ac(i__),e_Hp_ac(i__),r(i__),Delta(i__),alpha(i__)
c    ..............Just test output...........
        write(6,11)
     +  ID_1(i__),Ep_2440000_3(i__),Hp_dc(i__),e_Hp_dc(i__),
     +  Hp_ac(i__),e_Hp_ac(i__),r(i__),Delta(i__),alpha(i__)
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(i__),Ep_2440000_4(i__),RAdeg_3(i__),DEdeg_3(i__),
     +  BTmag_1(i__),VTmag_1(i__),Flag_3(i__),e_RAdeg_2(i__),
     +  e_DEdeg_2(i__),DE_RA_1(i__),theta_3(i__),sign_z(i__),
     +  e_incl(i__),e_vert(i__)
        if(ar__11(43:47) .EQ. '') BTmag_1(i__) = rNULL__
        if(ar__11(49:53) .EQ. '') VTmag_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,12)
     +  ID_2(i__),Ep_2440000_4(i__),RAdeg_3(i__),DEdeg_3(i__),
     +  BTmag_1(i__),VTmag_1(i__),Flag_3(i__),e_RAdeg_2(i__),
     +  e_DEdeg_2(i__),DE_RA_1(i__),theta_3(i__),sign_z(i__),
     +  e_incl(i__),e_vert(i__)
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(i__),Note1(i__),Note2(i__),Ntot(i__),Nline(i__),
     +  Text(i__)
c    ..............Just test output...........
        write(6,13)
     +  HIP_8(i__),Note1(i__),Note2(i__),Ntot(i__),Nline(i__),
     +  Text(i__)
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(i__),Note1_1(i__),Note2_1(i__),Ntot_1(i__),Nline_1(i__),
     +  Text_1(i__)
c    ..............Just test output...........
        write(6,14)
     +  HIP_9(i__),Note1_1(i__),Note2_1(i__),Ntot_1(i__),Nline_1(i__),
     +  Text_1(i__)
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(i__),Note1_2(i__),Note2_2(i__),Ntot_2(i__),
     +  Nline_2(i__),Text_2(i__)
c    ..............Just test output...........
        write(6,15)
     +  HIP_10(i__),Note1_2(i__),Note2_2(i__),Ntot_2(i__),
     +  Nline_2(i__),Text_2(i__)
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(i__),Ntot_3(i__),Nline_3(i__),nRef(i__)
c    ..............Just test output...........
        write(6,16)HIP_11(i__),Ntot_3(i__),Nline_3(i__),nRef(i__)
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(i__),Text_3(i__)
c    ..............Just test output...........
        write(6,17)nRef_1(i__),Text_3(i__)
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(i__),Ntot_4(i__),Nline_4(i__),Text_4(i__)
c    ..............Just test output...........
        write(6,18)dmRef_1(i__),Ntot_4(i__),Nline_4(i__),Text_4(i__)
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(i__),TYC(i__),Proxy_1(i__),RAhms_1(i__),
     +  DEdms_1(i__),Vmag_1(i__),r_Vmag_1(i__),RAdeg_4(i__),
     +  DEdeg_4(i__),AstroRef_1(i__),Plx_2(i__),pmRA_2(i__),
     +  pmDE_2(i__),e_RAdeg_3(i__),e_DEdeg_3(i__),e_Plx_2(i__),
     +  e_pmRA_2(i__),e_pmDE_2(i__),DE_RA_2(i__),Plx_RA_1(i__),
     +  Plx_DE_1(i__),pmRA_RA_1(i__),pmRA_DE_1(i__),pmRA_Plx_1(i__),
     +  pmDE_RA_1(i__),pmDE_DE_1(i__),pmDE_Plx_1(i__),
     +  pmDE_pmRA_1(i__),Nastro(i__),F2_1(i__),HIP_12(i__),
     +  BTmag_2(i__),e_BTmag_1(i__),VTmag_2(i__),e_VTmag_1(i__),
     +  r_BTmag(i__),B_V_1(i__),e_B_V_1(i__),Q(i__),Fs(i__),
     +  Source_3(i__),Nphoto(i__),VTscat(i__),VTmax(i__),VTmin(i__),
     +  Var(i__),VarFlag_1(i__),MultFlag_1(i__),morePhoto_3(i__),
     +  m_HIP_1(i__),PPM(i__),HD_1(i__),BD_1(i__),CoD_1(i__),
     +  CPD_1(i__),Remark(i__)
        if(ar__18(42:46) .EQ. '') Vmag_1(i__) = rNULL__
        if(ar__18(80:86) .EQ. '') Plx_2(i__) = rNULL__
        if(ar__18(88:95) .EQ. '') pmRA_2(i__) = rNULL__
        if(ar__18(97:104) .EQ. '') pmDE_2(i__) = rNULL__
        if(ar__18(106:111) .EQ. '') e_RAdeg_3(i__) = rNULL__
        if(ar__18(113:118) .EQ. '') e_DEdeg_3(i__) = rNULL__
        if(ar__18(120:125) .EQ. '') e_Plx_2(i__) = rNULL__
        if(ar__18(127:132) .EQ. '') e_pmRA_2(i__) = rNULL__
        if(ar__18(134:139) .EQ. '') e_pmDE_2(i__) = rNULL__
        if(ar__18(141:145) .EQ. '') DE_RA_2(i__) = rNULL__
        if(ar__18(147:151) .EQ. '') Plx_RA_1(i__) = rNULL__
        if(ar__18(153:157) .EQ. '') Plx_DE_1(i__) = rNULL__
        if(ar__18(159:163) .EQ. '') pmRA_RA_1(i__) = rNULL__
        if(ar__18(165:169) .EQ. '') pmRA_DE_1(i__) = rNULL__
        if(ar__18(171:175) .EQ. '') pmRA_Plx_1(i__) = rNULL__
        if(ar__18(177:181) .EQ. '') pmDE_RA_1(i__) = rNULL__
        if(ar__18(183:187) .EQ. '') pmDE_DE_1(i__) = rNULL__
        if(ar__18(189:193) .EQ. '') pmDE_Plx_1(i__) = rNULL__
        if(ar__18(195:199) .EQ. '') pmDE_pmRA_1(i__) = rNULL__
        if(ar__18(201:203) .EQ. '') Nastro(i__) = iNULL__
        if(ar__18(205:209) .EQ. '') F2_1(i__) = rNULL__
        if(ar__18(211:216) .EQ. '') HIP_12(i__) = iNULL__
        if(ar__18(218:223) .EQ. '') BTmag_2(i__) = rNULL__
        if(ar__18(225:229) .EQ. '') e_BTmag_1(i__) = rNULL__
        if(ar__18(231:236) .EQ. '') VTmag_2(i__) = rNULL__
        if(ar__18(238:242) .EQ. '') e_VTmag_1(i__) = rNULL__
        if(ar__18(246:251) .EQ. '') B_V_1(i__) = rNULL__
        if(ar__18(253:257) .EQ. '') e_B_V_1(i__) = rNULL__
        if(ar__18(261:261) .EQ. '') Q(i__) = iNULL__
        if(ar__18(263:266) .EQ. '') Fs(i__) = rNULL__
        if(ar__18(270:272) .EQ. '') Nphoto(i__) = iNULL__
        if(ar__18(274:278) .EQ. '') VTscat(i__) = rNULL__
        if(ar__18(280:284) .EQ. '') VTmax(i__) = rNULL__
        if(ar__18(286:290) .EQ. '') VTmin(i__) = rNULL__
        if(ar__18(303:308) .EQ. '') PPM(i__) = iNULL__
        if(ar__18(310:315) .EQ. '') HD_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,19)
     +  Catalog_1(i__),TYC(i__),Proxy_1(i__),RAhms_1(i__),
     +  DEdms_1(i__),Vmag_1(i__),r_Vmag_1(i__),RAdeg_4(i__),
     +  DEdeg_4(i__),AstroRef_1(i__),Plx_2(i__),pmRA_2(i__),
     +  pmDE_2(i__),e_RAdeg_3(i__),e_DEdeg_3(i__),e_Plx_2(i__),
     +  e_pmRA_2(i__),e_pmDE_2(i__),DE_RA_2(i__),Plx_RA_1(i__),
     +  Plx_DE_1(i__),pmRA_RA_1(i__),pmRA_DE_1(i__),pmRA_Plx_1(i__),
     +  pmDE_RA_1(i__),pmDE_DE_1(i__),pmDE_Plx_1(i__),
     +  pmDE_pmRA_1(i__),Nastro(i__),F2_1(i__),HIP_12(i__),
     +  BTmag_2(i__),e_BTmag_1(i__),VTmag_2(i__),e_VTmag_1(i__),
     +  r_BTmag(i__),B_V_1(i__),e_B_V_1(i__),Q(i__),Fs(i__),
     +  Source_3(i__),Nphoto(i__),VTscat(i__),VTmax(i__),VTmin(i__),
     +  Var(i__),VarFlag_1(i__),MultFlag_1(i__),morePhoto_3(i__),
     +  m_HIP_1(i__),PPM(i__),HD_1(i__),BD_1(i__),CoD_1(i__),
     +  CPD_1(i__),Remark(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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