FORTRAN Generation
(/./ftp/cats/J/AZh/72/291)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AZh/72/291 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/AZh/72/291   Variability at Frequencies 3.9 and 7.5 GHz   (Gorshkov+, 1995)
*================================================================================
*Variability at Frequencies 3.9 and 7.5 GHz: Radio Sources from the
*Zelenchukskii Survey with Fluxes > 200mJy.
*     Gorshkov, A.G., Konnikova, V. K.
*    <Astron. Zh. 72, 291 (1995)>
*    =1995AZh....72..291G
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! A complete sample of radio sources from the
                                 Zelenchukskii survey with fluxes S3.9>200mJy

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

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

      integer*4     RAh         ! (h) Right ascension (1950)
      integer*4     RAm         ! (min) Right ascension (1950)
      real*4        RAs         ! (s) Right ascension (1950)
      character*1   DE_         ! Declination sign
      integer*4     DEd         ! (deg) Declination (1950)
      integer*4     DEm         ! (arcmin) Declination (1950)
      integer*4     DEs         ! (arcsec) Declination (1950)
      character*1   Note        ! [ *]?  If * fluxes and spectral index were
*                                  changed by authors in comparison printed
*                                  version
      integer*4     F3_9        ! (mJy) Mean weighted fluxes at 3.9GHz (7.7cm)
      integer*4     F7_5        ! (mJy) Mean weighted fluxes at 7.5GHz (4cm)
      real*4        a           ! *? Spectral index (3.9-7.5GHz) for the mean
*                                     weighted fluxes.
      real*4        vF3_9       ! ? Relative variability amplitudes at 3.9GHz
*                                    for sources with probability of
*                                    variability greater than 0.95
      real*4        vF7_5       ! ? Relative variability amplitudes at 7.5GHz
*                                    for sources with probability of
*                                    variability greater than 0.95
      real*4        Pvar        ! ? Probability of variability 1-(n-1) at 7.5GHz
      character*2   Ident       ! *Optical identification
      real*4        Vmag        ! (mag) ? Stellar magnitude
      real*4        z           ! ? Redshift
      integer*4     Ref         ! [12-23]? References for the optical
*                                            identifications
*Note on a:
*   Spectral index was not calculated for extended sources.
*Note on Ident:
*    Q = known quasar
*    SO= stellar object
*    G = galaxy
*    EF= empty field on the Palomar Sky Survey image
*    L = Lacertae object
*    F = object too faint for establishment of morphological type
*    PR=populated region
*    St= star

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

C  Loading file 'table1.dat'	! A complete sample of radio sources from the
*                                 Zelenchukskii survey with fluxes S3.9>200mJy

C  Format for file interpretation

    1 format(
     +  I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,A1,1X,I5,2X,I4,2X,F5.2,
     +  2X,F4.2,2X,F4.2,2X,F4.2,2X,A2,2X,F4.1,2X,F5.3,2X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,180
        read(1,'(A79)')ar__
        read(ar__,1)
     +  RAh,RAm,RAs,DE_,DEd,DEm,DEs,Note,F3_9,F7_5,a,vF3_9,vF7_5,Pvar,
     +  Ident,Vmag,z,Ref
        if(ar__(36:40) .EQ. '') a = rNULL__
        if(ar__(43:46) .EQ. '') vF3_9 = rNULL__
        if(ar__(49:52) .EQ. '') vF7_5 = rNULL__
        if(ar__(55:58) .EQ. '') Pvar = rNULL__
        if(ar__(65:68) .EQ. '') Vmag = rNULL__
        if(ar__(71:75) .EQ. '') z = rNULL__
        if(ar__(78:79) .EQ. '') Ref = 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)
     +  RAh,RAm,RAs,DE_,DEd,DEm,DEs,Note,F3_9,F7_5,a,vF3_9,vF7_5,Pvar,
     +  Ident,Vmag,z,Ref
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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