FORTRAN Generation
(/./ftp/cats/J/ApJ/609/539)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/609/539 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/ApJ/609/539      Kinematics of parsec-scale radio jets     (Kellermann+, 2004)
*================================================================================
*Sub-milliarcsecond imaging of quasars and active galactic nuclei.
*III. Kinematics of parsec-scale radio jets.
*    Kellermann K.I., Lister M.L., Homan D.C., Vermeulen R.C., Cohen M.H.,
*    Ros E., Kadler M., Zensus J.A., Kovalev Y.Y.
*   <Astrophys. J., 609, 539-563 (2004)>
*   =2004ApJ...609..539K
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! General Source Properties

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

      character*8   Name        ! IAU name (based on B1950.0 position)
      character*1   f_Name      ! [def] Flag on Name (1)
      character*18  ONames      ! Other names
      character*1   Memb        ! Sub-sample member
      character*1   Class       ! Veron-Cetty & Veron (Cat. <VII/224>) catalog
*                                 optical classification (2)
      character*1   f_Class     ! [abc] Flag on Class (3)
      real*8        z           ! ? Redshift from Veron-Cetty & Veron
*                                 (Cat. <VII/224>) (4)
      character*6   RSpec       ! Radio spectrum description
      character*7   RMorph      ! Radio morphology classification (5)
      real*4        SVLBI       ! (Jy) Cleaned VLBI 15GHz (2cm) flux density
      character*1   E_ID        ! [YPN] EGRET gamma-ray identification (6)
*Note (1): Flag on Name is defined as follows:
*    d = The source 0218+357 is a gravitationally lensed AGN which we do not
*        include in our statistical analysis.
*    e = The complex source 0316+413 (3C 84) satisfies our selection criteria,
*        but was observed as part of a separate monitoring program
*        Walker et al. (2000ApJ...530..233W). It is included in our statistical
*        analysis for completeness.
*    f = Source has large or diffuse structures that are not adequately sampled
*        by our observations.
*Note (2): The optical classification has the following meaning:
*    Q = quasar;
*    B = BL Lac object;
*    G = active galaxy;
*    U = unidentified.
*Note (3): Flag on Class is defined as follows:
*    a = Source not in the Veron-Cetty & Veron (2001, Cat. <VII/224>) catalog.
*    b = Source classified as quasar in the Veron-Cetty & Veron (2001,
*        Cat. <VII/224>) catalog.
*    c = Source classified as galaxy in the Veron-Cetty & Veron (2001,
*        Cat. <VII/224>) catalog.
*Note (4): Redshift from Veron-Cetty & Veron (2001, Cat. <VII/224>).
*    (The redshifts for 0026+346, 0727-115 and 1155+251 are from Zensus 
*    et al. (2002AJ....124..662Z), that of 0202+149 is from Stickel et al.
*    (1996ApJ...468..556S), 0218+357 is from Lawrence
*    (1996IAUS..173..299L), 0754+100 is from Carangelo et al.
*    (2003A&A...412..651C) and that of 2209+236 is from Sowards-Emmerd,
*    Remani, and Michelson (2003, Cat. <J/ApJ/590/109>).)
*Note (5): The radio morphology classification is defined as follows:
*    CJ = core-jet;
*    2S = two-sided jet;
*     N = naked core;
*    LD = lobe-dominated;
*   CSS = compact-steep spectrum object.
*Note (6): According to Mattox et al. (2001, Cat. <J/ApJS/135/155>) and
*     Sowards-Emmerd, Remani, and Michelson (2003, Cat. <J/ApJ/590/109>),
*     the EGRET gamma-ray identification is defined as follows:
*    Y = highly probable identification;
*    P = probable identification;
*    N = no identification.

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

C  Declarations for 'table2.dat'	! Proper motions

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

      character*8   Name_1      ! IAU name
      character*1   f_Name_1    ! [s] Flag on Name (1)
      character*1   n_Name      ! [a] Note on Name (1)
      character*3   m_Name      ! Component identifier
      character*1   c1          ! [d] Indicates brightest component member
      integer*4     NEpoch      ! Number of epochs with measured position
      real*4        Rad         ! (mas) ? Mean radial position (2)
      integer*4     PA          ! (deg) [-178/179]? Mean structural position angle (2)
      real*4        muR         ! (mas/yr) Angular radial speed (total proper motion)
      real*4        e_muR       ! (mas/yr) The 1{sigma} uncertainty in muR
      real*4        beta        ! ? Radial speed in speed of light units (3)
      real*4        e_beta      ! ? The 1{sigma} uncertainty in beta
      real*8        t0          ! (yr) ? Extrapolated epoch of origin
      real*4        e_t0        ! (yr) ? The 1{sigma} uncertainty in t0
      character*1   Dcomp       ! Distinct component (4)
      character*1   QRate       ! Overall quality rating
*Note (1): Flags on Name, defined as follows:
*    s = Source is a member of the representative sample;
*    a = Motions for 0316+413 which are taken from unpublished data do not
*        appear in figure 10.
*        Motions for 0238-084 are from Vermeulen et al. (2003ASPC..300...43V).
*Note (2): Relative to the core.
*Note (3): Computed assuming H_0_=70km/s/Mpc, {Omega}_m_=0.3, and
*          {Omega}_{Lambda}_=0.7.
*Note (4): Does the component satisfy criterion "b" in Section 4.2, i.e.,
*          is it a well defined feature?

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

C  Loading file 'table1.dat'	! General Source Properties

C  Format for file interpretation

    1 format(
     +  A8,1X,A1,1X,A18,1X,A1,1X,A1,1X,A1,1X,F7.5,1X,A6,1X,A7,1X,F5.2,
     +  1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,135
        read(1,'(A67)')ar__
        read(ar__,1)
     +  Name,f_Name,ONames,Memb,Class,f_Class,z,RSpec,RMorph,SVLBI,
     +  E_ID
        if(ar__(37:43) .EQ. '') z = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name,f_Name,ONames,Memb,Class,f_Class,z,RSpec,RMorph,SVLBI,
     +  E_ID
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Proper motions

C  Format for file interpretation

    2 format(
     +  A8,1X,A1,A1,1X,A3,1X,A1,1X,I2,1X,F4.1,1X,I4,1X,F6.3,1X,F5.3,
     +  1X,F6.2,1X,F5.2,1X,F7.2,1X,F5.2,1X,A1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,282
        read(1,'(A75)')ar__1
        read(ar__1,2)
     +  Name_1,f_Name_1,n_Name,m_Name,c1,NEpoch,Rad,PA,muR,e_muR,beta,
     +  e_beta,t0,e_t0,Dcomp,QRate
        if(ar__1(22:25) .EQ. '') Rad = rNULL__
        if(ar__1(27:30) .EQ. '') PA = iNULL__
        if(ar__1(45:50) .EQ. '') beta = rNULL__
        if(ar__1(52:56) .EQ. '') e_beta = rNULL__
        if(ar__1(58:64) .EQ. '') t0 = rNULL__
        if(ar__1(66:70) .EQ. '') e_t0 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name_1,f_Name_1,n_Name,m_Name,c1,NEpoch,Rad,PA,muR,e_muR,beta,
     +  e_beta,t0,e_t0,Dcomp,QRate
c    .......End.of.Just test output...........
      end do
      close(1)

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