FORTRAN Generation
(/./ftp/cats/J/AJ/137/2949)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/137/2949 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.5, on 2013-May-25
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__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. J/AJ/137/2949             MY Cyg differential light curves       (Tucker+, 2009)
*================================================================================
*Orbital solutions and absolute elements of the eclipsing binary MY Cygni.
*    Tucker R.S., Sowell J.R., Williamon R.M., Coughlin J.L.
*   <Astron. J., 137, 2949-2955 (2009)>
*   =2009AJ....137.2949T
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! MY Cyg photometric observations

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

      real*8        HJDV        ! (d) V band Heliocentric Julian date
      real*4        DVmag       ! (mag) Differential Johnson V band magnitude
      real*8        HJDB        ! (d) B band Heliocentric Julian date
      real*4        DBmag       ! (mag) Differential Johnson B band magnitude
      real*8        HJDU        ! (d) U band Heliocentric Julian date
      real*4        DUmag       ! (mag) Differential Johnson U band magnitude

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

C  Declarations for 'table2.dat'	! *MY Cyg radial velocities

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

      real*8        HJD         ! (d) Heliocentric Julian date
      real*4        RVp         ! (km/s) Radial velocity of the primary
      real*4        RVs         ! (km/s) Radial velocity of the secondary

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

C  Declarations for 'table7.dat'	! MY Cyg photometric observations by Coughlin 2007
                              (undergraduate thesis, Emory Univ; same
                              comparison star as Williamon 1975AJ.....80..976W)

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

      real*8        HJDB_1      ! (d) ? B band Heliocentric Julian date
      real*4        DBmag_1     ! (mag) ? Differential B band magnitude
      real*8        HJDV_1      ! (d) ? V band Heliocentric Julian date
      real*4        DVmag_1     ! (mag) ? Differential V band magnitude
      real*8        HJDR        ! (d) ? R band Heliocentric Julian date
      real*4        DRmag       ! (mag) ? Differential R band magnitude
      real*8        HJDI        ! (d) ? I band Heliocentric Julian date
      real*4        DImag       ! (mag) ? Differential I band magnitude

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

C  Loading file 'table1.dat'	! MY Cyg photometric observations

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,file='table1.dat', status='old')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,1733
        read(1,'(A62)')ar__
        read(ar__,1)HJDV,DVmag,HJDB,DBmag,HJDU,DUmag
c    ..............Just test output...........
        write(6,1)HJDV,DVmag,HJDB,DBmag,HJDU,DUmag
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! *MY Cyg radial velocities

C  Format for file interpretation

    2 format(F11.3,1X,F6.1,1X,F6.1)

C  Effective file loading

      open(unit=1,file='table2.dat', status='old')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,28
        read(1,'(A25)')ar__1
        read(ar__1,2)HJD,RVp,RVs
c    ..............Just test output...........
        write(6,2)HJD,RVp,RVs
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! MY Cyg photometric observations by Coughlin 2007
*                              (undergraduate thesis, Emory Univ; same
*                              comparison star as Williamon 1975AJ.....80..976W)

C  Format for file interpretation

    3 format(
     +  F12.4,1X,F6.3,1X,F12.4,1X,F6.3,1X,F12.4,1X,F6.3,1X,F12.4,1X,
     +  F6.3)

C  Effective file loading

      open(unit=1,file='table7.dat', status='old')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,309
        read(1,'(A79)')ar__2
        read(ar__2,3)
     +  HJDB_1,DBmag_1,HJDV_1,DVmag_1,HJDR,DRmag,HJDI,DImag
        if(ar__2(1:12) .EQ. '') HJDB_1 = rNULL__
        if(ar__2(14:19) .EQ. '') DBmag_1 = rNULL__
        if(ar__2(21:32) .EQ. '') HJDV_1 = rNULL__
        if(ar__2(34:39) .EQ. '') DVmag_1 = rNULL__
        if(ar__2(41:52) .EQ. '') HJDR = rNULL__
        if(ar__2(54:59) .EQ. '') DRmag = rNULL__
        if(ar__2(61:72) .EQ. '') HJDI = rNULL__
        if(ar__2(74:79) .EQ. '') DImag = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  HJDB_1,DBmag_1,HJDV_1,DVmag_1,HJDR,DRmag,HJDI,DImag
c    .......End.of.Just test output...........
      end do
      close(1)

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