Conversion of standardized ReadMe file for
catalog into FORTRAN code for loading all data files into arrays.
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-22
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/140 Carbon stars from Baldone telescope (Alksne+ 1987)
*================================================================================
*Cool Carbon stars found with the Baldone Schmidt Telescope - Version 1987
* Alksne Z., Alksnis A., Ozolina V., Platais I.
* <Radioastrophys. Observatory Latvian Academy of Sciences, Riga (1987)>
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'data.dat' ! The Baldone Survey of Cool Carbon Stars
integer*4 nr__
parameter (nr__=318) ! Number of records
character*109 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension 1950
real*8 Dec (nr__) ! (deg) Declination 1950
integer*4 BC (nr__) ! [1/318] The number of the star in the Survey
integer*4 RAh (nr__) ! (h) Right Ascension 1950 (hours) (1)
integer*4 RAm (nr__) ! (min) Right Ascension 1950 (minutes) (1)
real*4 RAs (nr__) ! (s) Right Ascension 1950 (seconds) (1)
character*1 DE_ (nr__) ! Declination 1950 (sign) (1)
integer*4 DEd (nr__) ! (deg) Declination 1950 (degrees) (1)
integer*4 DEm (nr__) ! (arcmin) Declination 1950 (minutes) (1)
integer*4 DEs (nr__) ! (arcsec) Declination 1950 (seconds) (1)
real*4 GLON (nr__) ! (deg) Galactic longitude
real*4 GLAT (nr__) ! (deg) Galactic latitude
integer*4 Chart (nr__) ! [1/38] Reference for identification charts
* (see References below)
character*4 SpType (nr__) ! Spectral type: R or N or R-N (2)
character*1 l_Rmag (nr__) ! [(] Upper brightness limit (3)
real*4 Rmag (nr__) ! (mag) ? Photographic red magnitude (4)
real*4 Rmag2 (nr__) ! (mag) ? Upper limit for a range (4)
real*4 Vmag (nr__) ! (mag) ? Photographic V magnitude
real*4 Vmag2 (nr__) ! (mag) ? Upper limit for a range
character*1 l_Bmag (nr__) ! [(] Upper brightness limit (3)
real*4 Bmag (nr__) ! (mag) ? Photographic B magnitude
character*1 l_Bmag2 (nr__) ! [(] Upper brightness limit (3)
real*4 Bmag2 (nr__) ! (mag) ? Upper limit for a range
character*9 Name (nr__) ! Designation of the variable star or its
* preliminary name (SVS- Soviet Variable Star)
integer*4 Ref (nr__) ! [1/38]? Reference for the photometric data
* (see References below)
character*7 CrossId (nr__) ! Cross-identification (5)
*Note (1): position as measured on direct red plates taken with the
* Schmidt Telescope and reduced using the AGK3 catalogue.
*Note (2): the spectral classification (late R or early N) depends on
* density distribution of obtained spectra, thus the faint
* Carbon stars are denoted by C.
*Note (3): the symbol "(" indicates an upper limit of brightness in the
* case of extremely faint stars.
*Note (4): photographic red magnitude, effective wavelength 0.63 microns.
*Note (5): Cross-identification between this catalogue and the others.
* The number of the star is given in cols. 103-106, and the
* appropriate reference, coded according to the References
* below, in cols. 108-109.
C=============================================================================
C Loading file 'data.dat' ! The Baldone Survey of Cool Carbon Stars
C Format for file interpretation
1 format(
+ I3,2X,I2,1X,I2,1X,F4.1,2X,A1,I2,1X,I2,1X,I2,2X,F6.2,2X,F6.2,
+ 2X,I2,1X,A4,1X,A1,F4.1,1X,F4.1,2X,F4.1,1X,F4.1,1X,A1,F4.1,1X,
+ A1,F4.1,2X,A9,2X,I2,2X,A7)
C Effective file loading
open(unit=1,file='data.dat', status='old')
write(6,*) '....Loading file: data.dat'
do i__=1,318
read(1,'(A109)')ar__
read(ar__,1)
+ BC(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),GLON(i__),GLAT(i__),Chart(i__),SpType(i__),
+ l_Rmag(i__),Rmag(i__),Rmag2(i__),Vmag(i__),Vmag2(i__),
+ l_Bmag(i__),Bmag(i__),l_Bmag2(i__),Bmag2(i__),Name(i__),
+ Ref(i__),CrossId(i__)
if(ar__(54:57) .EQ. '') Rmag(i__) = rNULL__
if(ar__(59:62) .EQ. '') Rmag2(i__) = rNULL__
if(ar__(65:68) .EQ. '') Vmag(i__) = rNULL__
if(ar__(70:73) .EQ. '') Vmag2(i__) = rNULL__
if(ar__(76:79) .EQ. '') Bmag(i__) = rNULL__
if(ar__(82:85) .EQ. '') Bmag2(i__) = rNULL__
if(ar__(99:100) .EQ. '') Ref(i__) = iNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_(i__) = RAh(i__)
if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
Dec(i__) = DEd(i__)
if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,1)
+ BC(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),
+ DEs(i__),GLON(i__),GLAT(i__),Chart(i__),SpType(i__),
+ l_Rmag(i__),Rmag(i__),Rmag2(i__),Vmag(i__),Vmag2(i__),
+ l_Bmag(i__),Bmag(i__),l_Bmag2(i__),Bmag2(i__),Name(i__),
+ Ref(i__),CrossId(i__)
write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end