Conversion of standardized ReadMe file for
file /./ftp/cats/J/other/GCN/1070.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 2024-Apr-19
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/other/GCN/1070.1 HETE Trigger Information (Ricker+, 2001)
*================================================================================
*SGR1806-20: A Second, Intense Burst Localized by HETE.
* Ricker G., Lamb D., Woosley S., Vanderspek R., Crew G., Doty J.,
* Monnelly G., Villasenor J., Butler N., Cline T., Jernigan J.G., Levine A.,
* Martel F., Morgan E., Pizzichini G., Prigozhin G., Kawai N., Matsuoka M.,
* Shirasaki Y., Tamagawa T., Torii K., Sakamoto T., Yoshida A., Fenimore E.,
* Galassi M., Graziani C., Atteia J.-L., Boer M., Olive J.-F., Dezalay J.-P.,
* Hurley K.
* <GRB Circular Network, 1070, 1 (2001)>
* =2001GCN..1070....1R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hete.dat' ! HETE Trigger Information about GRBs
(and non-GRBs) detected by HETE (WXM and SXC)
integer*4 nr__
parameter (nr__=145) ! Number of records
character*703 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Trigger (nr__) ! ? Trigger number (in Simbad:
* HETE Trigger NNNNN)
character*1 n_Trigger (nr__) ! [n] 'n' for not available
character*8 Obs_Date (nr__) ! ("YY/MM/DD") Date of the start of the GRB
character*11 Obs_Time (nr__) ! ? Date the start of the GRB (hh:mm:ss.ss)
character*10 Notice (nr__) ! Type of notice (1)
character*4 Inst (nr__) ! Instrument that made the
* observation/detection (2)
real*8 RAdeg (nr__) ! (deg) ? Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) ? Declination in decimal degrees (J2000)
real*4 Error (nr__) ! (arcmin) ? Error uncertainty in that location
* (arcmin diameter, systematic+statistical)
character*1 Note (nr__) ! [on] Note on the absence of position (3)
character*632 Com (nr__) ! Comments describing the observation results
*Note (1): Type of notice: S/C_Alert, S/C_Update, S/C_Last, or GndAnalysis
*Note (2): Instrument: WXM, SXC, TBD if not known yet as in an S/C_Alert Notice
*Note (3): Note when no position as follows:
* n = not available
* o = none, if the s/c decided there was enough to trigger,
* but not enough events to produce a position
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hete2001.dat' ! Old HETE Trigger Information: 2001
integer*4 nr__1
parameter (nr__1=374) ! Number of records
character*703 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Trigger_1 (nr__1) ! ? Trigger number (in Simbad:
* HETE Trigger NNNNN)
character*1 n_Trigger_1(nr__1) ! [n] 'n' for not available
character*8 Obs_Date_1 (nr__1) ! ("YY/MM/DD") Date of the start of the GRB
character*11 Obs_Time_1 (nr__1) ! ? Date the start of the GRB (hh:mm:ss.ss)
character*10 Notice_1 (nr__1) ! Type of notice (1)
character*4 Inst_1 (nr__1) ! Instrument that made the
* observation/detection (2)
real*8 RAdeg_1 (nr__1) ! (deg) ? Right ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) ? Declination in decimal degrees (J2000)
real*4 Error_1 (nr__1) ! (arcmin) ? Error uncertainty in that location
* (arcmin diameter, systematic+statistical)
character*1 Note_1 (nr__1) ! [on] Note on the absence of position (3)
character*632 Com_1 (nr__1) ! Comments describing the observation results
*Note (1): Type of notice: S/C_Alert, S/C_Update, S/C_Last, or GndAnalysis
*Note (2): Instrument: WXM, SXC, TBD if not known yet as in an S/C_Alert Notice
*Note (3): Note when no position as follows:
* n = not available
* o = none, if the s/c decided there was enough to trigger,
* but not enough events to produce a position
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hete2002.dat' ! Old HETE Trigger Information: 2002
integer*4 nr__2
parameter (nr__2=419) ! Number of records
character*703 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Trigger_2 (nr__2) ! ? Trigger number (in Simbad:
* HETE Trigger NNNNN)
character*1 n_Trigger_2(nr__2) ! [n] 'n' for not available
character*8 Obs_Date_2 (nr__2) ! ("YY/MM/DD") Date of the start of the GRB
character*11 Obs_Time_2 (nr__2) ! ? Date the start of the GRB (hh:mm:ss.ss)
character*10 Notice_2 (nr__2) ! Type of notice (1)
character*4 Inst_2 (nr__2) ! Instrument that made the
* observation/detection (2)
real*8 RAdeg_2 (nr__2) ! (deg) ? Right ascension in decimal degrees (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) ? Declination in decimal degrees (J2000)
real*4 Error_2 (nr__2) ! (arcmin) ? Error uncertainty in that location
* (arcmin diameter, systematic+statistical)
character*1 Note_2 (nr__2) ! [on] Note on the absence of position (3)
character*632 Com_2 (nr__2) ! Comments describing the observation results
*Note (1): Type of notice: S/C_Alert, S/C_Update, S/C_Last, or GndAnalysis
*Note (2): Instrument: WXM, SXC, TBD if not known yet as in an S/C_Alert Notice
*Note (3): Note when no position as follows:
* n = not available
* o = none, if the s/c decided there was enough to trigger,
* but not enough events to produce a position
C=============================================================================
C Loading file 'hete.dat' ! HETE Trigger Information about GRBs
* (and non-GRBs) detected by HETE (WXM and SXC)
C Format for file interpretation
1 format(
+ I5,A1,5X,A8,1X,A11,1X,A10,1X,A4,1X,F7.3,F7.3,1X,F6.2,A1,1X,
+ A632)
C Effective file loading
open(unit=1,status='old',file=
+'hete.dat')
write(6,*) '....Loading file: hete.dat'
do i__=1,145
read(1,'(A703)')ar__
read(ar__,1)
+ Trigger(i__),n_Trigger(i__),Obs_Date(i__),Obs_Time(i__),
+ Notice(i__),Inst(i__),RAdeg(i__),DEdeg(i__),Error(i__),
+ Note(i__),Com(i__)
if(ar__(1:5) .EQ. '') Trigger(i__) = iNULL__
if(ar__(49:55) .EQ. '') RAdeg(i__) = rNULL__
if(ar__(56:62) .EQ. '') DEdeg(i__) = rNULL__
if(ar__(64:69) .EQ. '') Error(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Trigger(i__),n_Trigger(i__),Obs_Date(i__),Obs_Time(i__),
+ Notice(i__),Inst(i__),RAdeg(i__),DEdeg(i__),Error(i__),
+ Note(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hete2001.dat' ! Old HETE Trigger Information: 2001
C Format for file interpretation
2 format(
+ I5,A1,5X,A8,1X,A11,1X,A10,1X,A4,1X,F7.3,F7.3,1X,F6.2,A1,1X,
+ A632)
C Effective file loading
open(unit=1,status='old',file=
+'hete2001.dat')
write(6,*) '....Loading file: hete2001.dat'
do i__=1,374
read(1,'(A703)')ar__1
read(ar__1,2)
+ Trigger_1(i__),n_Trigger_1(i__),Obs_Date_1(i__),
+ Obs_Time_1(i__),Notice_1(i__),Inst_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),Error_1(i__),Note_1(i__),Com_1(i__)
if(ar__1(1:5) .EQ. '') Trigger_1(i__) = iNULL__
if(ar__1(49:55) .EQ. '') RAdeg_1(i__) = rNULL__
if(ar__1(56:62) .EQ. '') DEdeg_1(i__) = rNULL__
if(ar__1(64:69) .EQ. '') Error_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Trigger_1(i__),n_Trigger_1(i__),Obs_Date_1(i__),
+ Obs_Time_1(i__),Notice_1(i__),Inst_1(i__),RAdeg_1(i__),
+ DEdeg_1(i__),Error_1(i__),Note_1(i__),Com_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hete2002.dat' ! Old HETE Trigger Information: 2002
C Format for file interpretation
3 format(
+ I5,A1,5X,A8,1X,A11,1X,A10,1X,A4,1X,F7.3,F7.3,1X,F6.2,A1,1X,
+ A632)
C Effective file loading
open(unit=1,status='old',file=
+'hete2002.dat')
write(6,*) '....Loading file: hete2002.dat'
do i__=1,419
read(1,'(A703)')ar__2
read(ar__2,3)
+ Trigger_2(i__),n_Trigger_2(i__),Obs_Date_2(i__),
+ Obs_Time_2(i__),Notice_2(i__),Inst_2(i__),RAdeg_2(i__),
+ DEdeg_2(i__),Error_2(i__),Note_2(i__),Com_2(i__)
if(ar__2(1:5) .EQ. '') Trigger_2(i__) = iNULL__
if(ar__2(49:55) .EQ. '') RAdeg_2(i__) = rNULL__
if(ar__2(56:62) .EQ. '') DEdeg_2(i__) = rNULL__
if(ar__2(64:69) .EQ. '') Error_2(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Trigger_2(i__),n_Trigger_2(i__),Obs_Date_2(i__),
+ Obs_Time_2(i__),Notice_2(i__),Inst_2(i__),RAdeg_2(i__),
+ DEdeg_2(i__),Error_2(i__),Note_2(i__),Com_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end