FORTRAN Generation
(/./ftp/cats/J/AJ/109/359)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/109/359 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 2026-Oct-07
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/109/359                   H235 in NGC 752                    (Milone+ 1995)
*================================================================================
*Observations and Analysis of the Contact Binary H235 in the Open Cluster NGC 752
*     Milone E.F., Stagg C.R., Sugars B.A., McVean J.R., Schiller S.J.,
*     Kallrath  D.  H.  Bradstreet J.
*    <Astron. J. 109, 359 (1995)>
*    =1995AJ....109..359M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3'	! Photoelectric light curve data

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

      real*8        HJD        (nr__) ! (d) Heliocentric Julian Date
      real*4        Ph1        (nr__) ! *[0/1[ Phase of minima from ephemeris 1
      real*4        Ph3        (nr__) ! *[0/1[ Phase of minima from ephemeris 3
      real*4        dVmag      (nr__) ! (mag) Differential V magnitude
      real*4        dBmag      (nr__) ! (mag) Differential B magnitude
      real*4        dRmag      (nr__) ! (mag) Differential R magnitude
      real*4        dImag      (nr__) ! (mag) Differential I magnitude
      real*4        d_B_V      (nr__) ! (mag) Differential (B-V) color
      real*4        d_V_R      (nr__) ! (mag) Differential (V-R) color
      real*4        d_V_I      (nr__) ! (mag) ?Differential (V-I) color
*Note on Ph1:
*        Phase of minima calculated using ephemeris (1) from printed paper:
*        E0 = 2 447 815.338(3), P = 0.411 83(12)d
*Note on Ph3:
*        Phase of minima calculated using ephemeris (3) from printed paper:
*        E0 = 2 446 785.8660(23), P = 0.411 775 0(10)d

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

C  Loading file 'table3'	! Photoelectric light curve data

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3')
      write(6,*) '....Loading file: table3'
      do i__=1,867
        read(1,'(A80)')ar__
        read(ar__,1)
     +  HJD(i__),Ph1(i__),Ph3(i__),dVmag(i__),dBmag(i__),dRmag(i__),
     +  dImag(i__),d_B_V(i__),d_V_R(i__),d_V_I(i__)
        if(ar__(75:80) .EQ. '') d_V_I(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  HJD(i__),Ph1(i__),Ph3(i__),dVmag(i__),dBmag(i__),dRmag(i__),
     +  dImag(i__),d_B_V(i__),d_V_R(i__),d_V_I(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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