Conversion of standardized ReadMe file for
file /./ftp/cats/IX/60 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 2026-Oct-13
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. IX/60 Fermi superluminal sources (Xiao+, 2020)
*================================================================================
*Comparison between Fermi detected and non-Fermi detected superluminal sources
* Xiao H.B., Fan J.H., Yang J.H., Liu Y., Yuan Y.H., Tao J., Costantin D.,
* Zhang Y.T., Pei Z.Y., Zhang L.X., Yang W.X.
* <Science China Physics, Mechanics & Astronomy, 62, 129811 (2019)>
*Further evidence of superluminal active galactic nuclei as gamma-ray sources
* Xiao J.H., Fan R.. Rando J.T.. Zhu H.B., Hu L.J.
* <Astron. Nachrichten, 341, 462-470 (2020)>
*Efficient Fermi source identification with machine learning methods
* Xiao H.B., Cao H.T., Fan J.H., Costantin D., Luo G.Y., Pei Z.Y.
* <Astronomy and Computing, 32, 100387 (2020)>
* =2020yCat.9060....0X
* =2019SCPMA..62l9811X
* +2020AN....341..462X
* +2020A&C....3200387X
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat1-1.dat' ! Superluminal sources detected by Fermi-LAT
(data for 377 individual sources)
(Paper I, table 1)
integer*4 nr__
parameter (nr__=816) ! Number of records
character*127 ar__ ! Full-size record
character*18 Name ! Source name
character*5 Class ! Classification (G1)
real*8 z ! ? Redshift
real*4 delta ! ? Doppler factor
character*4 r_delta ! Reference for deltaR (G2)
real*4 Rmag ! (mag) ? R magnitude from BZCAT
real*4 Ext ! (mag) ? Galactic extinction, A_R_, from NED
integer*4 F1_4GHz ! ? Flux density at 1.4GHz from BZCAT
real*4 FX ! ? X-ray flux in the 0.1-2.4keV band
* from BZCAT
real*4 alphaph ! ? gamma-ray photom spectral index
real*4 Flux1000 ! ? gamma-ray photon flux in 0.1-100Gev range
character*1 n_Flux1000 ! [F] Note on Flux1000
real*8 pm ! (uas/yr) ? Proper motion
real*4 e_pm ! (uas/yr) ? rms uncertainty on pm
character*1 n_pm ! [e] Note on pm
character*3 Comp ! Component
character*6 r_pm ! Proper motion reference (G3)
real*4 beta ! ? Apparent velocity
real*8 e_beta ! ? rms uncertainty on beta
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat1-2.dat' ! Superluminal sources not detected by Fermi-LAT
(Paper I, table 2)
integer*4 nr__1
parameter (nr__1=383) ! Number of records
character*103 ar__1 ! Full-size record
character*10 Name_1 ! Name (HHMM+DDd)
character*6 Class_1 ! Classification (G1)
real*8 z_1 ! ? Redshift
real*4 delta_1 ! ? Doppler factor
character*4 r_delta_1 ! Reference for deltaR (G2)
real*4 Rmag_1 ! (mag) ? R magnitude from BZCAT
real*4 Ext_1 ! (mag) ? Calactic exteinction, A_R_, from NED
integer*4 F1_4GHz_1 ! ? Flux density at 1.4GHz from BZCAT
* (1.4/0.843GHz)
real*4 FX_1 ! ? X-ray flux in the 0.1-2.4keV band from
* BZCAT
real*8 pm_1 ! (uas/yr) ? Proper motion
real*4 e_pm_1 ! (uas/yr) ? rms uncertainty on pm
character*1 n_pm_1 ! [e] Note on pm
character*3 Comp_1 ! Component
character*8 r_pm_1 ! Proper motion reference (G3)
real*4 beta_1 ! ? Apparent velocity
real*4 e_beta_1 ! ? rms uncertainty on beta
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat2-1.dat' ! New Fermi detected superluminal sources from
4FGL (table1, Paper II) (Paper II, table 1)
integer*4 nr__2
parameter (nr__2=40) ! Number of records
character*73 ar__2 ! Full-size record
character*17 Name_2 ! Source name (4FGL JHHMM.m+DDMM)
character*8 OName ! Other name (HHMM+DDd)
character*5 Class_2 ! Classification (G1)
real*4 z_2 ! Redshift
real*4 deltaR ! ? Radio Doppler factor
character*3 r_deltaR ! Reference for deltaR (G2)
real*4 Flux1000_1 ! (ph/cm2/s) Integral photon flux from 1 to 100GeV
real*4 e_Flux1000 ! (ph/cm2/s) rms uncertainty on integral photon flux
* from 1 to 100GeV
real*4 alphaph_1 ! gamma-ray photon spectral index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat2-3.dat' ! 181 Fermi detected superluminal sources with
available Doppler factor (Paper II, table 3)
integer*4 nr__3
parameter (nr__3=181) ! Number of records
character*88 ar__3 ! Full-size record
character*17 Name_3 ! Source name
* (3FGL JHHMM.m+DDMM or 4FGl JHHMM.m+DDMM)
character*13 OName_1 ! Other name
character*5 Class_3 ! Classification (G1)
real*4 z_3 ! Redshift
real*4 deltaR_1 ! ? Radio Doppler factor
character*3 r_deltaR_1 ! Reference for deltaR (G2)
real*4 Flux1000_2 ! (ph/cm2/s) Integral photon flux from 1 to 100GeV
real*4 alphaph_2 ! gamma-ray photon spectral index
real*4 beta_2 ! Maximun apparent velocity, beta^max^*_app_c
real*8 deltaLg ! Estimated Doppler factor of this work
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat2-5.dat' ! 1505 blazars' luminosity-Doppler factor
(Paper II, table 5)
integer*4 nr__4
parameter (nr__4=1505) ! Number of records
character*56 ar__4 ! Full-size record
character*17 Name_4 ! Source name (4FGl JHHMM.m+DDMM)
character*3 Class_4 ! Classification (G1)
real*8 z_4 ! Redshift
real*4 Flux1000_3 ! (ph/cm2/s) Integral photon flux from 1 to 100GeV
real*4 alphaph_3 ! gamma-ray photon spectral index
real*4 deltaLg_1 ! Estimated Doppler factor
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat2-6.dat' ! 187 Fermi detected superluminal sources with
total and component flux from Lister et al.
(2018, Cat. J/ApJS/234/12) and Lister et al.
(2019, Cat. J/AJ/138/1874) (Paper II, table 6)
integer*4 nr__5
parameter (nr__5=187) ! Number of records
character*63 ar__5 ! Full-size record
character*17 Name_5 ! Source name
* (3FGL JHHMM.m+DDMM or 4FGl JHHMM.m+DDMM)
character*5 Class_5 ! Classification (G1)
real*4 z_5 ! Redshift
character*3 Comp_2 ! Component id
integer*4 IntF15GHzcomp ! (mJy) Component flux density at 15GHz
real*8 IntF15GHztot ! (mJy) Total flux density at 15GHz
real*4 alphaph_4 ! gamma-ray photon spectral index
real*4 Flux1000_4 ! (ph/cm2/s) Integral photon flux from 1 to 100GeV
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat3-7.dat' ! Identification results of the 1010 unassociated
sources with Best Model A (Paper III, table 7)
integer*4 nr__6
parameter (nr__6=1010) ! Number of records
character*59 ar__6 ! Full-size record
C Position composed of: GLON GLAT
character*17 Name_6 ! Source name (3FGL JHHMM.m+DDMM)
character*1 n_Name ! [c] Note on Name (1)
real*4 GLON ! (deg) Galactic longitude
real*4 GLAT ! (deg) Galactic latitude
real*4 SpIndex ! Spectral Index
real*8 VarIndex ! Variability Index
character*3 TypePred ! Prediction type (2)
real*8 Like ! Likelihood
*Note (1): c indicates that based on the region of the sky the source is is
* considered to be potentially confused with Galactic diffuse emission.
*Note (2): Prediction type: AGN or Non for cat3-7.dat, B or F for cat3-8.dat.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat3-8.dat' ! Identification results of the 573 BCUs with
Best Model B ((Paper III, table 8)
integer*4 nr__7
parameter (nr__7=573) ! Number of records
character*59 ar__7 ! Full-size record
C Position composed of: GLON GLAT
character*17 Name_7 ! Source name (3FGL JHHMM.m+DDMM)
character*1 n_Name_1 ! [c] Note on Name (1)
real*4 GLON_1 ! (deg) Galactic longitude
real*4 GLAT_1 ! (deg) Galactic latitude
real*4 SpIndex_1 ! Spectral Index
real*8 VarIndex_1 ! Variability Index
character*3 TypePred_1 ! Prediction type (2)
real*8 Like_1 ! Likelihood
*Note (1): c indicates that based on the region of the sky the source is is
* considered to be potentially confused with Galactic diffuse emission.
*Note (2): Prediction type: AGN or Non for cat3-7.dat, B or F for cat3-8.dat.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cat3-15.dat' ! 113 uncertain BCUs from Chiaro et al. (2016,
Cat. J/MNRAS/462/3180) and Lefaucheur and
Pita (2017, Cat. J/A+A/602/A86)
(Paper III, table 15)
integer*4 nr__8
parameter (nr__8=113) ! Number of records
character*35 ar__8 ! Full-size record
character*17 Name_8 ! Source name (3FGL JHHMM.m+DDMM)
character*3 C2016 ! [Unc ] Uncertain BCUs from Chiaro et al.
* (2016, Cat. J/MNRAS/462/3180)
character*3 L2017 ! [Unc ] Uncertain BCUs from Lefaucheur and Pita
* (2017, Cat. J/A+A/602/A86)
character*1 TW ! [FB] Prediction from this work
real*8 Like_2 ! Likelihood
C=============================================================================
C Loading file 'cat1-1.dat' ! Superluminal sources detected by Fermi-LAT
* (data for 377 individual sources)
* (Paper I, table 1)
C Format for file interpretation
1 format(
+ A18,1X,A5,1X,F8.6,1X,F5.2,A4,1X,F4.1,1X,F6.3,1X,I5,1X,F6.2,1X,
+ F6.4,1X,E8.3,A1,1X,F8.3,F6.2,A1,1X,A3,1X,A6,1X,F6.3,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'cat1-1.dat')
write(6,*) '....Loading file: cat1-1.dat'
do i__=1,816
read(1,'(A127)')ar__
read(ar__,1)
+ Name,Class,z,delta,r_delta,Rmag,Ext,F1_4GHz,FX,alphaph,
+ Flux1000,n_Flux1000,pm,e_pm,n_pm,Comp,r_pm,beta,e_beta
if(ar__(26:33) .EQ. '') z = rNULL__
if(ar__(35:39) .EQ. '') delta = rNULL__
if(ar__(45:48) .EQ. '') Rmag = rNULL__
if(ar__(50:55) .EQ. '') Ext = rNULL__
if(ar__(57:61) .EQ. '') F1_4GHz = iNULL__
if(ar__(63:68) .EQ. '') FX = rNULL__
if(ar__(70:75) .EQ. '') alphaph = rNULL__
if(ar__(77:84) .EQ. '') Flux1000 = rNULL__
if(ar__(87:94) .EQ. '') pm = rNULL__
if(ar__(95:100) .EQ. '') e_pm = rNULL__
if(ar__(114:119) .EQ. '') beta = rNULL__
if(ar__(121:127) .EQ. '') e_beta = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,Class,z,delta,r_delta,Rmag,Ext,F1_4GHz,FX,alphaph,
+ Flux1000,n_Flux1000,pm,e_pm,n_pm,Comp,r_pm,beta,e_beta
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat1-2.dat' ! Superluminal sources not detected by Fermi-LAT
* (Paper I, table 2)
C Format for file interpretation
2 format(
+ A10,1X,A6,1X,F8.6,1X,F5.2,1X,A4,1X,F4.1,1X,F5.3,1X,I5,1X,F5.2,
+ 1X,F7.2,1X,F6.2,A1,1X,A3,1X,A8,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'cat1-2.dat')
write(6,*) '....Loading file: cat1-2.dat'
do i__=1,383
read(1,'(A103)')ar__1
read(ar__1,2)
+ Name_1,Class_1,z_1,delta_1,r_delta_1,Rmag_1,Ext_1,F1_4GHz_1,
+ FX_1,pm_1,e_pm_1,n_pm_1,Comp_1,r_pm_1,beta_1,e_beta_1
if(ar__1(19:26) .EQ. '') z_1 = rNULL__
if(ar__1(28:32) .EQ. '') delta_1 = rNULL__
if(ar__1(39:42) .EQ. '') Rmag_1 = rNULL__
if(ar__1(44:48) .EQ. '') Ext_1 = rNULL__
if(ar__1(50:54) .EQ. '') F1_4GHz_1 = iNULL__
if(ar__1(56:60) .EQ. '') FX_1 = rNULL__
if(ar__1(62:68) .EQ. '') pm_1 = rNULL__
if(ar__1(70:75) .EQ. '') e_pm_1 = rNULL__
if(ar__1(91:96) .EQ. '') beta_1 = rNULL__
if(ar__1(98:103) .EQ. '') e_beta_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,Class_1,z_1,delta_1,r_delta_1,Rmag_1,Ext_1,F1_4GHz_1,
+ FX_1,pm_1,e_pm_1,n_pm_1,Comp_1,r_pm_1,beta_1,e_beta_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat2-1.dat' ! New Fermi detected superluminal sources from
* 4FGL (table1, Paper II) (Paper II, table 1)
C Format for file interpretation
3 format(
+ A17,1X,A8,1X,A5,1X,F5.3,1X,F5.2,1X,A3,1X,E8.3,1X,E8.3,1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'cat2-1.dat')
write(6,*) '....Loading file: cat2-1.dat'
do i__=1,40
read(1,'(A73)')ar__2
read(ar__2,3)
+ Name_2,OName,Class_2,z_2,deltaR,r_deltaR,Flux1000_1,
+ e_Flux1000,alphaph_1
if(ar__2(40:44) .EQ. '') deltaR = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,OName,Class_2,z_2,deltaR,r_deltaR,Flux1000_1,
+ e_Flux1000,alphaph_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat2-3.dat' ! 181 Fermi detected superluminal sources with
* available Doppler factor (Paper II, table 3)
C Format for file interpretation
4 format(
+ A17,1X,A13,1X,A5,2X,F5.3,3X,F5.2,1X,A3,1X,E8.3,1X,F6.4,2X,
+ F6.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'cat2-3.dat')
write(6,*) '....Loading file: cat2-3.dat'
do i__=1,181
read(1,'(A88)')ar__3
read(ar__3,4)
+ Name_3,OName_1,Class_3,z_3,deltaR_1,r_deltaR_1,Flux1000_2,
+ alphaph_2,beta_2,deltaLg
if(ar__3(48:52) .EQ. '') deltaR_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3,OName_1,Class_3,z_3,deltaR_1,r_deltaR_1,Flux1000_2,
+ alphaph_2,beta_2,deltaLg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat2-5.dat' ! 1505 blazars' luminosity-Doppler factor
* (Paper II, table 5)
C Format for file interpretation
5 format(A17,1X,A3,1X,F8.6,1X,E10.5,1X,F6.4,2X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'cat2-5.dat')
write(6,*) '....Loading file: cat2-5.dat'
do i__=1,1505
read(1,'(A56)')ar__4
read(ar__4,5)Name_4,Class_4,z_4,Flux1000_3,alphaph_3,deltaLg_1
c ..............Just test output...........
write(6,5)Name_4,Class_4,z_4,Flux1000_3,alphaph_3,deltaLg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat2-6.dat' ! 187 Fermi detected superluminal sources with
* total and component flux from Lister et al.
* (2018, Cat. J/ApJS/234/12) and Lister et al.
* (2019, Cat. J/AJ/138/1874) (Paper II, table 6)
C Format for file interpretation
6 format(A17,1X,A5,1X,F5.3,1X,A3,I4,1X,F8.2,1X,F6.4,2X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'cat2-6.dat')
write(6,*) '....Loading file: cat2-6.dat'
do i__=1,187
read(1,'(A63)')ar__5
read(ar__5,6)
+ Name_5,Class_5,z_5,Comp_2,IntF15GHzcomp,IntF15GHztot,
+ alphaph_4,Flux1000_4
c ..............Just test output...........
write(6,6)
+ Name_5,Class_5,z_5,Comp_2,IntF15GHzcomp,IntF15GHztot,
+ alphaph_4,Flux1000_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat3-7.dat' ! Identification results of the 1010 unassociated
* sources with Best Model A (Paper III, table 7)
C Format for file interpretation
7 format(A17,A1,1X,F6.2,1X,F6.2,1X,F5.3,1X,F8.3,1X,A3,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'cat3-7.dat')
write(6,*) '....Loading file: cat3-7.dat'
do i__=1,1010
read(1,'(A59)')ar__6
read(ar__6,7)
+ Name_6,n_Name,GLON,GLAT,SpIndex,VarIndex,TypePred,Like
c ..............Just test output...........
write(6,7)
+ Name_6,n_Name,GLON,GLAT,SpIndex,VarIndex,TypePred,Like
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat3-8.dat' ! Identification results of the 573 BCUs with
* Best Model B ((Paper III, table 8)
C Format for file interpretation
8 format(A17,A1,1X,F6.2,1X,F6.2,1X,F5.3,1X,F8.3,1X,A3,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'cat3-8.dat')
write(6,*) '....Loading file: cat3-8.dat'
do i__=1,573
read(1,'(A59)')ar__7
read(ar__7,8)
+ Name_7,n_Name_1,GLON_1,GLAT_1,SpIndex_1,VarIndex_1,TypePred_1,
+ Like_1
c ..............Just test output...........
write(6,8)
+ Name_7,n_Name_1,GLON_1,GLAT_1,SpIndex_1,VarIndex_1,TypePred_1,
+ Like_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cat3-15.dat' ! 113 uncertain BCUs from Chiaro et al. (2016,
* Cat. J/MNRAS/462/3180) and Lefaucheur and
* Pita (2017, Cat. J/A+A/602/A86)
* (Paper III, table 15)
C Format for file interpretation
9 format(A17,1X,A3,1X,A3,1X,A1,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'cat3-15.dat')
write(6,*) '....Loading file: cat3-15.dat'
do i__=1,113
read(1,'(A35)')ar__8
read(ar__8,9)Name_8,C2016,L2017,TW,Like_2
c ..............Just test output...........
write(6,9)Name_8,C2016,L2017,TW,Like_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end