FORTRAN Generation
(/./ftp/cats/J/A_A/406/L9)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/406/L9 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 2024-Mar-28
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/406/L9        TeV Gamma-Rays from 1ES1959+650          (Aharonian+, 2003)
*================================================================================
*Detection of TeV Gamma-Rays from the BL Lac 1ES1959+650
*in its low states and during a major outburst in 2002.
*   Aharonian F., Akhperjanian A., Beilicke M., Bernloehr K., Boerst H.-G.,
*   Bojahr H., Bolz O., Coarasa T., Contreras J.L., Cortina J.,
*   Denninghoff S., , Fonseca M.V., Girma M., Goetting N., Heinzelmann G.,
*   Hermann G., Heusler A., Hofmann W., Horns D., Jung I., Kankanyan R.,
*   Kestel M., Kohnle A., Konopelko A., Kornmeyer H., Kranich D., Lampeitl H.,
*   Lopez M., Lorenz E., Lucarelli F., Mang O., Meyer H., Mirzoyan R.,
*   Moralejo A., Ona-Wilhelmi E., Panter M., Plyasheshnikov A., Puehlhofer G.,
*   delosReyes R., Rhode W., Ripken J., Robrade J., Rowell G., Sahakian V.,
*   Samorski M., Schilling M., Siems M., Sobzynska D., Stamm W., Tluczykont M.,
*   Vitale V., Voelk H.J., Wiedner C.A., Wittek W.
*  <Astron. Astrophys. 406, L9 (2003)>
*  =2003A&A...406L...9A
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table4.dat'	! Diurnal integral flux values of 1ES1959+650 above
                            2TeV from 2000 to 2002 as measured with the HEGRA
                            stereoscopic system of 5 imaging atmospheric
                            Cherenkov telescopes.

      integer*4 nr__
      parameter (nr__=127)	! Number of records
      character*57 ar__   	! Full-size record

      real*8        MJD        (nr__) ! (d) Mean MJD value (JD-2400000.5)
*                                       of HEGRA observation
      real*4        Duration   (nr__) ! (h) Duration of HEGRA observation
      real*4        Flux       (nr__) ! (10-12/cm2/s) Integral flux above 2 TeV (1)
      real*4        e_Flux     (nr__) ! (10-12/cm2/s) Error on integral flux above 2 TeV (1)
      real*8        MJD1       (nr__) ! (d) Start MJD value of HEGRA observation
      real*8        MJD2       (nr__) ! (d) End MJD value of HEGRA observation
*Note (1):
*    A negative Flux value indicates an upper limit of the flux;
*    its error e_Flux is then set to 0.0

C=============================================================================

C  Loading file 'table4.dat'	! Diurnal integral flux values of 1ES1959+650 above
*                            2TeV from 2000 to 2002 as measured with the HEGRA
*                            stereoscopic system of 5 imaging atmospheric
*                            Cherenkov telescopes.

C  Format for file interpretation

    1 format(2X,F10.4,1X,F5.2,2X,F5.2,2X,F4.2,3X,F10.4,3X,F10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,127
        read(1,'(A57)')ar__
        read(ar__,1)
     +  MJD(i__),Duration(i__),Flux(i__),e_Flux(i__),MJD1(i__),
     +  MJD2(i__)
c    ..............Just test output...........
        write(6,1)
     +  MJD(i__),Duration(i__),Flux(i__),e_Flux(i__),MJD1(i__),
     +  MJD2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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