Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/153/212 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-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. J/AJ/153/212 Observations of binary stars at the WIYN telescope (Horch+, 2017)
*================================================================================
*Observations of binary stars with the Differential Speckle Survey Instrument.
*VII. Measures from 2010 September to 2012 February at the WIYN telescope.
* Horch E.P., Casetti-Dinescu D.I., Camarata M.A., Bidarian A.,
* van Altena W.F., Sherry W.H., Everett M.E., Howell S.B., Ciardi D.R.,
* Henry T.J., Nusdeo D.A., Winters J.G.
* <Astron. J., 153, 212 (2017)>
* =2017AJ....153..212H (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'stars.dat' ! Star list and Simbad positions
integer*4 nr__
parameter (nr__=1058) ! Number of records
character*36 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*1 Type ! [ml] Type of data (1)
character*10 Star ! Star name
integer*4 RAh ! (h) Simbad Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Simbad Minute of Right Ascension (J2000)
real*4 RAs ! (s) Simbad Second of Right Ascension (J2000)
character*1 DE_ ! Simbad Sign of the Declination (J2000)
integer*4 DEd ! (deg) Simbad Degree of Declination (J2000)
integer*4 DEm ! (arcmin) Simbad Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Simbad Arcsecond of Declination (J2000)
*Note (1): Type as follows:
* m = Double star speckle measures (table 2);
* l = 5{sigma} detection limits for unresolved stars (table 4).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Double star speckle measures
integer*4 nr__1
parameter (nr__1=2521) ! Number of records
character*80 ar__1 ! Full-size record
character*10 WDS ! Washington Double Star catalog (Cat. I/237)
* designation (HHMMm+DDMM, J2000)
character*11 OName ! Other designation
character*13 Disc ! Binary discover designation
integer*4 HIP ! [170/5000024]? Hipparcos catalogue number
* (Cat. I/239)
real*8 Date ! (yr) [10.7115/12.105] Besselian date of the
* observation (offset 2000)
real*4 PA ! (deg) [0/359.7] Position angle of the secondary star
* relative to the primary (1)
real*4 Sep ! (arcsec) [0.0101/1.7551] Separation of the two stars
real*4 Dmag ! (mag) [0/5.65]? Magnitude difference of the pair in
* the filter used
integer*4 lambda ! (nm) [88/880] Center wavelength of the filter used
integer*4 FWHM ! (nm) [40/50] Full width at half maximum of the filter
* transmission
character*3 fCat ! Catalog entry flag(s) (2)
*Note (1): North through East defining the positive sense of PA.
*Note (2): Catalog entry flag as follows:
* a = Quadrant ambiguous;
* b = Quadrant inconsistent with previous measures in the 4th Interferometric
* Catalog;
* c = These measures were previously published in Horch et al.
* (2011AJ....141..180H); this reanalysis should replace the previous values;
* d = These measures were previously published in Appourchaux et al.
* (2015A&A...582A..25A); this reanalysis should replace the previous values;
* K = This star is a member of the nearby K-dwarf sample discussed in the text.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! 5{sigma} detection limits for unresolved stars
integer*4 nr__2
parameter (nr__2=214) ! Number of records
character*90 ar__2 ! Full-size record
character*10 WDS_1 ! Washington Double Star catalog (Cat. I/237)
* designation (HHMMm+DDMM, J2000)
integer*4 HIP_1 ! [344/118077]? Hipparcos catalogue number
* (Cat. I/239)
real*8 Date_1 ! (yr) [10.7116/12.1049] Besselian date of the
* observation (offset 2000)
integer*4 lambdaA ! (nm) [692] Center wavelength of the 692 nm filter
* (Channel A) used
real*4 DlA_0_2 ! (mag) [2.59/4.76] 5{sigma} detection limit at
* 0.2 arcsec
real*4 DlA_1_0 ! (mag) [3.47/8.58] 5{sigma} detection limit at
* 1.0 arcsec
integer*4 lambdaB ! (nm) [562/880] Center wavelength of the filter
* (Channel B) used
real*4 DlB_0_2 ! (mag) [2.17/4.68] 5{sigma} detection limit at
* 0.2 arcsec
real*4 DlB_1_0 ! (mag) [3.08/7.84] 5{sigma} detection limit at
* 1.0 arcsec
character*36 Com ! Comment on source
C=============================================================================
C Loading file 'stars.dat' ! Star list and Simbad positions
C Format for file interpretation
1 format(A1,1X,A10,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'stars.dat')
write(6,*) '....Loading file: stars.dat'
do i__=1,1058
read(1,'(A36)')ar__
read(ar__,1)Type,Star,RAh,RAm,RAs,DE_,DEd,DEm,DEs
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(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)Type,Star,RAh,RAm,RAs,DE_,DEd,DEm,DEs
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Double star speckle measures
C Format for file interpretation
2 format(
+ A10,1X,A11,1X,A13,1X,I6,1X,F7.4,1X,F5.1,1X,F6.4,1X,F4.2,1X,I3,
+ 1X,I2,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,2521
read(1,'(A80)')ar__1
read(ar__1,2)
+ WDS,OName,Disc,HIP,Date,PA,Sep,Dmag,lambda,FWHM,fCat
if(ar__1(38:43) .EQ. '') HIP = iNULL__
if(ar__1(66:69) .EQ. '') Dmag = rNULL__
c ..............Just test output...........
write(6,2)
+ WDS,OName,Disc,HIP,Date,PA,Sep,Dmag,lambda,FWHM,fCat
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! 5{sigma} detection limits for unresolved stars
C Format for file interpretation
3 format(
+ A10,1X,I6,1X,F7.4,1X,I3,1X,F4.2,1X,F4.2,1X,I3,1X,F4.2,1X,F4.2,
+ 1X,A36)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,214
read(1,'(A90)')ar__2
read(ar__2,3)
+ WDS_1,HIP_1,Date_1,lambdaA,DlA_0_2,DlA_1_0,lambdaB,DlB_0_2,
+ DlB_1_0,Com
if(ar__2(12:17) .EQ. '') HIP_1 = iNULL__
c ..............Just test output...........
write(6,3)
+ WDS_1,HIP_1,Date_1,lambdaA,DlA_0_2,DlA_1_0,lambdaB,DlB_0_2,
+ DlB_1_0,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end