Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/555/A65 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/A+A/555/A65 X-ray variability of M31 central field (Hofmann+, 2013)
*================================================================================
*X-ray source variability study of the M 31 central field using Chandra HRC-I.
* Hofmann F., Pietsch W., Henze M., Haberl F., Strum R., Della Valle M.,
* Hartmann D.H., Hatzidimitriou D.
* <Astron. Astrophys. 555, A65 (2013)>
* =2013A&A...555A..65H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Source catalogue with general source information
integer*4 nr__
parameter (nr__=318) ! Number of records
character*396 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 SRCID ! [1/318] Source identification number
real*8 RAdeg ! (deg) Right ascension (J2000)
real*8 DEdeg ! (deg) Declination (J2000)
real*8 ePos ! (arcsec) Positional error (1)
real*4 ctN ! (ct) Net counts
real*8 e_ctN ! (ct) Error of net counts
real*4 ctBkg ! (ct) Net background counts
real*8 CRN ! (ct/s) Net count rate
real*4 e_CRN ! (ct/s) Error of net count rate
real*8 CRBkg ! (ct/s) Background count rate
real*4 e_CRBkg ! (ct/s) Error of background count rate
real*4 ExpTime ! (s) Exposure time
real*4 Signi ! Highest detection significance
integer*4 Res ! Detection-image resolution (binning) (2)
real*4 Flux ! (mW/m2) Flux in 0.2-10keV band, in erg/s/cm^2^ (3)
real*4 e_Flux ! (mW/m2) Flux error
real*4 Lum ! (10-7W) Luminosity, assuming a distance of 780kpc
real*4 e_Lum ! (10-7W) Luminosity error
real*4 DistNN ! (arcsec) Distance to nearest neighbour
real*4 svar ! Source variability parameter
* (Fmax-Fmin)/sqrt({sigma}max^2^+{sigma}min^2^)
real*4 fvar ! Highest source flux ratio (Fmax/Fmin)
real*8 ndet ! [0/1] Number of detections (8)
real*8 ssvar ! Highest svar between observations
real*4 Fmax ! (mW/m2) Maximum flux
real*4 e_Fmax ! (mW/m2) Error of maximum flux
real*4 Fmin ! (mW/m2) Minimum flux
real*4 e_Fmin ! (mW/m2) Error of minimum flux
character*2 lcC ! [HV OB] Light curve classification:
* HV (highly variable) or OB (outburst)
character*28 nova ! Nova name(s) (4)
character*12 glc ! ? Globular-cluster name (5)
character*8 Class ! Source classification (6)
character*58 Com ! Comments and references (7)
*Note (1): 3 x (1-{sigma} statistical error plus systematic errors)
*Note (2): Detection-image resolution as:
* 1 = 1x1 binning
* 2 = 2x2 binning
* 3 = 4x4 binning
* 4 = 8x8 binning
* 5 = 16x16 binning
*Note (3): 0.2-10keV flux calculated for a power-law spectrum with a photon
* index of 1.7 and Galactic foreground absorption (column density of
* N_H_=6.6x10^20^cm^-2^).
*Note (4): Nova identifier (M31N YYYY-MMa) gives year (YYYY) and month (MM)
* of outburst plus a letter to distinguish multiple novae per month.
*Note (5): Globular cluster name followed by quality flag (see Revised Bologna
* Catalog of Globular Clusters: Galleti et al., 2007, Cat. J/A+A/471/127)
*Note (6): Classes (candidates are marked with <>) are:
* XRB = X-ray binary
* GlC = globular cluster
* SSS = super-soft source
* nova = nova
* SNR = supernova remnant
* AGN = active galactic nucleus
* fgStar = foreground star
*Note (7): Abbreviations and Codes used:
* TPC2006 = Trudolyubov et al., 2006ApJ...645..277T
* HHP2009 = Haberl et al., 2009ATel.2262....1H
* GGS2007 = Galache et al., 2007ATel.1171....1G
* HPG2006 = Haberl et al., 2006ATel..881....1H
* BGG2012 = Barnard et al., 2012ApJ...756...32B
* B90 = Braun, 1990ApJS...72..761B
* HPH2010 = Henze et al., 2010A&A...523A..89H
* GGT2009 = Galache et al., 2009ATel.1978....1G
* KGG2009 = Kong et al., 2009ATel.2074....1K
* HHP2012 = Kong et al., 2012ATel.4164....1H
* HPH2012a = Henze et al., 2012ATel.3890....1H
* HPH2012b = Henze et al., 2012ATel.3921....1H
* HPH2012c = Henze et al., 2013A&A...549A.120H
* GGT2008 = Galache et al., 2008ATel.1698....1G
* TP2008 = Trudolyubov & Priedhorsky, 2008ApJ...676.1218T
* HPH2009 = Henze et al., 2009ATel.2356....1H
* KGM2001 = Kong et al., 2001ATel...76....1K
* PBG2008 = Pietsch et al., 2008ATel.1674....1P
* HHP2007 = Harbel et al., 2007ATel.1296....1H
* TBP2001 = Trudolyubov et al., 2001ApJ...563L.119T
* TBP2002 = Trudolyubov et al., 2002ApJ...581L..27T
* PFH2008 = Pietsch et al., 2008ATel.1647....1P
* HPH2009 = Henze et al., 2009A&A...500..769H
* WGM2005 = Williams et al., 2005ApJ...632.1086W
* BK2009 = Barnard & Kolb, 2009MNRAS.397L..92B
* BKO2004 = Barnard et al., 2004A&A...423..147B
* PHH2012 = Pietsch et al., 2010ATel.3038....1P
* no lc = no light curve available
* AP = activity period
* SOB = short outburst
* rec = recurrent
* XT[SRCID] = cross-talk with SRCID
* LOB = long outburst
* ULX = ultra-luminous X-ray source
* M31* = nucleus of M 31
* radio B90 = correlation with radio source in B90 catalogue
*Note (8): ndet is the fraction of the number of data points in the light
* curve having a count>1-{sigma}.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'xcorr.dat' ! X-ray correlations
integer*4 nr__1
parameter (nr__1=1735) ! Number of records
character*101 ar__1 ! Full-size record
integer*4 SRCID_1 ! Source identification number in the catalogue
character*9 Cat ! X-ray correlations Catalogue (1)
character*85 XID ! Identification and description in "Cat" (1)
*Note (1): The correlation contains the Catalogue "Cat", identification and
* description in the Catalogue. References and column name in reference
* catalogues are:
* KGP2002 = Kong et al., 2002 Cat. J/ApJ/577/38, Notes
* Ka2002 = Kaaret, 2002, Cat. J/ApJ/578/114, CXOM31
* PFH2005 = Pietsch et al., 2005, Cat. J/A+A/434/483, Class
* SHL2001 = Supper et al., 2001, Cat. J/A+A/373/63, n_SII
* SHP97 = Supper et al., 1997, Cat. J/A+A/317/328, n_ROSAT
* SPH2008 = Stiele et al., 2008, Cat. J/A+A/480/599, Class
* SPH2011 = Stiele et al., 2011, Cat. J/A+A/534/A55, Class
* VG2007 = Voss & Gilfanov, 2007, Cat. J/A+A/468/49, Type
* WGK2004 = Williams et al., 2004, Cat. J/ApJ/609/735, Type
* WNG2006 = Williams et al., 2006, Cat. J/ApJ/643/356, F
* HPH2011 = Henze et al., 2011, Cat. J/A+A/533/A52, M31N
* PHS2007 = Pietsch et al., 2007, Cat. J/A+A/465/375, M31N
* CSC = Evans et al., 2010ApJS..189...37E (Cat. IX/45), msid
* D2002 = Di Stefano et al., 2002ApJ...570..618D, opt_ident
* DKG2004 = Di Stefano et al., 2004, Cat. J/ApJ/610/247, category
* O2006 = Orio, 2006ApJ...643..844O, type
* OBT2001 = Osborne et al., 2001A&A...378..800O, Comments
* SBK2009 = Shaw Greening et al., 2009, Cat. J/A&A...495..733S, type
* TF91 = Trinchieri & Fabbiano, 1991ApJ...382...82T, INDEX
* TP2004 = Trudolyubov & Priedhorsky, 2004ApJ...616..821T, optical_id
* TP2006 = Trudolyubov et al., 2006ApJ...645..277T, class
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Light curve table containing flux and errors
integer*4 nr__2
parameter (nr__2=19716) ! Number of records
character*62 ar__2 ! Full-size record
integer*4 SRCID_2 ! [1/318] Sequential number
integer*4 ObsID ! [267/13281] HRC-I observation identifier
character*10 Date ! ("YYYY-MM-DD") Date of the observation
real*8 MJD ! (d) Modified Julian Date
character*1 l_Flux ! [<] when 3-{sigma} upper limit in Flux for
* sources with detection-significance <1{sigma}
real*4 Flux_1 ! (mW/m2) ? Flux (1)
real*4 B_Flux ! (mW/m2) ? Upper value of flux with 1{sigma} error (1)
real*4 b_Flux_1 ! (mW/m2) ? Lower value of flux with 1{sigma} error (1)
*Note (1): When no values, the source was outside the field of view of the
* observation.
C=============================================================================
C Loading file 'table1.dat' ! Source catalogue with general source information
C Format for file interpretation
1 format(
+ I3,1X,F9.5,1X,F9.5,1X,F9.5,1X,E12.6,1X,F9.5,1X,E12.6,1X,F8.6,
+ 1X,E12.6,1X,F8.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,I1,1X,E9.2,1X,
+ E9.2,1X,E9.2,1X,E9.2,1X,E9.2,1X,E9.2,1X,E9.2,1X,F8.6,1X,F10.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,A2,1X,A28,1X,A12,1X,A8,
+ 1X,A58)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,318
read(1,'(A396)')ar__
read(ar__,1)
+ SRCID,RAdeg,DEdeg,ePos,ctN,e_ctN,ctBkg,CRN,e_CRN,CRBkg,
+ e_CRBkg,ExpTime,Signi,Res,Flux,e_Flux,Lum,e_Lum,DistNN,svar,
+ fvar,ndet,ssvar,Fmax,e_Fmax,Fmin,e_Fmin,lcC,nova,glc,Class,Com
c ..............Just test output...........
write(6,1)
+ SRCID,RAdeg,DEdeg,ePos,ctN,e_ctN,ctBkg,CRN,e_CRN,CRBkg,
+ e_CRBkg,ExpTime,Signi,Res,Flux,e_Flux,Lum,e_Lum,DistNN,svar,
+ fvar,ndet,ssvar,Fmax,e_Fmax,Fmin,e_Fmin,lcC,nova,glc,Class,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'xcorr.dat' ! X-ray correlations
C Format for file interpretation
2 format(I3,1X,A9,3X,A85)
C Effective file loading
open(unit=1,status='old',file=
+'xcorr.dat')
write(6,*) '....Loading file: xcorr.dat'
do i__=1,1735
read(1,'(A101)')ar__1
read(ar__1,2)SRCID_1,Cat,XID
c ..............Just test output...........
write(6,2)SRCID_1,Cat,XID
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Light curve table containing flux and errors
C Format for file interpretation
3 format(I3,1X,I5,1X,A10,1X,F11.5,1X,A1,E8.3,2X,E8.3,2X,E8.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,19716
read(1,'(A62)')ar__2
read(ar__2,3)
+ SRCID_2,ObsID,Date,MJD,l_Flux,Flux_1,B_Flux,b_Flux_1
if(ar__2(35:42) .EQ. '') Flux_1 = rNULL__
if(ar__2(45:52) .EQ. '') B_Flux = rNULL__
if(ar__2(55:62) .EQ. '') b_Flux_1 = rNULL__
c ..............Just test output...........
write(6,3)
+ SRCID_2,ObsID,Date,MJD,l_Flux,Flux_1,B_Flux,b_Flux_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end