FORTRAN Generation
(/./ftp/cats/J/A_A/619/A73)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/619/A73 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-12
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/A+A/619/A73 Solar Twins age-chromospheric activity   (Lorenzo-Oliveira+, 2018)
*================================================================================
*The Solar Twin Planet Search. VI. The age-chromospheric activity relation.
*    Lorenzo-Oliveira D., Freitas F., Melendez J., Bedell M., Ramirez I.,
*    Bean J., Asplund M., Spina L., Dreizler S., Alves-Brito A., Casagrande L.
*    <Astron. Astrophys. 619, A73 (2018)>
*    =2018A&A...619A..73L        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! All the solar twins observed by HARPS, their
                                 respective programs and S values compiled from
                                 the literature

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

      integer*4     HIP        (nr__) ! HIP number
      character*1   n_HIP      (nr__) ! [*] Note on HIP (G1)
      character*106 ESO        (nr__) ! ESO programs
      real*4        S          (nr__) ! ? CaII H and K activity S-index
      real*4        e_S        (nr__) ! ? rms uncertainty on S
      character*9   r_S        (nr__) ! Reference for S value (1)
*Note (1): References as follows:
*           1 = Ramirez et al. (2014A&A...572A..48R)
*           2 = Jenkins et al. (2011, Cat. J/A+A/531/A8)
*           3 = Jenkins et al. (2006, Cat. J/MNRAS/372/163)
*           4 = Wright et al. (2004, Cat. J/ApJS/152/261)
*           5 = Henry et al. (1996, Cat. J/AJ/111/439)
*           6 = Duncan et al. (1991ApJS...76..383D, Cat. III/159)
*           7 = Melendez et al. (2009ApJ...704L..66M)

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

C  Declarations for 'table2.dat'	! Parameters of all sample solar twins

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

      integer*4     HIP_1      (nr__1) ! HIP number
      character*1   n_HIP_1    (nr__1) ! [*] Note on HIP (G1)
      real*4        Vmag       (nr__1) ! (mag) V magnitude
      real*4        e_Vmag     (nr__1) ! (mag) ? rms uncertainty on Vmag
      real*4        B_V        (nr__1) ! (mag) ? B-V colour index
      real*4        e_B_V      (nr__1) ! (mag) ? rms uncertainty on B-V
      integer*4     Teff       (nr__1) ! (K) Effective temperature
      integer*4     e_Teff     (nr__1) ! (K) ? rms uncertainty on Teff
      real*4        logg       (nr__1) ! ([cm/s2]) ? Surface gravity
      real*4        e_logg     (nr__1) ! ([cm/s2]) ? rms uncertainty on logg
      real*4        v_Fe_H_    (nr__1) ! ([-]) ? Metallicity
      real*4        e__Fe_H_   (nr__1) ! ([-]) ? rms uncertainty on [Fe/H]
      real*4        SMW        (nr__1) ! Mount Wildson CaII H and K activity
*                                          S-index, S_MW_
      real*4        e_SMW      (nr__1) ! rms uncertainty on SMW
      real*4        logR_HK_Teff(nr__1) ! logR'HK activity index (from Teff)
      real*4        e_logR_HK_Teff(nr__1) ! rms uncertainty on logR'HK(Teff)
      real*4        Age        (nr__1) ! (Gyr) ? Age
      real*4        E_Age      (nr__1) ! (Gyr) ? Error on Age (upper limit)
      real*4        e_Age_1    (nr__1) ! (Gyr) ? Error on Age (lower limit)
      integer*4     NObs       (nr__1) ! Number of observations
      real*4        Time       (nr__1) ! (yr) Time-span

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

C  Loading file 'table1.dat'	! All the solar twins observed by HARPS, their
*                                 respective programs and S values compiled from
*                                 the literature

C  Format for file interpretation

    1 format(I6,A1,1X,A106,1X,F5.3,1X,F5.3,1X,A9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,82
        read(1,'(A136)')ar__
        read(ar__,1)
     +  HIP(i__),n_HIP(i__),ESO(i__),S(i__),e_S(i__),r_S(i__)
        if(ar__(116:120) .EQ. '') S(i__) = rNULL__
        if(ar__(122:126) .EQ. '') e_S(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  HIP(i__),n_HIP(i__),ESO(i__),S(i__),e_S(i__),r_S(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Parameters of all sample solar twins

C  Format for file interpretation

    2 format(
     +  I6,A1,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,I4,1X,I2,1X,F5.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F3.1,
     +  1X,F3.1,1X,F3.1,1X,I4,1X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,82
        read(1,'(A111)')ar__1
        read(ar__1,2)
     +  HIP_1(i__),n_HIP_1(i__),Vmag(i__),e_Vmag(i__),B_V(i__),
     +  e_B_V(i__),Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),
     +  v_Fe_H_(i__),e__Fe_H_(i__),SMW(i__),e_SMW(i__),
     +  logR_HK_Teff(i__),e_logR_HK_Teff(i__),Age(i__),E_Age(i__),
     +  e_Age_1(i__),NObs(i__),Time(i__)
        if(ar__1(15:19) .EQ. '') e_Vmag(i__) = rNULL__
        if(ar__1(21:25) .EQ. '') B_V(i__) = rNULL__
        if(ar__1(27:31) .EQ. '') e_B_V(i__) = rNULL__
        if(ar__1(38:39) .EQ. '') e_Teff(i__) = iNULL__
        if(ar__1(41:45) .EQ. '') logg(i__) = rNULL__
        if(ar__1(47:51) .EQ. '') e_logg(i__) = rNULL__
        if(ar__1(53:58) .EQ. '') v_Fe_H_(i__) = rNULL__
        if(ar__1(60:64) .EQ. '') e__Fe_H_(i__) = rNULL__
        if(ar__1(91:93) .EQ. '') Age(i__) = rNULL__
        if(ar__1(95:97) .EQ. '') E_Age(i__) = rNULL__
        if(ar__1(99:101) .EQ. '') e_Age_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  HIP_1(i__),n_HIP_1(i__),Vmag(i__),e_Vmag(i__),B_V(i__),
     +  e_B_V(i__),Teff(i__),e_Teff(i__),logg(i__),e_logg(i__),
     +  v_Fe_H_(i__),e__Fe_H_(i__),SMW(i__),e_SMW(i__),
     +  logR_HK_Teff(i__),e_logR_HK_Teff(i__),Age(i__),E_Age(i__),
     +  e_Age_1(i__),NObs(i__),Time(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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