FORTRAN Generation
(/./ftp/cats/J/ApJ/893/L20)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/893/L20 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-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        ! Name (4FGLJHHMM.m+DDMM)
      character*11  Oname       ! Alternative name
      real*8        Eh          ! (GeV) Peak energy of high energy component
      real*4        logLumh     ! ([10-7W]) Peak luminosity of high energy component
*                                       (erg/s)
      real*8        Bem         ! (gauss) Magnetic field strength
      real*8        E_Bem       ! (gauss) Upper uncertainty on
*                                       magnetic field strength
      real*8        e_Bem_1     ! (gauss) Lower uncertainty on
*                                       magnetic field strength
      real*4        logPjmin    ! ([10-7W]) Minimum jet power, erg/s
      real*4        E_logPjmin  ! ([10-7W]) Upper uncertainty on minimum jet power
*                                       (erg/s)
      real*4        e_logPjmin_1 ! ([10-7W]) Lower uncertainty on minimum jet power
*                                       (erg/s)
      real*4        logMBH      ! ([Msun]) ? Black hole mass, solar units
      real*4        e_logMBH    ! ([Msun]) ? Uncertainty on black hole mass
      real*4        logLumd     ! ([10-7W]) ? Accretion disk luminosity (erg/s)
      real*4        e_logLumd   ! ([10-7W]) ? Uncertainty on accretion disk luminosity
*                                       (erg/s)
      real*4        logPBZ      ! ([10-7W]) Blandford-Znajek jet power (erg/s)
      integer*4     f_logPBZ    ! [1/2] Type of Blandford-Znajek jet power
*                                       estimate (1)
      real*4        e_logPBZ    ! ([10-7W]) Uncertainty on Blandford-Znajek jet power
*                                       (erg/s)
      real*4        logEpk      ! ([TeV]) Peak energy of muon neutrinos
      real*4        E_logEpk    ! ([TeV]) Upper uncertainty on peak energy of
*                                       muon neutrinos
      real*4        e_logEpk_1  ! ([TeV]) Lower uncertainty on peak energy of
*                                       muon neutrinos
      real*4        logFpk      ! ([TeV/cm2/s]) Peak flux of muon neutrinos
      real*4        E_logFpk    ! ([TeV/cm2/s]) Upper uncertainty on peak flux of
*                                       muon neutrinos
      real*4        e_logFpk_1  ! ([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,Oname,Eh,logLumh,Bem,E_Bem,e_Bem_1,logPjmin,E_logPjmin,
     +  e_logPjmin_1,logMBH,e_logMBH,logLumd,e_logLumd,logPBZ,
     +  f_logPBZ,e_logPBZ,logEpk,E_logEpk,e_logEpk_1,logFpk,E_logFpk,
     +  e_logFpk_1
        if(ar__(85:88) .EQ. '') logMBH = rNULL__
        if(ar__(90:93) .EQ. '') e_logMBH = rNULL__
        if(ar__(95:99) .EQ. '') logLumd = rNULL__
        if(ar__(101:104) .EQ. '') e_logLumd = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name,Oname,Eh,logLumh,Bem,E_Bem,e_Bem_1,logPjmin,E_logPjmin,
     +  e_logPjmin_1,logMBH,e_logMBH,logLumd,e_logLumd,logPBZ,
     +  f_logPBZ,e_logPBZ,logEpk,E_logEpk,e_logEpk_1,logFpk,E_logFpk,
     +  e_logFpk_1
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end