FORTRAN Generation
(/./ftp/cats/J/ApJS/230/7)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/230/7 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-18
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/ApJS/230/7     Spectral flux densities from 50MHz to 50GHz     (Perley+, 2017)
*================================================================================
*An accurate flux density scale from 50 MHz to 50 GHz.
*    Perley R.A., Butler B.J.
*   <Astrophys. J. Suppl. Ser., 230, 7-7 (2017)>
*   =2017ApJS..230....7P    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Source list

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

      character*10  Name       (nr__) ! Source name
      character*21  OName      (nr__) ! Alternate name
      character*1   Band4      (nr__) ! [*] *: Band 4 used for this source (1)
      character*1   BandP      (nr__) ! [*] *: Band P used for this source (1)
      character*1   BandL      (nr__) ! [*] *: Band L used for this source (1)
      character*1   BandS      (nr__) ! [*] *: Band S used for this source (1)
      character*1   BandC      (nr__) ! [*] *: Band C used for this source (1)
      character*1   BandX      (nr__) ! [*] *: Band X used for this source (1)
      character*1   BandKu     (nr__) ! [*] *: Band Ku used for this source (1)
      character*1   BandK      (nr__) ! [*] *: Band K used for this source (1)
      character*1   BandKa     (nr__) ! [*] *: Band Ka used for this source (1)
      character*1   BandQ      (nr__) ! [*] *: Band Q used for this source (1)
      real*8        LAS        (nr__) ! (arcsec) [0.6/3000] Largest Angular Size
*Note (1): Band code as in table 1:
*        -----------------------------------------------------------------
*         Band     lambda    Frequency span
*                             (MHz)
*        -----------------------------------------------------------------
*         4         4m           73-74.6
*         P        90cm         224-480
*         L        20cm        1000-2000
*         S        10cm        2000-4000
*         C         5cm        4000-8000
*         X         3cm        8000-12000
*         Ku        2cm       12000-18000
*         K        13mm       18000-26500
*         Ka        9mm       26500-40000
*         Q         7mm       40000-50000
*        -----------------------------------------------------------------
*        The VLA's low-frequency system now comprises a single receiver
*        covering 50-480MHz and two feeds, covering 50-80, and 224-480MHz,
*        respectively. The 90cm observations described here are taken with
*        this new system. The 4m observations described here are from the
*        narrowband "legacy" system (Kassim+ 2007ApJS..172..686K), which
*        is now disabled.

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

C  Declarations for 'table5.dat'	! Derived spectral flux densities

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

      real*4        Freq       (nr__1) ! (GHz) [0.07/48.1] Frequency
      real*4        J0133      (nr__1) ! (Jy) Flux density of J0133-3629 at Freq (1)
      real*4        v3C48      (nr__1) ! (Jy) Flux density of J0137+3309 (3C48) at Freq (1)
      integer*4     ForA       (nr__1) ! (Jy) Flux density of J0322-3712 (Fornax A) at Freq (1)
      real*4        v3C123     (nr__1) ! (Jy) Flux density of J0437+2940 (3C123) at Freq (1)
      real*4        J0444      (nr__1) ! (Jy) Flux density of J0444-2809 at Freq (1)
      real*4        v3C138     (nr__1) ! (Jy) Flux density of J0521+1638 (3C138) at Freq (1)
      real*4        PicA       (nr__1) ! (Jy) Flux density of J0519-4546 (Pictor A) at Freq (1)
      integer*4     v3C144     (nr__1) ! (Jy) Flux density of J0534+2200 (3C144) at Freq (1)
      real*4        v3C147     (nr__1) ! (Jy) Flux density of J0542+4951 (3C147) at Freq (1)
      real*8        v3C196     (nr__1) ! (Jy) Flux density of J0813+4813 (3C196) at Freq (1)
      real*4        v3C218     (nr__1) ! (Jy) Flux density of J0918-1205 (3C218) at Freq (1)
      integer*4     v3C274     (nr__1) ! (Jy) Flux density of J1230+1223 (3C274) at Freq (1)
      real*4        v3C286     (nr__1) ! (Jy) Flux density of J1331+3030 (3C286) at Freq (1)
      real*8        v3C295     (nr__1) ! (Jy) Flux density of J1411+5212 (3C295) at Freq (1)
      real*4        v3C348     (nr__1) ! (Jy) Flux density of J1651+0459 (3C348) at Freq (1)
      real*4        v3C353     (nr__1) ! (Jy) Flux density of J1720-0058 (3C353) at Freq (1)
      real*4        v3C380     (nr__1) ! (Jy) Flux density of J1829+4844 (3C380) at Freq (1)
      integer*4     v3C405     (nr__1) ! (Jy) Flux density of J1959+4044 (3C405) at Freq (1)
      real*4        v3C444     (nr__1) ! (Jy) Flux density of J2214-1701 (3C444) at Freq (1)
      integer*4     v3C461     (nr__1) ! (Jy) Flux density of J2323+5848 (3C461) at Freq (1)
      real*4        Err        (nr__1) ! (%) [0.2/4.7] Percent post-fit standard deviation
*Note (1): A value of zero indicates no observation at that frequency.

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

C  Declarations for 'table10.dat'	! Flux calibrator suitability for VLA

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

      character*10  Name_1     (nr__2) ! Source name
      character*6   Var        (nr__2) ! Variability criterion (1)
      real*4        MaxFreq    (nr__2) ! (GHz) [0.2/30]? Maximum frequency (if null = All) (2)
      character*1   f_MaxFreq  (nr__2) ! [A] A = All
      character*1   Range1     (nr__2) ! [0] Lower baseline range
      character*1   l_Range2   (nr__2) ! Upper limit flag on Range2
      real*8        Range2     (nr__2) ! [0.05/5000] Upper baseline range
*                                  (kilowavelenghts) to meet confusion,
*                                  structure, and sensitive criteria
      character*47  Comm       (nr__2) ! Comments
*Note (1): Variability criterion--maximum frequency in GHz for sources with
*          strong nuclei.
*Note (2): Frequency below which the primary beam resolution criterion is met.

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

C  Loading file 'table2.dat'	! Source list

C  Format for file interpretation

    1 format(
     +  A10,1X,A21,1X,A1,1X,A1,1X,A1,1X,A1,1X,A1,1X,A1,1X,A1,1X,A1,1X,
     +  A1,1X,A1,1X,F7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,20
        read(1,'(A60)')ar__
        read(ar__,1)
     +  Name(i__),OName(i__),Band4(i__),BandP(i__),BandL(i__),
     +  BandS(i__),BandC(i__),BandX(i__),BandKu(i__),BandK(i__),
     +  BandKa(i__),BandQ(i__),LAS(i__)
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),OName(i__),Band4(i__),BandP(i__),BandL(i__),
     +  BandS(i__),BandC(i__),BandX(i__),BandKu(i__),BandK(i__),
     +  BandKa(i__),BandQ(i__),LAS(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Derived spectral flux densities

C  Format for file interpretation

    2 format(
     +  F6.3,1X,F5.2,1X,F6.3,1X,I3,1X,F6.2,1X,F5.2,1X,F6.3,1X,F5.1,1X,
     +  I4,1X,F5.2,1X,F7.3,1X,F6.2,1X,I4,1X,F5.2,1X,F7.3,1X,F6.2,1X,
     +  F5.1,1X,F6.2,1X,I5,1X,F6.3,1X,I5,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,34
        read(1,'(A137)')ar__1
        read(ar__1,2)
     +  Freq(i__),J0133(i__),v3C48(i__),ForA(i__),v3C123(i__),
     +  J0444(i__),v3C138(i__),PicA(i__),v3C144(i__),v3C147(i__),
     +  v3C196(i__),v3C218(i__),v3C274(i__),v3C286(i__),v3C295(i__),
     +  v3C348(i__),v3C353(i__),v3C380(i__),v3C405(i__),v3C444(i__),
     +  v3C461(i__),Err(i__)
c    ..............Just test output...........
        write(6,2)
     +  Freq(i__),J0133(i__),v3C48(i__),ForA(i__),v3C123(i__),
     +  J0444(i__),v3C138(i__),PicA(i__),v3C144(i__),v3C147(i__),
     +  v3C196(i__),v3C218(i__),v3C274(i__),v3C286(i__),v3C295(i__),
     +  v3C348(i__),v3C353(i__),v3C380(i__),v3C405(i__),v3C444(i__),
     +  v3C461(i__),Err(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! Flux calibrator suitability for VLA

C  Format for file interpretation

    3 format(A10,1X,A6,1X,F4.1,A1,1X,A1,1X,A1,F7.2,1X,A47)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,20
        read(1,'(A82)')ar__2
        read(ar__2,3)
     +  Name_1(i__),Var(i__),MaxFreq(i__),f_MaxFreq(i__),Range1(i__),
     +  l_Range2(i__),Range2(i__),Comm(i__)
        if(ar__2(19:22) .EQ. '') MaxFreq(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Name_1(i__),Var(i__),MaxFreq(i__),f_MaxFreq(i__),Range1(i__),
     +  l_Range2(i__),Range2(i__),Comm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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