Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-22
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJS/179/360 Thermonuclear X-ray bursts observed by RXTE (Galloway+, 2008)
*================================================================================
*Thermonuclear (type-I) X-ray bursts observed by the Rossi X-ray Timing Explorer.
* Galloway D.K., Muno M.P., Hartman J.M., Psaltis D., Chakrabarty D.
* <Astrophys. J. Suppl. Ser. 179, 360 (2008)>
* =2008ApJS..179..360G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sources.dat' ! Sources list
integer*4 nr__
parameter (nr__=48) ! Number of records
character*60 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension J2000
real*8 Dec ! (deg) Declination J2000
character*17 Name ! Object identification
integer*4 RAh ! (h) Simbad right ascension (J2000)
integer*4 RAm ! (min) Simbad right ascension (J2000)
real*4 RAs ! (s) Simbad right ascension (J2000)
character*1 DE_ ! Simbad declination sign (J2000)
integer*4 DEd ! (deg) Simbad declination (J2000)
integer*4 DEm ! (arcmin) Simbad declination (J2000)
integer*4 DEs ! (arcsec) Simbad declination (J2000)
character*18 SName ! Simbad name
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Table of burst parameters
integer*4 nr__1
parameter (nr__1=1187) ! Number of records
character*301 ar__1 ! Full-size record
character*17 Name_1 ! Object identification
integer*4 ID ! Burst number
character*5 f_ID ! [a-h, ] Note on the burst (1)
character*15 OID ! RXTE observational identification,
* (NNNNN-NN-NN-NNA)
character*10 Obs_date ! ("YYYY/MM/DD") UT year of start time (2)
integer*4 Obs_h ! (h) UT hour of start time (2)
integer*4 Obs_m ! (min) UT minute of start time (2)
integer*4 Obs_s ! (s) UT second of start time (2)
real*8 MJD ! (d) Modified Julian Date of start time (2)
character*1 l_Flux ! Limit flag on Flux
real*8 Flux ! (pW/m2) ? Peak flux (3)
real*4 e_Flux ! (pW/m2) ? Uncertainty in Flux (3)
character*1 l_Fluence ! Limit flag on Fluence
real*8 Fluence ! (nJ/m2) ? Burst fluence (5)
real*4 e_Fluence ! (nJ/m2) ? Uncertainty in Fluence (5)
character*2 Exp ! Presence of radius expansion? (6)
character*1 l_Rtime ! Limit flag on Rtime
real*4 Rtime ! (s) ? Rise time
real*4 e_Rtime ! (s) ? Uncertainty in Rtime
integer*4 Prate ! (ct/s) ? Peak count rate per PCU (proportional
* counter unit), including persistent
* emission and background
real*8 tau1 ! (s) ? Exponential decay constant {tau}_1_
real*4 tau2 ! (s) ? Exponential decay constant {tau}_2_
real*4 Scale ! (s) ? Burst time scale
real*4 e_Scale ! (s) ? Uncertainty in Scale
character*1 l_Nfluence ! Limit flag on Nfluence
real*8 Nfluence ! ? Normalized fluence (7)
real*4 e_Nfluence ! ? Uncertainty in Nfluence
real*8 Pflux ! (pW/m2) ? Persistent 2.5-25 keV flux
* level prior to burst (3)
real*4 e_Pflux ! (pW/m2) ? Uncertainty in Pflux(3)
real*8 NPflux ! (pW/m2) ? Persistent 2.5-25 keV flux level
* prior to burst normalized to F_Edd_ (3)
real*8 e_NPflux ! (pW/m2) ? Uncertainty in NPflux (3)
real*4 Scol ! ? Soft color prior to burst (4)
real*4 e_Scol ! ? Uncertainty in Scolor
real*4 Hcol ! ? Hard color prior to burst (4)
real*4 e_Hcol ! ? Uncertainty in Hcolor
character*1 l_SZ ! Limit flag on SZ
real*8 SZ ! ? Position on color-color diagram
real*4 e_SZ ! ? Uncertainty in SZ
real*4 Delt ! (h) ? Inferred recurrence time
real*4 cbol ! ? Correction used to estimate bolometric
* flux
real*4 e_cbol ! ? Uncertainty in cbol
real*8 alpha ! ? Ratio of integrated persistent flux to
* burst fluence
real*8 e_alpha ! ? Uncertainty in alpha
character*17 Refs ! Previously published analyses,
* in refs.dat file
*Note (1): Flag as follows:
* a = The burst was observed during a slew, and thus offset from the
* source position
* b = The observation was offset from the source position. In cases (a)
* and (b) we scaled the flux and fluence by the mean collimator
* response appropriate for the position of the source in the field
* of view, as described in appendix B
* c = The origin of the burst is uncertain; the burst may have been
* from another source in the field-of-view (we rescaled the flux
* and fluence, if necessary, based upon the assumed origin)
* d = Buffer overruns (or some other instrumental effect) caused gaps
* in the high time-resolution data
* e = The burst was so faint that only the peak flux could be measured,
* and not the fluence or other parameters
* f = An extremely faint burst or possibly problems with the background
* subtraction, resulting in no fit results
* g = The full burst profile was not observed, so that the event can be
* considered an unconfirmed burst candidate. Typically in these
* cases the initial burst rise is missed, so that the measured
* peak flux and fluence are lower limits only
* h = High-time resolution datamodes did not cover the burst.
*Note (2): We neglect corrections to give the time at the solar-system
* barycenter.
*Note (3): In units of 10^-9^erg/s/cm^-2^.
*Note (4): The hard and soft X-ray colors are the ratios of the
* background-subtracted PCA counts in the (8.6-18.0)/(5.0-8.6)keV and
* (3.6-5.0)/(2.2-3.6)keV energy bands, respectively, calculated from 64s
* integrations when the source intensity was above 100ct/s, and 256s
* integrations otherwise, and corrected for gain variations over the
* mission lifetime by comparison to the colors of the Crab nebula
*Note (5): In units of 10^-6^erg/cm^2^.
*Note (6): Presence of radius expansion flag as follows:
* Y = yes
* N = no
* Y? = possible
* blank = unknown
*Note (7): Normalized by the mean peak flux of the PRE bursts F_Edd_,
* where available (U_b_ in van Paradijs et al., 1988MNRAS.233..437V).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__2
parameter (nr__2=51) ! Number of records
character*49 ar__2 ! Full-size record
integer*4 Ref ! Reference number
character*19 BibCode ! BibCode
character*26 Aut ! Author's name
C=============================================================================
C Loading file 'sources.dat' ! Sources list
C Format for file interpretation
1 format(A17,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,3X,A18)
C Effective file loading
open(unit=1,file='sources.dat', status='old')
write(6,*) '....Loading file: sources.dat'
do i__=1,48
read(1,'(A60)')ar__
read(ar__,1)Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,SName
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAs .GE. 0) RA_=RA_+RAs/240.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DEs .GE. 0) Dec=Dec+DEs/3600.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,SName
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Table of burst parameters
C Format for file interpretation
2 format(
+ A17,1X,I3,1X,A5,1X,A15,1X,A10,1X,I2,1X,I2,1X,I2,1X,F11.5,1X,
+ A1,F8.3,1X,F5.2,1X,A1,F7.4,1X,F6.4,1X,A2,1X,A1,F6.2,1X,F4.2,
+ 1X,I5,1X,F7.3,1X,F6.2,1X,F6.2,1X,F6.2,1X,A1,F7.3,1X,F6.3,1X,
+ F7.4,1X,F6.4,1X,F7.5,1X,F7.5,1X,F5.3,1X,F5.3,1X,F6.4,1X,F6.4,
+ 1X,A1,F7.4,1X,F6.4,1X,F6.3,1X,F5.3,1X,F5.3,1X,F8.3,1X,F7.3,1X,
+ A17)
C Effective file loading
open(unit=1,file='table6.dat', status='old')
write(6,*) '....Loading file: table6.dat'
do i__=1,1187
read(1,'(A301)')ar__1
read(ar__1,2)
+ Name_1,ID,f_ID,OID,Obs_date,Obs_h,Obs_m,Obs_s,MJD,l_Flux,Flux,
+ e_Flux,l_Fluence,Fluence,e_Fluence,Exp,l_Rtime,Rtime,e_Rtime,
+ Prate,tau1,tau2,Scale,e_Scale,l_Nfluence,Nfluence,e_Nfluence,
+ Pflux,e_Pflux,NPflux,e_NPflux,Scol,e_Scol,Hcol,e_Hcol,l_SZ,SZ,
+ e_SZ,Delt,cbol,e_cbol,alpha,e_alpha,Refs
if(ar__1(78:85) .EQ. '') Flux = rNULL__
if(ar__1(87:91) .EQ. '') e_Flux = rNULL__
if(ar__1(94:100) .EQ. '') Fluence = rNULL__
if(ar__1(102:107) .EQ. '') e_Fluence = rNULL__
if(ar__1(113:118) .EQ. '') Rtime = rNULL__
if(ar__1(120:123) .EQ. '') e_Rtime = rNULL__
if(ar__1(125:129) .EQ. '') Prate = iNULL__
if(ar__1(131:137) .EQ. '') tau1 = rNULL__
if(ar__1(139:144) .EQ. '') tau2 = rNULL__
if(ar__1(146:151) .EQ. '') Scale = rNULL__
if(ar__1(153:158) .EQ. '') e_Scale = rNULL__
if(ar__1(161:167) .EQ. '') Nfluence = rNULL__
if(ar__1(169:174) .EQ. '') e_Nfluence = rNULL__
if(ar__1(176:182) .EQ. '') Pflux = rNULL__
if(ar__1(184:189) .EQ. '') e_Pflux = rNULL__
if(ar__1(191:197) .EQ. '') NPflux = rNULL__
if(ar__1(199:205) .EQ. '') e_NPflux = rNULL__
if(ar__1(207:211) .EQ. '') Scol = rNULL__
if(ar__1(213:217) .EQ. '') e_Scol = rNULL__
if(ar__1(219:224) .EQ. '') Hcol = rNULL__
if(ar__1(226:231) .EQ. '') e_Hcol = rNULL__
if(ar__1(234:240) .EQ. '') SZ = rNULL__
if(ar__1(242:247) .EQ. '') e_SZ = rNULL__
if(ar__1(249:254) .EQ. '') Delt = rNULL__
if(ar__1(256:260) .EQ. '') cbol = rNULL__
if(ar__1(262:266) .EQ. '') e_cbol = rNULL__
if(ar__1(268:275) .EQ. '') alpha = rNULL__
if(ar__1(277:283) .EQ. '') e_alpha = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,ID,f_ID,OID,Obs_date,Obs_h,Obs_m,Obs_s,MJD,l_Flux,Flux,
+ e_Flux,l_Fluence,Fluence,e_Fluence,Exp,l_Rtime,Rtime,e_Rtime,
+ Prate,tau1,tau2,Scale,e_Scale,l_Nfluence,Nfluence,e_Nfluence,
+ Pflux,e_Pflux,NPflux,e_NPflux,Scol,e_Scol,Hcol,e_Hcol,l_SZ,SZ,
+ e_SZ,Delt,cbol,e_cbol,alpha,e_alpha,Refs
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
3 format(I2,1X,A19,1X,A26)
C Effective file loading
open(unit=1,file='refs.dat', status='old')
write(6,*) '....Loading file: refs.dat'
do i__=1,51
read(1,'(A49)')ar__2
read(ar__2,3)Ref,BibCode,Aut
c ..............Just test output...........
write(6,3)Ref,BibCode,Aut
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end