FORTRAN Generation
(/./ftp/cats/J/ApJ/481/95)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/481/95 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-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/481/95        Radio identification of EGRET sources (Mattox+ 1997)
*================================================================================
*The identification of EGRET sources with flat-spectrum radio sources.
*       Mattox J.R., Schachter J., Molnar L., Hartman R.C., Patnaik A.R.
*      <Astrophys. J. 481, 95 (1997)>
*      =1997ApJ...481...95M      (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Potential radio identifications of sources in
                                  the second EGRET catalog

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

C  Position composed of: GLON GLAT
      character*10  v2EG       (nr__) ! J2000 position name of the EGRET source
      character*1   n_2EG      (nr__) ! [cfj] Note (1)
      real*4        GLON       (nr__) ! (deg) Galactic longitude
      real*4        GLAT       (nr__) ! (deg) Galactic latitude
      character*1   n_GLAT     (nr__) ! [km] Note on GLAT (2)
      real*4        VI         (nr__) ! ? Variability index (3)
      character*8   ID         (nr__) ! ? EGRET catalogue identification
*                                     (OtherName in EGRET catalog description)
      character*1   u_ID       (nr__) ! [?] ?: lower confidence level,
*                                         else high-confidence level
      character*12  RadioID    (nr__) ! Possible radio counterpart (4)
      character*1   n_RadioID  (nr__) ! [abdginxy] Note (5)
      character*12  Names      (nr__) ! Other name of counterpart
      character*1   n_Names    (nr__) ! [hlew] Note (6)
      integer*4     S5GHz      (nr__) ! (mJy) ? Density flux at 4.85GHz
      real*4        SI         (nr__) ! ? Radio spectral index
      real*4        r0         (nr__) ! (arcmin) ? Characteristic angle between confusing
*                                     sources that are at least as bright and
*                                     flat as this radio source
      real*8        eta        (nr__) ! ? a priori probability {eta} = p(id) that the
*                                      source is a {gamma}-ray source
      real*4        r          (nr__) ! (arcmin) ? Angle between the EGRET position estimate
*                                     and the radio position
      real*4        r95        (nr__) ! (arcmin) ? Radius of the 95% confidence contour in
*                                     the direction of the radio source
      real*8        Contour    (nr__) ! (%) ? Position contour at the radio position
      real*4        LR         (nr__) ! ? Likelihood ratio indicating the strength of
*                                     the indication for the determination
      real*4        p_id_r     (nr__) ! ? a posteriori probability that the
*                                     identification is correct
*Note (1): c: Verbunt et al. (1996A&A...311L...9V) identify this source with
*              millisecond pulsar J0218+4232.
*          f: Identification analysis not performed, radio survey source confused
*          j: Not (completely) in GB or PMN survey.
*Note (2): k: Using the position of Hunter et al. (1994ApJ...436..216H), which
*              differs significantly from that of Thompson et al.
*              (1995, Cat. <J/ApJS/101/259>) because of the use of a different
*              Galactic diffuse model.
*          m: Position estimate (and error ellipse a=35', b=20', {phi}=8{deg})
*              obtained from VP7.2 with E>1000Mev data, the analysis used for the
*              identification of Thompson et al. (1995, Cat. <J/ApJS/101/259>)
*              with 2022-077.
*Note (3): Variability index, -logP, where P is the probability that the source
*           is invariant based on a {chi}^2^ test of flux for all available
*           viewing periods. A grater variability index corresponds to a higher
*           probability that the source is variable
*Note (4): B1950 position are used for Green Bank sources and original
*           J2000 position names are used for PMN sources
*Note (5): a: Planetary nebula.
*          b: Shows OVV behavior.
*          d: A gamma-ray hare was reported by Lundgren et al.
*              (1995IAUC.6258....1L) during the 1995 August 8-22 EGRET exposure.
*              An optical spectrum (S. Lundgren 1997, private communication)
*              found only a BL Lac-type continuum spectrum.
*          g: This potential identification indicated by Mukherjee et al.
*              (1995ApJ...445..189M) and Nolan et al. (1996ApJ...459..100N).
*          i: Dingus et al. (1996ApJ...467..589D) note this possible
*              identification.
*          n: Fich (1993ApJS...86..475F) identified this extended radio source
*              with an H II region. Therefore, it is unlikely to be the EGRET
*              source.
*          x: A {gamma}-flare of this source (NRAO 190) was observed by
*              McGlynn et al., 1997ApJ...481..625M
*          y: This potential identification noted by Thompson et al.
*              (1996ApJS..107..227T , their table2)
*Note (6): e: Classified as a galaxy by Stickel & Kuhr 1994. However, S. Wagner
*              (1996, private communication) finds optical variability on a
*              ~ 1week timescale.
*          h: This source is very likely to be the EGRET source because it is
*              detected at TeV energies (Buckley et al., 1996ApJ...472L...9B).
*          l: Classified as a galaxy by Stickel & Kuhr
*              (1994, Cat. <J/A+AS/105/211>)
*          w: 0219+428 = 3C 66

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

C  Declarations for 'table2.dat'	! Potential radio identifications of EGRET sources
                                 in the supplement to the second EGRET catalog

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

C  Position composed of: GLON GLAT
      character*10  v2EG_1     (nr__1) ! J2000 position name of the EGRET source
      character*1   n_2EG_1    (nr__1) ! [cfj] Note (1)
      real*4        GLON_1     (nr__1) ! (deg) Galactic longitude
      real*4        GLAT_1     (nr__1) ! (deg) Galactic latitude
      character*1   n_GLAT_1   (nr__1) ! [km] Note on GLAT (2)
      real*4        VI_1       (nr__1) ! ? Variability index (3)
      character*8   ID_1       (nr__1) ! ? EGRET catalogue identification
*                                     (OtherName in EGRET catalog description)
      character*1   u_ID_1     (nr__1) ! [?] ?: lower confidence level,
*                                         else high-confidence level
      character*12  RadioID_1  (nr__1) ! Possible radio counterpart (4)
      character*1   n_RadioID_1(nr__1) ! [abdginxy] Note (5)
      character*12  Names_1    (nr__1) ! Other name of counterpart
      character*1   n_Names_1  (nr__1) ! [hlew] Note (6)
      integer*4     S5GHz_1    (nr__1) ! (mJy) ? Density flux at 4.85GHz
      real*4        SI_1       (nr__1) ! ? Radio spectral index
      real*4        r0_1       (nr__1) ! (arcmin) ? Characteristic angle between confusing
*                                     sources that are at least as bright and
*                                     flat as this radio source
      real*8        eta_1      (nr__1) ! ? a priori probability {eta} = p(id) that the
*                                      source is a {gamma}-ray source
      real*4        r_1        (nr__1) ! (arcmin) ? Angle between the EGRET position estimate
*                                     and the radio position
      real*4        r95_1      (nr__1) ! (arcmin) ? Radius of the 95% confidence contour in
*                                     the direction of the radio source
      real*8        Contour_1  (nr__1) ! (%) ? Position contour at the radio position
      real*4        LR_1       (nr__1) ! ? Likelihood ratio indicating the strength of
*                                     the indication for the determination
      real*4        p_id_r_1   (nr__1) ! ? a posteriori probability that the
*                                     identification is correct
*Note (1): c: Verbunt et al. (1996A&A...311L...9V) identify this source with
*              millisecond pulsar J0218+4232.
*          f: Identification analysis not performed, radio survey source confused
*          j: Not (completely) in GB or PMN survey.
*Note (2): k: Using the position of Hunter et al. (1994ApJ...436..216H), which
*              differs significantly from that of Thompson et al.
*              (1995, Cat. <J/ApJS/101/259>) because of the use of a different
*              Galactic diffuse model.
*          m: Position estimate (and error ellipse a=35', b=20', {phi}=8{deg})
*              obtained from VP7.2 with E>1000Mev data, the analysis used for the
*              identification of Thompson et al. (1995, Cat. <J/ApJS/101/259>)
*              with 2022-077.
*Note (3): Variability index, -logP, where P is the probability that the source
*           is invariant based on a {chi}^2^ test of flux for all available
*           viewing periods. A grater variability index corresponds to a higher
*           probability that the source is variable
*Note (4): B1950 position are used for Green Bank sources and original
*           J2000 position names are used for PMN sources
*Note (5): a: Planetary nebula.
*          b: Shows OVV behavior.
*          d: A gamma-ray hare was reported by Lundgren et al.
*              (1995IAUC.6258....1L) during the 1995 August 8-22 EGRET exposure.
*              An optical spectrum (S. Lundgren 1997, private communication)
*              found only a BL Lac-type continuum spectrum.
*          g: This potential identification indicated by Mukherjee et al.
*              (1995ApJ...445..189M) and Nolan et al. (1996ApJ...459..100N).
*          i: Dingus et al. (1996ApJ...467..589D) note this possible
*              identification.
*          n: Fich (1993ApJS...86..475F) identified this extended radio source
*              with an H II region. Therefore, it is unlikely to be the EGRET
*              source.
*          x: A {gamma}-flare of this source (NRAO 190) was observed by
*              McGlynn et al., 1997ApJ...481..625M
*          y: This potential identification noted by Thompson et al.
*              (1996ApJS..107..227T , their table2)
*Note (6): e: Classified as a galaxy by Stickel & Kuhr 1994. However, S. Wagner
*              (1996, private communication) finds optical variability on a
*              ~ 1week timescale.
*          h: This source is very likely to be the EGRET source because it is
*              detected at TeV energies (Buckley et al., 1996ApJ...472L...9B).
*          l: Classified as a galaxy by Stickel & Kuhr
*              (1994, Cat. <J/A+AS/105/211>)
*          w: 0219+428 = 3C 66

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

C  Loading file 'table1.dat'	! Potential radio identifications of sources in
*                                  the second EGRET catalog

C  Format for file interpretation

    1 format(
     +  A10,A1,1X,F6.2,1X,F6.2,A1,F5.2,1X,A8,A1,1X,A12,A1,1X,A12,A1,
     +  I5,1X,F4.1,1X,F6.1,1X,F8.6,1X,F5.1,1X,F5.1,1X,F7.3,1X,E15.3,
     +  1X,E8.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,154
        read(1,'(A141)')ar__
        read(ar__,1)
     +  v2EG(i__),n_2EG(i__),GLON(i__),GLAT(i__),n_GLAT(i__),VI(i__),
     +  ID(i__),u_ID(i__),RadioID(i__),n_RadioID(i__),Names(i__),
     +  n_Names(i__),S5GHz(i__),SI(i__),r0(i__),eta(i__),r(i__),
     +  r95(i__),Contour(i__),LR(i__),p_id_r(i__)
        if(ar__(27:31) .EQ. '') VI(i__) = rNULL__
        if(ar__(70:74) .EQ. '') S5GHz(i__) = iNULL__
        if(ar__(76:79) .EQ. '') SI(i__) = rNULL__
        if(ar__(81:86) .EQ. '') r0(i__) = rNULL__
        if(ar__(88:95) .EQ. '') eta(i__) = rNULL__
        if(ar__(97:101) .EQ. '') r(i__) = rNULL__
        if(ar__(103:107) .EQ. '') r95(i__) = rNULL__
        if(ar__(109:115) .EQ. '') Contour(i__) = rNULL__
        if(ar__(117:131) .EQ. '') LR(i__) = rNULL__
        if(ar__(133:140) .EQ. '') p_id_r(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  v2EG(i__),n_2EG(i__),GLON(i__),GLAT(i__),n_GLAT(i__),VI(i__),
     +  ID(i__),u_ID(i__),RadioID(i__),n_RadioID(i__),Names(i__),
     +  n_Names(i__),S5GHz(i__),SI(i__),r0(i__),eta(i__),r(i__),
     +  r95(i__),Contour(i__),LR(i__),p_id_r(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Potential radio identifications of EGRET sources
*                                 in the supplement to the second EGRET catalog

C  Format for file interpretation

    2 format(
     +  A10,A1,1X,F6.2,1X,F6.2,A1,F5.2,1X,A8,A1,1X,A12,A1,1X,A12,A1,
     +  I5,1X,F4.1,1X,F6.1,1X,F8.6,1X,F5.1,1X,F5.1,1X,F7.3,1X,E15.3,
     +  1X,E8.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,60
        read(1,'(A140)')ar__1
        read(ar__1,2)
     +  v2EG_1(i__),n_2EG_1(i__),GLON_1(i__),GLAT_1(i__),
     +  n_GLAT_1(i__),VI_1(i__),ID_1(i__),u_ID_1(i__),RadioID_1(i__),
     +  n_RadioID_1(i__),Names_1(i__),n_Names_1(i__),S5GHz_1(i__),
     +  SI_1(i__),r0_1(i__),eta_1(i__),r_1(i__),r95_1(i__),
     +  Contour_1(i__),LR_1(i__),p_id_r_1(i__)
        if(ar__1(27:31) .EQ. '') VI_1(i__) = rNULL__
        if(ar__1(70:74) .EQ. '') S5GHz_1(i__) = iNULL__
        if(ar__1(76:79) .EQ. '') SI_1(i__) = rNULL__
        if(ar__1(81:86) .EQ. '') r0_1(i__) = rNULL__
        if(ar__1(88:95) .EQ. '') eta_1(i__) = rNULL__
        if(ar__1(97:101) .EQ. '') r_1(i__) = rNULL__
        if(ar__1(103:107) .EQ. '') r95_1(i__) = rNULL__
        if(ar__1(109:115) .EQ. '') Contour_1(i__) = rNULL__
        if(ar__1(117:131) .EQ. '') LR_1(i__) = rNULL__
        if(ar__1(133:140) .EQ. '') p_id_r_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  v2EG_1(i__),n_2EG_1(i__),GLON_1(i__),GLAT_1(i__),
     +  n_GLAT_1(i__),VI_1(i__),ID_1(i__),u_ID_1(i__),RadioID_1(i__),
     +  n_RadioID_1(i__),Names_1(i__),n_Names_1(i__),S5GHz_1(i__),
     +  SI_1(i__),r0_1(i__),eta_1(i__),r_1(i__),r95_1(i__),
     +  Contour_1(i__),LR_1(i__),p_id_r_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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