Conversion of standardized ReadMe file for
file /./ftp/cats/II/36 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.81 (2015-09-23), on 2024-Apr-19
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. II/36 UBV and uvby-beta observations of Orion OB1 Association (Warren+ 1977)
*================================================================================
*Catalogue of Individual UBV and uvby-beta Observations in the Region
*of the Orion OB 1 Association.
* Warren W. H. Jr., Hesser J. E.
* <Astrophys. J. Suppl. 34, 115 (1977)>
* =1977ApJS...34..115W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ubv.dat' ! UBV data
integer*4 nr__
parameter (nr__=357) ! Number of records
character*42 ar__ ! Full-size record
integer*4 Brun (nr__) ! Number in Brun's catalogue (1935POLyo...1...12B)
integer*4 WH (nr__) ! ? Number assigned by Warren and Hesser
* (1977ApJS...34..115W)
real*4 Vmag (nr__) ! (mag) Observed visual magnitude (1)
character*1 u_Vmag (nr__) ! [ :] certain nightly V are discordant (2)
real*4 B_V (nr__) ! (mag) Observed B-V (3)
character*1 u_B_V (nr__) ! [ :] if uncertain nightly B-V
real*4 U_B (nr__) ! (mag) ? observed U-B color (4)
character*2 u_U_B (nr__) ! [ :] if uncertain nightly U-B; two colons
* if nightly U-B very uncertain.
character*8 DATE_UBV (nr__) ! ("YY.MM.DD") [0-9.]! Date of observation in form YY.MM.DD
* (e.g. 68.12.04 = 4 December 1968).
*Note (1):
* Byte 15 is used only when a night's observations averaged
* to yield a 5 in the thousandths column.
*Note (2):
* The following stars each have a single discordant value which careful
* inspection of the original data fails to explain; while the existence of
* these discrepancies may be indicative of variability, it is more likely
* that they simply reflect an undetected error at the telescope; therefore,
* they should not be used in forming means; Brun 37 (68.12.05) Brun 62
* (68.12.06) Brun 202 (68.12.04) and Brun 244 (68.12.04).
* Star Brun 490 appears variable in V, which Walker's (1969) data weakly
* suggest too.
*Note (3):
* Byte 23 is used only when a night's observations averaged to yield a
* 5 in the thousandths column.
*Note (4):
* Byte 31 is used only when a night's observations averaged to yield
* a 5 in the thousandths column.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'uvbybeta.dat' ! uvby-beta data
integer*4 nr__1
parameter (nr__1=1595) ! Number of records
character*70 ar__1 ! Full-size record
integer*4 WH_1 (nr__1) ! Number assigned by Warren and Hesser
* (1977ApJS...34..115W)
character*3 m_WH (nr__1) ! Component identification for multiple systems
character*9 Ident (nr__1) ! Other identification (1)
real*4 Vmag_1 (nr__1) ! (mag) ? V magnitude transformed from y.
real*4 b_y (nr__1) ! (mag) ? b-y color index
real*4 m1 (nr__1) ! (mag) ? m1 ((v-b) - (b-y)) color index
real*4 c1 (nr__1) ! (mag) ? c1 ((u-v) - (v-b)) color index
character*8 DATE_uvby (nr__1) ! ("YY.MM.DD") [0-9.]? UT date of uvby observation
real*4 beta (nr__1) ! (mag) ? beta index
character*8 DATE_beta (nr__1) ! ("YY.MM.DD") [0-9.]? UT date of {beta} observation
*Note (1):
* HD = Henry Draper Catalogue
* BD = Bonner Durchmusterung
* P = Parenago 1954 (Cat. II/171)
C=============================================================================
C Loading file 'ubv.dat' ! UBV data
C Format for file interpretation
1 format(I4,1X,I3,1X,F6.3,A1,1X,F6.3,A1,1X,F6.3,A2,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'ubv.dat')
write(6,*) '....Loading file: ubv.dat'
do i__=1,357
read(1,'(A42)')ar__
read(ar__,1)
+ Brun(i__),WH(i__),Vmag(i__),u_Vmag(i__),B_V(i__),u_B_V(i__),
+ U_B(i__),u_U_B(i__),DATE_UBV(i__)
if(ar__(6:8) .EQ. '') WH(i__) = iNULL__
if(ar__(26:31) .EQ. '') U_B(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Brun(i__),WH(i__),Vmag(i__),u_Vmag(i__),B_V(i__),u_B_V(i__),
+ U_B(i__),u_U_B(i__),DATE_UBV(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'uvbybeta.dat' ! uvby-beta data
C Format for file interpretation
2 format(
+ I5,A3,1X,A9,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,A8,1X,F5.3,1X,
+ A8)
C Effective file loading
open(unit=1,status='old',file=
+'uvbybeta.dat')
write(6,*) '....Loading file: uvbybeta.dat'
do i__=1,1595
read(1,'(A70)')ar__1
read(ar__1,2)
+ WH_1(i__),m_WH(i__),Ident(i__),Vmag_1(i__),b_y(i__),m1(i__),
+ c1(i__),DATE_uvby(i__),beta(i__),DATE_beta(i__)
if(ar__1(20:25) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__1(27:32) .EQ. '') b_y(i__) = rNULL__
if(ar__1(34:39) .EQ. '') m1(i__) = rNULL__
if(ar__1(41:46) .EQ. '') c1(i__) = rNULL__
if(ar__1(57:61) .EQ. '') beta(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ WH_1(i__),m_WH(i__),Ident(i__),Vmag_1(i__),b_y(i__),m1(i__),
+ c1(i__),DATE_uvby(i__),beta(i__),DATE_beta(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end