Conversion of standardized ReadMe file for
file /./ftp/cats/III/144 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-Apr-16
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/144 Radial Velocities for High Proper Motion Stars (Fouts+, 1986)
*================================================================================
*New Subdwarfs. V. Radial Velocities for 889 High-Proper-Motion Stars Measured
*with the Mount Wilson 100-inch Reflector
* Fouts G., Sandage A.
* <Astron. J. 91, 1189 (1986)>
* =1986AJ.....91.1189F
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! The catalog ordered by B1950.0 right ascension
integer*4 nr__
parameter (nr__=878) ! Number of records
character*104 ar__ ! Full-size record
C B1950.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension B1950.0
real*8 DEdeg ! (deg) Declination B1950.0
C ---------------------------------- ! (position vector(s) in degrees)
character*1 flag ! *Flag
character*8 ID ! *Star designation
integer*4 RAh ! (h) Right ascension, RA, B1950.0, hours
integer*4 RAm ! (min) RA, B1950.0, minutes
integer*4 RAs ! (s) RA, B1950.0, seconds
character*1 DE_ ! Declination sign
integer*4 DEd ! (deg) Declination, Dec, B1950.0 degrees
integer*4 DEm ! (arcmin) Dec, B1950.0, minutes
integer*4 DEs ! (arcsec) Dec, B1950.0, seconds
real*4 Glon ! (deg) Galactic longitude
real*4 Glat ! (deg) Galactic latitude
real*4 Vmag ! (mag) *V magnitude
real*4 B_V ! (mag) *B-V color
real*4 U_B ! (mag) *U-B color
real*4 pmRA ! (arcsec/a) *Annual proper motion, RA, B1950.0
real*4 pmDE ! (arcsec/a) *Annual proper motion, Dec, B1950.0
real*4 phi ! (deg) *Parallactic angle
real*4 pmGlon ! (arcsec/a) *Annual proper mot., long., B1950.0
real*4 pmGlat ! (arcsec/a) *Annual proper motion, lat., B1950.0
real*4 RV ! (km/s) Radial velocity
integer*4 o_RV ! Number of observations
*Note on flag:
* An asterisk (*) flag to indicate stars for which prior radial-velocity
* data were reported in Sandage (1969ApJ...158.1115S), Sandage
* (1981AJ.....86.1643S), Wilson (1953, Cat. III/21), and Eggen
* (1964RGOB...84..111E, "A catalogue of high-velocity stars").
*Note on ID:
* Star designation. Lowell Observatory designations (Giclas, Burnham,
* and Thomas 1971) are given if they exist. Otherwise, designations from
* the Luyten proper-motion catalogs (Luyten and Hughes 1980, NLTT and
* supplement, Cat. I/98) the Durchmusterungen (Bonner Durchmusterung,
* Argelander 1859-62, Kuestner 1903, Schoenfeld 1886, Cat. I/119 and
* I/122; Cordoba Durchmusterung, Thome 1892-1932, Cat. I/114; Cape
* Photographic Durchmusterung Gill and Kapteyn 1895-1900, Cat. I/108),
* or the Selected Area catalogs are given.
*Note on Vmag, B-V, and U-B:
* UBV data taken from Sandage and Kowal (1986).
*Note on pmRA and pmDE:
* The adopted proper motions have been averaged from the Lowell and
* Luyten catalogs, the latter proper motions in r and theta being
* converted to proper motions in right ascension and declination via
* the following relations:
* pmRA = pm(r)*sin(theta)
* pmDE = pm(r)*cos(theta)
*Note on pmGlon and pmGlat:
* The galactic proper motions, pmGlon and pmGlat, are calculated from
* the equatorial motions and the parallactic angle (phi):
* pmGlon = pmRA*cos(phi) + pmDE*sin(phi)
* pmGlat = -pmRA*sin(phi) + pm_DE*cos(phi)
*Note on phi:
* The angle between the equatorial and galactic coordinate systems
* (see Smart 1968, "Stellar Kinematics", Cambridge University Press,
* pp. 17-19).
C=============================================================================
C Loading file 'catalog.dat' ! The catalog ordered by B1950.0 right ascension
C Format for file interpretation
1 format(
+ A1,A8,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,F6.2,1X,F6.2,
+ 1X,F5.2,1X,F5.2,1X,F5.2,1X,F6.3,1X,F6.3,1X,F5.1,1X,F6.3,1X,
+ F6.3,1X,F6.1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,878
read(1,'(A104)')ar__
read(ar__,1)
+ flag,ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Glon,Glat,Vmag,B_V,U_B,
+ pmRA,pmDE,phi,pmGlon,pmGlat,RV,o_RV
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ flag,ID,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Glon,Glat,Vmag,B_V,U_B,
+ pmRA,pmDE,phi,pmGlon,pmGlat,RV,o_RV
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end