Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/893/L20 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-Jun-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. J/ApJ/893/L20 Proton synchrotron gamma-rays in blazars (Liodakis+, 2020)
*================================================================================
*Proton synchrotron gamma-rays and the energy crisis in blazars.
* Liodakis I., Petropoulou M.
* <Astrophys. J., 893, L20 (2020)>
* =2020ApJ...893L..20L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Parameter estimates for the sources in our sample
integer*4 nr__
parameter (nr__=145) ! Number of records
character*150 ar__ ! Full-size record
character*16 Name (nr__) ! Name (4FGLJHHMM.m+DDMM)
character*11 Oname (nr__) ! Alternative name
real*8 Eh (nr__) ! (GeV) Peak energy of high energy component
real*4 logLumh (nr__) ! ([10-7W]) Peak luminosity of high energy component
* (erg/s)
real*8 Bem (nr__) ! (gauss) Magnetic field strength
real*8 E_Bem (nr__) ! (gauss) Upper uncertainty on
* magnetic field strength
real*8 e_Bem_1 (nr__) ! (gauss) Lower uncertainty on
* magnetic field strength
real*4 logPjmin (nr__) ! ([10-7W]) Minimum jet power, erg/s
real*4 E_logPjmin (nr__) ! ([10-7W]) Upper uncertainty on minimum jet power
* (erg/s)
real*4 e_logPjmin_1(nr__) ! ([10-7W]) Lower uncertainty on minimum jet power
* (erg/s)
real*4 logMBH (nr__) ! ([Msun]) ? Black hole mass, solar units
real*4 e_logMBH (nr__) ! ([Msun]) ? Uncertainty on black hole mass
real*4 logLumd (nr__) ! ([10-7W]) ? Accretion disk luminosity (erg/s)
real*4 e_logLumd (nr__) ! ([10-7W]) ? Uncertainty on accretion disk luminosity
* (erg/s)
real*4 logPBZ (nr__) ! ([10-7W]) Blandford-Znajek jet power (erg/s)
integer*4 f_logPBZ (nr__) ! [1/2] Type of Blandford-Znajek jet power
* estimate (1)
real*4 e_logPBZ (nr__) ! ([10-7W]) Uncertainty on Blandford-Znajek jet power
* (erg/s)
real*4 logEpk (nr__) ! ([TeV]) Peak energy of muon neutrinos
real*4 E_logEpk (nr__) ! ([TeV]) Upper uncertainty on peak energy of
* muon neutrinos
real*4 e_logEpk_1 (nr__) ! ([TeV]) Lower uncertainty on peak energy of
* muon neutrinos
real*4 logFpk (nr__) ! ([TeV/cm2/s]) Peak flux of muon neutrinos
real*4 E_logFpk (nr__) ! ([TeV/cm2/s]) Upper uncertainty on peak flux of
* muon neutrinos
real*4 e_logFpk_1 (nr__) ! ([TeV/cm2/s]) Lower uncertainty on peak flux of
* muon neutrinos
*Note (1): Type of Blandford-Znajek jet power estimate as follows:
* 1 = estimate derived from core-shift measurements
* 2 = estimate derived without core-shift measurement
C=============================================================================
C Loading file 'table1.dat' ! Parameter estimates for the sources in our sample
C Format for file interpretation
1 format(
+ A16,1X,A11,1X,F7.4,1X,F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,
+ F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,I1,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,145
read(1,'(A150)')ar__
read(ar__,1)
+ Name(i__),Oname(i__),Eh(i__),logLumh(i__),Bem(i__),E_Bem(i__),
+ e_Bem_1(i__),logPjmin(i__),E_logPjmin(i__),e_logPjmin_1(i__),
+ logMBH(i__),e_logMBH(i__),logLumd(i__),e_logLumd(i__),
+ logPBZ(i__),f_logPBZ(i__),e_logPBZ(i__),logEpk(i__),
+ E_logEpk(i__),e_logEpk_1(i__),logFpk(i__),E_logFpk(i__),
+ e_logFpk_1(i__)
if(ar__(85:88) .EQ. '') logMBH(i__) = rNULL__
if(ar__(90:93) .EQ. '') e_logMBH(i__) = rNULL__
if(ar__(95:99) .EQ. '') logLumd(i__) = rNULL__
if(ar__(101:104) .EQ. '') e_logLumd(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),Oname(i__),Eh(i__),logLumh(i__),Bem(i__),E_Bem(i__),
+ e_Bem_1(i__),logPjmin(i__),E_logPjmin(i__),e_logPjmin_1(i__),
+ logMBH(i__),e_logMBH(i__),logLumd(i__),e_logLumd(i__),
+ logPBZ(i__),f_logPBZ(i__),e_logPBZ(i__),logEpk(i__),
+ E_logEpk(i__),e_logEpk_1(i__),logFpk(i__),E_logFpk(i__),
+ e_logFpk_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end