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

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

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



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Mar-29
C=============================================================================
*  This code was generated from the ReadMe file documenting a catalogue
*  according to the "Standard for Documentation of Astronomical Catalogues"
*  currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
*  (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
*  Please report problems or questions to   
C=============================================================================

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

C=============================================================================
Cat. III/239     Radial Velocities with Astrometric Data      (Kharchenko+ 2004)
*================================================================================
*The catalogue of radial velocities of galactic stars with high
*precision astrometric data (CRVAD)
*   Kharchenko N.V., Piskunov A.E., Scholz R.-D.
*  <Astron. Nachr. 325, 439 (2004)>
*  =2004AN....325..439K
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'catalog.dat'	! *The Catalogue

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

C  J2000.0 position composed of: RAhour DEdeg (Epoch=1991.25)
      real*8        RAhour      ! (h) Right Ascension J2000.0, epoch 1991.25
      real*8        DEdeg       ! (deg) Declination J2000.0, epoch 1991.25
      real*4        e_RAhour    ! (mas) Standard error in RA*cos(DEdeg)
      real*4        e_DEdeg     ! (mas) Standard error in DEdeg
      real*8        Plx         ! (mas) ?=9999.99 Trigonometric parallax
      real*8        e_Plx       ! (mas) ?=999.99 Standard error in Plx
      real*8        pmRA        ! (mas/yr) Proper Motion in RA*cos(DEdeg)
      real*8        pmDE        ! (mas/yr) Proper Motion in DE
      real*8        e_pmRA      ! (mas/yr) Standard error in pmRA
      real*8        e_pmDE      ! (mas/yr) Standard error in pmDE
      real*8        Bmag        ! (mag) ?=99.999 B magnitude in Johnson system
      real*8        Vmag        ! (mag) ?=99.999 V magnitude in Johnson system
      real*4        e_Bmag      ! (mag) ?=9.999 Standard error on B magnitude
      real*4        e_Vmag      ! (mag) ?=9.999 Standard error on V magnitude
      real*4        Scat        ! (mag) ?=9.999 Scatter on magnitude
      character*1   v1          ! *[GN ] Known variability from GCVS/NSV
      character*1   v2          ! *[UVW ] Variability from Tycho-1
      character*1   v3          ! *[CDMPRU ] Variability type
      character*1   v4          ! *[VYIXR ] Variability from CMC11
      character*2   d12         ! *[A-S ] CCDM component identifier (Cat. I/274)
      character*1   d3          ! *[A-S ] Component identifier
      character*1   d4          ! *[DRSYZ ] Duplicity from Tycho-1
      character*1   d5          ! *[CGOVX ] Double/Multiple Systems flag
      character*1   d6          ! [D ] Duplicity flag from PPM
      character*12  SpType      ! MK Spectral type (from ASCC-2.5)
      integer*4     HIP         ! ?=0 Hipparcos number (Cat. I/239)
      integer*4     HD          ! ?=0 HD number (Cat. III/135)
      integer*4     ASCC        ! ASCC-2.5 (Cat. I/280) number
      integer*4     GCRV        ! GCRV (Cat. III/213) number
      character*13  CCDM        ! *CCDM number of multiple star and its components
      real*4        mag         ! (mag) ?=99.99 Visual or photographic GCRV magnitude
      character*1   n_mag       ! *[* ] indicates a photographic mag
      character*1   u_mag       ! *[V: ] variability or uncertainty on mag
      character*11  SpGCRV      ! ? MK Spectral type from GCRV (Cat. III/213)
      real*8        RV          ! (km/s) Average Radial Velocity
      real*4        e_RV        ! (km/s) ?=-9.9 Mean standard error in RV
      character*1   n_RV        ! *[Gge ] Note on RV
      character*1   u_RV        ! *[:* ]
      character*1   q_RV        ! *[A-E I] Quality index of the RV
      character*1   f_RV        ! *[SOCR ] Flag of star type
      integer*4     o_RV        ! ?=0 Number of observations
      integer*4     N_mult      ! [1,2] Number of matches for given GCRV entry
*Note on v1:
*    Byte 292 from Tycho-1 file tyc_main.dat
*    This flag is set when the variability is known:
*    G : variable known in General Catalog of Variable Stars
*    N : variable known in New Suspected Variables catalog
*Note on v2:
*    Byte 294 from Tycho-1 file tyc_main.dat
*    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 v3:
*    Byte 322 from Hipparcos file hip_main.dat
*    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 v4:
*    Flag from CMC11 catalogue
*    V      : variable star
*    Y or I : if the CAMC did not measure the mag, an approximate value
*             from other sources is given to a precision of 0.1 mag
*    X or R : star fainter than 12 with unknown mag is given as 12.
*Note on d12:
*    Bytes 149, 150 from Tycho-2 file catalog.dat.
*    The CCDM component identifiers for double or multiple Hipparcos stars
*    contributing to this Tycho-2 entry.
*Note on d3:
*    Byte 41 from Hipparcos file hip_dm_com.dat
*Note on d4:
*    Byte 296 from Tycho-1 file tyc_main.dat
*    Unresolved duplicity status from Tycho data analysis:
*    D     : duplicity clearly indicated (B_Tmag and V_Tmag 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 d5:
*    Byte 347 from Hipparcos file hip_main.dat
*    Indicates that further details are given in the Hipparcos 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 CCDM:
*    CCDM index (Cat. I/274) is copied from the GCRV (Cat. III/213);
*Note on n_mag:
*    This flag is "*" if mag is photographic or B magnitude;
*Note on u_mag:
*    This flag has the following meanings:
*      V for variable star,
*      ":" if the data is uncertain;
*Note on n_RV:
*    This flag has the following meanings:
*      e: for emission radial velocity
*      G: RV of the gravity center for double or multiple system
*      g: RV of the gravity center for a sub-system
*Note on u_RV:
*    This flag has the following meanings:
*      ":" uncertainty on RV,
*      "*" approximate value of RV;
*Note on q_RV:
*    from A to E decreasing quality, I for insufficient data;
*Note on f_RV:
*    This flag has the following meanings:
*          S = STND : UAI's standard star
*          O = ORB  : systemic velocity
*          C = CEPH : cepheid star
*          R = RRLY : RR Lyrae star

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

C  Declarations for 'rv_std.dat'	! Candidates to RV standards

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

C  J2000.0 position composed of: RAhour DEdeg (Epoch=1991.25)
      real*8        RAhour_1    ! (h) Right Ascension J2000.0, epoch 1991.25
      real*8        DEdeg_1     ! (deg) Declination J2000.0, epoch 1991.25
      real*4        e_RAhour_1  ! (mas) Standard error in RA*cos(DEdeg)
      real*4        e_DEdeg_1   ! (mas) Standard error in DEdeg
      real*8        Plx_1       ! (mas) ?=9999.99 Trigonometric parallax
      real*8        e_Plx_1     ! (mas) ?=999.99 Standard error in Plx
      real*8        pmRA_1      ! (mas/yr) Proper Motion in RA*cos(DE)
      real*8        pmDE_1      ! (mas/yr) Proper Motion in DE
      real*8        e_pm        ! (mas/yr) Standard error in pm
      real*8        Vmag_1      ! (mag) V magnitude in Johnson system
      real*4        e_Vmag_1    ! (mag) Standard error on V magnitude
      real*8        B_V         ! (mag) (B-V) colour in Johnson system
      real*4        e_B_V       ! (mag) Standard error on (B-V) colour
      character*12  Sp          ! ? Spectral type and luminosity class
      real*8        RV_1        ! (km/s) Mean Radial Velocity
      real*4        e_RV_1      ! (km/s) *?=-9.9 Mean standard error in RV
      character*1   q_RV_1      ! *[AB ] Quality index
      integer*4     o_RV_1      ! Number of observations
      integer*4     HIP_1       ! ?=0 Hipparcos number (Cat. I/239)
      integer*4     HD_1        ! ?=0 HD number (Cat. III/135)
      integer*4     ASCC_1      ! ASCC-2.5 number (Cat. I/280)
      integer*4     GCRV_1      ! GCRV number  (Cat. III/213)
*Note on e_RV, q_RV:
*    the accuracy of radial velocity is described by a standard error
*    or quality index A or B, that correspond to an error of 0.74 and
*    1.78 km/s, respectively.

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

C  Loading file 'catalog.dat'	! *The Catalogue

C  Format for file interpretation

    1 format(
     +  F12.9,1X,F12.8,F6.1,F6.1,F8.2,F7.2,F9.2,F9.2,F7.2,F7.2,F7.3,
     +  F7.3,F6.3,F6.3,F6.3,1X,A1,A1,A1,A1,A2,A1,A1,A1,A1,1X,A12,1X,
     +  I6,1X,I6,1X,I7,1X,I5,1X,A13,1X,F5.2,1X,A1,A1,1X,A11,1X,F7.2,
     +  F5.1,1X,A1,A1,A1,A1,I4,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,34553
        read(1,'(A227)')ar__
        read(ar__,1)
     +  RAhour,DEdeg,e_RAhour,e_DEdeg,Plx,e_Plx,pmRA,pmDE,e_pmRA,
     +  e_pmDE,Bmag,Vmag,e_Bmag,e_Vmag,Scat,v1,v2,v3,v4,d12,d3,d4,d5,
     +  d6,SpType,HIP,HD,ASCC,GCRV,CCDM,mag,n_mag,u_mag,SpGCRV,RV,
     +  e_RV,n_RV,u_RV,q_RV,f_RV,o_RV,N_mult
c    ..............Just test output...........
        write(6,1)
     +  RAhour,DEdeg,e_RAhour,e_DEdeg,Plx,e_Plx,pmRA,pmDE,e_pmRA,
     +  e_pmDE,Bmag,Vmag,e_Bmag,e_Vmag,Scat,v1,v2,v3,v4,d12,d3,d4,d5,
     +  d6,SpType,HIP,HD,ASCC,GCRV,CCDM,mag,n_mag,u_mag,SpGCRV,RV,
     +  e_RV,n_RV,u_RV,q_RV,f_RV,o_RV,N_mult
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'rv_std.dat'	! Candidates to RV standards

C  Format for file interpretation

    2 format(
     +  F12.9,1X,F12.8,F6.1,F6.1,F8.2,F7.2,F9.2,F9.2,F7.2,F7.3,F6.3,
     +  F7.3,F6.3,1X,A12,1X,F7.2,F5.1,1X,A1,I4,1X,I6,1X,I6,1X,I7,1X,I5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'rv_std.dat')
      write(6,*) '....Loading file: rv_std.dat'
      do i__=1,3967
        read(1,'(A163)')ar__1
        read(ar__1,2)
     +  RAhour_1,DEdeg_1,e_RAhour_1,e_DEdeg_1,Plx_1,e_Plx_1,pmRA_1,
     +  pmDE_1,e_pm,Vmag_1,e_Vmag_1,B_V,e_B_V,Sp,RV_1,e_RV_1,q_RV_1,
     +  o_RV_1,HIP_1,HD_1,ASCC_1,GCRV_1
c    ..............Just test output...........
        write(6,2)
     +  RAhour_1,DEdeg_1,e_RAhour_1,e_DEdeg_1,Plx_1,e_Plx_1,pmRA_1,
     +  pmDE_1,e_pm,Vmag_1,e_Vmag_1,B_V,e_B_V,Sp,RV_1,e_RV_1,q_RV_1,
     +  o_RV_1,HIP_1,HD_1,ASCC_1,GCRV_1
c    .......End.of.Just test output...........
      end do
      close(1)

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