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

Conversion of standardized ReadMe file for file /./ftp/cats/V/137D into FORTRAN code for reading data files line by line.

Note that special values are assigned to unknown or unspecified numbers (also called NULL numbers); when necessary, the coordinate components making up the right ascension and declination are converted into floating-point numbers representing these angles in degrees.



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2017-Dec-11
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         ! Hipparcos identifier
      character*6   Comp        ! Component(s) (1)
      character*41  Classes     ! SIMBAD classifications (comma separated) (2)
      integer*4     Gr          ! ? Cluster or Association membership(s) (3)
      real*8        RAdeg       ! (deg) Right ascension (ICRS, Epoch=J1991.25)
      real*8        DEdeg       ! (deg) Declination (ICRS, Epoch=J1991.25)
      real*4        Plx         ! (mas) Trigonometric parallax
      real*8        pmRA        ! (mas/yr) Proper motion in RA*cos(DEdeg)
      real*8        pmDE        ! (mas/yr) Proper motion in Declination
      real*4        e_RAdeg     ! (mas) Standard error on RA*cos(DEdeg)
      real*4        e_DEdeg     ! (mas) Standard error on DE
      real*4        e_Plx       ! (mas) Standard error on Plx
      real*4        e_pmRA      ! (mas/yr) Standard error on pmRA
      real*4        e_pmDE      ! (mas/yr) Standard error on pmDE
      integer*4     r_HIP       ! [1/2] Reference for Hipparcos astrometry (4)
      integer*4     r_pm        ! [1/3] Reference for proper motion (5)
      real*8        GLon        ! (deg) Galactic longitude (6)
      real*8        GLat        ! (deg) Galactic latitude (6)
      real*8        Dist        ! (pc) ? Heliocentric distance (7)
      real*4        e_Dist      ! (%) ? Distance error expressed as percentage (8)
      real*8        pmGLon      ! (mas/yr) Proper motion in GLon*cos(GLat)
      real*8        pmGLat      ! (mas/yr) Proper motion in GLat
      real*4        X           ! (pc) ? Heliocentric distance towards Gal. center
      real*8        Y           ! (pc) ? Heliocentric distance towards Gal. rotation
      real*4        Z           ! (pc) ? Heliocentric distance towards N. Gal. Pole
      integer*4     RGal        ! (pc) ? Galactocentric distance (9)
      real*4        vT          ! (km/s) ? Transverse velocity
      character*26  SpType      ! Spectral type (MK, HD, or other)
      integer*4     Tc          ! ]0/140[? Temperature class codified (10)
      integer*4     Lc          ! [1/6]? Luminosity class codified (11)
      real*8        RV          ! (km/s) ? Radial velocity
      real*4        e_RV        ! (km/s) ? Standard error on RV (12)
      character*1   q_RV        ! [ABCD] Quality flag on RV (13)
      real*4        v_Fe_H_     ! ([Sun]) ? Iron abundance
      real*4        e__Fe_H_    ! ([Sun]) ? Standard error on [Fe/H] (14)
      character*1   q__Fe_H_    ! [ABC] Quality flag on [Fe/H] (15)
      real*4        age         ! (Gyr) ? Age, in billions of years
      real*4        b_age       ! (Gyr) ? Lower confidence limit on age
      real*4        B_age_1     ! (Gyr) ? Upper confidence limit on age
      real*4        U           ! (km/s) ? Heliocentric velocity towards Gal. center
      real*4        V           ! (km/s) ? Heliocentric velocity towards Gal. rotation
      real*4        W           ! (km/s) ? Heliocentric velocity towards N. Gal. Pole
      real*4        UVW         ! (km/s) ? Total heliocentric velocity
      real*8        Dmin        ! (pc) ? Minimum distance at solar encounter (16)
      integer*4     Tmin        ! (kyr) ? Timing of Dmin, in thousands of years (17)
      real*4        e           ! ? Total orbital eccentricity (18)
      real*4        phi         ! (deg) [-180/180]? Pericenter position angle (18)
      integer*4     a           ! (pc) ? Semi-major axis of orbit (18)
      integer*4     b           ! (pc) ? Semi-minor axis of orbit (18)
      integer*4     c           ! (pc) ? Focus-to-center distance of orbit (18)
      integer*4     L           ! (deg) ? Semilatus rectum of orbit (18)
      integer*4     Rmin        ! (pc) ? Orbital radius at pericenter (18)
      integer*4     Rmax        ! (pc) ? Orbital radius at apocenter (18)
      integer*4     Npl         ! ? Number of exoplanets (known in April 2012)
      character*5   Mpl         ! 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       ! (pc) Hipparcos identifier
      real*8        Hpmag       ! (mag) Median magnitude in Hipparcos system
      real*4        e_Hpmag     ! (mag) Standard error on Hpmag
      character*1   m_Hpmag     ! [A-Z*-] Reference flag for Hpmag (20)
      real*4        Hpmax       ! (mag) ? Hpmag at maximum (5th percentile)
      real*4        Hpmin       ! (mag) ? Hpmag at minimum (95th percentile)
      real*4        Per         ! (d) ? Variability period
      character*1   Hvar        ! (mag) [CDMPRU] Variability type (21)
      real*4        Umag        ! (mag) ? Apparent magnitude in Johnson U
      real*4        Bmag        ! (mag) ? Apparent magnitude in Johnson B
      real*4        Vmag        ! (mag) Apparent magnitude in Johnson V
      real*4        Rmag        ! (mag) ? Apparent magnitude in R
      real*4        Imag        ! (mag) ? Apparent magnitude in I
      real*4        Jmag        ! (mag) ? J selected default magnitude (22)
      real*4        Hmag        ! (mag) ? H selected default magnitude (22)
      real*4        Kmag        ! (mag) ? K selected default magnitude (22)
      real*4        e_Jmag      ! (mag) ? J total magnitude uncertainty (22)
      real*4        e_Hmag      ! (mag) ? H total magnitude uncertainty (22)
      real*4        e_Kmag      ! (mag) ? K total magnitude uncertainty (22)
      character*17  v2MASS      ! 2MASS source designation (22)
      character*3   q2M         ! JHK photometric quality flag (22)
      real*4        B_V         ! (mag) ? Johnson B-V color index
      real*4        V_I         ! (mag) ? Color index in Cousins' system
      real*4        e_B_V       ! (mag) ? Standard error on B-V
      real*4        e_V_I       ! (mag) ? Standard error on V-I
      real*4        HpMag_1     ! (mag) ? Absolute Magnitude of Hpmag (23)
      real*4        UMag_1      ! (mag) ? Absolute Magnitude in Johnson U (23)
      real*4        BMag_1      ! (mag) ? Absolute Magnitude in Johnson B (23)
      real*4        VMag_1      ! (mag) ? Absolute Magnitude in Johnson V (23)
      real*4        RMag_1      ! (mag) ? Absolute Magnitude in R (23)
      real*4        IMag_1      ! (mag) ? Absolute Magnitude in I (23)
      real*4        JMag_1      ! (mag) ? Absolute Magnitude J (23)
      real*4        HMag_1      ! (mag) ? Absolute Magnitude H (23)
      real*4        KMag_1      ! (mag) ? Absolute Magnitude K (23)
      real*8        Lum         ! (Lsun) ? Stellar luminosity (23)
      real*4        magmin      ! (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       ! Hipparcos identifier
      integer*4     HD          ! ? Henry Draper catalog identifier
      character*3   Cst         ! Constellation membership
      integer*4     Atlas       ! Millennium Star Atlas page number
      character*10  Coords      ! RA, DE in compact format; n = + ; s = -
      character*48  Name        ! Star name(s)
      character*106 GrpName     ! Cluster or Association name(s)
      character*16  CompID      ! Reference for Comp (main.dat)
      character*25  rSpType     ! Reference for SpType (main.dat) (24)
      character*15  r_RV        ! Reference for RV (main.dat)
      character*127 r__Fe_H_    ! 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     ! [1,703] 3-digit shorthand code (25)
      character*19  BibCode     ! 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   ! Group name
      character*21  AltName     ! Alternative group name
      integer*4     N           ! Number of Hiipparcos candidate stars
      real*4        n_1         ! [3.5/16] Separation parameter (27)
      character*13  Stream      ! Kinematic group (28)
      real*4        RAdeg_1     ! (deg) Group mean right ascension (ICRS)
      real*4        DEdeg_1     ! (deg) Group mean declination (ICRS)
      real*4        Dist_1      ! (pc) Group mean parallax distance
      real*4        pmRA_1      ! (mas/yr) Group mean proper motion in RA*cos(DEdeg)
      real*4        pmDE_1      ! (mas/yr) Mean proper motion in declination
      real*4        RV_1        ! (km/s) Group mean radial velocity
      real*4        s_RAdeg     ! (deg) Group width on RA*cos(DEdeg) (1-{sigma})
      real*4        s_DEdeg     ! (deg) Group width on DE (1-{sigma})
      real*4        s_Dist      ! (pc) Group width on Dist (1-{sigma})
      real*4        s_pmRA      ! (mas/yr) Group width on pmRA (1-{sigma})
      real*4        s_pmDE      ! (mas/yr) Group width on pmDE (1-{sigma})
      real*4        s_RV        ! (km/s) Group width on radial velocity
      real*4        e_RAdeg_1   ! (deg) Standard error on RA*cos(DEdeg)
      real*4        e_DEdeg_1   ! (deg) Standard error on DE
      real*4        e_Dist_1    ! (pc) Standard error on Dist
      real*4        e_pmRA_1    ! (mas/yr) Standard error on pmRA
      real*4        e_pmDE_1    ! (mas/yr) Standard error on pmDE
      real*4        e_RV_1      ! (km/s) Standard error on radial velocity
      real*4        rad         ! (pc) Physical radius (1-{sigma})
      real*4        arad        ! (deg) Angular radius (1-{sigma})
      real*4        X_1         ! (pc) Heliocentric distance towards Gal. center
      real*4        Y_1         ! (pc) Heliocentric distance towards Gal. rotation
      real*4        Z_1         ! (pc) Heliocentric distance towards N. Gal. Pole
      real*4        U_1         ! (km/s) Heliocentric velocity towards Gal. center
      real*4        V_1         ! (km/s) Heliocentric velocity towards Gal. rotation
      real*4        W_1         ! (km/s) Heliocentric velocity towards N. Gal. Pole
      real*4        s_X         ! (pc) Group width on X (1-{sigma})
      real*4        s_Y         ! (pc) Group width on Y (1-{sigma})
      real*4        s_Z         ! (pc) Group width on Z (1-{sigma})
      real*4        s_U         ! (km/s) Group width on U (1-{sigma})
      real*4        s_V         ! (km/s) Group width on V (1-{sigma})
      real*4        s_W         ! (km/s) Group width on W (1-{sigma})
      real*4        e_X         ! (pc) Standard error on X
      real*4        e_Y         ! (pc) Standard error on Y
      real*4        e_Z         ! (pc) Standard error on Z
      real*4        e_U         ! (km/s) Standard error on U
      real*4        e_V         ! (km/s) Standard error on V
      real*4        e_W         ! (km/s) Standard error on W
      real*4        v_Fe_H__1   ! ([Sun]) ? Iron abundance
      real*4        e__Fe_H__1  ! ([Sun]) ? Standard error on [Fe/H]
      integer*4     n__Fe_H_    ! ? Number of Hipparcos member stars with [Fe/H]
      real*4        E_B_V_1     ! (mag) ? Color excess in B-V, CDS Cat. B/ocl
      real*4        E_B_V_2     ! (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,Comp,Classes,Gr,RAdeg,DEdeg,Plx,pmRA,pmDE,e_RAdeg,e_DEdeg,
     +  e_Plx,e_pmRA,e_pmDE,r_HIP,r_pm,GLon,GLat,Dist,e_Dist,pmGLon,
     +  pmGLat,X,Y,Z,RGal,vT,SpType,Tc,Lc,RV,e_RV,q_RV,v_Fe_H_,
     +  e__Fe_H_,q__Fe_H_,age,b_age,B_age_1,U,V,W,UVW,Dmin,Tmin,e,phi,
     +  a,b,c,L,Rmin,Rmax,Npl,Mpl
        if(ar__(57:57) .EQ. '') Gr = iNULL__
        if(ar__(173:179) .EQ. '') Dist = rNULL__
        if(ar__(181:184) .EQ. '') e_Dist = rNULL__
        if(ar__(204:209) .EQ. '') X = rNULL__
        if(ar__(211:217) .EQ. '') Y = rNULL__
        if(ar__(219:224) .EQ. '') Z = rNULL__
        if(ar__(226:229) .EQ. '') RGal = iNULL__
        if(ar__(231:235) .EQ. '') vT = rNULL__
        if(ar__(264:266) .EQ. '') Tc = iNULL__
        if(ar__(268:268) .EQ. '') Lc = iNULL__
        if(ar__(270:276) .EQ. '') RV = rNULL__
        if(ar__(278:283) .EQ. '') e_RV = rNULL__
        if(ar__(287:291) .EQ. '') v_Fe_H_ = rNULL__
        if(ar__(293:296) .EQ. '') e__Fe_H_ = rNULL__
        if(ar__(300:303) .EQ. '') age = rNULL__
        if(ar__(305:308) .EQ. '') b_age = rNULL__
        if(ar__(310:313) .EQ. '') B_age_1 = rNULL__
        if(ar__(315:320) .EQ. '') U = rNULL__
        if(ar__(322:327) .EQ. '') V = rNULL__
        if(ar__(329:334) .EQ. '') W = rNULL__
        if(ar__(336:340) .EQ. '') UVW = rNULL__
        if(ar__(342:348) .EQ. '') Dmin = rNULL__
        if(ar__(350:354) .EQ. '') Tmin = iNULL__
        if(ar__(356:361) .EQ. '') e = rNULL__
        if(ar__(363:368) .EQ. '') phi = rNULL__
        if(ar__(370:376) .EQ. '') a = iNULL__
        if(ar__(378:382) .EQ. '') b = iNULL__
        if(ar__(384:390) .EQ. '') c = iNULL__
        if(ar__(392:396) .EQ. '') L = iNULL__
        if(ar__(398:401) .EQ. '') Rmin = iNULL__
        if(ar__(403:409) .EQ. '') Rmax = iNULL__
        if(ar__(411:411) .EQ. '') Npl = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  HIP,Comp,Classes,Gr,RAdeg,DEdeg,Plx,pmRA,pmDE,e_RAdeg,e_DEdeg,
     +  e_Plx,e_pmRA,e_pmDE,r_HIP,r_pm,GLon,GLat,Dist,e_Dist,pmGLon,
     +  pmGLat,X,Y,Z,RGal,vT,SpType,Tc,Lc,RV,e_RV,q_RV,v_Fe_H_,
     +  e__Fe_H_,q__Fe_H_,age,b_age,B_age_1,U,V,W,UVW,Dmin,Tmin,e,phi,
     +  a,b,c,L,Rmin,Rmax,Npl,Mpl
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,Hpmag,e_Hpmag,m_Hpmag,Hpmax,Hpmin,Per,Hvar,Umag,Bmag,
     +  Vmag,Rmag,Imag,Jmag,Hmag,Kmag,e_Jmag,e_Hmag,e_Kmag,v2MASS,q2M,
     +  B_V,V_I,e_B_V,e_V_I,HpMag_1,UMag_1,BMag_1,VMag_1,RMag_1,
     +  IMag_1,JMag_1,HMag_1,KMag_1,Lum,magmin
        if(ar__1(25:29) .EQ. '') Hpmax = rNULL__
        if(ar__1(31:35) .EQ. '') Hpmin = rNULL__
        if(ar__1(37:42) .EQ. '') Per = rNULL__
        if(ar__1(46:51) .EQ. '') Umag = rNULL__
        if(ar__1(53:58) .EQ. '') Bmag = rNULL__
        if(ar__1(66:70) .EQ. '') Rmag = rNULL__
        if(ar__1(72:77) .EQ. '') Imag = rNULL__
        if(ar__1(79:84) .EQ. '') Jmag = rNULL__
        if(ar__1(86:91) .EQ. '') Hmag = rNULL__
        if(ar__1(93:98) .EQ. '') Kmag = rNULL__
        if(ar__1(100:104) .EQ. '') e_Jmag = rNULL__
        if(ar__1(106:110) .EQ. '') e_Hmag = rNULL__
        if(ar__1(112:116) .EQ. '') e_Kmag = rNULL__
        if(ar__1(140:145) .EQ. '') B_V = rNULL__
        if(ar__1(147:151) .EQ. '') V_I = rNULL__
        if(ar__1(153:157) .EQ. '') e_B_V = rNULL__
        if(ar__1(159:162) .EQ. '') e_V_I = rNULL__
        if(ar__1(164:169) .EQ. '') HpMag_1 = rNULL__
        if(ar__1(171:176) .EQ. '') UMag_1 = rNULL__
        if(ar__1(178:182) .EQ. '') BMag_1 = rNULL__
        if(ar__1(184:188) .EQ. '') VMag_1 = rNULL__
        if(ar__1(190:194) .EQ. '') RMag_1 = rNULL__
        if(ar__1(196:200) .EQ. '') IMag_1 = rNULL__
        if(ar__1(202:206) .EQ. '') JMag_1 = rNULL__
        if(ar__1(208:212) .EQ. '') HMag_1 = rNULL__
        if(ar__1(214:219) .EQ. '') KMag_1 = rNULL__
        if(ar__1(221:229) .EQ. '') Lum = rNULL__
        if(ar__1(231:235) .EQ. '') magmin = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  HIP_1,Hpmag,e_Hpmag,m_Hpmag,Hpmax,Hpmin,Per,Hvar,Umag,Bmag,
     +  Vmag,Rmag,Imag,Jmag,Hmag,Kmag,e_Jmag,e_Hmag,e_Kmag,v2MASS,q2M,
     +  B_V,V_I,e_B_V,e_V_I,HpMag_1,UMag_1,BMag_1,VMag_1,RMag_1,
     +  IMag_1,JMag_1,HMag_1,KMag_1,Lum,magmin
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,HD,Cst,Atlas,Coords,Name,GrpName,CompID,rSpType,r_RV,
     +  r__Fe_H_
        if(ar__2(8:13) .EQ. '') HD = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  HIP_2,HD,Cst,Atlas,Coords,Name,GrpName,CompID,rSpType,r_RV,
     +  r__Fe_H_
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,BibCode
c    ..............Just test output...........
        write(6,4)KeyCode,BibCode
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,AltName,N,n_1,Stream,RAdeg_1,DEdeg_1,Dist_1,pmRA_1,
     +  pmDE_1,RV_1,s_RAdeg,s_DEdeg,s_Dist,s_pmRA,s_pmDE,s_RV,
     +  e_RAdeg_1,e_DEdeg_1,e_Dist_1,e_pmRA_1,e_pmDE_1,e_RV_1,rad,
     +  arad,X_1,Y_1,Z_1,U_1,V_1,W_1,s_X,s_Y,s_Z,s_U,s_V,s_W,e_X,e_Y,
     +  e_Z,e_U,e_V,e_W,v_Fe_H__1,e__Fe_H__1,n__Fe_H_,E_B_V_1,E_B_V_2
        if(ar__4(295:299) .EQ. '') v_Fe_H__1 = rNULL__
        if(ar__4(301:304) .EQ. '') e__Fe_H__1 = rNULL__
        if(ar__4(306:308) .EQ. '') n__Fe_H_ = iNULL__
        if(ar__4(310:313) .EQ. '') E_B_V_1 = rNULL__
        if(ar__4(315:318) .EQ. '') E_B_V_2 = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  GrpName_1,AltName,N,n_1,Stream,RAdeg_1,DEdeg_1,Dist_1,pmRA_1,
     +  pmDE_1,RV_1,s_RAdeg,s_DEdeg,s_Dist,s_pmRA,s_pmDE,s_RV,
     +  e_RAdeg_1,e_DEdeg_1,e_Dist_1,e_pmRA_1,e_pmDE_1,e_RV_1,rad,
     +  arad,X_1,Y_1,Z_1,U_1,V_1,W_1,s_X,s_Y,s_Z,s_U,s_V,s_W,e_X,e_Y,
     +  e_Z,e_U,e_V,e_W,v_Fe_H__1,e__Fe_H__1,n__Fe_H_,E_B_V_1,E_B_V_2
c    .......End.of.Just test output...........
      end do
      close(1)

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