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. J/AJ/105/499 OB Associations and Spiral Structure in M33 (Regan+ 1993)
*================================================================================
*OB Associations and Spiral Structure in the Southern Arm of M33
* Regan M.W., Wilson C.D.
* <Astron. J. 105, 499 (1993)>
* =1993AJ....105..499R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! The observations
integer*4 nr__
parameter (nr__=2018) ! Number of records
character*87 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension 1950
real*8 Dec ! (deg) Declination 1950
integer*4 ID ! [1/2018]+ Star ID
integer*4 RAh ! (h) Right Ascension (1950)
integer*4 RAm ! (min) Right Ascension (1950)
real*4 RAs ! (s) Right Ascension (1950)
integer*4 X ! (pix) X Coordinate
real*4 V ! (mag) V band Magnitude
real*4 B ! (mag) B band Magnitude
real*4 U ! (mag) ? U band Magnitude (blank if undetermined)
real*4 B_V ! (mag) B-V Magnitude
real*4 U_B ! (mag) ? U-B Magnitude (blank if undetermined)
integer*4 ass ! [22/58]? Association number containing star
character*1 DE_ ! Sign of Declination (1950)
integer*4 DEd ! (deg) [30] Declination (1950)
integer*4 DEm ! (arcmin) Declination (1950)
integer*4 DEs ! (arcsec) Declination (1950)
integer*4 Y ! (pix) Y Coordinate
real*4 e_V ! (mag) V uncertainty
real*4 e_B ! (mag) B uncertainty
real*4 e_U ! (mag) ? U uncertainty (blank if undetermined)
C=============================================================================
C Loading file 'table1.dat' ! The observations
C Format for file interpretation
1 format(
+ 1X,I4,2X,I1,1X,I2,1X,F4.1,1X,I4,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,I4,2X,F4.2,2X,F4.2,
+ 2X,F4.2)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,2018
read(1,'(A87)')ar__
read(ar__,1)
+ ID,RAh,RAm,RAs,X,V,B,U,B_V,U_B,ass,DE_,DEd,DEm,DEs,Y,e_V,e_B,
+ e_U
if(ar__(35:39) .EQ. '') U = rNULL__
if(ar__(47:51) .EQ. '') U_B = rNULL__
if(ar__(53:54) .EQ. '') ass = iNULL__
if(ar__(84:87) .EQ. '') e_U = rNULL__
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(DEs .GE. 0) Dec=Dec+DEs/3600.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)
+ ID,RAh,RAm,RAs,X,V,B,U,B_V,U_B,ass,DE_,DEd,DEm,DEs,Y,e_V,e_B,
+ e_U
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end