FORTRAN Generation
()

Conversion of standardized ReadMe file for catalog 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.5, on 2014-Jul-25
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__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  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

      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

      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,file='catalog.dat', status='old')
      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,file='rv_std.dat', status='old')
      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