Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/151/153 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. J/AJ/151/153 Speckle interferometry at SOAR in 2015 (Tokovinin+, 2016)
*================================================================================
*Speckle interferometry at SOAR in 2015.
* Tokovinin A., Mason B.D., Hartkopf W.I., Mendez R.A., Horch E.P.
* <Astron. J., 151, 153-153 (2016)>
* =2016AJ....151..153T (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Measurements of double stars at SOAR
integer*4 nr__
parameter (nr__=1303) ! Number of records
character*111 ar__ ! Full-size record
character*10 WDS (nr__) ! WDS identification (from cat. B/wds)
character*9 Disc (nr__) ! Discoverer designation (as adopted in the WDS)
character*5 Comp (nr__) ! Components
integer*4 HIP (nr__) ! ? Hipparcos number (cat. I/239)
character*1 m_HIP (nr__) ! [BC] Multiplicity index on HIP
character*10 OName (nr__) ! Other name
real*8 Epoch (nr__) ! (yr) [15.1/15.92] Besselian epoch of observation
character*1 Flt (nr__) ! [HIVy] Filter used: V, H, I=near-infrared I
* filter (788/132nm), or y=Stromgren y filter
* (543/22nm)
integer*4 No (nr__) ! [1/6] Number of observations used
real*4 PA (nr__) ! (deg) Position angle {theta}
real*4 Err (nr__) ! (mas) [0/28] Internal measurement error in
* tangential direction {rho}{sigma}_{theta}_
real*4 Sep (nr__) ! (arcsec) [0.012/3.37] Separation {rho}
real*4 e_Sep (nr__) ! (mas) [0/40] The 1{sigma} uncertainty in Sep
* {sigma}{rho}
real*4 Dmag (nr__) ! (mag) [0/6.8] Magnitude difference of pair {Delta}m
character*1 f_Dmag (nr__) ! [*:q] Flag on Dmag (1)
real*4 v_O_C_t (nr__) ! (deg) [-75.5/265.2]? The (O-C) in PA [O-C]_{theta}_
real*4 v_O_C_r (nr__) ! (arcsec) [-1.1/0.32]? The (O-C) in Sep [O-C]_{rho}_
character*8 Ref (nr__) ! Reference code to the orbit adopted in the
* Sixth Catalog of Orbits of Visual Binary Stars
* (Hartkopf et al. 2001AJ....122.3472H) also
* available at
* http://ad.usno.navy.mil/wds/orb6/wdsref.html);
* in refs.dat file
*Note (1): The flags are defined as follows:
* * = Dmag and the true quadrant are determined from the resolved
* long-exposure image;
* : = Data are noisy and {Delta}m is likely over-estimated (see Tokovinin et
* al. 2010, Cat. J/AJ/139/743);
* q = True quadrant is established from the shift-and-add (SAA or "lucky")
* image.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Unresolved stars
integer*4 nr__1
parameter (nr__1=360) ! Number of records
character*82 ar__1 ! Full-size record
character*10 WDS_1 (nr__1) ! WDS identification (from cat. B/wds) (1)
character*10 Disc_1 (nr__1) ! Discoverer designation (or other name for
* objects without discoverer designations)
character*5 Comp_1 (nr__1) ! Components
integer*4 HIP_1 (nr__1) ! [93/118319]? Hipparcos number (cat. I/239)
character*1 m_HIP_1 (nr__1) ! [BC] Multiplicity index on HIP
character*10 OName_1 (nr__1) ! Other name (mostly HD identifier)
real*8 Epoch_1 (nr__1) ! (yr) [15.1/15.92] Epoch of observation
character*1 Flt_1 (nr__1) ! [IVy] Filter used (V, I=near-infrared I filter
* (788/132nm), or y=Stromgren y filter
* (543/22nm)
integer*4 No_1 (nr__1) ! [1/4] Number of observations used
real*4 MSep (nr__1) ! (arcsec) [0.02/0.06] Mininimum resolvable separation
* {rho}_min_ (2)
real*4 Dm0_15 (nr__1) ! (mag) [1.9/5.7] The 5{sigma} magnitude difference
* {Delta}m at 0.15'' (3)
real*4 Dm1 (nr__1) ! (mag) [2/7.6] The 5{sigma} magnitude difference
* {Delta}m at 1'' (3)
character*1 u_Dm1 (nr__1) ! [:] Flag for noisy data that are mostly
* associated with faint stars (4)
*Note (1): For stars that do not have entries in the WDS, fictitious WDS-style
* codes based on the J2000 position are listed.
*Note (2): The estimated resolution limit equal to the diffraction radius
* {lambda}/D.
*Note (3): The 5{sigma} detection limit {Delta}m_5_ at 0.15'' and 1'' separations
* determined by the procedure described in Tokovinin et al. 2010
* (Cat. J/AJ/139/743). When two or more data cubes are processed, the largest
* {Delta}m value is listed.
*Note (4): In such cases, the quoted {Delta}m might be too large (optimistic);
* however, the knowledge that these stars were observed and no companions
* were found is still useful for statistics.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Newly resolved pairs
integer*4 nr__2
parameter (nr__2=27) ! Number of records
character*89 ar__2 ! Full-size record
character*10 WDS_2 (nr__2) ! WDS identification (from cat. B/wds)
character*8 Disc_2 (nr__2) ! Discoverer designation (1)
character*5 Comp_2 (nr__2) ! Components
integer*4 HIP_2 (nr__2) ! ? Hipparcos number (cat. I/239)
character*1 m_HIP_2 (nr__2) ! [B] Multiplicity index on HIP
character*11 OName_2 (nr__2) ! Other name
real*8 Epoch_2 (nr__2) ! [15.16/15.92] Epoch of observation
character*1 Flt_2 (nr__2) ! [Iy] Filter used: I=near-infrared I filter
* (788/132nm), or y=Stromgren y filter (543/22nm)
real*4 PA_1 (nr__2) ! (deg) [2.4/322.6] Position angle {theta}
real*4 Sep_1 (nr__2) ! (arcsec) [0.021/1.78] Separation {rho}
character*5 Dmag_1 (nr__2) ! (mag) Magnitude difference of pair {Delta}m
character*5 SpT (nr__2) ! Spectral type (2)
real*4 plx (nr__2) ! (mas) [3.5/29]? Hipparcos parallax (van Leeuwen, 2007,
* Cat. I/311) {pi}_HIP2_
*Note (1): For some multiple systems, we used existing discover codes and simply
* added new component designations.
*Note (2): As given in SIMBAD or estimated from absolute magnitude.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__3
parameter (nr__3=132) ! Number of records
character*158 ar__3 ! Full-size record
character*8 Ref_1 (nr__3) ! Reference code
character*19 BibCode (nr__3) ! Bibliographic code
character*80 Aut (nr__3) ! Author's name
character*48 Com (nr__3) ! Comments
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Notes
integer*4 nr__4
parameter (nr__4=26) ! Number of records
character*1213 ar__4 ! Full-size record
character*10 WDS_3 (nr__4) ! WDS identification (from cat. B/wds)
character*1202 Note (nr__4) ! Notes
C=============================================================================
C Loading file 'table2.dat' ! Measurements of double stars at SOAR
C Format for file interpretation
1 format(
+ A10,1X,A9,A5,5X,I6,A1,1X,A10,1X,F7.4,1X,A1,1X,I1,1X,F5.1,1X,
+ F4.1,1X,F6.4,1X,F4.1,1X,F3.1,1X,A1,1X,F5.1,1X,F6.3,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,1303
read(1,'(A111)')ar__
read(ar__,1)
+ WDS(i__),Disc(i__),Comp(i__),HIP(i__),m_HIP(i__),OName(i__),
+ Epoch(i__),Flt(i__),No(i__),PA(i__),Err(i__),Sep(i__),
+ e_Sep(i__),Dmag(i__),f_Dmag(i__),v_O_C_t(i__),v_O_C_r(i__),
+ Ref(i__)
if(ar__(31:36) .EQ. '') HIP(i__) = iNULL__
if(ar__(91:95) .EQ. '') v_O_C_t(i__) = rNULL__
if(ar__(97:102) .EQ. '') v_O_C_r(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ WDS(i__),Disc(i__),Comp(i__),HIP(i__),m_HIP(i__),OName(i__),
+ Epoch(i__),Flt(i__),No(i__),PA(i__),Err(i__),Sep(i__),
+ e_Sep(i__),Dmag(i__),f_Dmag(i__),v_O_C_t(i__),v_O_C_r(i__),
+ Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Unresolved stars
C Format for file interpretation
2 format(
+ A10,1X,A10,1X,A5,5X,I6,A1,1X,A10,1X,F7.4,1X,A1,1X,I1,1X,F5.3,
+ 1X,F5.2,1X,F5.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,360
read(1,'(A82)')ar__1
read(ar__1,2)
+ WDS_1(i__),Disc_1(i__),Comp_1(i__),HIP_1(i__),m_HIP_1(i__),
+ OName_1(i__),Epoch_1(i__),Flt_1(i__),No_1(i__),MSep(i__),
+ Dm0_15(i__),Dm1(i__),u_Dm1(i__)
if(ar__1(33:38) .EQ. '') HIP_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ WDS_1(i__),Disc_1(i__),Comp_1(i__),HIP_1(i__),m_HIP_1(i__),
+ OName_1(i__),Epoch_1(i__),Flt_1(i__),No_1(i__),MSep(i__),
+ Dm0_15(i__),Dm1(i__),u_Dm1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Newly resolved pairs
C Format for file interpretation
3 format(
+ A10,1X,A8,1X,A5,5X,I6,A1,1X,A11,1X,F7.4,1X,A1,1X,F5.1,1X,F6.4,
+ 1X,A5,1X,A5,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,27
read(1,'(A89)')ar__2
read(ar__2,3)
+ WDS_2(i__),Disc_2(i__),Comp_2(i__),HIP_2(i__),m_HIP_2(i__),
+ OName_2(i__),Epoch_2(i__),Flt_2(i__),PA_1(i__),Sep_1(i__),
+ Dmag_1(i__),SpT(i__),plx(i__)
if(ar__2(31:36) .EQ. '') HIP_2(i__) = iNULL__
if(ar__2(86:89) .EQ. '') plx(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ WDS_2(i__),Disc_2(i__),Comp_2(i__),HIP_2(i__),m_HIP_2(i__),
+ OName_2(i__),Epoch_2(i__),Flt_2(i__),PA_1(i__),Sep_1(i__),
+ Dmag_1(i__),SpT(i__),plx(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
4 format(A8,1X,A19,1X,A80,1X,A48)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,132
read(1,'(A158)')ar__3
read(ar__3,4)Ref_1(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,4)Ref_1(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'notes.dat' ! Notes
C Format for file interpretation
5 format(A10,1X,A1202)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,26
read(1,'(A1213)')ar__4
read(ar__4,5)WDS_3(i__),Note(i__)
c ..............Just test output...........
write(6,5)WDS_3(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end