Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/690/A111 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 2026-Aug-12
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/A+A/690/A111 IceCube neutrinos-blazars association (Kouch+, 2024)
*================================================================================
*Association of the IceCube neutrinos with blazars in the CGRaBS sample.
* Kouch P.M., Lindfors E., Hovatta T., Liodakis I., Koljonen K.I.I.,
* Nilsson K., Kiehlmann S., Max-Moerbeck W., Readhead A.C.S., Reeves R.A.,
* Pearson T.J., Jormanainen J., Fallah Ramazani V., Graham M.J.
* <Astron. Astrophys. 690, A111 (2024)>
* =2024A&A...690A.111K (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'neutrino.dat' ! IceCat1+: high-energy neutrino events used
integer*4 nr__
parameter (nr__=283) ! Number of records
character*120 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*9 IC (nr__) ! IC ID of the high-energy neutrino event (1)
real*8 MJD (nr__) ! (d) Modified Julian Date (MJD) of the arrival time
real*4 E (nr__) ! (TeV) Energy estimate
real*4 Sness (nr__) ! Signalness
real*4 RAdeg (nr__) ! (deg) Right Ascension (J2000)
real*4 E_RAdeg (nr__) ! (deg) Error on Right Ascension (upper value)
real*4 e_RAdeg_1 (nr__) ! (deg) Error on Right Ascension (lower value)
real*4 DEdeg (nr__) ! (deg) Declination (J2000)
real*4 E_DEdeg (nr__) ! (deg) Error on Declination (upper value)
real*4 e_DEdeg_1 (nr__) ! (deg) Error on Declination (lower value)
real*8 OmMin (nr__) ! (deg+2) Omega in the published error region scenario (2)
real*8 OmMax (nr__) ! (deg+2) Omega in the enlarged error region scenario (2)
real*8 wMin (nr__) ! Weight in the published error region scenario (2)
real*8 wMax (nr__) ! Weight in the enlarged error region scenario (2)
character*8 Ref (nr__) ! IceCube paper reference for the event
*Note (1): The high-energy neutrino events (283) are taken from two references.
* 267 are from Ice-Cat1 (Abbasi et al., 2023ApJS..269...25A,
* [10.3847/1538-4365/acfa95]).
* 16 are from off-line events analyzed in (Abbasi et al. 2022ApJ...928...50A,
* [10.3847/1538-4357/ac4d29]).
* The IC ID of the latter events ends in an 'X'.
*Note (2): The published and enlarged error region scenarios are described in
* Section 3.3 of the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'blazars.dat' ! CGRaBS blazars: monitored by OVRO
integer*4 nr__1
parameter (nr__1=1157) ! Number of records
character*104 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*10 Blz (nr__1) ! CGRaBS blazar name (JHHMM+DDMM)
real*8 RAdeg_1 (nr__1) ! (deg) Blazar right ascension (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Blazar declination (J2000)
character*1 BZCc (nr__1) ! Blazar 5BZCat type (1)
character*3 FGLc (nr__1) ! Blazar 4FGL class (2)
character*71 Alias (nr__1) ! Aliases of the blazar separated by '/'
*Note (1): The 5BZCat blazar types are 'B' for BL Lac objects, 'Q' for Flast
* Spectrum Radio Quasars, 'U' is for uncertain, and '-' means that the
* source is not in the 5BZCat catalog.
*Note (2): The 4FGL classes are low-, intermediate-, and high-synchrotron peaked
* sources, respectively, denoted by LSP, ISP, and HSP. Most of the
* sources do not have a 4FGL class and are marked with '-'.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'assoc.dat' ! All blazar-neutrino spatial associations
integer*4 nr__2
parameter (nr__2=275) ! Number of records
character*148 ar__2 ! Full-size record
character*9 IC_1 (nr__2) ! IC ID of spatially-associated neutrino events
character*10 Blz_1 (nr__2) ! CGRaBS name of spatially-associated blazars
real*8 dBN (nr__2) ! (deg) Distance between the neutrino and blazar
real*8 dMin (nr__2) ! (deg) Distance between the neutrino and published
* edge (1)
real*8 dMax (nr__2) ! (deg) Distance between the neutrino and enlarged
* edge (1)
real*8 Sr (nr__2) ! (Jy) Unweighted mean radio flux density (1)
real*8 wSrMin (nr__2) ! (Jy) Weighted mean radio flux density (published) (1)
real*8 wSrMax (nr__2) ! (Jy) Weighted mean radio flux density (enlarged) (1)
character*8 So (nr__2) ! (mJy) Unweighted mean optical flux density (1)
character*7 wSoMin (nr__2) ! (mJy) Weighted mean optical flux density
* (published) (1)
character*7 wSoMax (nr__2) ! (mJy) Weighted mean optical flux density
* (enlarged) (1)
real*8 AIr (nr__2) ! Unweighted radio activity index (1)
real*8 wAIrMin (nr__2) ! Weighted radio activity index (published) (1)
real*8 wAIrMax (nr__2) ! Weighted radio activity index (enlarged) (1)
character*7 AIo (nr__2) ! Unweighted optical activity index (1)
character*7 wAIoMin (nr__2) ! Weighted optical activity index (published) (1)
character*7 wAIoMax (nr__2) ! Weighted optical activity index (enlarged) (1)
*Note (1): 'Published' refers to the published error region scenario and
* 'enlarged' refers to the enlarged error region scenario. These are
* explained in Section 3.3 of the paper.
* The weighting scheme is explained in Section 3.2 of the paper.
C=============================================================================
C Loading file 'neutrino.dat' ! IceCat1+: high-energy neutrino events used
C Format for file interpretation
1 format(
+ A9,1X,F11.5,1X,F6.1,1X,F5.3,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F4.2,1X,F5.2,1X,F10.5,1X,F10.5,1X,F7.5,1X,F7.5,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'neutrino.dat')
write(6,*) '....Loading file: neutrino.dat'
do i__=1,283
read(1,'(A120)')ar__
read(ar__,1)
+ IC(i__),MJD(i__),E(i__),Sness(i__),RAdeg(i__),E_RAdeg(i__),
+ e_RAdeg_1(i__),DEdeg(i__),E_DEdeg(i__),e_DEdeg_1(i__),
+ OmMin(i__),OmMax(i__),wMin(i__),wMax(i__),Ref(i__)
c ..............Just test output...........
write(6,1)
+ IC(i__),MJD(i__),E(i__),Sness(i__),RAdeg(i__),E_RAdeg(i__),
+ e_RAdeg_1(i__),DEdeg(i__),E_DEdeg(i__),e_DEdeg_1(i__),
+ OmMin(i__),OmMax(i__),wMin(i__),wMax(i__),Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'blazars.dat' ! CGRaBS blazars: monitored by OVRO
C Format for file interpretation
2 format(A10,1X,F7.3,1X,F7.3,1X,A1,1X,A3,1X,A71)
C Effective file loading
open(unit=1,status='old',file=
+'blazars.dat')
write(6,*) '....Loading file: blazars.dat'
do i__=1,1157
read(1,'(A104)')ar__1
read(ar__1,2)
+ Blz(i__),RAdeg_1(i__),DEdeg_1(i__),BZCc(i__),FGLc(i__),
+ Alias(i__)
c ..............Just test output...........
write(6,2)
+ Blz(i__),RAdeg_1(i__),DEdeg_1(i__),BZCc(i__),FGLc(i__),
+ Alias(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'assoc.dat' ! All blazar-neutrino spatial associations
C Format for file interpretation
3 format(
+ A9,1X,A10,1X,F9.5,1X,F9.5,1X,F9.5,1X,F8.5,1X,F7.5,1X,F7.5,1X,
+ A8,1X,A7,1X,A7,1X,F7.5,1X,F7.5,1X,F7.5,1X,A7,1X,A7,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'assoc.dat')
write(6,*) '....Loading file: assoc.dat'
do i__=1,275
read(1,'(A148)')ar__2
read(ar__2,3)
+ IC_1(i__),Blz_1(i__),dBN(i__),dMin(i__),dMax(i__),Sr(i__),
+ wSrMin(i__),wSrMax(i__),So(i__),wSoMin(i__),wSoMax(i__),
+ AIr(i__),wAIrMin(i__),wAIrMax(i__),AIo(i__),wAIoMin(i__),
+ wAIoMax(i__)
c ..............Just test output...........
write(6,3)
+ IC_1(i__),Blz_1(i__),dBN(i__),dMin(i__),dMax(i__),Sr(i__),
+ wSrMin(i__),wSrMax(i__),So(i__),wSoMin(i__),wSoMax(i__),
+ AIr(i__),wAIrMin(i__),wAIrMax(i__),AIo(i__),wAIoMin(i__),
+ wAIoMax(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end