Conversion of standardized ReadMe file for
file /./ftp/cats/I/112A 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-Jun-20
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. 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
real*8 RA_ ! (deg) Right Ascension 1950
real*8 Dec ! (deg) Declination 1950
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 Declarations for 'synonyms.dat' ! Synonyms of G stars measured 2 or more times
integer*4 nr__1
parameter (nr__1=626) ! Number of records
character*19 ar__1 ! Full-size record
character*8 Gsyn ! G number in X-id column of "south.dat" file
character*8 G_1 ! The main G number
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'pos158.dat' ! Accurate positions for field#158 (from B. Skiff)
integer*4 nr__2
parameter (nr__2=69) ! Number of records
character*41 ar__2 ! Full-size record
real*8 RA__1 ! (deg) Right Ascension J2000
real*8 Dec_1 ! (deg) Declination J2000
character*8 G_2 ! Original G number assigned (GFFF-NNN)
integer*4 RAh_1 ! (h) Right Ascension J2000 (hours)
integer*4 RAm_1 ! (min) Right Ascension J2000 (minutes)
real*4 RAs_1 ! (s) Right Ascension J2000 (seconds)
character*1 DE__1 ! Declination J2000 (sign)
integer*4 DEd_1 ! (deg) Declination J2000 (degrees)
integer*4 DEm_1 ! (arcmin) Declination J2000 (minutes)
real*4 DEs ! (arcsec) Declination J2000 (seconds)
character*1 s ! [CMSTUb] Source of position (1)
real*4 Vmag ! (mag) Approximate visual magnitude
*Note (1): source of position as follows:
* b = USNO-B1.0 (2003AJ....125..984M, I/284)
* C = Carlsberg meridian circle, CMC14 (I/304)
* M = 2MASS (2006AJ....131.1163S, II/246)
* S = SDSS-DR7 (2009ApJS..182..543A, II/294)
* T = Tycho-2 (2000A&A...357..367H, I/259)
* U = UCAC2 (2004AJ....127.3043Z, I/289) + UCAC3 (I/315)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'acro.dat' ! References of used acronyms
integer*4 nr__3
parameter (nr__3=68) ! Number of records
character*332 ar__3 ! Full-size record
character*17 Acronym ! Designation of catalog or list (1)
character*1 Sep ! [=:] original list (=) or supplement (:) (1)
character*311 Text ! Text of Reference
*Note (1):
* The references to other proper-motion catalogs and lists are
* identified by character prefixes, followed by the numerical
* designations of the objects in a catalog or list, according to
* standard astronomical usage. The prefixes are identified with their
* corresponding references. (Note: not all of the identifiers listed
* actually appear in the catalog. The key is repeated from the Lowell
* Proper Motion Survey, Northern Hemisphere (Giclas, Burnham, and Thomas
* 1971) and some identifiers may only appear in that catalog).
*
* When the column 'Sep' (byte 20) contains the equal sign (=), the
* bibliographical sources that provided information for the catalog
* consist of references for the star designations in other catalogs and
* lists, plus references for general information. The general references
* are combined with some supplementary identification references in
* the source publications and that order has not been changed in the
* following sections. The information given below appears in the source
* reference and is repeated here for users who do not have ready access
* to that publication.
*
* When the column 'Sep' (byte 20) contains the colon (:), the
* supplementary references include sources of some identifiers used
* in the catalog, plus general reference sources consulted during
* preparation of the catalog.
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,file='south.dat', status='old')
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__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAs .GE. 0) RA_=RA_+RAs/240.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
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)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'synonyms.dat' ! Synonyms of G stars measured 2 or more times
C Format for file interpretation
2 format(A8,3X,A8)
C Effective file loading
open(unit=1,file='synonyms.dat', status='old')
write(6,*) '....Loading file: synonyms.dat'
do i__=1,626
read(1,'(A19)')ar__1
read(ar__1,2)Gsyn,G_1
c ..............Just test output...........
write(6,2)Gsyn,G_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'pos158.dat' ! Accurate positions for field#158 (from B. Skiff)
C Format for file interpretation
3 format(
+ A8,3X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A1,1X,F4.1)
C Effective file loading
open(unit=1,file='pos158.dat', status='old')
write(6,*) '....Loading file: pos158.dat'
do i__=1,69
read(1,'(A41)')ar__2
read(ar__2,3)
+ G_2,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs,s,Vmag
c Derive coordinates RA__1 and Dec_1 from input data
c (RA__1 and Dec_1 are set to rNULL_ when unknown)
RA__1 = RAh_1
if(RAh_1 .GE. 0) RA__1=RAh_1*15.
if(RAm_1 .GE. 0) RA__1=RA__1+RAm_1/4.
if(RAs_1 .GE. 0) RA__1=RA__1+RAs_1/240.
Dec_1 = DEd_1
if(DEm_1 .GE. 0) Dec_1=Dec_1+DEm_1/60.
if(DEs .GE. 0) Dec_1=Dec_1+DEs/3600.
if(DE__1.EQ.'-'.AND.Dec_1.GE.0) Dec_1=-Dec_1
c ..............Just test output...........
write(6,3)
+ G_2,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs,s,Vmag
write(6,'(6H Pos: 2F8.4)') RA__1,Dec_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'acro.dat' ! References of used acronyms
C Format for file interpretation
4 format(A17,2X,A1,1X,A311)
C Effective file loading
open(unit=1,file='acro.dat', status='old')
write(6,*) '....Loading file: acro.dat'
do i__=1,68
read(1,'(A332)')ar__3
read(ar__3,4)Acronym,Sep,Text
c ..............Just test output...........
write(6,4)Acronym,Sep,Text
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end