FORTRAN Generation
(/./ftp/cats/J/MNRAS/416/2477)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/416/2477 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 2024-Mar-29
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/MNRAS/416/2477    STEREO observations of variable stars    (Wraight+, 2011)
*================================================================================
*STEREO observations of stars and the search for exoplanets.
*     Wraight K.T., White G.J., Bewsher D., Norton A.J..
*    <Mon. Not. R. Astron. Soc. 416, 2477 (2011)>
*    =2011MNRAS.416.2477W
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'tablea1.dat'	! Statistics for 263 EBs and the significance of
                              their detection by different algorithms

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg      (nr__) ! (deg) Right Ascension (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination (J2000)
      real*4        PED        (nr__) ! (mag) Primary Eclipse Depth
      real*8        Per        (nr__) ! (d) Period
      real*8        Wmag       (nr__) ! (mag) Weighted mean magnitude (4)
      real*8        BLS        (nr__) ! Box-Least-Squares variability index 
*                                 (Kovacs et al., 2002A&A...391..369K)
      real*8        SID        (nr__) ! Stetson variability index
      real*8        LS         (nr__) ! Peak Lomb-Scargle power (1982ApJ...263..835S)
      real*4        ecosw      (nr__) ! ?=- Minimum eccentricity |e*cos({omega})|
      real*4        e_ecosw    (nr__) ! ?=- rms uncertainty on |e*cos({omega})|
      character*11  SpType     (nr__) ! MK spectral type (1)
      character*23  SimID      (nr__) ! ID of closest star in SIMBAD at RA/DE
      character*24  SimID2     (nr__) ! ID of alternative nearby star (2)
      integer*4     EB         (nr__) ! [0,1] 1 = known eclipsing binary flag
      real*8        MJD        (nr__) ! (d) mid-eclipse Modified Julian Date
      character*1   Note       (nr__) ! [*] Note for CCDM J03295+2107 (3)
*Note (1): The spectral type given is of the star in SIMBAD closest to the
*  RA/DE coordinates (See also Note (2)).
*Note (2): Due to the large (70") pixel scale of the Heliospheric Imager it is
*  not always possible to be confident of which star is responsible for the
*  variability detected. When a nearby star may potentially be the source of the
*  variability, its identity is recorded here. If either of these stars is
*  recorded as an eclipsing binary in SIMBAD, or the references imply or state
*  they are, then the Known EB flag is set to 1.
*Note (3): For CCDM J03295+2107AB (row 11), has since been found from an
*  analysis including data from the STEREO-Behind spacecraft (STEREO/HI-1B) data
*  to be an eccentric eclipsing binary with a period of 16.6291 days and a
*  minimum eccentricity (|e*cos(w)|) of 0.34+/-0.01. The original published
*  version has about one third of this period, as the analysis instead identified
*  by mistake the time between primary and secondary eclipses. The additional
*  data from STEREO/HI-1B also allowed a refinement of the mid-eclipse time to
*  54212.12009 MJD.
*  Tablea1.dat contains these new data, and then differs from published version
*  in the paper.
*Note (4): the filter bandpass is 630-730nm, with some contributions
*  around 400 and 950nm)

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

C  Loading file 'tablea1.dat'	! Statistics for 263 EBs and the significance of
*                              their detection by different algorithms

C  Format for file interpretation

    1 format(
     +  F10.6,1X,F11.7,1X,F5.3,1X,F8.5,1X,F8.5,1X,F7.5,1X,F14.7,1X,
     +  F8.4,1X,F5.3,1X,F5.3,1X,A11,1X,A23,1X,A24,1X,I1,1X,F12.6,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea1.dat')
      write(6,*) '....Loading file: tablea1.dat'
      do i__=1,263
        read(1,'(A167)')ar__
        read(ar__,1)
     +  RAdeg(i__),DEdeg(i__),PED(i__),Per(i__),Wmag(i__),BLS(i__),
     +  SID(i__),LS(i__),ecosw(i__),e_ecosw(i__),SpType(i__),
     +  SimID(i__),SimID2(i__),EB(i__),MJD(i__),Note(i__)
        if (idig(ar__(80:84)).EQ.0) ecosw(i__) =  rNULL__
        if (idig(ar__(86:90)).EQ.0) e_ecosw(i__) =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  RAdeg(i__),DEdeg(i__),PED(i__),Per(i__),Wmag(i__),BLS(i__),
     +  SID(i__),LS(i__),ecosw(i__),e_ecosw(i__),SpType(i__),
     +  SimID(i__),SimID2(i__),EB(i__),MJD(i__),Note(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end