FORTRAN Generation
()

Conversion of standardized ReadMe file for catalog 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.5, on 2013-May-23
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__
ter  (rNULL__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. J/A+A/505/859       M dwarfs radial velocities             (Zechmeister+, 2009)
*================================================================================
*The M dwarf planet search programme at the ESO VLT + UVES.
*A search for terrestrial planets in the habitable zone of M dwarfs.
*     Zechmeister M., Kuerster M., Endl M.
*    <Astron. Astrophys. 505, 859 (2009)>
*    =2009A&A...505..859Z
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'targets.dat'	! Targets

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

      character*8   Name       (nr__) ! Star name
      character*8   OName      (nr__) ! Other name
      integer*4     N          (nr__) ! Number of measurements
      integer*4     T          (nr__) ! (d) Time baseline
      real*4        sigma      (nr__) ! (m/s) Radial velocity scatter
      real*4        v_sigma_   (nr__) ! (m/s) Mean radial velocity measurement error
      character*9   Com        (nr__) ! Comment
      character*7   SpType     (nr__) ! MK spectral type
      real*4        Vmag       (nr__) ! (mag) V magnitude
      real*4        d          (nr__) ! (pc) Heliocentric distance (1)
      real*4        Mass       (nr__) ! (solMass) Stellar mass (2)
      real*8        pmRA       (nr__) ! (mas/yr) Proper motion along RA
      real*4        e_pmRA     (nr__) ! (mas/yr) rms uncertainty on pmRA
      real*8        pmDE       (nr__) ! (mas/yr) Proper motion along DE
      real*4        e_pmDE     (nr__) ! (mas/yr) rms uncertainty on pmDE
      real*4        plx        (nr__) ! (mas) Parallax (1)
      real*4        e_plx      (nr__) ! (mas) rms uncertainty on plx
      real*4        SA         (nr__) ! (m/s/yr) Secular acceleration (d(RV(t))/dt)
      real*4        e_SA       (nr__) ! (m/s/yr) rms uncertainty on SA
      character*11  FileName   (nr__) ! Name of the file with radial velocity
*                                     measurements in "rv" subdirectory
*Note (1): From Hipparcos, van Leuuwen, 2007A&A...474..653V, Cat. <I/311>
*Note (2): Stellar mass derived from the K-band mass-luminosity relation by
*           Delfosse (2000A&A...364..217D)

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

C  Declarations for 'rv/*'	! Individual RV and error bars

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

      real*8        BJD        (nr__1) ! (d) Barycentric Julian date
      real*8        RV         (nr__1) ! (m/s) Differential radial velocity (1)
      real*4        e_RV       (nr__1) ! (m/s) Uncertainty on RV
*Note (1): RV is not corrected for the secular acceleration. See the paper
*     or the summary table (targets.dat) for the correction values derived
*     from Hipparcos parallax and proper motion.

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

C  Loading file 'targets.dat'	! Targets

C  Format for file interpretation

    1 format(
     +  A8,A8,1X,I2,1X,I4,1X,F6.1,1X,F4.1,1X,A9,1X,A7,1X,F5.2,1X,F6.2,
     +  1X,F4.2,1X,F8.2,1X,F4.2,1X,F8.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,
     +  F5.3,1X,F5.3,1X,A11)

C  Effective file loading

      open(unit=1,file='targets.dat', status='old')
      write(6,*) '....Loading file: targets.dat'
      do i__=1,41
        read(1,'(A136)')ar__
        read(ar__,1)
     +  Name(i__),OName(i__),N(i__),T(i__),sigma(i__),v_sigma_(i__),
     +  Com(i__),SpType(i__),Vmag(i__),d(i__),Mass(i__),pmRA(i__),
     +  e_pmRA(i__),pmDE(i__),e_pmDE(i__),plx(i__),e_plx(i__),SA(i__),
     +  e_SA(i__),FileName(i__)
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),OName(i__),N(i__),T(i__),sigma(i__),v_sigma_(i__),
     +  Com(i__),SpType(i__),Vmag(i__),d(i__),Mass(i__),pmRA(i__),
     +  e_pmRA(i__),pmDE(i__),e_pmDE(i__),plx(i__),e_plx(i__),SA(i__),
     +  e_SA(i__),FileName(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'rv/*'	! Individual RV and error bars

C  Format for file interpretation

    2 format(F14.6,1X,F8.2,1X,F6.2)

C  Effective file loading

      open(unit=1,file='rv/*', status='old')
      write(6,*) '....Loading file: rv/*'
      do i__=1,41
        read(1,'(A30)')ar__1
        read(ar__1,2)BJD(i__),RV(i__),e_RV(i__)
c    ..............Just test output...........
        write(6,2)BJD(i__),RV(i__),e_RV(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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