FORTRAN Generation
(/./ftp/cats/V/137D)

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

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

C=============================================================================
Cat. V/137D        Extended Hipparcos Compilation (XHIP)          (Anderson+, 2012)
*================================================================================
*XHIP: An Extended Hipparcos Compilation
*     Anderson E., Francis C.
*    <Astron. Letters 38 (2012)>
*   =2012AstL...38..331A  (http://arxiv.org/abs/1108.4971)
*   =2012PAZh...38..374A
*   =2012yCat.5137....0A
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'main.dat'	! Astrometry, spectrography, space motions, exoplanets

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

C  Position composed of: RAdeg DEdeg (Epoch=J1991.25)
      integer*4     HIP        (nr__) ! Hipparcos identifier
      character*6   Comp       (nr__) ! Component(s) (1)
      character*41  Classes    (nr__) ! SIMBAD classifications (comma separated) (2)
      integer*4     Gr         (nr__) ! ? Cluster or Association membership(s) (3)
      real*8        RAdeg      (nr__) ! (deg) Right ascension (ICRS, Epoch=J1991.25)
      real*8        DEdeg      (nr__) ! (deg) Declination (ICRS, Epoch=J1991.25)
      real*4        Plx        (nr__) ! (mas) Trigonometric parallax
      real*8        pmRA       (nr__) ! (mas/yr) Proper motion in RA*cos(DEdeg)
      real*8        pmDE       (nr__) ! (mas/yr) Proper motion in Declination
      real*4        e_RAdeg    (nr__) ! (mas) Standard error on RA*cos(DEdeg)
      real*4        e_DEdeg    (nr__) ! (mas) Standard error on DE
      real*4        e_Plx      (nr__) ! (mas) Standard error on Plx
      real*4        e_pmRA     (nr__) ! (mas/yr) Standard error on pmRA
      real*4        e_pmDE     (nr__) ! (mas/yr) Standard error on pmDE
      integer*4     r_HIP      (nr__) ! [1/2] Reference for Hipparcos astrometry (4)
      integer*4     r_pm       (nr__) ! [1/3] Reference for proper motion (5)
      real*8        GLon       (nr__) ! (deg) Galactic longitude (6)
      real*8        GLat       (nr__) ! (deg) Galactic latitude (6)
      real*8        Dist       (nr__) ! (pc) ? Heliocentric distance (7)
      real*4        e_Dist     (nr__) ! (%) ? Distance error expressed as percentage (8)
      real*8        pmGLon     (nr__) ! (mas/yr) Proper motion in GLon*cos(GLat)
      real*8        pmGLat     (nr__) ! (mas/yr) Proper motion in GLat
      real*4        X          (nr__) ! (pc) ? Heliocentric distance towards Gal. center
      real*8        Y          (nr__) ! (pc) ? Heliocentric distance towards Gal. rotation
      real*4        Z          (nr__) ! (pc) ? Heliocentric distance towards N. Gal. Pole
      integer*4     RGal       (nr__) ! (pc) ? Galactocentric distance (9)
      real*4        vT         (nr__) ! (km/s) ? Transverse velocity
      character*26  SpType     (nr__) ! Spectral type (MK, HD, or other)
      integer*4     Tc         (nr__) ! ]0/140[? Temperature class codified (10)
      integer*4     Lc         (nr__) ! [1/6]? Luminosity class codified (11)
      real*8        RV         (nr__) ! (km/s) ? Radial velocity
      real*4        e_RV       (nr__) ! (km/s) ? Standard error on RV (12)
      character*1   q_RV       (nr__) ! [ABCD] Quality flag on RV (13)
      real*4        v_Fe_H_    (nr__) ! ([Sun]) ? Iron abundance
      real*4        e__Fe_H_   (nr__) ! ([Sun]) ? Standard error on [Fe/H] (14)
      character*1   q__Fe_H_   (nr__) ! [ABC] Quality flag on [Fe/H] (15)
      real*4        age        (nr__) ! (Gyr) ? Age, in billions of years
      real*4        b_age      (nr__) ! (Gyr) ? Lower confidence limit on age
      real*4        B_age_1    (nr__) ! (Gyr) ? Upper confidence limit on age
      real*4        U          (nr__) ! (km/s) ? Heliocentric velocity towards Gal. center
      real*4        V          (nr__) ! (km/s) ? Heliocentric velocity towards Gal. rotation
      real*4        W          (nr__) ! (km/s) ? Heliocentric velocity towards N. Gal. Pole
      real*4        UVW        (nr__) ! (km/s) ? Total heliocentric velocity
      real*8        Dmin       (nr__) ! (pc) ? Minimum distance at solar encounter (16)
      integer*4     Tmin       (nr__) ! (kyr) ? Timing of Dmin, in thousands of years (17)
      real*4        e          (nr__) ! ? Total orbital eccentricity (18)
      real*4        phi        (nr__) ! (deg) [-180/180]? Pericenter position angle (18)
      integer*4     a          (nr__) ! (pc) ? Semi-major axis of orbit (18)
      integer*4     b          (nr__) ! (pc) ? Semi-minor axis of orbit (18)
      integer*4     c          (nr__) ! (pc) ? Focus-to-center distance of orbit (18)
      integer*4     L          (nr__) ! (deg) ? Semilatus rectum of orbit (18)
      integer*4     Rmin       (nr__) ! (pc) ? Orbital radius at pericenter (18)
      integer*4     Rmax       (nr__) ! (pc) ? Orbital radius at apocenter (18)
      integer*4     Npl        (nr__) ! ? Number of exoplanets (known in April 2012)
      character*5   Mpl        (nr__) ! Exoplanet discovery method(s) (19)
*Note (1): from CCDM (Cat. I/274) or WDS (Cat. B/wds)
*Note (2): object types, see
*     http://simbad.u-strasbg.fr/simbad/sim-display?data=otypes
*Note (3): Group IDs are in biblio.dat (GrpName)
*     see also: groups.dat and Cat. B/Ocl
*Note (4): flag indicates astrometry source:
*     1 = HIP (Cat. I/239) for 3,504 problematic cases in HIP2
*     2 = HIP2 (Cat. I/311) for all remaining cases
*Note (5): flag indicates the following:
*     1 = From HIP2 (Cat. I/311)
*     2 = From Tycho-2 (Cat. I/259)
*     3 = HIP2 & Tycho-2 combined for improved error bounds
*Note (6): transformation from ICRS defined by:
*     NGP RAdeg = 192.85948,
*     NGP DEdeg = +27.12825,
*     GLon of ascending node @ celestial equator = 32.93192 deg
*Note (7): null where parallax errors <20%
*     Dist = 1000/Plx * (1+1.2*(e_Plx/Plx)^2^)
*Note (8): null where Dist is fitted to cluster
*Note (9): adopts Sgr A* (X,Y,Z) as (7400,-7.2,-6.0) pc
*Note (10): Temperature classes are:
*      O=10 B=20 A=30 F=40 G=50 K=60 M=70
*      L=80 T=90 S=100 C=110 R=120 N=130
*      Subclasses 0-9 are additionally summed.
*Note (11): I=1 II=2 III=3 IV=4 V=5 VI=6
*Note (12): 999 = error not available though RV is available
*Note (13): quality flag means:
*       A = Most likely to be within stated error bounds
*       B = May have small systematic errors
*       C = May have larger systematic errors
*       D = No error available or other serious problems
*Note (14): 9.99 = error not available though [Fe/H] is available
*Note (15): quality flag means:
*       A = Calibrated
*       B = Calibrated, but no error bounds assigned
*       C = Uncalibrated and no error bounds assigned
*Note (16): based on linear approximation to space motion
*Note (17): Past times are negative; future times positive
*Note (18): LSR (U,V,W) = (14.1,14.6,6.9) km/s are
*      adopted.  Note (9) above with regard to Rgal also applies.
*Note (19): method is codified:
*     RA = "radial velocity or astrometric methods"
*     I  = "imaging"
*     T  = "timing"
*     X  = "transit"

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

C  Declarations for 'photo.dat'	! Photometry

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

      integer*4     HIP_1      (nr__1) ! (pc) Hipparcos identifier
      real*8        Hpmag      (nr__1) ! (mag) Median magnitude in Hipparcos system
      real*4        e_Hpmag    (nr__1) ! (mag) Standard error on Hpmag
      character*1   m_Hpmag    (nr__1) ! [A-Z*-] Reference flag for Hpmag (20)
      real*4        Hpmax      (nr__1) ! (mag) ? Hpmag at maximum (5th percentile)
      real*4        Hpmin      (nr__1) ! (mag) ? Hpmag at minimum (95th percentile)
      real*4        Per        (nr__1) ! (d) ? Variability period
      character*1   Hvar       (nr__1) ! (mag) [CDMPRU] Variability type (21)
      real*4        Umag       (nr__1) ! (mag) ? Apparent magnitude in Johnson U
      real*4        Bmag       (nr__1) ! (mag) ? Apparent magnitude in Johnson B
      real*4        Vmag       (nr__1) ! (mag) Apparent magnitude in Johnson V
      real*4        Rmag       (nr__1) ! (mag) ? Apparent magnitude in R
      real*4        Imag       (nr__1) ! (mag) ? Apparent magnitude in I
      real*4        Jmag       (nr__1) ! (mag) ? J selected default magnitude (22)
      real*4        Hmag       (nr__1) ! (mag) ? H selected default magnitude (22)
      real*4        Kmag       (nr__1) ! (mag) ? K selected default magnitude (22)
      real*4        e_Jmag     (nr__1) ! (mag) ? J total magnitude uncertainty (22)
      real*4        e_Hmag     (nr__1) ! (mag) ? H total magnitude uncertainty (22)
      real*4        e_Kmag     (nr__1) ! (mag) ? K total magnitude uncertainty (22)
      character*17  v2MASS     (nr__1) ! 2MASS source designation (22)
      character*3   q2M        (nr__1) ! JHK photometric quality flag (22)
      real*4        B_V        (nr__1) ! (mag) ? Johnson B-V color index
      real*4        V_I        (nr__1) ! (mag) ? Color index in Cousins' system
      real*4        e_B_V      (nr__1) ! (mag) ? Standard error on B-V
      real*4        e_V_I      (nr__1) ! (mag) ? Standard error on V-I
      real*4        HpMag_1    (nr__1) ! (mag) ? Absolute Magnitude of Hpmag (23)
      real*4        UMag_1     (nr__1) ! (mag) ? Absolute Magnitude in Johnson U (23)
      real*4        BMag_1     (nr__1) ! (mag) ? Absolute Magnitude in Johnson B (23)
      real*4        VMag_1     (nr__1) ! (mag) ? Absolute Magnitude in Johnson V (23)
      real*4        RMag_1     (nr__1) ! (mag) ? Absolute Magnitude in R (23)
      real*4        IMag_1     (nr__1) ! (mag) ? Absolute Magnitude in I (23)
      real*4        JMag_1     (nr__1) ! (mag) ? Absolute Magnitude J (23)
      real*4        HMag_1     (nr__1) ! (mag) ? Absolute Magnitude H (23)
      real*4        KMag_1     (nr__1) ! (mag) ? Absolute Magnitude K (23)
      real*8        Lum        (nr__1) ! (Lsun) ? Stellar luminosity (23)
      real*4        magmin     (nr__1) ! (mag) ? Apparent magnitude V at Tmin (main.dat)
*Note (20): 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 (21): Hipparcos-defined type of variability
* blank = 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 (22): See Cat. II/246
*Note (23): does not factor absorption
*      by the interstellar medium; may also contain integrated components.

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

C  Declarations for 'biblio.dat'	! References

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

      integer*4     HIP_2      (nr__2) ! Hipparcos identifier
      integer*4     HD         (nr__2) ! ? Henry Draper catalog identifier
      character*3   Cst        (nr__2) ! Constellation membership
      integer*4     Atlas      (nr__2) ! Millennium Star Atlas page number
      character*10  Coords     (nr__2) ! RA, DE in compact format; n = + ; s = -
      character*48  Name       (nr__2) ! Star name(s)
      character*106 GrpName    (nr__2) ! Cluster or Association name(s)
      character*16  CompID     (nr__2) ! Reference for Comp (main.dat)
      character*25  rSpType    (nr__2) ! Reference for SpType (main.dat) (24)
      character*15  r_RV       (nr__2) ! Reference for RV (main.dat)
      character*127 r__Fe_H_   (nr__2) ! Reference(s) for [Fe/H] (main.dat)
*Note (24): BibCodes "1993BICDS..43....5T" & "2003AJ....125..359W" are
*     augmented (after a colon :) with the content of the r_Sp or r_SpType
*     columns from those sources respectively.

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

C  Declarations for 'refs.dat'	! Reference key for r_RV & r_[Fe/H]

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

      integer*4     KeyCode    (nr__3) ! [1,703] 3-digit shorthand code (25)
      character*19  BibCode    (nr__3) ! Corresponding 19-digit bibliographic code (26)
*Note (25): Range of key code:
*  * [001-176]: Calibrated [Fe/H] sources
*  * [201-495]: Uncalibrated [Fe/H] sources
*  * [501-703]: Radial velocity sources
*Note (26): Radial velocities obtained through SIMBAD having no source
*     reference are indicated with the pseudo-code "2010.........SIMBAD"

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

C  Declarations for 'groups.dat'	! Statistics on open clusters and stellar associations

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

C  Position composed of: RAdeg DEdeg
      character*21  GrpName_1  (nr__4) ! Group name
      character*21  AltName    (nr__4) ! Alternative group name
      integer*4     N          (nr__4) ! Number of Hiipparcos candidate stars
      real*4        n_1        (nr__4) ! [3.5/16] Separation parameter (27)
      character*13  Stream     (nr__4) ! Kinematic group (28)
      real*4        RAdeg_1    (nr__4) ! (deg) Group mean right ascension (ICRS)
      real*4        DEdeg_1    (nr__4) ! (deg) Group mean declination (ICRS)
      real*4        Dist_1     (nr__4) ! (pc) Group mean parallax distance
      real*4        pmRA_1     (nr__4) ! (mas/yr) Group mean proper motion in RA*cos(DEdeg)
      real*4        pmDE_1     (nr__4) ! (mas/yr) Mean proper motion in declination
      real*4        RV_1       (nr__4) ! (km/s) Group mean radial velocity
      real*4        s_RAdeg    (nr__4) ! (deg) Group width on RA*cos(DEdeg) (1-{sigma})
      real*4        s_DEdeg    (nr__4) ! (deg) Group width on DE (1-{sigma})
      real*4        s_Dist     (nr__4) ! (pc) Group width on Dist (1-{sigma})
      real*4        s_pmRA     (nr__4) ! (mas/yr) Group width on pmRA (1-{sigma})
      real*4        s_pmDE     (nr__4) ! (mas/yr) Group width on pmDE (1-{sigma})
      real*4        s_RV       (nr__4) ! (km/s) Group width on radial velocity
      real*4        e_RAdeg_1  (nr__4) ! (deg) Standard error on RA*cos(DEdeg)
      real*4        e_DEdeg_1  (nr__4) ! (deg) Standard error on DE
      real*4        e_Dist_1   (nr__4) ! (pc) Standard error on Dist
      real*4        e_pmRA_1   (nr__4) ! (mas/yr) Standard error on pmRA
      real*4        e_pmDE_1   (nr__4) ! (mas/yr) Standard error on pmDE
      real*4        e_RV_1     (nr__4) ! (km/s) Standard error on radial velocity
      real*4        rad        (nr__4) ! (pc) Physical radius (1-{sigma})
      real*4        arad       (nr__4) ! (deg) Angular radius (1-{sigma})
      real*4        X_1        (nr__4) ! (pc) Heliocentric distance towards Gal. center
      real*4        Y_1        (nr__4) ! (pc) Heliocentric distance towards Gal. rotation
      real*4        Z_1        (nr__4) ! (pc) Heliocentric distance towards N. Gal. Pole
      real*4        U_1        (nr__4) ! (km/s) Heliocentric velocity towards Gal. center
      real*4        V_1        (nr__4) ! (km/s) Heliocentric velocity towards Gal. rotation
      real*4        W_1        (nr__4) ! (km/s) Heliocentric velocity towards N. Gal. Pole
      real*4        s_X        (nr__4) ! (pc) Group width on X (1-{sigma})
      real*4        s_Y        (nr__4) ! (pc) Group width on Y (1-{sigma})
      real*4        s_Z        (nr__4) ! (pc) Group width on Z (1-{sigma})
      real*4        s_U        (nr__4) ! (km/s) Group width on U (1-{sigma})
      real*4        s_V        (nr__4) ! (km/s) Group width on V (1-{sigma})
      real*4        s_W        (nr__4) ! (km/s) Group width on W (1-{sigma})
      real*4        e_X        (nr__4) ! (pc) Standard error on X
      real*4        e_Y        (nr__4) ! (pc) Standard error on Y
      real*4        e_Z        (nr__4) ! (pc) Standard error on Z
      real*4        e_U        (nr__4) ! (km/s) Standard error on U
      real*4        e_V        (nr__4) ! (km/s) Standard error on V
      real*4        e_W        (nr__4) ! (km/s) Standard error on W
      real*4        v_Fe_H__1  (nr__4) ! ([Sun]) ? Iron abundance
      real*4        e__Fe_H__1 (nr__4) ! ([Sun]) ? Standard error on [Fe/H]
      integer*4     n__Fe_H_   (nr__4) ! ? Number of Hipparcos member stars with [Fe/H]
      real*4        E_B_V_1    (nr__4) ! (mag) ? Color excess in B-V, CDS Cat. B/ocl
      real*4        E_B_V_2    (nr__4) ! (mag) ? Color excess in B-V, CDS Cat. J/A+A/477/165
*Note (27): the separation parameter n is a measure of the concentration of
*     a cluster, and of its separation from the surrounding star field (low
*     values show poor dynamical separation from surrounding stars). A value
*     of 5 is usual for a cluster, while 3.5 is typical for an
*     association.
*Note (28): see 2012MNRAS.422.1283F (arXiv:1202.1375)

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

C  Loading file 'main.dat'	! Astrometry, spectrography, space motions, exoplanets

C  Format for file interpretation

    1 format(
     +  I6,1X,A6,1X,A41,1X,I1,1X,F12.8,1X,F12.8,1X,F6.2,1X,F8.2,1X,
     +  F8.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,I1,1X,I1,1X,
     +  F12.8,1X,F12.8,1X,F7.2,1X,F4.1,1X,F8.2,1X,F8.2,1X,F6.1,1X,
     +  F7.1,1X,F6.1,1X,I4,1X,F5.1,1X,A26,1X,I3,1X,I1,1X,F7.2,1X,F6.2,
     +  1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F6.1,
     +  1X,F6.1,1X,F6.1,1X,F5.1,1X,F7.2,1X,I5,1X,F6.4,1X,F6.1,1X,I7,
     +  1X,I5,1X,I7,1X,I5,1X,I4,1X,I7,1X,I1,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'main.dat')
      write(6,*) '....Loading file: main.dat'
      do i__=1,117955
        read(1,'(A417)')ar__
        read(ar__,1)
     +  HIP(i__),Comp(i__),Classes(i__),Gr(i__),RAdeg(i__),DEdeg(i__),
     +  Plx(i__),pmRA(i__),pmDE(i__),e_RAdeg(i__),e_DEdeg(i__),
     +  e_Plx(i__),e_pmRA(i__),e_pmDE(i__),r_HIP(i__),r_pm(i__),
     +  GLon(i__),GLat(i__),Dist(i__),e_Dist(i__),pmGLon(i__),
     +  pmGLat(i__),X(i__),Y(i__),Z(i__),RGal(i__),vT(i__),
     +  SpType(i__),Tc(i__),Lc(i__),RV(i__),e_RV(i__),q_RV(i__),
     +  v_Fe_H_(i__),e__Fe_H_(i__),q__Fe_H_(i__),age(i__),b_age(i__),
     +  B_age_1(i__),U(i__),V(i__),W(i__),UVW(i__),Dmin(i__),
     +  Tmin(i__),e(i__),phi(i__),a(i__),b(i__),c(i__),L(i__),
     +  Rmin(i__),Rmax(i__),Npl(i__),Mpl(i__)
        if(ar__(57:57) .EQ. '') Gr(i__) = iNULL__
        if(ar__(173:179) .EQ. '') Dist(i__) = rNULL__
        if(ar__(181:184) .EQ. '') e_Dist(i__) = rNULL__
        if(ar__(204:209) .EQ. '') X(i__) = rNULL__
        if(ar__(211:217) .EQ. '') Y(i__) = rNULL__
        if(ar__(219:224) .EQ. '') Z(i__) = rNULL__
        if(ar__(226:229) .EQ. '') RGal(i__) = iNULL__
        if(ar__(231:235) .EQ. '') vT(i__) = rNULL__
        if(ar__(264:266) .EQ. '') Tc(i__) = iNULL__
        if(ar__(268:268) .EQ. '') Lc(i__) = iNULL__
        if(ar__(270:276) .EQ. '') RV(i__) = rNULL__
        if(ar__(278:283) .EQ. '') e_RV(i__) = rNULL__
        if(ar__(287:291) .EQ. '') v_Fe_H_(i__) = rNULL__
        if(ar__(293:296) .EQ. '') e__Fe_H_(i__) = rNULL__
        if(ar__(300:303) .EQ. '') age(i__) = rNULL__
        if(ar__(305:308) .EQ. '') b_age(i__) = rNULL__
        if(ar__(310:313) .EQ. '') B_age_1(i__) = rNULL__
        if(ar__(315:320) .EQ. '') U(i__) = rNULL__
        if(ar__(322:327) .EQ. '') V(i__) = rNULL__
        if(ar__(329:334) .EQ. '') W(i__) = rNULL__
        if(ar__(336:340) .EQ. '') UVW(i__) = rNULL__
        if(ar__(342:348) .EQ. '') Dmin(i__) = rNULL__
        if(ar__(350:354) .EQ. '') Tmin(i__) = iNULL__
        if(ar__(356:361) .EQ. '') e(i__) = rNULL__
        if(ar__(363:368) .EQ. '') phi(i__) = rNULL__
        if(ar__(370:376) .EQ. '') a(i__) = iNULL__
        if(ar__(378:382) .EQ. '') b(i__) = iNULL__
        if(ar__(384:390) .EQ. '') c(i__) = iNULL__
        if(ar__(392:396) .EQ. '') L(i__) = iNULL__
        if(ar__(398:401) .EQ. '') Rmin(i__) = iNULL__
        if(ar__(403:409) .EQ. '') Rmax(i__) = iNULL__
        if(ar__(411:411) .EQ. '') Npl(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  HIP(i__),Comp(i__),Classes(i__),Gr(i__),RAdeg(i__),DEdeg(i__),
     +  Plx(i__),pmRA(i__),pmDE(i__),e_RAdeg(i__),e_DEdeg(i__),
     +  e_Plx(i__),e_pmRA(i__),e_pmDE(i__),r_HIP(i__),r_pm(i__),
     +  GLon(i__),GLat(i__),Dist(i__),e_Dist(i__),pmGLon(i__),
     +  pmGLat(i__),X(i__),Y(i__),Z(i__),RGal(i__),vT(i__),
     +  SpType(i__),Tc(i__),Lc(i__),RV(i__),e_RV(i__),q_RV(i__),
     +  v_Fe_H_(i__),e__Fe_H_(i__),q__Fe_H_(i__),age(i__),b_age(i__),
     +  B_age_1(i__),U(i__),V(i__),W(i__),UVW(i__),Dmin(i__),
     +  Tmin(i__),e(i__),phi(i__),a(i__),b(i__),c(i__),L(i__),
     +  Rmin(i__),Rmax(i__),Npl(i__),Mpl(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'photo.dat'	! Photometry

C  Format for file interpretation

    2 format(
     +  I6,1X,F7.4,1X,F6.4,1X,A1,1X,F5.2,1X,F5.2,1X,F6.2,1X,A1,1X,
     +  F6.3,1X,F6.3,1X,F5.2,1X,F5.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,
     +  1X,F5.3,1X,F5.3,1X,F5.3,1X,A17,1X,A3,1X,F6.3,1X,F5.2,1X,F5.3,
     +  1X,F4.2,1X,F6.3,1X,F6.3,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
     +  F5.2,1X,F5.2,1X,F6.2,1X,F9.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'photo.dat')
      write(6,*) '....Loading file: photo.dat'
      do i__=1,117955
        read(1,'(A235)')ar__1
        read(ar__1,2)
     +  HIP_1(i__),Hpmag(i__),e_Hpmag(i__),m_Hpmag(i__),Hpmax(i__),
     +  Hpmin(i__),Per(i__),Hvar(i__),Umag(i__),Bmag(i__),Vmag(i__),
     +  Rmag(i__),Imag(i__),Jmag(i__),Hmag(i__),Kmag(i__),e_Jmag(i__),
     +  e_Hmag(i__),e_Kmag(i__),v2MASS(i__),q2M(i__),B_V(i__),
     +  V_I(i__),e_B_V(i__),e_V_I(i__),HpMag_1(i__),UMag_1(i__),
     +  BMag_1(i__),VMag_1(i__),RMag_1(i__),IMag_1(i__),JMag_1(i__),
     +  HMag_1(i__),KMag_1(i__),Lum(i__),magmin(i__)
        if(ar__1(25:29) .EQ. '') Hpmax(i__) = rNULL__
        if(ar__1(31:35) .EQ. '') Hpmin(i__) = rNULL__
        if(ar__1(37:42) .EQ. '') Per(i__) = rNULL__
        if(ar__1(46:51) .EQ. '') Umag(i__) = rNULL__
        if(ar__1(53:58) .EQ. '') Bmag(i__) = rNULL__
        if(ar__1(66:70) .EQ. '') Rmag(i__) = rNULL__
        if(ar__1(72:77) .EQ. '') Imag(i__) = rNULL__
        if(ar__1(79:84) .EQ. '') Jmag(i__) = rNULL__
        if(ar__1(86:91) .EQ. '') Hmag(i__) = rNULL__
        if(ar__1(93:98) .EQ. '') Kmag(i__) = rNULL__
        if(ar__1(100:104) .EQ. '') e_Jmag(i__) = rNULL__
        if(ar__1(106:110) .EQ. '') e_Hmag(i__) = rNULL__
        if(ar__1(112:116) .EQ. '') e_Kmag(i__) = rNULL__
        if(ar__1(140:145) .EQ. '') B_V(i__) = rNULL__
        if(ar__1(147:151) .EQ. '') V_I(i__) = rNULL__
        if(ar__1(153:157) .EQ. '') e_B_V(i__) = rNULL__
        if(ar__1(159:162) .EQ. '') e_V_I(i__) = rNULL__
        if(ar__1(164:169) .EQ. '') HpMag_1(i__) = rNULL__
        if(ar__1(171:176) .EQ. '') UMag_1(i__) = rNULL__
        if(ar__1(178:182) .EQ. '') BMag_1(i__) = rNULL__
        if(ar__1(184:188) .EQ. '') VMag_1(i__) = rNULL__
        if(ar__1(190:194) .EQ. '') RMag_1(i__) = rNULL__
        if(ar__1(196:200) .EQ. '') IMag_1(i__) = rNULL__
        if(ar__1(202:206) .EQ. '') JMag_1(i__) = rNULL__
        if(ar__1(208:212) .EQ. '') HMag_1(i__) = rNULL__
        if(ar__1(214:219) .EQ. '') KMag_1(i__) = rNULL__
        if(ar__1(221:229) .EQ. '') Lum(i__) = rNULL__
        if(ar__1(231:235) .EQ. '') magmin(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  HIP_1(i__),Hpmag(i__),e_Hpmag(i__),m_Hpmag(i__),Hpmax(i__),
     +  Hpmin(i__),Per(i__),Hvar(i__),Umag(i__),Bmag(i__),Vmag(i__),
     +  Rmag(i__),Imag(i__),Jmag(i__),Hmag(i__),Kmag(i__),e_Jmag(i__),
     +  e_Hmag(i__),e_Kmag(i__),v2MASS(i__),q2M(i__),B_V(i__),
     +  V_I(i__),e_B_V(i__),e_V_I(i__),HpMag_1(i__),UMag_1(i__),
     +  BMag_1(i__),VMag_1(i__),RMag_1(i__),IMag_1(i__),JMag_1(i__),
     +  HMag_1(i__),KMag_1(i__),Lum(i__),magmin(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'biblio.dat'	! References

C  Format for file interpretation

    3 format(
     +  I6,1X,I6,1X,A3,1X,I4,1X,A10,1X,A48,1X,A106,1X,A16,1X,A25,1X,
     +  A15,1X,A127)

C  Effective file loading

      open(unit=1,status='old',file=
     +'biblio.dat')
      write(6,*) '....Loading file: biblio.dat'
      do i__=1,117955
        read(1,'(A376)')ar__2
        read(ar__2,3)
     +  HIP_2(i__),HD(i__),Cst(i__),Atlas(i__),Coords(i__),Name(i__),
     +  GrpName(i__),CompID(i__),rSpType(i__),r_RV(i__),r__Fe_H_(i__)
        if(ar__2(8:13) .EQ. '') HD(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  HIP_2(i__),HD(i__),Cst(i__),Atlas(i__),Coords(i__),Name(i__),
     +  GrpName(i__),CompID(i__),rSpType(i__),r_RV(i__),r__Fe_H_(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! Reference key for r_RV & r_[Fe/H]

C  Format for file interpretation

    4 format(I3,1X,A19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,674
        read(1,'(A23)')ar__3
        read(ar__3,4)KeyCode(i__),BibCode(i__)
c    ..............Just test output...........
        write(6,4)KeyCode(i__),BibCode(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'groups.dat'	! Statistics on open clusters and stellar associations

C  Format for file interpretation

    5 format(
     +  A21,1X,A21,1X,I3,1X,F4.1,1X,A13,1X,F6.2,1X,F6.2,1X,F6.1,1X,
     +  F6.1,1X,F6.1,1X,F5.1,1X,F6.2,1X,F6.2,1X,F6.1,1X,F5.1,1X,F5.1,
     +  1X,F4.1,1X,F5.2,1X,F5.2,1X,F6.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,
     +  F5.1,1X,F5.1,1X,F6.1,1X,F6.1,1X,F6.1,1X,F5.1,1X,F5.1,1X,F5.1,
     +  1X,F5.1,1X,F5.1,1X,F5.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F5.1,1X,
     +  F5.1,1X,F5.1,1X,F3.1,1X,F3.1,1X,F3.1,1X,F5.2,1X,F4.2,1X,I3,1X,
     +  F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'groups.dat')
      write(6,*) '....Loading file: groups.dat'
      do i__=1,87
        read(1,'(A318)')ar__4
        read(ar__4,5)
     +  GrpName_1(i__),AltName(i__),N(i__),n_1(i__),Stream(i__),
     +  RAdeg_1(i__),DEdeg_1(i__),Dist_1(i__),pmRA_1(i__),pmDE_1(i__),
     +  RV_1(i__),s_RAdeg(i__),s_DEdeg(i__),s_Dist(i__),s_pmRA(i__),
     +  s_pmDE(i__),s_RV(i__),e_RAdeg_1(i__),e_DEdeg_1(i__),
     +  e_Dist_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),e_RV_1(i__),
     +  rad(i__),arad(i__),X_1(i__),Y_1(i__),Z_1(i__),U_1(i__),
     +  V_1(i__),W_1(i__),s_X(i__),s_Y(i__),s_Z(i__),s_U(i__),
     +  s_V(i__),s_W(i__),e_X(i__),e_Y(i__),e_Z(i__),e_U(i__),
     +  e_V(i__),e_W(i__),v_Fe_H__1(i__),e__Fe_H__1(i__),
     +  n__Fe_H_(i__),E_B_V_1(i__),E_B_V_2(i__)
        if(ar__4(295:299) .EQ. '') v_Fe_H__1(i__) = rNULL__
        if(ar__4(301:304) .EQ. '') e__Fe_H__1(i__) = rNULL__
        if(ar__4(306:308) .EQ. '') n__Fe_H_(i__) = iNULL__
        if(ar__4(310:313) .EQ. '') E_B_V_1(i__) = rNULL__
        if(ar__4(315:318) .EQ. '') E_B_V_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  GrpName_1(i__),AltName(i__),N(i__),n_1(i__),Stream(i__),
     +  RAdeg_1(i__),DEdeg_1(i__),Dist_1(i__),pmRA_1(i__),pmDE_1(i__),
     +  RV_1(i__),s_RAdeg(i__),s_DEdeg(i__),s_Dist(i__),s_pmRA(i__),
     +  s_pmDE(i__),s_RV(i__),e_RAdeg_1(i__),e_DEdeg_1(i__),
     +  e_Dist_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),e_RV_1(i__),
     +  rad(i__),arad(i__),X_1(i__),Y_1(i__),Z_1(i__),U_1(i__),
     +  V_1(i__),W_1(i__),s_X(i__),s_Y(i__),s_Z(i__),s_U(i__),
     +  s_V(i__),s_W(i__),e_X(i__),e_Y(i__),e_Z(i__),e_U(i__),
     +  e_V(i__),e_W(i__),v_Fe_H__1(i__),e__Fe_H__1(i__),
     +  n__Fe_H_(i__),E_B_V_1(i__),E_B_V_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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