Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/196/1 into FORTRAN code for reading data files line by line.
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 ! Date of burst (YYMMDD)
character*1 m_Burst ! [ABC- ] Flag when multiple bursts;
* added by CDS from figures
character*3 f_Burst ! [a-w, ] Localisation problems (1)
character*8 Time ! ("h:m:s") Universal Time of the burst at Earth
integer*4 Trig ! [5591,8121] BATSE trigger number
* (<Trigger NNNN> in Simbad)
real*4 RAdeg ! (deg) BATSE J2000 right ascension of the center of
* the error circle
real*4 DEdeg ! (deg) BATSE J2000 declination of the center of the
* error circle (J2000)
real*4 Sig ! (deg) Total 1{sigma} statistical BATSE error
* circle radius (2)
real*8 RA1deg ! (deg) Right ascension of the center of the first
* IPN annulus (J2000) in the heliocentric frame
real*8 DE1deg ! (deg) Declination of the center of the first IPN
* annulus (J2000) in the heliocentric frame
real*4 Rad1 ! (deg) Angular radius of the first IPN annulus, in
* the heliocentric frame
real*4 DelR1 ! (deg) Half-width of the first IPN annulus
* ({delta}R_IPN1_) (3)
real*8 RA2deg ! (deg) ? Right ascension of the center of the second
* IPN annulus (J2000) in the heliocentric frame
real*8 DE2deg ! (deg) ? Declination of the center of the second IPN
* annulus (J2000) in the heliocentric frame
real*4 Rad2 ! (deg) ? Angular radius of the second IPN annulus, in
* the heliocentric frame
real*4 DelR2 ! (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 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Burst_1 ! Date of the burst (YYMMDD)
character*1 f_Burst_1 ! [o] o: alternative error box for trigger 7647
integer*4 Trig_1 ! [6380,8116] BATSE trigger number
integer*4 RAh ! (h) Hour of right ascension (J2000)
integer*4 RAm ! (min) Minute of right ascension (J2000)
real*4 RAs ! (s) Second of right ascension (J2000)
character*1 DE_ ! Sign of declination (J2000)
integer*4 DEd ! (deg) Degree of declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of declination (J2000)
integer*4 loc ! [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,m_Burst,f_Burst,Time,Trig,RAdeg,DEdeg,Sig,RA1deg,DE1deg,
+ Rad1,DelR1,RA2deg,DE2deg,Rad2,DelR2
if(ar__(76:82) .EQ. '') RA2deg = rNULL__
if(ar__(84:90) .EQ. '') DE2deg = rNULL__
if(ar__(92:97) .EQ. '') Rad2 = rNULL__
if(ar__(99:104) .EQ. '') DelR2 = rNULL__
c ..............Just test output...........
write(6,1)
+ Burst,m_Burst,f_Burst,Time,Trig,RAdeg,DEdeg,Sig,RA1deg,DE1deg,
+ Rad1,DelR1,RA2deg,DE2deg,Rad2,DelR2
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,f_Burst_1,Trig_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,loc
RAdeg_1 = rNULL__
DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh*15.
if(RAm .GT. -180) RAdeg_1=RAdeg_1+RAm/4.
if(RAs .GT. -180) RAdeg_1=RAdeg_1+RAs/240.
if(DEd .GE. 0) DEdeg_1=DEd
if(DEm .GE. 0) DEdeg_1=DEdeg_1+DEm/60.
if(DEs .GE. 0) DEdeg_1=DEdeg_1+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,2)
+ Burst_1,f_Burst_1,Trig_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,loc
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end