FORTRAN Generation
(/./ftp/cats/J/AJ/165/268)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/165/268 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 2026-Oct-10
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/165/268  Radial velocity and rotation period of HIP 33069b  (Vowell+, 2023)
*================================================================================
*HIP 33609 b; An Eccentric Brown Dwarf Transiting a V=7.3 Rapidly Rotating B Star
*    Vowell N., Rodriguez J.E., Quinn S.N., Zhou G., Vanderburg A., Mann A.W.,
*    Hooton M.J., Stassun K.G., Howard S., Bieryla A., Latham D.W.,
*    Howell S.B., Guillot T., Ziegler C., Collins K.A., Carmichael T.W.,
*    Jenkins J.M., Shporer A., ABE L., Bendjoya P., Bush J.L., Buttu M.,
*    Collins K.I., Eastman J.D., Fields M.J., Gasparetto T., Gunther M.N.,
*    Kostov V.B., Kraus A.L., Lester K.V., Levine A.M., Littlefield C.,
*    Marie-Sainte W., Mekarnia D., Osborn H.P., Rapetti D., Ricker G.R.,
*    Seager S., Sefako R., Srdoc G., Suarez O., Torres G., Triaud A.H.M.J.,
*    Vanderspek R., Winn J.N.
*   <Astron. J., 165, 268 (2023)>
*   =2023AJ....165..268V
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Radial velocity measurements for HIP 33609 system

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

      real*8        BJD         ! [2458872/2459521e] Barycentric Julian Date at TBD
      real*8        RVel        ! (m/s) [26993/33897] Radial velocity
      real*4        e_RVel      ! (m/s) [591/1782] Error on RVel

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

C  Declarations for 'fig6.dat'	! Rotation periods of candidate members of MELANGE-6
                             as a function of BP-RP color

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     TIC         ! TESS Input Catalog identifier
      real*8        RAdeg       ! (deg) [91/118] Right Ascension (J2000)
      real*8        DEdeg       ! (deg) [-57/-38] Declination in decimal degrees (J2000)
      real*4        Gmag        ! (mag) [5.35/20.5] Gaia DR3 G band magnitude
      real*4        BpRp        ! (mag) [-0.13/4.02] Gaia DR3 (Bp-Rp) color
      real*4        Voff        ! (km/s) [0/5] Relative velocity offset
      real*8        Prot        ! (d) [0.19/24.2]? Rotational period
      integer*4     Q           ! ? Quality flag from Rampalli+, 2021, J/ApJ/921/167
*                             (1)
*Note (1): Flags as follows:
*    0 = A clear periodic modulation in the light curve and a strong,
*        sharp peak in the periodogram.
*    1 = Evidence of modulation in the light curve, but the structure in
*        the first half of the light curve makes it difficult to
*        confidently claim a period.
*    2 = Evidence for a long-term periodic trend likely due to systematics.
*    3 = Systematics completely dominate the light curve structure.
*   -1 = stars that have identical light curves, meaning they fall on the
*        same TESS pixel and so even if there is rotation we wouldn't
*        be able to distinguish which star it belongs to.
*   -2 = An eclipsing binary where, even though at least one of the stars
*        may be rotating, our period search kept picking up the period
*        of the binary.
*  -99 = no data or star is beyond TESS rotation period detection limits.

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

C  Loading file 'table2.dat'	! Radial velocity measurements for HIP 33609 system

C  Format for file interpretation

    1 format(F13.5,1X,F7.1,1X,F6.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,39
        read(1,'(A28)')ar__
        read(ar__,1)BJD,RVel,e_RVel
c    ..............Just test output...........
        write(6,1)BJD,RVel,e_RVel
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig6.dat'	! Rotation periods of candidate members of MELANGE-6
*                             as a function of BP-RP color

C  Format for file interpretation

    2 format(
     +  I9,1X,F11.7,1X,F11.7,1X,F6.3,1X,F6.3,1X,F5.3,1X,F12.9,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig6.dat')
      write(6,*) '....Loading file: fig6.dat'
      do i__=1,283
        read(1,'(A70)')ar__1
        read(ar__1,2)TIC,RAdeg,DEdeg,Gmag,BpRp,Voff,Prot,Q
        if(ar__1(55:66) .EQ. '') Prot = rNULL__
        if(ar__1(68:70) .EQ. '') Q = iNULL__
c    ..............Just test output...........
        write(6,2)TIC,RAdeg,DEdeg,Gmag,BpRp,Voff,Prot,Q
c    .......End.of.Just test output...........
      end do
      close(1)

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