FORTRAN Generation
(/./ftp/cats/J/ApJ/703/1672)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/703/1672 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-Apr-20
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/ApJ/703/1672   Far-IR and H{alpha} fluxes in galaxies   (Kennicutt+, 2009)
*================================================================================
*Dust-corrected star formation rates of galaxies.
*I. Combinations of H{alpha} and infrared tracers.
*    Kennicutt R.C., Hao C.-N., Calzetti D., Moustakas J., Dale D.A., Bendo G.,
*    Engelbracht C.W., Johnson B.D., Lee J.C.
*   <Astrophys. J., 703, 1672-1695 (2009)>
*   =2009ApJ...703.1672K
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Integrated measurements of SINGS galaxies

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

      character*8   Name       (nr__) ! Galaxy name
      real*4        Dist       (nr__) ! (Mpc) Distance from SINGS (Kennicutt et al.
*                                        2003PASP..115..928K)
      real*4        log_Ha_N   (nr__) ! ([mW/m2]) H{alpha}+[NII] flux (in erg/s/cm^2^)
      real*4        e_log_Ha_N (nr__) ! ([mW/m2]) ? 1{sigma} uncertainty in log(Ha+NII)
      character*1   n_log_Ha_N (nr__) ! [ab] why e_log(Ha+NII) is missing (1)
      real*4        Ha_Hb      (nr__) ! H{alpha}/H{beta} flux ratio
      real*4        e_Ha_Hb    (nr__) ! The 1{sigma} uncertainty in Ha/Hb
      character*1   f_Ha_Hb    (nr__) ! [c] Assumed value of Ha/Hb (2)
      real*4        NII_Ha     (nr__) ! NII/H{alpha} flux ratio
      real*4        e_NII_Ha   (nr__) ! 1{sigma} uncertainty in NII/Ha
      real*4        S25        (nr__) ! (Jy) ? 25 micron flux density
      real*4        e_S25      (nr__) ! (Jy) ? 1{sigma} uncertainty in S25
      real*4        S60        (nr__) ! (Jy) ? 60 micron flux density
      real*4        e_S60      (nr__) ! (Jy) ? 1{sigma} uncertainty in S60
      real*4        S100       (nr__) ! (Jy) ? 100 micron flux density
      real*4        e_S100     (nr__) ! (Jy) ? 1{sigma} uncertainty in S100
      character*9   r_Ha_Hb    (nr__) ! Reference(s) on Ha/Hb flux ratio (see
*                                     refs.dat file; 'a' = assumed value
      character*11  r_NII_Ha   (nr__) ! Reference(s) on NII/Ha flux ratio
*                                     (see refs.dat file)
      character*7   Ref        (nr__) ! Reference(s) on IRAS fluxes (see
*                                     refs.dat file)
*Note (1): Flag as follows:
*     a = No error was given in the original paper.
*     b = This value is from R. C. Kennicutt et al. 2009 (in preparation).
*         No error was given.
*Note (2):
*     c= They are low metallicity galaxies. The H{alpha}/H{beta} ratio is
*        assumed to be 2.86 so that the estimated extinction is zero under
*        the Case B assumption made in this paper. The error is the median
*        error of the remaining galaxies.

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

C  Declarations for 'table2.dat'	! Integrated measurements of Moustakas & Kennicutt
                             (2006, Cat. J/ApJS/164/81) sample

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

      character*14  Name_1     (nr__1) ! Galaxy name
      real*4        Dist_1     (nr__1) ! (Mpc) Distance from Moustakas & Kennicutt
*                                 (2006, Cat. J/ApJS/164/81)
      real*8        F_Ha       (nr__1) ! (aW/m2) H{alpha} flux (in 10^-15^arg/cm^2^/s = aW/m^2^)
      real*4        e_F_Ha     (nr__1) ! (aW/m2) 1{sigma} uncertainty in F(Ha)
      real*8        F_Hb       (nr__1) ! (aW/m2) H{beta} flux
      real*4        e_F_Hb     (nr__1) ! (aW/m2) 1{sigma} uncertainty in F(Hb)
      real*4        S25_1      (nr__1) ! (Jy) IRAS 25 micron flux density
      real*4        e_S25_1    (nr__1) ! (Jy) 1{sigma} uncertainty in S25
      real*4        S60_1      (nr__1) ! (Jy) IRAS 60 micron flux density
      real*4        e_S60_1    (nr__1) ! (Jy) 1{sigma} uncertainty in S60
      real*4        S100_1     (nr__1) ! (Jy) IRAS 100 micron flux density
      real*4        e_S100_1   (nr__1) ! (Jy) 1{sigma} uncertainty in S100
      real*4        S1_4       (nr__1) ! (mJy) ? The 1.4GHz flux density
      real*4        e_S1_4     (nr__1) ! (mJy) ? The 1{sigma} uncertainty in S1.4 (1)
      integer*4     Ref_1      (nr__1) ! Reference on IRAS fluxes (see refs.dat file)
*Note (1): Source with a 1.4GHz flux density but no uncertainty are taken
*     from Condon (1987, Cat. J/ApJS/65/485). No uncertainties were given in
*     that paper.

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

C  Declarations for 'refs.dat'	! References

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

      integer*4     Ref_2      (nr__2) ! Reference number
      character*19  BibCode    (nr__2) ! Bibcode
      character*21  Aut        (nr__2) ! Auhtor's name(s)
      character*47  Comm       (nr__2) ! Comment

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

C  Loading file 'table1.dat'	! Integrated measurements of SINGS galaxies

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,58
        read(1,'(A123)')ar__
        read(ar__,1)
     +  Name(i__),Dist(i__),log_Ha_N(i__),e_log_Ha_N(i__),
     +  n_log_Ha_N(i__),Ha_Hb(i__),e_Ha_Hb(i__),f_Ha_Hb(i__),
     +  NII_Ha(i__),e_NII_Ha(i__),S25(i__),e_S25(i__),S60(i__),
     +  e_S60(i__),S100(i__),e_S100(i__),r_Ha_Hb(i__),r_NII_Ha(i__),
     +  Ref(i__)
        if(ar__(24:27) .EQ. '') e_log_Ha_N(i__) = rNULL__
        if(ar__(55:59) .EQ. '') S25(i__) = rNULL__
        if(ar__(61:65) .EQ. '') e_S25(i__) = rNULL__
        if(ar__(67:72) .EQ. '') S60(i__) = rNULL__
        if(ar__(74:78) .EQ. '') e_S60(i__) = rNULL__
        if(ar__(80:85) .EQ. '') S100(i__) = rNULL__
        if(ar__(87:92) .EQ. '') e_S100(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),Dist(i__),log_Ha_N(i__),e_log_Ha_N(i__),
     +  n_log_Ha_N(i__),Ha_Hb(i__),e_Ha_Hb(i__),f_Ha_Hb(i__),
     +  NII_Ha(i__),e_NII_Ha(i__),S25(i__),e_S25(i__),S60(i__),
     +  e_S60(i__),S100(i__),e_S100(i__),r_Ha_Hb(i__),r_NII_Ha(i__),
     +  Ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Integrated measurements of Moustakas & Kennicutt
*                             (2006, Cat. J/ApJS/164/81) sample

C  Format for file interpretation

    2 format(
     +  A14,1X,F5.1,1X,F7.1,1X,F5.1,1X,F8.2,1X,F5.2,1X,F5.2,1X,F5.3,
     +  1X,F5.2,1X,F5.3,1X,F6.2,1X,F5.3,1X,F6.1,1X,F5.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,113
        read(1,'(A103)')ar__1
        read(ar__1,2)
     +  Name_1(i__),Dist_1(i__),F_Ha(i__),e_F_Ha(i__),F_Hb(i__),
     +  e_F_Hb(i__),S25_1(i__),e_S25_1(i__),S60_1(i__),e_S60_1(i__),
     +  S100_1(i__),e_S100_1(i__),S1_4(i__),e_S1_4(i__),Ref_1(i__)
        if(ar__1(88:93) .EQ. '') S1_4(i__) = rNULL__
        if(ar__1(95:99) .EQ. '') e_S1_4(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name_1(i__),Dist_1(i__),F_Ha(i__),e_F_Ha(i__),F_Hb(i__),
     +  e_F_Hb(i__),S25_1(i__),e_S25_1(i__),S60_1(i__),e_S60_1(i__),
     +  S100_1(i__),e_S100_1(i__),S1_4(i__),e_S1_4(i__),Ref_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    3 format(I3,1X,A19,1X,A21,2X,A47)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,29
        read(1,'(A94)')ar__2
        read(ar__2,3)Ref_2(i__),BibCode(i__),Aut(i__),Comm(i__)
c    ..............Just test output...........
        write(6,3)Ref_2(i__),BibCode(i__),Aut(i__),Comm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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