Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/473/1512 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-Oct-06
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/MNRAS/473/1512 AMI 15.7GHz GRB catalogue (Anderson+, 2018)
*================================================================================
*The Arcminute Microkelvin Imager catalogue of gamma-ray burst afterglows
*at 15.7 GHz.
* Anderson G.E., Staley T.D., van der Horst A.J., Fender R.P., Rowlinson A.,
* Mooley K.P., Broderick J.W., Wijers R.A.M.J., Rumsey C., Titterington D.J.
* <Mon. Not. R. Astron. Soc., 473, 1512-1536 (2018)>
* =2018MNRAS.473.1512A (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! The AMI 15.7GHz GRB catalogue: GRBs that have
been detected or possibly detected with AMI
integer*4 nr__
parameter (nr__=227) ! Number of records
character*83 ar__ ! Full-size record
character*7 GRB (nr__) ! Name of GRB (1)
character*2 Flag (nr__) ! Radio detection flag (2)
character*9 Tel (nr__) ! Best Swift telescope position (3)
character*10 Date (nr__) ! ("date") Date of AMI observation (4)
character*1 n_Date (nr__) ! [C] C for Concat (4)
real*8 Start (nr__) ! (d) ? Start of observation in Modified Julian
* Date (MJD) format
real*4 Int (nr__) ! (h) ? Length of AMI observation
real*4 DayPburst (nr__) ! (d) ? Number of days post-burst since the
* start of the AMI observation
real*4 MinPburst (nr__) ! (min) ? Number of minutes post-burst since the
* start of the AMI observation (5)
character*1 Det (nr__) ! [*] * indicates the observations where the
* listed peak flux is a radio detection (6)
real*4 Fpeak (nr__) ! (mJy/beam) ? Peak flux at 15.7GHz of the
* radio afterglow (7)
real*4 e_Fpeak (nr__) ! (mJy/beam) ? Error in the peak flux at 15.7GHz
* measurement (1 sigma) (8)
real*4 Signi (nr__) ! Significance of the flux measured in
* units of sigma above the local RMS
real*4 RMS (nr__) ! (mJy/beam) The global RMS of the middle quarter
* of the AMI image
*Note (1): All GRBs in this table are classed as long GRBs. All GRB observations
* from GRB 130907A onwards followed the updated ALARRM strategy (see
* Section 2.1). GRB 130702A and GRB 140320C were discovered by Fermi
* and Integral, respectively.
*Note (2): Radio detection abbreviations as follows:
* A = new radio GRB discovered with AMI
* AC = new radio GRB discovered with AMI that was only detected in the
* concatenated image
* R = radio afterglow first detected with another radio telescope
* P = possible new candidate radio GRB discovered with AMI
*Note (3): Best Swift telescope position used to search for a radio counterpart
* BAT = Swift Burst Alert Telescope
* XRT = Swift X-ray Telescope
* UVOT = Swift Ultraviolet/Optical Telescope.
* For GRB 150413A, the best optical position provided by the MASTER II
* robotic telescope (Ivanov et al. 2015, GCN, 17689) was used instead.
*Note (4): Date of the AMI observation in yyyy-mm-dd.
* C = Concat = concatenation of all epochs with the same pointing (in the case of
* GRB 120320A, two observations were taken with different pointings so
* the concatenated image is the same as the 2012-04-05 epoch).
*Note (5): This time is only provided for those AMI observations where the start
* time was <=0.01 days post-burst.
*Note (6): * indicates those GRBs for which the listed AMI peak flux is the
* measured flux of a radio source (signal-to-noise >=4) that is coincident with
* the best known Swift position of the GRB (i.e. within 3 positional sigma).
* All other listed fluxes are derived from a forced Gaussian fit at the best
* known Swift position.
*Note (7): Peak flux density as reported by PySE.
*Note (8): The 1 sigma error bar is the flux error output by PySE
* added in quadrature to the AMI 5% calibration error.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! The AMI 15.7GHz GRB catalogue: GRBs that were not
detected with AMI, have a possible steady source
association or a concatenated detection
integer*4 nr__1
parameter (nr__1=783) ! Number of records
character*91 ar__1 ! Full-size record
character*7 GRB_1 (nr__1) ! Name of GRB (1)
character*1 Type (nr__1) ! Short (S) or long (L) duration GRB
character*2 Flag_1 (nr__1) ! Radio detection flag (2)
character*6 Tel_1 (nr__1) ! Best Swift telescope position (3)
character*10 Date_1 (nr__1) ! ("date") Date of AMI observation (4)
character*1 n_Date_1 (nr__1) ! [C] C for Concat (4)
real*8 Start_1 (nr__1) ! (d) ? Start of observation in Modified Julian
* Date (MJD) format
real*4 Int_1 (nr__1) ! (h) ? Length of AMI observation
real*4 DayPburst_1(nr__1) ! (d) ? Number of days post-burst since the start
* of the AMI observation
real*4 MinPburst_1(nr__1) ! (min) ? Number of minutes post-burst since the
* start of the AMI observation (5)
character*1 Det_1 (nr__1) ! [*] * indicates the observations where the
* listed peak flux is a radio detection (6)
real*4 Fpeak_1 (nr__1) ! (mJy/beam) Peak flux 15.7GHz of the radio afterglow (7)
real*8 e_Fpeak_1 (nr__1) ! (mJy/beam) Error in the peak flux 15.7GHz measurement
* (1 sigma) (8)
real*4 Signi_1 (nr__1) ! Significance of the flux measured in units
* of sigma above the local RMS
real*4 RMS_1 (nr__1) ! (mJy/beam) The global RMS of the middle quarter of
* the AMI image
*Note (1): All GRB observations from GRB 130813A onwards followed the updated
* ALARRM strategy (see Section 2.1). GRB 140320B was discovered by Integral and
* GRB 140508A, GRB 140606B, GRB 140623A, and GRB 140801A were discovered
* by Fermi.
*Note (2): Radio detection abbreviations
* R = radio afterglow first detected with another radio telescope
* S = AMI detected coincident radio source is likely steady
* C = coincident radio source detected in AMI concatenated image
* CS = coincident radio source detected in AMI concatenated image that is
* confirmed to be a steady source
* N = no detection
*Note (3): Best Swift telescope position used to search for a radio counterpart
* BAT = Swift Burst Alert Telescope
* XRT = Swift X-ray Telescope
* UVOT = Swift Ultraviolet/Optical Telescope.
* For GRB 130521A, the best optical position provided by Skynet/PROMPT
* (James et al., 2013, GCN, 14713) was used instead.
*Note (4): Date of the AMI observation in yyyy-mm-dd.
* C = Concat = concatenation of all epochs with the same pointing
*Note (5): This time is only provided for those AMI observations where the
* start time was <=0.01 days post-burst
*Note (6): * indicates those GRBs for which the listed AMI peak flux is the
* measured flux of a radio source (signal-to-noise >=4) that is coincident with
* the best known Swift position of the GRB (i.e. within 3 positional sigma).
* All other listed fluxes are derived from a forced Gaussian fit at the best
* known Swift position.
*Note (7): Peak flux density as reported by PySE. The PySE fitting algorithm
* failed on the 2012-04-05 and 2012-04-08 AMI observations of GRB 120311A and
* all of the GRB 121128A observations. The forced fits were instead conducted
* using MIRIAD and the reported significance is just the forced fitted peak flux
* divided by the global RMS.
*Note (8): The 1 sigma error bar is the flux error output by PySE added in
* quadrature to the AMI 5% calibration error.
C=============================================================================
C Loading file 'table1.dat' ! The AMI 15.7GHz GRB catalogue: GRBs that have
* been detected or possibly detected with AMI
C Format for file interpretation
1 format(
+ A7,1X,A2,1X,A9,1X,A10,A1,F8.2,1X,F3.1,1X,F5.2,1X,F5.2,3X,A1,
+ 1X,F5.2,1X,F6.2,1X,F4.1,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,227
read(1,'(A83)')ar__
read(ar__,1)
+ GRB(i__),Flag(i__),Tel(i__),Date(i__),n_Date(i__),Start(i__),
+ Int(i__),DayPburst(i__),MinPburst(i__),Det(i__),Fpeak(i__),
+ e_Fpeak(i__),Signi(i__),RMS(i__)
if(ar__(33:40) .EQ. '') Start(i__) = rNULL__
if(ar__(42:44) .EQ. '') Int(i__) = rNULL__
if(ar__(46:50) .EQ. '') DayPburst(i__) = rNULL__
if(ar__(52:56) .EQ. '') MinPburst(i__) = rNULL__
if(ar__(62:66) .EQ. '') Fpeak(i__) = rNULL__
if(ar__(68:73) .EQ. '') e_Fpeak(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ GRB(i__),Flag(i__),Tel(i__),Date(i__),n_Date(i__),Start(i__),
+ Int(i__),DayPburst(i__),MinPburst(i__),Det(i__),Fpeak(i__),
+ e_Fpeak(i__),Signi(i__),RMS(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! The AMI 15.7GHz GRB catalogue: GRBs that were not
* detected with AMI, have a possible steady source
* association or a concatenated detection
C Format for file interpretation
2 format(
+ A7,1X,A1,1X,A2,2X,A6,1X,A10,A1,F8.2,1X,F4.1,1X,F6.2,3X,F5.2,
+ 2X,A1,1X,F6.2,1X,F8.2,2X,F4.1,2X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,783
read(1,'(A91)')ar__1
read(ar__1,2)
+ GRB_1(i__),Type(i__),Flag_1(i__),Tel_1(i__),Date_1(i__),
+ n_Date_1(i__),Start_1(i__),Int_1(i__),DayPburst_1(i__),
+ MinPburst_1(i__),Det_1(i__),Fpeak_1(i__),e_Fpeak_1(i__),
+ Signi_1(i__),RMS_1(i__)
if(ar__1(33:40) .EQ. '') Start_1(i__) = rNULL__
if(ar__1(42:45) .EQ. '') Int_1(i__) = rNULL__
if(ar__1(47:52) .EQ. '') DayPburst_1(i__) = rNULL__
if(ar__1(56:60) .EQ. '') MinPburst_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ GRB_1(i__),Type(i__),Flag_1(i__),Tel_1(i__),Date_1(i__),
+ n_Date_1(i__),Start_1(i__),Int_1(i__),DayPburst_1(i__),
+ MinPburst_1(i__),Det_1(i__),Fpeak_1(i__),e_Fpeak_1(i__),
+ Signi_1(i__),RMS_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end