FORTRAN Generation
(/./ftp/cats/J/ApJS/131/95)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/131/95 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-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/ApJS/131/95       VSOP 5 GHz continuum survey              (Fomalont+, 2000)
*================================================================================
*The VSOP 5 GHz continuum survey: the prelaunch VLBA observations.
*    Fomalont E.B., Frey S., Paragi Z., Gurvits L.I., Scott W.K., Taylor A.R.,
*    Edwards P.G., Hirabayashi H.
*   <Astrophys. J. Suppl. Ser. 131, 95 (2000)>
*   =2000ApJS..131...95F
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! The 5 GHz VLBA Prelaunch (VLBApls) Source List

      integer*4 nr__
      parameter (nr__=374)	! Number of records
      character*98 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)

      character*10  Name        ! IAU source name (based on J2000 position)
      character*8   CName       ! Common name
      integer*4     RAh         ! (h) Right Ascension (J2000)
      integer*4     RAm         ! (min) Right Ascension (J2000)
      real*8        RAs         ! (s) Right Ascension (J2000)
      character*1   f_RAs       ! [*] Indicates a positional accuracy of 10 mas
      character*1   DE_         ! Sign of the Declination (J2000)
      integer*4     DEd         ! (deg) Declination (J2000)
      integer*4     DEm         ! (arcmin) Declination (J2000)
      real*4        DEs         ! (arcsec) Declination (J2000)
      real*4        S5GHz       ! (Jy) Flux density at 5GHz from finding catalogs
      integer*4     GLAT        ! (deg) Galactic latitude
      character*1   OCount      ! Optical counterpart (1)
      character*1   l_z         ! Limit flag on the redshift
      real*4        z           ! ? Redshift (2)
      character*1   v5GHz       ! [gc] Other 5 GHz observations (3)
      character*1   v1_6GHz     ! [l] Indicates a 1.6 GHz VLBI image in the CJ
*                                   survey (Polandis, 1995, Cat. <J/ApJS/98/1>)
      character*1   Image       ! [nNv] Other 2.3 and/or 8.4 GHz
*                                   observations (4)
      character*1   v15GHzU     ! [uU] Single (u) or multiple (U) epoch 15 GHz
*                                   VLBI images exist in the USNO database
      character*1   v15GHzK     ! [uU] Single (u) or multiple (U) epoch 15 GHz
*                                   VLBI images exist in Kellermann et al.
*                                   (1998AJ....115.1295K)
      character*1   q_Name      ! Quality flag for the object (5)
*Note (1): Optical counter part:
*    Q = Quasar
*    B = BL Lac object
*    E = Empty field or unidentified optical counterpart.
*    G = Radio galaxy or other type of AGN other than B and Q
*        (e.g. Seyfert galaxy).
*Note (2): Redshift and optical identification are according to Veron-Cetty &
*     Veron (1998, see Cat. <VII/215>) for Q and B objects and some G unless
*     otherwise stated in the comments (see last column and Section 4 of the
*     paper). For several sources identified as G, redshift value is from NED.
*Note (3): 5GHz observations:
*    c = Caltech-Jodrell Bank Flat Spectrum (CJF) survey
*         (Taylor et al., 1996ApJS..107...37T)
*    g = Other published images in
*         Gurvits, Kellermann & Frey (1999, Cat. <J/A+A/342/378>).
*Note (4): 2.3 and/or 8.4 GHz observations:
*    n = Single epoch observation in the USNO Radio Reference Frame Image
*        Database.
*    N = Multiple epoch observations in the USNO Radio Reference Frame Image
*        Database.
*    v = VLBA Calibrators database
*Note (5): Quality flags:
*    No entry = Image is of good quality and presented in Figure 2, with the
*               source model in Table 2.
*    "=" = Source was very resolved and only a correlated flux density versus
*          u-v spacing plot is shown in Figure 3.
*    x = A comment is also given for this source in Section 4 of the paper.
*    # = Source was not detected and is probably fully resolved with VLBA
*        resolution.

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

C  Declarations for 'table2.dat'	! Source Component Parameters

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

      character*10  Name_1      ! IAU source name
      real*4        FluxDen     ! (Jy) Total flux density of the component
      real*4        EWPos       ! (mas) East-west offset from reference position
      real*4        NSPos       ! (mas) North-south offset from reference position
      character*1   l_MajDiam   ! Limit flag on MajDiam
      real*4        MajDiam     ! (mas) The FWHM angular size of major diameter
      real*4        e_MajDiam   ! (mas) ? Uncertainty in MajDiam
      character*1   l_MinDiam   ! Limit flag on MinDiam
      real*4        MinDiam     ! (mas) The FWHM angular size of minor diameter
      real*4        e_MinDiam   ! (mas) ? Uncertainty in MinDiam
      integer*4     PA          ! (deg) ? Position angle of the major axis

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

C  Loading file 'table1.dat'	! The 5 GHz VLBA Prelaunch (VLBApls) Source List

C  Format for file interpretation

    1 format(
     +  A10,2X,A8,2X,I2,1X,I2,1X,F7.4,1X,A1,2X,A1,I2,1X,I2,1X,F6.3,2X,
     +  F5.1,3X,I3,3X,A1,3X,A1,F5.3,3X,A1,2X,A1,2X,A1,2X,A1,2X,A1,2X,
     +  A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,374
        read(1,'(A98)')ar__
        read(ar__,1)
     +  Name,CName,RAh,RAm,RAs,f_RAs,DE_,DEd,DEm,DEs,S5GHz,GLAT,
     +  OCount,l_z,z,v5GHz,v1_6GHz,Image,v15GHzU,v15GHzK,q_Name
        if(ar__(74:78) .EQ. '') z = rNULL__
        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)
     +  Name,CName,RAh,RAm,RAs,f_RAs,DE_,DEd,DEm,DEs,S5GHz,GLAT,
     +  OCount,l_z,z,v5GHz,v1_6GHz,Image,v15GHzU,v15GHzK,q_Name
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Source Component Parameters

C  Format for file interpretation

    2 format(
     +  A10,3X,F6.3,2X,F5.1,2X,F5.1,2X,A1,F4.1,2X,F3.1,3X,A1,F4.1,2X,
     +  F3.1,2X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,930
        read(1,'(A63)')ar__1
        read(ar__1,2)
     +  Name_1,FluxDen,EWPos,NSPos,l_MajDiam,MajDiam,e_MajDiam,
     +  l_MinDiam,MinDiam,e_MinDiam,PA
        if(ar__1(43:45) .EQ. '') e_MajDiam = rNULL__
        if(ar__1(56:58) .EQ. '') e_MinDiam = rNULL__
        if(ar__1(61:63) .EQ. '') PA = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name_1,FluxDen,EWPos,NSPos,l_MajDiam,MajDiam,e_MajDiam,
     +  l_MinDiam,MinDiam,e_MinDiam,PA
c    .......End.of.Just test output...........
      end do
      close(1)

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