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