FORTRAN Generation
(/./ftp/cats/J/A_A/382/92)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/382/92 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-20
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/A+A/382/92           Sco OB2 association                      (Shatsky+, 2002)
*================================================================================
*The mass ratio distribution of B-type visual binaries in the Sco OB2 association
*    Shatsky N., Tokovinin A.
*   <Astron. Astrophys. 382, 92 (2002)>
*   =2002A&A...382...92S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Target star photometry

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

      integer*4     HD          ! HD number of target star
      character*8   SpType      ! MK spectral classification (from Simbad)
      real*4        Vmag        ! (mag) V magnitude (from Simbad)
      real*4        Av          ! (mag) Interstellar extinction in the visible (1)
      real*4        Kmag        ! (mag) K-band (2.2um) magnitude (2)
      real*4        e_Kmag      ! (mag) ? Error of Kmag (3)
      real*4        J_K         ! (mag) ? (J-K) color index (absent for HD 143275)
      real*4        e_J_K       ! (mag) ? Error of (J-K) (4)
      integer*4     Nobs        ! Number of observations (pointings)
      integer*4     Nphys       ! Number of physical companions around
      integer*4     Nopt        ! Number of optical companions around
      integer*4     Nsky        ! Number of stars in sky-offset field
      character*3   Flag        ! [ANTcjs+] flags (5)
*Note (1):
*    Taken from deGeus et al. (1989A&A...216...44D) or, if `A' is given in
*    the field Flag, estimated from (B-V) color excess.
*Note (2):
*    For non-observed objects the theoretical K magnitude is given as
*    computed from SpType, Vmag and Av
*Note (3):
*    Error is a deviation of observation results from the mean if Nobs>1 or
*    difference between magnitudes by two halves of data cube if Nobs=1.
*    Not given for non-photometric condition observations (`c' or `+' in
*    Flag) or for objects which weren't observed explicitly (`N' in Flag)
*Note (4): Error is quadratically added e_Kmag and e_Jmag. See note (3) also.
*Note (5): Flags meaning:
*    A - extinction is estimated from color excess E(B-V)
*    N - not observed explicitly
*    T - object studied in Tokovinin et al., 1999A&A...346..481T
*    c - non-photometric conditions
*    + - "likely" photometric conditions: K and J-K don't deviate much from 
*        expected but night is not photometric
*    j - joint magnitudes of the very close pair are given
*    s - a few pixels of detector in the core of image were near saturation limit

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

C  Declarations for 'table2.dat'	! Secondaries photometry and positions

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

      integer*4     HD_1        ! HD number of target star
      character*2   ID          ! Component identifier (1)
      real*4        Kmag_1      ! (mag) K-band magnitude as K(primary) + delta(K)
      real*4        e_Kmag_1    ! (mag) ? Error: deviation of delta(K) from mean
      real*4        J_K_1       ! (mag) ? (J-K) color as (J-K)(primary) + delta(J-K)
      real*4        e_J_K_1     ! (mag) ? Error: deviation of delta(J-K) from mean
      character*2   Flag_1      ! [c+!?: ] flags (2)
      real*4        Rho         ! (arcsec) Separation from primary
      real*4        e_Rho       ! (mas) Error of Rho (in milliarcseconds)
      real*4        PA          ! (deg) position angle
      real*4        e_PA        ! (deg) Error of PA
      integer*4     Nobs_1      ! Number of observations (pointings)
      character*2   Status      ! [PO? ] `P'-physical,`O'-optical (`?'-uncertain)
*Note (1):
*    Taken from the literature for known companions, assigned as `P' or `Q'
*    for new likely physical companions, or assigned numerically as
*    `1',`2',... in increasing order of separation for opticals. Followed
*    by "?" for (two) new companions with uncertain detection.
*Note (2): Flags meaning:
*    c - non-photometric conditions
*    + - "likely" photometric conditions: K and J-K of primaries don't 
*        deviate much from expected but night is not photometric
*    ! - new physical companion
*    : - pair is too close for reliable measurements

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

C  Declarations for 'table3.dat'	! Background sources data

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

      integer*4     HD_2        ! HD number of target star
      character*1   ID_1        ! [abcdefghi] background star identifier
      real*4        Jmag        ! (mag) ? J-band (1.2mu) magnitude
      real*4        e_Jmag      ! (mag) ? error: spread of different observations
      real*4        Kmag_2      ! (mag) K-band (2.2mu) magnitude
      real*4        e_Kmag_2    ! (mag) ? error: spread of different observations
      character*1   Flag_2      ! [c+] flag (1)
      real*4        delRA       ! (arcsec) Coordinate difference with target in
*                                     right ascension direction
      real*4        delDE       ! (arcsec) Coordinate difference  with target in
*                                     declination direction
*Note (1): Flags meaning:
*    c: non-photometric conditions
*    +: "likely" photometric conditions: K and J-K of primaries don't 
*       deviate much from expected but the night is not photometric

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

C  Loading file 'table1.dat'	! Target star photometry

C  Format for file interpretation

    1 format(
     +  I6,1X,A8,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
     +  I1,1X,I1,1X,I1,1X,I1,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,115
        read(1,'(A58)')ar__
        read(ar__,1)
     +  HD,SpType,Vmag,Av,Kmag,e_Kmag,J_K,e_J_K,Nobs,Nphys,Nopt,Nsky,
     +  Flag
        if(ar__(32:35) .EQ. '') e_Kmag = rNULL__
        if(ar__(37:41) .EQ. '') J_K = rNULL__
        if(ar__(43:46) .EQ. '') e_J_K = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  HD,SpType,Vmag,Av,Kmag,e_Kmag,J_K,e_J_K,Nobs,Nphys,Nopt,Nsky,
     +  Flag
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Secondaries photometry and positions

C  Format for file interpretation

    2 format(
     +  I6,1X,A2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A2,1X,F6.3,1X,
     +  F5.3,1X,F5.1,1X,F3.1,1X,I1,1X,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,96
        read(1,'(A63)')ar__1
        read(ar__1,2)
     +  HD_1,ID,Kmag_1,e_Kmag_1,J_K_1,e_J_K_1,Flag_1,Rho,e_Rho,PA,
     +  e_PA,Nobs_1,Status
        if(ar__1(17:20) .EQ. '') e_Kmag_1 = rNULL__
        if(ar__1(22:26) .EQ. '') J_K_1 = rNULL__
        if(ar__1(28:31) .EQ. '') e_J_K_1 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  HD_1,ID,Kmag_1,e_Kmag_1,J_K_1,e_J_K_1,Flag_1,Rho,e_Rho,PA,
     +  e_PA,Nobs_1,Status
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Background sources data

C  Format for file interpretation

    3 format(
     +  I6,1X,A1,1X,F4.1,1X,F3.1,1X,F4.1,1X,F3.1,1X,A1,1X,F4.1,1X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,97
        read(1,'(A38)')ar__2
        read(ar__2,3)
     +  HD_2,ID_1,Jmag,e_Jmag,Kmag_2,e_Kmag_2,Flag_2,delRA,delDE
        if(ar__2(10:13) .EQ. '') Jmag = rNULL__
        if(ar__2(15:17) .EQ. '') e_Jmag = rNULL__
        if(ar__2(24:26) .EQ. '') e_Kmag_2 = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  HD_2,ID_1,Jmag,e_Jmag,Kmag_2,e_Kmag_2,Flag_2,delRA,delDE
c    .......End.of.Just test output...........
      end do
      close(1)

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