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.81 (2015-09-23), on 2024-May-07
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. I/112A Lowell Proper Motion Survey - Southern Hemisphere (Giclas+ 1979)
*================================================================================
*Lowell Proper Motion Survey - Southern Hemisphere Catalog 1979
* Giclas H.L., Burnham Jr. R., Thomas N.G.
* <Lowell Observatory Bulletin No. 164, (1978)>
* =1978LowOB...8...89G
* =1979LowOB...8..145G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'south.dat' ! The catalog of Southern G-numbered stars
integer*4 nr__
parameter (nr__=2827) ! Number of records
character*98 ar__ ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg ! (deg) Right Ascension 1950
real*8 DEdeg ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*8 G ! Original G number assigned (1)
character*1 n_G ! [*] Note in the printed catalog (2)
integer*4 RAh ! (h) Right ascension hours (Equinox 1950)
integer*4 RAm ! (min) Right ascension minutes
integer*4 RAs ! (s) Right ascension seconds
character*1 DE_ ! Declination sign (Equinox 1950)
integer*4 DEd ! (deg) Declination degrees
real*4 DEm ! (arcmin) Declination arcminutes
real*4 pm ! (arcsec/yr) Annual total proper motion (3)
integer*4 pmPA ! (deg) Proper motion position angle (3)
real*4 Pmag ! (mag) Estimated Photographic (blue) magnitude
integer*4 col ! [-1/4]? Color class (4)
character*44 X_id ! Other identifications (5)
integer*4 LOB ! Identification chart reference (6)
*Note (1):
* byte 1 of the field always contains a G and byte 5 always contains a
* dash (-). The number in bytes 2-4 is the plate number, while that in
* bytes 6-8 is the star number on the plate. Original G number assigned
* to the star on the plate region where it was first found. Stars that
* were subsequently found on later plates were assigned new G numbers;
* however, in this catalog all duplicate observations have been combined
* and averaged with the original G numbers retained. In cases where a
* close companion has been found on a later plate, the source of a
* finding chart prepared from that plate may be found in the
* identification chart field (LOB). When the catalog entry is a mean
* formed from two or more observations, the additional stars making
* up the mean are listed first in the X-id column.
*Note (2):
* An asterisk (*) indicates that a note concerning the star appears in
* the Bulletin in which it was first published. The notes are included
* in the "notes.htx" file (a note about a peculiar entry is introduced
* by the "\entry" text)
*Note (3):
* The value "pm" represents the annual proper motion in seconds of arc.
* The components of the proper motion are:
*
* {mu}_{alpha}_*cos({delta}) = pm*sin(pmPA) ;
* {mu}_{delta}_ = pm*cos(pmPA)
*Note (4):
* The estimated color class is on a scale from -2 (the bluest) to +4
* (the very reddest) stars found. Almost all stars of color classes -2
* and -1 have been subsequently identified from photometry or spectra
* as white dwarfs.
*
* Because there is such a wide variety of different types of objects
* among these motions, it is impossible to define a narrow relationship
* between these color classes and the conventional B-V and U-B colors.
* The table below is given with the warning that the true color of any
* one individual star may deviate from the mean value by as much as half
* a magnitude for stars redder than color class 0.
* ---------------------------------
* Color Class B-V U-B
* ---------------------------------
* -1 +0.11 -0.74
* 0 +0.26 -0.57
* +1 +0.80 +0.33
* +2 +1.30 +1.09
* +3 +1.61 +1.18
* ---------------------------------
*Note (5):
* Catalog entries are, in many cases, the combined means of measurements
* on two or more plate regions; in these cases, the additional G
* number(s) appear as the first entry(ies) in the Reference field.
* If a G number in the references carries an asterisk, but the field
* following the identification field (byte 9) does not, then the note on
* that star appears only in the Lowell Observatory Bulletin where the
* latter number was published. Following the G numbers, references to
* other proper-motion lists and catalogs are given, as identified in the
* reference key in the appendix of this document.
*Note (6):
* The Lowell Bulletin number in which an identification chart (18x18')
* is given for the star. The summary catalog was compiled from
* collections of data appearing in earlier Lowell Observatory Bulletins.
* There are no finding charts in the source reference, but they can be
* found in the original references given in the table of the Note on
* n_ID.
C=============================================================================
C Loading file 'south.dat' ! The catalog of Southern G-numbered stars
C Format for file interpretation
1 format(
+ A8,A1,1X,I2,1X,I2,1X,I2,2X,A1,I2,1X,F4.1,1X,F5.2,2X,I3,2X,
+ F4.1,2X,I2,2X,A44,I3)
C Effective file loading
open(unit=1,status='old',file=
+'I/112A/south.dat')
write(6,*) '....Loading file: south.dat'
do i__=1,2827
read(1,'(A98)')ar__
read(ar__,1)
+ G,n_G,RAh,RAm,RAs,DE_,DEd,DEm,pm,pmPA,Pmag,col,X_id,LOB
if(ar__(48:49) .EQ. '') col = iNULL__
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(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ G,n_G,RAh,RAm,RAs,DE_,DEd,DEm,pm,pmPA,Pmag,col,X_id,LOB
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end