FORTRAN Generation
(/./ftp/cats/J/AJ/133/2559)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/133/2559 into FORTRAN code for reading data files line by line.

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/AJ/133/2559    Compact radio sources in spiral galaxies. II.   (Maddox+, 2007)
*================================================================================
*A study of compact radio sources in nearby face-on spiral galaxies.
*II. Multiwavelength analyses of sources in M 51.
*    Maddox L.A., Cowan J.J., Kilgard R.E., Schinnerer E., Stockdale C.J.
*   <Astron. J., 133, 2559-2568 (2007)>
*   =2007AJ....133.2559M
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Positions, Flux Densities and Spectral Indices
                             of Point Sources in M51

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     v_MCK2007_  ! Source number
      character*1   n__MCK2007_ ! [b] Radio coincidence (1)
      integer*4     RAh         ! (h) Hour of Right Ascension (J2000)
      integer*4     RAm         ! (min) Minute of Right Ascension (J2000)
      real*4        RAs         ! (s) Second of Right Ascension (J2000)
      character*1   DE_         ! Sign of the Declination (J2000)
      integer*4     DEd         ! (deg) Degree of Declination (J2000)
      integer*4     DEm         ! (arcmin) Arcminute of Declination (J2000)
      real*4        DEs         ! (arcsec) Arcsecond of Declination (J2000)
      integer*4     S20cm       ! (uJy) ? 20cm (1.5GHz) flux density in {mu}Jy/beam
      integer*4     e_S20cm     ! (uJy) ? Error in S20cm in {mu}Jy/beam
      integer*4     S6cm        ! (uJy) ? 6cm (5GHz) flux density in {mu}Jy/beam
      integer*4     e_S6cm      ! (uJy) ? Error in S6cm in {mu}Jy/beam
      character*1   l_alpha     ! Limit flag on spectral index
      real*4        alpha       ! Spectral Index (2)
      real*4        e_alpha     ! ? Error in spectral index
      character*16  CXOU        ! CXOU designation (JHHMMSS.s+DDMMSS) of
*                                   X-ray counterpart
      integer*4     LX          ! (10+30W) ? X-ray luminosity in 0.3-2keV band of
*                                    X-ray counterpart
      character*7   XType       ! Source classification based on X-ray colours
      character*1   n_CXOU      ! [a] X-ray coincidence (1)
*Note (1): Note as follows:
*    a = Coincident with SN 1994I
*    b = Coincident with optical HII region
*Note (2): Where S_{nu}_{propto}{nu}^+{alpha}^.

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

C  Loading file 'table2.dat'	! Positions, Flux Densities and Spectral Indices
*                             of Point Sources in M51

C  Format for file interpretation

    1 format(
     +  I3,1X,A1,1X,I2,1X,I2,1X,F5.2,A1,I2,1X,I2,1X,F5.2,1X,I4,1X,I2,
     +  1X,I4,1X,I2,1X,A1,F6.3,1X,F5.3,1X,A16,1X,I1,1X,A7,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,107
        read(1,'(A88)')ar__
        read(ar__,1)
     +  v_MCK2007_,n__MCK2007_,RAh,RAm,RAs,DE_,DEd,DEm,DEs,S20cm,
     +  e_S20cm,S6cm,e_S6cm,l_alpha,alpha,e_alpha,CXOU,LX,XType,n_CXOU
        if(ar__(31:34) .EQ. '') S20cm = iNULL__
        if(ar__(36:37) .EQ. '') e_S20cm = iNULL__
        if(ar__(39:42) .EQ. '') S6cm = iNULL__
        if(ar__(44:45) .EQ. '') e_S6cm = iNULL__
        if(ar__(55:59) .EQ. '') e_alpha = rNULL__
        if(ar__(78:78) .EQ. '') LX = iNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  v_MCK2007_,n__MCK2007_,RAh,RAm,RAs,DE_,DEd,DEm,DEs,S20cm,
     +  e_S20cm,S6cm,e_S6cm,l_alpha,alpha,e_alpha,CXOU,LX,XType,n_CXOU
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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