Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/229/31 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-May-10
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/ApJS/229/31 IPN supplement to the 2nd Fermi GBM catalog (Hurley+, 2017)
*================================================================================
*The InterPlanetary Network supplement to the second Fermi GBM catalog of
*cosmic Gamma-Ray bursts.
* Hurley K., Aptekar R.L., Golenetskii S.V., Frederiks D.D., Svinkin D.S.,
* Pal'shin V.D., Briggs M.S., Meegan C., Connaughton V., Goldsten J.,
* Boynton W., Fellows C., Harshman K., Mitrofanov I.G., Golovin D.V.,
* Kozyrev A.S., Litvak M.L., Sanin A.B., Rau A., von Kienlin A., Zhang X.,
* Yamaoka K., Fukazawa Y., Ohno M., Tashiro M., Terada Y., Barthelmy S.,
* Cline T., Gehrels N., Cummings J., Krimm H.A., Smith D.M., Del Monte E.,
* Feroci M., Marisaldi M.
* <Astrophys. J. Suppl. Ser., 229, 31-31 (2017)>
* =2017ApJS..229...31H (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! IPN/GBM gamma-ray bursts
integer*4 nr__
parameter (nr__=431) ! Number of records
character*71 ar__ ! Full-size record
character*11 obs_date (nr__) ! ("Y/M/D") UT date of observation
character*8 obs_time (nr__) ! ("h:m:s") Universal Time of observation (1)
character*6 n_Fermi (nr__) ! Note on Fermi (2)
character*9 Fermi (nr__) ! Fermi GBM designation (YYMMDDddd)
character*30 Inst (nr__) ! Observation Instrument(s) (3)
character*1 f_Inst (nr__) ! Flag on Inst (4)
character*1 f_Fermi (nr__) ! Flag on Fermi; column added by CDS (5)
*Note (1): Universal time is the trigger time of a near-Earth spacecraft .
*Note (2):
* UNCERT = Four events were listed as UNCERTAIN in the GBM catalog -- GBM:
* UNCERT10102526 UNCERT11092348 UNCERT11123162 UNCERT12011729; however,
* they are likely to be valid cosmic events.
*Note (3): Instruments observing GBM as follows:
* AGI = Astro-rivelatore Gamma a Immagini LEggero (AGILE);
* INT = International Gamma-Ray Laboratory;
* Kon = Konus-Wind;
* MES = Mercury Surface, Space Environment, Geochemistry, and Ranging
* mission;
* MO = Mars Odyssey;
* RHE = Ramaty High Energy Solar Spectroscopic Imager;
* Suz = Suzaku;
* Swi = Swift.
*Note (4): Flags on Observation Instrument(s) as follows:
* d = Burst was outside the coded field of view of the BAT, and not
* localized by it;
* e = Burst was localized by Swift-BAT; IPN triangulation cannot improve
* on this localization;
* f = Burst was localized by INTEGRAL-IBIS; IPN triangulation cannot
* improve on this localization.
*Note (5):
* n = 101127093 has been added by CDS in this table to match the object listed
* in tables 5 and 6.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! IPN annuli for 165 bursts
integer*4 nr__1
parameter (nr__1=165) ! Number of records
character*79 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*9 Fermi_1 (nr__1) ! Fermi Gamma-Ray Burst Designation (YYMMDDddd)
real*8 RAdeg (nr__1) ! (deg) Right Ascension in decimal degrees (J2000),
* 1st IPN annulus
real*8 DEdeg (nr__1) ! (deg) Declination in decimal degrees (J2000),
* 1st IPN annulus
real*8 RIPN1 (nr__1) ! (deg) [5.4/90] Angular radius of the 1st IPN annulus
real*4 dRIPN1 (nr__1) ! (deg) [0.04/6.3] Half-width dRIPN1 of the
* 1st IPN annulus
real*8 RA2deg (nr__1) ! (deg) ? Right Ascension in decimal degrees (J2000),
* 2nd IPN annulus
real*8 DE2deg (nr__1) ! (deg) ? Declination in decimal degrees (J2000),
* 2nd IPN annulus
real*8 RIPN2 (nr__1) ! (deg) ? Angular radius of the 2nd IPN annulus
real*8 dRIPN2 (nr__1) ! (deg) ? Half-width dRIPN2 of the 2nd IPN annulus
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! IPN error box coordinates for 164 bursts
integer*4 nr__2
parameter (nr__2=164) ! Number of records
character*115 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*9 Fermi_2 (nr__2) ! Fermi Gamma-Ray Burst Designation (YYMMDDddd)
real*8 RAdeg_1 (nr__2) ! (deg) Right Ascension in decimal degrees (J2000),
* error box center
real*8 DEdeg_1 (nr__2) ! (deg) Declination in decimal degrees (J2000),
* error box center
real*8 RA1deg (nr__2) ! (deg) Right Ascension in decimal degrees (J2000),
* 1st error box corner
real*8 DE1deg (nr__2) ! (deg) Declination in decimal degrees (J2000),
* 1st error box corner
real*8 RA2deg_1 (nr__2) ! (deg) Right Ascension in decimal degrees (J2000),
* 2nd error box corner
real*8 DE2deg_1 (nr__2) ! (deg) Declination in decimal degrees (J2000),
* 2nd error box corner
real*8 RA3deg (nr__2) ! (deg) Right Ascension in decimal degrees (J2000),
* 3rd error box corner
real*8 DE3deg (nr__2) ! (deg) Declination in decimal degrees (J2000),
* 3rd error box corner
real*8 RA4deg (nr__2) ! (deg) ? Right Ascension in decimal degrees (J2000),
* 4th error box corner
real*8 DE4deg (nr__2) ! (deg) ? Declination in decimal degrees (J2000),
* 4th error box corner
real*4 aGBM (nr__2) ! (deg2) [35/13800] Area of the 3{sigma} GBM error
* contour, in square degrees
real*8 aIPN (nr__2) ! (deg2) [0.02/380] Area of the 3{sigma} IPN error box,
* in square degrees (1)
real*8 Ratio (nr__2) ! [1/37900] Ratio of the GBM area to the IPN area
*Note (1): If there is a single IPN annulus, this is the area of the entire
* annulus; if there are two annuli, this is the area of the error box
* formed by their intersection.
C=============================================================================
C Loading file 'table2.dat' ! IPN/GBM gamma-ray bursts
C Format for file interpretation
1 format(A11,1X,A8,1X,A6,A9,1X,A30,1X,A1,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,431
read(1,'(A71)')ar__
read(ar__,1)
+ obs_date(i__),obs_time(i__),n_Fermi(i__),Fermi(i__),Inst(i__),
+ f_Inst(i__),f_Fermi(i__)
c ..............Just test output...........
write(6,1)
+ obs_date(i__),obs_time(i__),n_Fermi(i__),Fermi(i__),Inst(i__),
+ f_Inst(i__),f_Fermi(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! IPN annuli for 165 bursts
C Format for file interpretation
2 format(
+ 3X,A9,1X,F8.4,1X,F8.4,1X,F7.4,1X,F6.4,1X,F8.4,1X,F8.4,1X,F7.4,
+ 1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,165
read(1,'(A79)')ar__1
read(ar__1,2)
+ Fermi_1(i__),RAdeg(i__),DEdeg(i__),RIPN1(i__),dRIPN1(i__),
+ RA2deg(i__),DE2deg(i__),RIPN2(i__),dRIPN2(i__)
if(ar__1(47:54) .EQ. '') RA2deg(i__) = rNULL__
if(ar__1(56:63) .EQ. '') DE2deg(i__) = rNULL__
if(ar__1(65:71) .EQ. '') RIPN2(i__) = rNULL__
if(ar__1(73:79) .EQ. '') dRIPN2(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Fermi_1(i__),RAdeg(i__),DEdeg(i__),RIPN1(i__),dRIPN1(i__),
+ RA2deg(i__),DE2deg(i__),RIPN2(i__),dRIPN2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! IPN error box coordinates for 164 bursts
C Format for file interpretation
3 format(
+ 3X,A9,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,
+ 1X,F7.3,1X,F7.3,1X,F7.3,1X,F6.0,1X,F7.3,1X,F7.1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,164
read(1,'(A115)')ar__2
read(ar__2,3)
+ Fermi_2(i__),RAdeg_1(i__),DEdeg_1(i__),RA1deg(i__),
+ DE1deg(i__),RA2deg_1(i__),DE2deg_1(i__),RA3deg(i__),
+ DE3deg(i__),RA4deg(i__),DE4deg(i__),aGBM(i__),aIPN(i__),
+ Ratio(i__)
if(ar__2(78:84) .EQ. '') RA4deg(i__) = rNULL__
if(ar__2(86:92) .EQ. '') DE4deg(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Fermi_2(i__),RAdeg_1(i__),DEdeg_1(i__),RA1deg(i__),
+ DE1deg(i__),RA2deg_1(i__),DE2deg_1(i__),RA3deg(i__),
+ DE3deg(i__),RA4deg(i__),DE4deg(i__),aGBM(i__),aIPN(i__),
+ Ratio(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end