FORTRAN Generation
(/./ftp/cats/J/AJ/131/2687)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/131/2687 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-Mar-28
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/131/2687       Speckle interferometry at USNO. XI       (Mason+, 2006)
*================================================================================
*Speckle interferometry at the US naval observatory. XI.
*    Mason B.D., Hartkopf W.I., Wycoff G.L., Rafferty T.J.
*   <Astron. J., 131, 2687-2694 (2006)>
*   =2006AJ....131.2687M
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! Speckle Interferometric Measurements of Double Stars

      integer*4 nr__
      parameter (nr__=702)	! Number of records
      character*57 ar__   	! Full-size record

      character*10  WDS         ! WDS name (HHMMm+DDMM, J2000)
      character*7   Disc        ! Discoverer designation
      character*4   Comp        ! Additional system components
      real*8        Epoch       ! (yr) Epoch of observation
      real*4        PA          ! (deg) Position angle
      character*1   u_PA        ! [:] Uncertainty flag for PA
      real*4        Sep         ! (arcsec) Separation
      character*1   u_Sep       ! [:] Uncertainty flag for Sep
      integer*4     o_Sep       ! Number of observations contained in the mean
      character*7   Note        ! Object notes (1)
*Note (1): Object notes as follows:
*      C = Confirming Observation
*      L = Linear elements determined. See Hartkopf et al.
*          (2007AJ....134.1671M)
*      a = Large change in separation and/or angle of position
*      b = Much better measure than the uncertain one of
*          Mason et al. (2001, Cat. J/AJ/122/1586)
*      N = 50-184 : Number of years since last measure

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! Binaries not found

      integer*4 nr__1
      parameter (nr__1=7)	! Number of records
      character*44 ar__1  	! Full-size record

      character*10  WDS_1       ! WDS name (HHMMm+DDMM, J2000)
      character*8   Disc_1      ! Discoverer designation
      integer*4     Date        ! (yr) Date
      integer*4     PA_1        ! (deg) Published position angle
      real*4        Sep_1       ! (arcsec) Published separation
      real*4        mag1        ! (mag) Published primary magnitude
      real*4        mag2        ! (mag) Published secondary magnitude
      character*1   Note_1      ! [1] Inaccurate position (1)
*Note (1):  1 = Only arcminute coordinates are available for this pair,
*     resolved by Herschel (1833, MmRAS, 6, 1), although he could
*     "not define it. Night not good; telescopes give double images."
*     Probably an erroneous detection.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! Speckle Interferometric Measurements and Residuals
                             to Systems with Orbits

      integer*4 nr__2
      parameter (nr__2=116)	! Number of records
      character*127 ar__2  	! Full-size record

      character*10  WDS_2       ! WDS name (HHMMm+DDMM, J2000)
      character*8   Disc_2      ! Discoverer designation
      character*4   Comp_1      ! Additional system components
      real*8        Epoch_1     ! (yr) ? Epoch of observation
      real*4        PA_2        ! (deg) ? Position angle
      real*4        Sep_2       ! (arcsec) ? Separation
      integer*4     o_Sep_1     ! ? Number of observations contained in the mean
      real*4        d_PA        ! (deg) The O-C position angle residual
      real*4        d_Sep       ! (arcsec) The O-C separation residual
      character*3   n_Ref       ! Notes on Ref (1)
      character*63  Ref         ! Observation reference
*Note (1): Notes as follows:
*      * = System used in characterizing errors.
*      1 = This orbit might possibly be fit better with linear elements.
*          See Hartkopf et al. (2007AJ....134.1671M).
*      2 = This measure was inconsistent with previous measures, so not
*          included in Mason et al. (2004, Cat. <J/AJ/127/539>). However,
*          available data are deemed insufficient for a new orbital
*          calculation at this time.
*      3 = This system was expected to show significant motion over the
*          calendar year, so multiple observations have been obtained.
*      4 = The new orbit is listed in Table 4, ephemerides based on these
*          elements are listed in Table 6, and the orbit is illustrated
*          in Figure 2.
*      5 = This orbit calculates the opposite quadrant for the position
*          angle, however, it is clearly in this quadrant with an obvious
*          {Delta}m of ~1.5.

C=============================================================================

C  Loading file 'table1.dat'	! Speckle Interferometric Measurements of Double Stars

C  Format for file interpretation

    1 format(
     +  A10,2X,A7,1X,A4,1X,F8.3,1X,F5.1,A1,1X,F5.2,A1,1X,I1,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,702
        read(1,'(A57)')ar__
        read(ar__,1)WDS,Disc,Comp,Epoch,PA,u_PA,Sep,u_Sep,o_Sep,Note
c    ..............Just test output...........
        write(6,1)WDS,Disc,Comp,Epoch,PA,u_PA,Sep,u_Sep,o_Sep,Note
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table2.dat'	! Binaries not found

C  Format for file interpretation

    2 format(A10,1X,A8,1X,I4,1X,I3,1X,F3.1,1X,F4.1,1X,F4.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,7
        read(1,'(A44)')ar__1
        read(ar__1,2)WDS_1,Disc_1,Date,PA_1,Sep_1,mag1,mag2,Note_1
c    ..............Just test output...........
        write(6,2)WDS_1,Disc_1,Date,PA_1,Sep_1,mag1,mag2,Note_1
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3.dat'	! Speckle Interferometric Measurements and Residuals
*                             to Systems with Orbits

C  Format for file interpretation

    3 format(
     +  A10,2X,A8,1X,A4,1X,F8.3,1X,F5.1,1X,F4.2,1X,I2,1X,F4.1,1X,F5.2,
     +  1X,A3,1X,A63)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,116
        read(1,'(A127)')ar__2
        read(ar__2,3)
     +  WDS_2,Disc_2,Comp_1,Epoch_1,PA_2,Sep_2,o_Sep_1,d_PA,d_Sep,
     +  n_Ref,Ref
        if(ar__2(27:34) .EQ. '') Epoch_1 = rNULL__
        if(ar__2(36:40) .EQ. '') PA_2 = rNULL__
        if(ar__2(42:45) .EQ. '') Sep_2 = rNULL__
        if(ar__2(47:48) .EQ. '') o_Sep_1 = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  WDS_2,Disc_2,Comp_1,Epoch_1,PA_2,Sep_2,o_Sep_1,d_PA,d_Sep,
     +  n_Ref,Ref
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end