Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/135/155 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-Apr-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/135/155 Radio identification of 3EG EGRET sources (Mattox+, 2001)
*================================================================================
*A quantitative evaluation of potential radio identifications
*for 3EG EGRET sources.
* Mattox J.R., Hartman R.C., Reimer O.
* <Astrophys. J. Suppl. Ser. 135, 155 (2001)>
* =2001ApJS..135..155M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Complete EGRET Identification Results
integer*4 nr__
parameter (nr__=366) ! Number of records
character*141 ar__ ! Full-size record
character*10 v3EG ! The J2000 3EG EGRET position name
character*1 f_BName ! [aA] Flag on BName (1)
character*21 BName ! 3EG Blazar identification
character*1 n_BName ! [bde] Note on BName (2)
character*15 RName ! Position name of potential radio counterpart
character*1 f_RName ! [acfg] Flag on RName (3)
character*14 OName ! Other or common name for RName
integer*4 S4_85GHz ! (mJy) ? Flux density at 4.85GHz from catalog of
* Becker et al., <VIII/13>
real*4 Sp_Index ! ? Radio spectral index (4)
real*4 R0 ! (arcmin) ? Characteristic angle (5)
real*4 eta ! ? Probability radio source is a {gamma}-ray
* source (see equ. 6)
real*4 Rad ! (arcmin) ? Angle between the EGRET position estimate
* and the radio position
real*4 R95_ ! (arcmin) ? The 95% confidence contour radius in the
* direction of the radio source
character*1 l_Cont ! Limit flag on Cont
real*4 Cont ! (%) ? EGRET position confidence contour at the
* radio position (6)
real*4 LR ! ? Likelihood ratio (7)
real*4 Prob ! ? The a posteriori probability that the
* identification is correct as determined with
* equ. 4
*Note (1):
* A = Indicates that H-EG3 designated this as a high-confidence AGN
* identification.
* a = Indicates that H-EG3 designated this as a
* ``lower-confidence potential'' AGN identification.
*Note (2):
* b = The poor resolution of the single dish PMN survey does not
* provide for definitive results for radio sources near the LMC.
* d = The Cen A region is not resolved in the single dish PMN survey.
* Therefore, this study is not useful for judging the claim that EGRET
* source J1324-4314 is Cen A.
* e = EGRET source J1627-2419 is indicated by H-3EG as a possible artifact
* due to improper modeling of diffuse {gamma}-ray emission of the
* {rho} Oph cloud.
*Note (3):
* a = B0220+4245 is an extended radio source, making it an unlikely
* counterpart.
* c = EGRET sources J1308+8744 and J1621+8203 lie outside of a 5 GHz
* survey region.
* f = Radio source J1631-4015 corresponds to a planetary nebula.
* g = Radio source B2226+6122 corresponds to a Galactic H II region.
*Note (4): S({nu}) {prop.to} {nu}^{alpha}^.
*Note (5): Angle between potentially confusing radio sources which are at least
* as bright and flat as this radio source (equation 3 of
* Mattox et al., 1997, Cat. <J/ApJ/481/95>).
*Note (6): See equation 6 of Mattox et al. 1997, Cat. <J/ApJ/481/95>.
* Ranges from 0% at the EGRET position estimate to 100% far away.
*Note (7): Indicating the strength of the indication for the identification
* (see equation 5).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! EGRET Identifications with a High Probability of
Being Correct
integer*4 nr__1
parameter (nr__1=46) ! Number of records
character*91 ar__1 ! Full-size record
character*12 RName_1 ! Radio name, as in table1.
character*8 Name ! Radio source name based on 1950 position
real*4 Omag ! (mag) ? Optical magnitude from NED
integer*4 S4_85GHz_1 ! (mJy) Flux density at 4.85GHz from catalog of
* Becker et al., <VIII/13>
real*4 Sp_Index_1 ! Radio spectral index (2)
character*1 l_SVLBI ! Limit flag on a SVLBI
integer*4 SVLBI ! (mJy) ? VLBI flux density at 2.29 GHz
integer*4 e_SVLBI ! (mJy) ? Uncertainty in SVLBI
real*8 Prob_1 ! The a posteriori probability of a
* correct identification (3)
real*4 Flux ! (10-6ct/cm2/s) The {gamma}-ray flux (E>100MeV) in
* units of 10^-6^cm^-2^s^-1^
real*4 e_Flux ! (10-6ct/cm2/s) Uncertainty in Flux
character*1 l_z ! Limit flag on z
real*4 z ! ? Redshift
character*12 n_Name ! Further notes on Name
*Note (2): S({nu}) {prop.to} {nu}^{alpha}^.
*Note (3): See equation 10 of Mattox et al. 1997, Cat. <J/ApJ/481/95>)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Plausible EGRET Identifications
integer*4 nr__2
parameter (nr__2=39) ! Number of records
character*96 ar__2 ! 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*12 RName_2 ! Radio name, as in table1.
character*9 Name_1 ! Radio source name, based on 1950 position
character*1 f_Name ! [gh] Flag on Name (2)
real*4 Omag_1 ! (mag) ? Optical magnitude from NED
integer*4 RAh ! (h) Right Ascension (J2000)
integer*4 RAm ! (min) Right Ascension (J2000)
real*8 RAs ! (s) Right Ascension (J2000)
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)
integer*4 S4_85GHz_2 ! (mJy) Flux density at 4.85GHz from catalog of
* Becker et al., <VIII/13>
real*4 Sp_Index_2 ! ? Radio spectral index (3)
real*4 Prob_2 ! The a posteriori probability of a correct
* identification (eq. 4).
character*1 f_Prob ! [f] Flag on Prob (4)
character*12 Morph ! The ``Morphological Type'' from NED
character*1 f_Morph ! [aA] Flag on Morph (5)
*Note (2):
* g = Both 1930-397 and 1933-400 are potential identifications of
* 3EG J1935-4022.
* h = Both 1936-155 and 1938-155 are potential identifications of
* 3EG J1937-1529.
*Note (3): S({nu}) {prop.to} {nu}^{alpha}^.
*Note (4):
* f = The high confidence identification of 1219+285 by H-3EG is based on
* the E>1 GeV position of Lamb & Macomb (1997). Using their position
* result, we find Prob=0.04 for 1219+285.
*Note (5):
* A = Indicates that H-EG3 designated this source as a high-confidence
* EGRET identification.
* a = Indicates that H-EG3 designated this source as a
* ``lower-confidence potential'' EGRET identification.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Elliptical fits to the 95% confidence position
contours of all 3EG EGRET sources
integer*4 nr__3
parameter (nr__3=271) ! Number of records
character*64 ar__3 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*10 v3EG_1 ! The 3EG EGRET {gamma}-ray position name
real*4 RAdeg_1 ! (deg) Right Ascension expressed in degrees (J2000)
real*4 DEdeg_1 ! (deg) Declination expressed in degrees (J2000)
real*4 GLON ! (deg) Galactic longitude
real*4 GLAT ! (deg) Galactic latitude
real*4 MajAxis ! (arcmin) Semi-major axis length
real*4 MinAxis ! (arcmin) Semi-minor axis length
integer*4 PA ! (deg) Position angle of MajAxis
character*2 f_PA ! Flag on PA (1)
*Note (1):
* a = A poor elliptical fit, a fit to a compromised map, or an anomalous
* error region. The elliptical fits should only be used with caution for
* these sources.
* b = The 68% contour is not centered in the 95% contour, causing an offset
* from the position of H-3EG.
C=============================================================================
C Loading file 'table1.dat' ! Complete EGRET Identification Results
C Format for file interpretation
1 format(
+ A10,1X,A1,1X,A21,1X,A1,1X,A15,1X,A1,1X,A14,1X,I6,1X,F4.1,1X,
+ F6.1,1X,E10.4,1X,F5.1,1X,F5.1,1X,A1,F6.3,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,366
read(1,'(A141)')ar__
read(ar__,1)
+ v3EG,f_BName,BName,n_BName,RName,f_RName,OName,S4_85GHz,
+ Sp_Index,R0,eta,Rad,R95_,l_Cont,Cont,LR,Prob
if(ar__(71:76) .EQ. '') S4_85GHz = iNULL__
if(ar__(78:81) .EQ. '') Sp_Index = rNULL__
if(ar__(83:88) .EQ. '') R0 = rNULL__
if(ar__(90:99) .EQ. '') eta = rNULL__
if(ar__(101:105) .EQ. '') Rad = rNULL__
if(ar__(107:111) .EQ. '') R95_ = rNULL__
if(ar__(114:119) .EQ. '') Cont = rNULL__
if(ar__(121:130) .EQ. '') LR = rNULL__
if(ar__(132:141) .EQ. '') Prob = rNULL__
c ..............Just test output...........
write(6,1)
+ v3EG,f_BName,BName,n_BName,RName,f_RName,OName,S4_85GHz,
+ Sp_Index,R0,eta,Rad,R95_,l_Cont,Cont,LR,Prob
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! EGRET Identifications with a High Probability of
* Being Correct
C Format for file interpretation
2 format(
+ A12,4X,A8,1X,F4.1,1X,I5,1X,F4.1,1X,A1,I4,1X,I3,1X,F7.5,1X,
+ F5.2,1X,F4.2,1X,A1,F6.4,1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,46
read(1,'(A91)')ar__1
read(ar__1,2)
+ RName_1,Name,Omag,S4_85GHz_1,Sp_Index_1,l_SVLBI,SVLBI,e_SVLBI,
+ Prob_1,Flux,e_Flux,l_z,z,n_Name
if(ar__1(26:29) .EQ. '') Omag = rNULL__
if(ar__1(43:46) .EQ. '') SVLBI = iNULL__
if(ar__1(48:50) .EQ. '') e_SVLBI = iNULL__
if(ar__1(72:77) .EQ. '') z = rNULL__
c ..............Just test output...........
write(6,2)
+ RName_1,Name,Omag,S4_85GHz_1,Sp_Index_1,l_SVLBI,SVLBI,e_SVLBI,
+ Prob_1,Flux,e_Flux,l_z,z,n_Name
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Plausible EGRET Identifications
C Format for file interpretation
3 format(
+ A12,4X,A9,1X,A1,1X,F5.2,1X,I2,1X,I2,1X,F7.4,2X,A1,I2,1X,I2,1X,
+ F6.3,2X,I4,1X,F4.1,1X,F5.3,1X,A1,1X,A12,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,39
read(1,'(A96)')ar__2
read(ar__2,3)
+ RName_2,Name_1,f_Name,Omag_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,
+ S4_85GHz_2,Sp_Index_2,Prob_2,f_Prob,Morph,f_Morph
if(ar__2(29:33) .EQ. '') Omag_1 = rNULL__
if(ar__2(70:73) .EQ. '') Sp_Index_2 = 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,3)
+ RName_2,Name_1,f_Name,Omag_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,
+ S4_85GHz_2,Sp_Index_2,Prob_2,f_Prob,Morph,f_Morph
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Elliptical fits to the 95% confidence position
* contours of all 3EG EGRET sources
C Format for file interpretation
4 format(
+ A10,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F5.1,1X,F5.1,1X,I3,1X,
+ A2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,271
read(1,'(A64)')ar__3
read(ar__3,4)
+ v3EG_1,RAdeg_1,DEdeg_1,GLON,GLAT,MajAxis,MinAxis,PA,f_PA
c ..............Just test output...........
write(6,4)
+ v3EG_1,RAdeg_1,DEdeg_1,GLON,GLAT,MajAxis,MinAxis,PA,f_PA
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end