Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/179/360 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 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/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
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
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,status='old',file=
+'sources.dat')
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
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,SName
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
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,status='old',file=
+'table6.dat')
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,status='old',file=
+'refs.dat')
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