FORTRAN Generation
(/./ftp/cats/J/ApJ/814/1)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/814/1 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-15
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/814/1 GRB 120326A, 100418A & 100901A multi-wavelength obs. (Laskar+, 2015)
*================================================================================
*Energy injection in gamma-ray burst afterglows.
*    Laskar T., Berger E., Margutti R., Perley D., Zauderer B.A., Sari R.,
*    Fong W.-F.
*   <Astrophys. J., 814, 1 (2015)>
*   =2015ApJ...814....1L    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Swift/UVOT Observations of GRB 120326A

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

      real*8        t          (nr__) ! (d) [0.001/19.5] Time since discovery; t-t_0_
      character*5   Filt       (nr__) ! Filter used in the observation (White, u, b,
*                                   v, uvm2, uvw1 or uvw2)
      real*4        Freq       (nr__) ! (Hz) [5.5e+14/1.5e+15] Frequency of Filter peak
      real*8        Flux       (nr__) ! (uJy) [0.2/150] Flux density (1)
      real*4        e_Flux     (nr__) ! (uJy) [0.5/15] Uncertainty in Flux
      integer*4     Det_       (nr__) ! [0/1] Detection flag; 1 = Detection
*Note (1): In cases of non-detections, we report the formal flux density
*          measurement from aperture photometry.

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

C  Declarations for 'table2.dat'	! Optical Observations of GRB 120326A

      integer*4 nr__1
      parameter (nr__1=138)	! Number of records
      character*112 ar__1  	! Full-size record

      real*8        t_1        (nr__1) ! (d) [0.0002/21.4] Time since discovery; t-t_0_
      character*12  Obs        (nr__1) ! Observatory identifier
      character*11  T_I        (nr__1) ! Telescope/Instrument description
      character*2   Filt_1     (nr__1) ! Filter used in the observation (CR, JHK, RI,
*                                   RcIc, g'r'i'z', bv or Y')
      real*4        Freq_1     (nr__1) ! (Hz) [1.3e+14/7e+14] Effective frequency
      real*8        Flux_1     (nr__1) ! (uJy) [1.6/4040] Flux density (1)
      real*8        e_Flux_1   (nr__1) ! (uJy) [0.09/1350] Uncertainty in Flux
      integer*4     Det__1     (nr__1) ! [0/1] Detection flag; 1 = Detection
      character*25  Aut        (nr__1) ! First author's name(s) of the reference (2)
      character*19  BibCode    (nr__1) ! Bibcode of the reference (2)
*Note (1): In cases of non-detections, we report the 3{sigma} upper limit
*          on the flux density.
*Note (2): We note that the photometry reported in Table 1 of
*          Urata et al. (2014ApJ...789..146U) yields an unphysical spectral index
*          of {Beta}~-4.5 for GRB 120326A. We assume that these numbers have been
*          scaled by the same factors as reported in Figure 4 of their paper, and
*          divide the gr'Rizy light curves by a factor of 0.5, 1, 1, 2, 4, and 6,
*          respectively. We further assume that the data have not been corrected
*          for galactic extinction.

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

C  Declarations for 'table3.dat'	! Millimeter and Radio Observations of GRB 120326A

      integer*4 nr__2
      parameter (nr__2=34)	! Number of records
      character*57 ar__2  	! Full-size record

      character*14  Date       (nr__2) ! ("Y/M/D") UT date of the observation
      real*8        t_2        (nr__2) ! (d) [0.5/122.1] Time since discovery; t-t_0_
      character*5   Obs_1      (nr__2) ! Observatory identifier (CARMA, SMA or VLA/B,
*                                  VLA/C) (1)
      real*4        Freq_2     (nr__2) ! (GHz) [5/222] Frequency of the observation
      real*4        Int        (nr__2) ! (min) [4.5/386] Integration time
      character*1   l_Flux     (nr__2) ! Limit flag on Flux
      real*4        Flux_2     (nr__2) ! (mJy) [0.05/3.4] Flux density
      real*4        e_Flux_2   (nr__2) ! (mJy) [0.01/0.9] The 1{sigma} statistical
*                                  uncertainty in Flux
      integer*4     Det__2     (nr__2) ! [0/1] Detection flag; 1 = Detection
*Note (1): The letter following "VLA" indicates the array configuration. The
*          PI of the SMA data set is Y. Urata.

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

C  Declarations for 'table6.dat'	! Swift/UVOT Observations of GRB 100418A

      integer*4 nr__3
      parameter (nr__3=64)	! Number of records
      character*40 ar__3  	! Full-size record

      real*8        t_3        (nr__3) ! (d) [0.001/34.4] Time since discovery; t-t_0_
      character*5   Filt_2     (nr__3) ! Filter used in the observation White, u, b,
*                                   v, uvm2, uvw1 or uvw2)
      real*4        Freq_3     (nr__3) ! (Hz) [5.5e+14/1.5e+15] Frequency of Filter peak
      real*8        Flux_3     (nr__3) ! (uJy) [1.3/291] Flux density (1)
      real*4        e_Flux_3   (nr__3) ! (uJy) [0.3/76] Uncertainty in Flux
      integer*4     Det__3     (nr__3) ! [0/1] Detection flag; 1 = Detection
*Note (1): In cases of non-detections, we report the formal flux density
*          measurement from aperture photometry.

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

C  Declarations for 'table7.dat'	! Swift/UVOT Observations of GRB 100901A

      integer*4 nr__4
      parameter (nr__4=125)	! Number of records
      character*40 ar__4  	! Full-size record

      real*8        t_4        (nr__4) ! (d) [0.001/34.4] Time since discovery; t-t_0_
      character*5   Filt_3     (nr__4) ! Filter used in the observation White, u, b,
*                                   v, uvm2, uvw1 or uvw2)
      real*4        Freq_4     (nr__4) ! (Hz) [5.5e+14/1.5e+15] Frequency of Filter peak
      real*8        Flux_4     (nr__4) ! (uJy) [1.3/291] Flux density (1)
      real*4        e_Flux_4   (nr__4) ! (uJy) [0.3/76] Uncertainty in Flux
      integer*4     Det__4     (nr__4) ! [0/1] Detection flag; 1 = Detection
*Note (1): In cases of non-detections, we report the formal flux density
*          measurement from aperture photometry.

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

C  Declarations for 'table8.dat'	! Optical Observations of GRB 100901A

      integer*4 nr__5
      parameter (nr__5=263)	! Number of records
      character*112 ar__5  	! Full-size record

      real*8        t_5        (nr__5) ! (d) [0.0002/21.4] Time since discovery; t-t_0_
      character*12  Obs_2      (nr__5) ! Observatory identifier
      character*11  T_I_1      (nr__5) ! Telescope/Instrument description
      character*2   Filt_4     (nr__5) ! Filter used in the observation (CR, JHK, RI,
*                                   RcIc, g'r'i'z', bv or Y')
      real*4        Freq_5     (nr__5) ! (Hz) [1.3e+14/7e+14] Effective frequency
      real*8        Flux_5     (nr__5) ! (uJy) [1.6/4040] Flux density (1)
      real*8        e_Flux_5   (nr__5) ! (uJy) [0.09/1350] Uncertainty in Flux
      integer*4     Det__5     (nr__5) ! [0/1] Detection flag; 1 = Detection
      character*25  Aut_1      (nr__5) ! First author's name(s) of the reference (2)
      character*19  BibCode_1  (nr__5) ! Bibcode of the reference (2)
*Note (1): In cases of non-detections, we report the 3{sigma} upper limit
*          on the flux density.
*Note (2): We note that the photometry reported in Table 1 of
*          Urata et al. (2014ApJ...789..146U) yields an unphysical spectral index
*          of {Beta}~-4.5 for GRB 120326A. We assume that these numbers have been
*          scaled by the same factors as reported in Figure 4 of their paper, and
*          divide the gr'Rizy light curves by a factor of 0.5, 1, 1, 2, 4, and 6,
*          respectively. We further assume that the data have not been corrected
*          for galactic extinction.

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

C  Loading file 'table1.dat'	! Swift/UVOT Observations of GRB 120326A

C  Format for file interpretation

    1 format(F8.5,1X,A5,1X,E8.2,1X,F7.3,1X,F6.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,52
        read(1,'(A40)')ar__
        read(ar__,1)
     +  t(i__),Filt(i__),Freq(i__),Flux(i__),e_Flux(i__),Det_(i__)
c    ..............Just test output...........
        write(6,1)
     +  t(i__),Filt(i__),Freq(i__),Flux(i__),e_Flux(i__),Det_(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Optical Observations of GRB 120326A

C  Format for file interpretation

    2 format(
     +  F9.6,1X,A12,1X,A11,1X,A2,1X,E8.1,1X,F7.2,1X,F9.4,1X,I1,1X,A25,
     +  1X,A19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,138
        read(1,'(A112)')ar__1
        read(ar__1,2)
     +  t_1(i__),Obs(i__),T_I(i__),Filt_1(i__),Freq_1(i__),
     +  Flux_1(i__),e_Flux_1(i__),Det__1(i__),Aut(i__),BibCode(i__)
c    ..............Just test output...........
        write(6,2)
     +  t_1(i__),Obs(i__),T_I(i__),Filt_1(i__),Freq_1(i__),
     +  Flux_1(i__),e_Flux_1(i__),Det__1(i__),Aut(i__),BibCode(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Millimeter and Radio Observations of GRB 120326A

C  Format for file interpretation

    3 format(
     +  A14,1X,F7.3,1X,A5,1X,F5.1,1X,F6.2,1X,A1,1X,F5.3,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,34
        read(1,'(A57)')ar__2
        read(ar__2,3)
     +  Date(i__),t_2(i__),Obs_1(i__),Freq_2(i__),Int(i__),
     +  l_Flux(i__),Flux_2(i__),e_Flux_2(i__),Det__2(i__)
c    ..............Just test output...........
        write(6,3)
     +  Date(i__),t_2(i__),Obs_1(i__),Freq_2(i__),Int(i__),
     +  l_Flux(i__),Flux_2(i__),e_Flux_2(i__),Det__2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Swift/UVOT Observations of GRB 100418A

C  Format for file interpretation

    4 format(F8.5,1X,A5,1X,E8.2,1X,F7.3,1X,F6.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,64
        read(1,'(A40)')ar__3
        read(ar__3,4)
     +  t_3(i__),Filt_2(i__),Freq_3(i__),Flux_3(i__),e_Flux_3(i__),
     +  Det__3(i__)
c    ..............Just test output...........
        write(6,4)
     +  t_3(i__),Filt_2(i__),Freq_3(i__),Flux_3(i__),e_Flux_3(i__),
     +  Det__3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Swift/UVOT Observations of GRB 100901A

C  Format for file interpretation

    5 format(F8.5,1X,A5,1X,E8.2,1X,F7.3,1X,F6.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,125
        read(1,'(A40)')ar__4
        read(ar__4,5)
     +  t_4(i__),Filt_3(i__),Freq_4(i__),Flux_4(i__),e_Flux_4(i__),
     +  Det__4(i__)
c    ..............Just test output...........
        write(6,5)
     +  t_4(i__),Filt_3(i__),Freq_4(i__),Flux_4(i__),e_Flux_4(i__),
     +  Det__4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Optical Observations of GRB 100901A

C  Format for file interpretation

    6 format(
     +  F9.6,1X,A12,1X,A11,1X,A2,1X,E8.1,1X,F7.2,1X,F9.4,1X,I1,1X,A25,
     +  1X,A19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,263
        read(1,'(A112)')ar__5
        read(ar__5,6)
     +  t_5(i__),Obs_2(i__),T_I_1(i__),Filt_4(i__),Freq_5(i__),
     +  Flux_5(i__),e_Flux_5(i__),Det__5(i__),Aut_1(i__),
     +  BibCode_1(i__)
c    ..............Just test output...........
        write(6,6)
     +  t_5(i__),Obs_2(i__),T_I_1(i__),Filt_4(i__),Freq_5(i__),
     +  Flux_5(i__),e_Flux_5(i__),Det__5(i__),Aut_1(i__),
     +  BibCode_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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