Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/196/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-Oct-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/196/1 The IPN supplement to the BATSE 5B catalog (Hurley+, 2011)
*================================================================================
*The Interplanetary Network supplement to the Burst And Transient Source
*Experiment 5B catalog of cosmic gamma-ray bursts.
* Hurley K., Briggs M.S., Kippen R.M., Kouveliotou C., Fishman G., Meegan C.,
* Cline T., Trombka J., McClanahan T., Boynton W., Starr R., McNutt R.,
* Boer M.
* <Astrophys. J. Suppl. Ser., 196, 1 (2011)>
* =2011ApJS..196....1H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Interplanetary Network (IPN) annuli
integer*4 nr__
parameter (nr__=344) ! Number of records
character*104 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Burst (nr__) ! Date of burst (YYMMDD)
character*1 m_Burst (nr__) ! [ABC- ] Flag when multiple bursts;
* added by CDS from figures
character*3 f_Burst (nr__) ! [a-w, ] Localisation problems (1)
character*8 Time (nr__) ! ("h:m:s") Universal Time of the burst at Earth
integer*4 Trig (nr__) ! [5591,8121] BATSE trigger number
* (<Trigger NNNN> in Simbad)
real*4 RAdeg (nr__) ! (deg) BATSE J2000 right ascension of the center of
* the error circle
real*4 DEdeg (nr__) ! (deg) BATSE J2000 declination of the center of the
* error circle (J2000)
real*4 Sig (nr__) ! (deg) Total 1{sigma} statistical BATSE error
* circle radius (2)
real*8 RA1deg (nr__) ! (deg) Right ascension of the center of the first
* IPN annulus (J2000) in the heliocentric frame
real*8 DE1deg (nr__) ! (deg) Declination of the center of the first IPN
* annulus (J2000) in the heliocentric frame
real*4 Rad1 (nr__) ! (deg) Angular radius of the first IPN annulus, in
* the heliocentric frame
real*4 DelR1 (nr__) ! (deg) Half-width of the first IPN annulus
* ({delta}R_IPN1_) (3)
real*8 RA2deg (nr__) ! (deg) ? Right ascension of the center of the second
* IPN annulus (J2000) in the heliocentric frame
real*8 DE2deg (nr__) ! (deg) ? Declination of the center of the second IPN
* annulus (J2000) in the heliocentric frame
real*4 Rad2 (nr__) ! (deg) ? Angular radius of the second IPN annulus, in
* the heliocentric frame
real*4 DelR2 (nr__) ! (deg) ? Half-width of the second IPN annulus
*Note (1): Flag as follows:
* a = Degenerate localization.
* b = Localized by the Rossi X-Ray Timing Explorer (RXTE) All-Sky Monitor
* (Smith et al. 1999ApJ...526..683S).
* c = Localized by the Advanced Satellite for Cosmology and Astrophysics (ASCA)
* (Murakami et al. 1997IAUC.6732....1M).
* d = Localized by BeppoSAX (Galama et al. 1997ApJ...486L...5G).
* e = Localized by BeppoSAX (Piro et al. 1997IAUC.6656....1P).
* f = Localized by the RXTE Proportional Counter Array (Marshall et al.
* 1997IAUC.6683....1C) and by ASCA (Murakami et al. 1997IAUC.6722....1M).
* g = Localized by BeppoSAX (Antonelli et al. 1997IAUC.6792....1A).
* h = Localized by BeppoSAX (Piro et al. 1997IAUC.6797....1P).
* i = Localized by BeppoSAX (in't Zand et al. 1998IAUC.6854....1I).
* j = Localized by BeppoSAX (Celidonio et al. 1998IAUC.6851....1C).
* k = Localized by BeppoSAX (in't Zand et al. 1998IAUC.6805....1I).
* l = Localized by BeppoSAX (Nicastro et al. 1998IAUC.6912....2N).
* m = Localized by BeppoSAX (Galama et al. 1998GCN...127....1G).
* n = Localized by BeppoSAX (Heise et al. 1999IAUC.7099....1H).
* o = Localized by the RXTE All-Sky Monitor (Smith et al. 1999GCN...187....1S).
* p = Localized by the RXTE Proportional Counter Array (Takeshima, T. &
* Marshall, F. 1999, GCN/RXTE PCA Burst Position Notice, 07 May).
* q = Localized by BeppoSAX (Dadina et al. 1999IAUC.7160....1D).
* r = Localized by BeppoSAX (Piro et al. 1999GCN...409....1P).
* s = Localized by BeppoSAX (Piro et al. 1999GCN...500....1P).
* t = Localized by BeppoSAX (in't Zand et al. 2000ApJ...545..266I).
* u = Localized by BeppoSAX (Smith et al. 1999IAUC.7301....2S).
* v = Localized by the RXTE Proportional Counter Array (Takeshima, T.,
* Markwardt, C., Marshall, F., Giblin, T., & Kippen, R. 1999, GCN Circ.,
* 478) and by the Chandra X-Ray Observatory (Piro et al.
* 1999GCN...393....1G).
* w = Localized by the RXTE All-Sky Monitor (Smith et al. 2000GCN...665....1S).
*Note (2): The approximate total 1{sigma} radius is obtained by adding
* 1.6{deg} in quadrature, but see Briggs et al. 1999ApJS..122..503B
* for an improved error model.
*Note (3): The 3{sigma} confidence annulus is given by R_IPN1_+/-{delta}R_IPN1_.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Interplanetary Network (IPN) error boxes
integer*4 nr__1
parameter (nr__1=405) ! Number of records
character*39 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Burst_1 (nr__1) ! Date of the burst (YYMMDD)
character*1 f_Burst_1 (nr__1) ! [o] o: alternative error box for trigger 7647
integer*4 Trig_1 (nr__1) ! [6380,8116] BATSE trigger number
integer*4 RAh (nr__1) ! (h) Hour of right ascension (J2000)
integer*4 RAm (nr__1) ! (min) Minute of right ascension (J2000)
real*4 RAs (nr__1) ! (s) Second of right ascension (J2000)
character*1 DE_ (nr__1) ! Sign of declination (J2000)
integer*4 DEd (nr__1) ! (deg) Degree of declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Arcminute of declination (J2000)
real*4 DEs (nr__1) ! (arcsec) Arcsecond of declination (J2000)
integer*4 loc (nr__1) ! [1/4] 1=center or 4=corner position (1)
*Note (1): This flag indicates whether the row contains:
* 1 = the coordinates of the center of the error box
* 4 = the coordinates of one of the four corners
C=============================================================================
C Loading file 'table2.dat' ! Interplanetary Network (IPN) annuli
C Format for file interpretation
1 format(
+ I6,A1,A3,1X,A8,1X,I4,1X,F6.2,1X,F6.2,1X,F5.2,1X,F7.3,1X,F7.3,
+ 1X,F6.3,1X,F6.3,1X,F7.3,1X,F7.3,1X,F6.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,344
read(1,'(A104)')ar__
read(ar__,1)
+ Burst(i__),m_Burst(i__),f_Burst(i__),Time(i__),Trig(i__),
+ RAdeg(i__),DEdeg(i__),Sig(i__),RA1deg(i__),DE1deg(i__),
+ Rad1(i__),DelR1(i__),RA2deg(i__),DE2deg(i__),Rad2(i__),
+ DelR2(i__)
if(ar__(76:82) .EQ. '') RA2deg(i__) = rNULL__
if(ar__(84:90) .EQ. '') DE2deg(i__) = rNULL__
if(ar__(92:97) .EQ. '') Rad2(i__) = rNULL__
if(ar__(99:104) .EQ. '') DelR2(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Burst(i__),m_Burst(i__),f_Burst(i__),Time(i__),Trig(i__),
+ RAdeg(i__),DEdeg(i__),Sig(i__),RA1deg(i__),DE1deg(i__),
+ Rad1(i__),DelR1(i__),RA2deg(i__),DE2deg(i__),Rad2(i__),
+ DelR2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Interplanetary Network (IPN) error boxes
C Format for file interpretation
2 format(
+ I6,A1,1X,I4,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,405
read(1,'(A39)')ar__1
read(ar__1,2)
+ Burst_1(i__),f_Burst_1(i__),Trig_1(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),loc(i__)
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg_1(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg_1(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ Burst_1(i__),f_Burst_1(i__),Trig_1(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),loc(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end