Conversion of standardized ReadMe file for
file /./ftp/cats/II/235 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. II/235 The Photometric Magnitude Difference Catalog (Worley+ 2001)
*================================================================================
*The 2001 US Naval Observatory Double Star CD-ROM. IV. The Photometric
*Magnitude Difference Catalog
* Worley C.E., Mason B.D., Wycoff G.L.
* <Astron. J. 122, 3482 (2001)>
* =2001AJ....122.3482W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'deltam.dat' ! Delta-M Catalog
integer*4 nr__
parameter (nr__=19589) ! Number of records
character*86 ar__ ! Full-size record
C J2000 position composed of: RAh DE- DEd DEm
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh ! (h) *Right Ascension J2000 (hours)
integer*4 RAdm ! (0.1min) *Right Ascension J2000 (minutes)
character*1 DE_ ! *Declination J2000 (sign)
integer*4 DEd ! (deg) *Declination J2000 (degrees)
integer*4 DEm ! (arcmin) *Declination J2000 (minutes)
character*7 DiscName ! *Discoverer Code & Number
character*5 Comp ! *Component Identification
real*4 dm ! (mag) Magnitude Difference
integer*4 o_dm ! ?Number of observations used to
* compute dm
character*8 r_dm ! *Reference Code
character*1 Meth ! *Observational Method
character*43 Notes ! *?Notes
*Note on RAh, RAdm, DE-, DEd, DEm:
* The positions given represent our best estimates of these values. Where
* possible, these are based on the WDS data, with proper motion incorporated.
*Note on DiscName:
* The discoverer, identified by a one-to-three letter code, and the
* discoverer's number, if assigned. The reference list ("refs.dat")
* should be consulted first to identify individuals.
*Note on Comp:
* Components, when the object has more than two. The WDS 2001 scheme is
* maintained.
*Note on r_dm:
* The NewCode from refs.dat file. The first three letters corresponds to the
* WDS discoverer while the next four numbers indicate the date of publication.
*Note on Meth:
* The various methods are detailed below, and in "method.txt" file;
* statistics on the different methods can be found in file "table1.txt".
* There are a plethora of interferometric techniques which are distinguished
* due to the familiarity of the authors with same. For specifics, consult
* the appropriate reference cited in the references.
* A = Adaptive Optics.
* B = Long-baseline interferometer visibility amplitude fitting.
* C = Visual measures with a double-image micrometer.
* D = speckle interferometric comparison of DVA peak heights of stars with
* known delta-m and approximately the same brightness.
* E = speckle interferometric visibility curve-fitting.
* F = speckle interferometric fork.
* G = Visual estimate made with an objective grating.
* H = Wickes and Dicke interferometer - rotating prism and high-frequency
* chopping wheel.
* I = Integrated image taken with speckle camera. No autocorrelation is
* done, and the scaled peak heights are compared.
* J = CCD.
* K = Focal-grating photometer.
* L = speckle interferometric relative intensity of autocorrelation.
* M = IR speckle interferometric shift-and-add.
* N = Photocell with diaphragm in focal plane.
* O = Lunar Occultation.
* P = Photographic plates.
* Q = Tycho photometry of Hipparcos Stars
* R = Area scanner.
* S = speckle interferometric modified shift-and-add.
* T = speckle interferometric triple correlation or bispectrum.
* V = Photoelectric photometry.
* W = Wedge photometer.
* Z = Visual photometry with a double-image or polarizing photometer.
*Note on Notes:
* This may contain specifics of the filter used (both effective wavelength and
* FWHM in nanometers), errors (i.e., ``SIGMA = ''), colors, and indications of
* uncertainty, variability when a component is not in the WDS, or when a
* reference is to an additional publication.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! Reference file
integer*4 nr__1
parameter (nr__1=4195) ! Number of records
character*338 ar__1 ! Full-size record
character*8 NewCode ! *New reference code used in Orbit,
* Speckle and Delta-M Catalogs
character*3 OldCode ! Discoverer Code
integer*4 RefNo ! *?Reference number of publication
character*37 Obs ! Observer(s)
character*284 Publ ! Publication reference for measures
*Note on NewCode:
* The first three letters corresponds to the WDS discoverer while the next four
* numbers indicate the date of publication. New Code = XYZ1999x
*Note on RefNo:
* 99 = Unpublished, manuscript or reference not found.
C=============================================================================
C Loading file 'deltam.dat' ! Delta-M Catalog
C Format for file interpretation
1 format(I2,I3,A1,I2,I2,1X,A7,A5,F5.2,2X,I2,1X,A8,1X,A1,A43)
C Effective file loading
open(unit=1,status='old',file=
+'deltam.dat')
write(6,*) '....Loading file: deltam.dat'
do i__=1,19589
read(1,'(A86)')ar__
read(ar__,1)
+ RAh,RAdm,DE_,DEd,DEm,DiscName,Comp,dm,o_dm,r_dm,Meth,Notes
if(ar__(31:32) .EQ. '') o_dm = iNULL__
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(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ RAh,RAdm,DE_,DEd,DEm,DiscName,Comp,dm,o_dm,r_dm,Meth,Notes
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! Reference file
C Format for file interpretation
2 format(A8,2X,A3,I2,2X,A37,A284)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,4195
read(1,'(A338)')ar__1
read(ar__1,2)NewCode,OldCode,RefNo,Obs,Publ
if(ar__1(14:15) .EQ. '') RefNo = iNULL__
c ..............Just test output...........
write(6,2)NewCode,OldCode,RefNo,Obs,Publ
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end