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 2013-May-18
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