FORTRAN Generation
(/./ftp/cats/J/A_AS/65/485)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_AS/65/485 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-Apr-19
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+AS/65/485       22MHz flux densities of radio sources    (Roger+, 1986)
*================================================================================
*Spectral flux densities of radio sources at 22 MHz.
*    Roger R.S., Costain C.H., Stewart D.I.
*   <Astron. Astrophys. Suppl. Ser. 65, 485 (1986)>
*   =1986A&AS...65..485R (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Flux densities at 22.25MHz

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

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

      character*8   Bname       ! Name based on 1950 position
      integer*4     RAh         ! (h) Right Ascension B1950 (hours) (1)
      integer*4     RAm         ! (min) Right Ascension B1950 (minutes) (1)
      real*4        RAs         ! (s) [0,60[? Right Ascension B1950 (seconds) (1)
      character*1   DE_         ! Declination B1950 (sign) (1)
      integer*4     DEd         ! (deg) Declination B1950 (degrees) (1)
      integer*4     DEm         ! (arcmin) Declination B1950 (minutes) (1)
      integer*4     DEs         ! (arcsec) [0,60]? Declination B1950 (seconds) (1)
      character*18  Source      ! Source name
      integer*4     S22MHz      ! (Jy) Flux density
      integer*4     e_S22MHz    ! (Jy) ? Mean error on S22MHz
      integer*4     o_S22MHz    ! ? Number of observations
      character*19  Conf        ! Name of confusing sources
      character*16  Bkg         ! Name of background complexity
      character*7   Notes       ! Notes (2)
*Note (1): the position was added by Trushkin (1998) with CATS
*Note (2): the symbols have the following meaning:
*       C: a 26.3 MHz flux density measured at Clark Lake is given by
*          Viner and Erickson (1975AJ.....80..931V)
*   G1, G2, G3, G4: Low frequency flux densities measured at
*          Grakovo are given by Braude et al. (1978Ap&SS..54...37B,
*          1979Ap&SS..64...73B, 1981Ap&SS..74..409B, 1985Ap&SS.111..237B)
*          respectively.
*       1: See notes for this source in Roger et al. (1973AJ.....78.1030R).
*       2: Comprises 3C66A and 3C66B and possibly a broad emission component.
*       3: Confusing sources surrounding 3C196 or emission from
*          broad region centered on the source: near Abell 637.
*       4: Transits 0.8 min earlier than 4C14.31, near right ascension
*          of Abell 795.
*       5: Evidence of E-W broadening.
*       6: Confused region; other sources not catalogued. Alternative
*          interpretation with source on broad background feature
*          gives a flux density of 128 +/- 20 Jy.
*       7: Broad emission centred on 3C298 (~50' EW).
*       8: Extended source, corrected for broadening.
*       9: Right ascension 45 sec later than listed in 4C (lobe shift ?).
*      10: Possibly extended.
*      11: Primary calibration source.
*      12: Confused region.
*      13: Epoch 1966.5.

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

C  Loading file 'table1.dat'	! Flux densities at 22.25MHz

C  Format for file interpretation

    1 format(
     +  A8,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,2X,A18,1X,I5,1X,
     +  I4,2X,I2,1X,A19,A16,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,395
        read(1,'(A108)')ar__
        read(ar__,1)
     +  Bname,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Source,S22MHz,e_S22MHz,
     +  o_S22MHz,Conf,Bkg,Notes
        if(ar__(17:20) .EQ. '') RAs = rNULL__
        if(ar__(29:30) .EQ. '') DEs = iNULL__
        if(ar__(58:61) .EQ. '') e_S22MHz = iNULL__
        if(ar__(64:65) .EQ. '') o_S22MHz = 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)
     +  Bname,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Source,S22MHz,e_S22MHz,
     +  o_S22MHz,Conf,Bkg,Notes
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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