FORTRAN Generation
(/./ftp/cats/J/MNRAS/383/479)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/383/479 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/383/479 Radio obs. of optically obscured galaxies (Magliocchetti+, 2008)
*================================================================================
*The radio properties of optically obscured Spitzer sources.
*    Magliocchetti M., Andreani P., Zwaan M.A.
*   <Mon. Not. R. Astron. Soc., 383, 479-496 (2008)>
*   =2008MNRAS.383..479M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Radio counterparts to R>25.5 Spitzer sources

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

      character*16  Name       (nr__) ! Galaxy name (JHHMMSS.s+DDMMSS)
      real*4        F24        (nr__) ! (mJy) MIPS flux at 24um from Fadda et al.
*                                   (2006, Cat. J/AJ/131/2859)
      character*1   l_F3_6     (nr__) ! Limit flag on F3.6
      real*8        F3_6       (nr__) ! (uJy) IRAC flux at 3.6um (1)
      character*1   l_F4_5     (nr__) ! Limit flag on F4.5
      real*8        F4_5       (nr__) ! (uJy) IRAC flux at 4.5um (1)
      character*1   l_F5_8     (nr__) ! Limit flag on F5.8
      real*8        F5_8       (nr__) ! (uJy) IRAC flux at 5.8um (1)
      character*1   l_F8       (nr__) ! Limit flag on F8
      real*8        F8         (nr__) ! (uJy) IRAC flux at 8um (1)
      character*1   l_F1_4G    (nr__) ! Limit flag on F1.4G
      real*4        F1_4G      (nr__) ! (mJy) Flux at 1.4GHz
      character*1   r_F1_4G    (nr__) ! [M] Reference for F1.4G (2)
      real*4        F610M      (nr__) ! (mJy) ?=- Flux at 610MHz from Garn et al.
*                                      (2007, Cat. J/MNRAS/376/1251)
      character*1   n_F1_4G    (nr__) ! [*] Note on radio flux (3)
      character*5   Class      (nr__) ! Class (4)
      integer*4     zt         (nr__) ! [1/5] Likely redshift range (5)
      character*1   l_alpha    (nr__) ! Limit flag on alpha
      real*4        alpha      (nr__) ! ?=- Radio spectral index
*                                      (S({nu}){prop.to}{nu}^-{alpha}^)
*Note (1): IRAC data from Lacy et al. (2005ApJS..161...41L)
*Note (2): All radio data at 1.4 GHz are taken from Condon et al. (2003,
*     Cat. J/AJ/125/2411), except for those sources labelled with M, where
*     the radio information comes from Morganti et al. (2004, Cat.
*     J/A+A/424/371)
*Note (3): In all cases but J171143.9+600741 and J172005.0+592430 (marked
*     with *) the reported radio fluxes correspond to the single radio
*     source closest to the centre of 24-um emission
*Note (4): Class as follows:
*    SB0   = star forming, F24>=0.6; F8/F24<0.1
*    LI    = local interloper
*    AGN0? = AGN0 ?
*    AGN0  = AGN, F24>=0.6; F8/F24>=0.1
*    AGN1  = AGN, F24<0.6; F4.5/F3.6>=4.5/3.6
*    UNCL  = unclassified, F24<0.6; F4.5/F3.6<4.5/3.6 or no IRAC information
*Note (5): Likely redshift range as follows:
*      1 = 0<=z<=1
*      2 = 1<=z<=1.5
*      3 = 1.5<=z<=3
*      4 = pure {beta}>1 Power-Law (PL) AGN SED; no redshift assignment
*          is possible in these cases
*      5 = no redshift estimate was possible in that case

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

C  Loading file 'table1.dat'	! Radio counterparts to R>25.5 Spitzer sources

C  Format for file interpretation

    1 format(
     +  A16,1X,F4.2,1X,A1,F7.2,1X,A1,F7.2,1X,A1,F7.2,1X,A1,F7.2,1X,A1,
     +  F5.2,A1,1X,F5.2,A1,1X,A5,3X,I1,1X,A1,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,96
        read(1,'(A89)')ar__
        read(ar__,1)
     +  Name(i__),F24(i__),l_F3_6(i__),F3_6(i__),l_F4_5(i__),
     +  F4_5(i__),l_F5_8(i__),F5_8(i__),l_F8(i__),F8(i__),
     +  l_F1_4G(i__),F1_4G(i__),r_F1_4G(i__),F610M(i__),n_F1_4G(i__),
     +  Class(i__),zt(i__),l_alpha(i__),alpha(i__)
        if (idig(ar__(67:71)).EQ.0) F610M(i__) =  rNULL__
        if (idig(ar__(85:89)).EQ.0) alpha(i__) =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),F24(i__),l_F3_6(i__),F3_6(i__),l_F4_5(i__),
     +  F4_5(i__),l_F5_8(i__),F5_8(i__),l_F8(i__),F8(i__),
     +  l_F1_4G(i__),F1_4G(i__),r_F1_4G(i__),F610M(i__),n_F1_4G(i__),
     +  Class(i__),zt(i__),l_alpha(i__),alpha(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