Conversion of standardized ReadMe file for
file /./ftp/cats/IX/56 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-Aug-12
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. IX/56 3XMM-DR7s serendipitous source catalogue from stacks (Traulsen+, 2019)
*================================================================================
*The XMM-Newton serendipitous survey.
*VIII. The first XMM-Newton serendipitous source catalogue from overlapping
*observations.
* Traulsen I., Schwope A.D., Lamer G., Ballet J., Carrera F., Coriat M.,
* Freyberg M.J., Michel L., Motch C., Rosen S.R., Webb N., Ceballos M.T.,
* Koliopanos F., Kurpas J., Page M., Watson M.G.
* <Astron. Astrophys. 624, A77 (2019)>
* =2019A&A...624A..77T
* =2019yCat.9056....0T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '3xmmdr7s.dat' ! Source catalogue
integer*4 nr__
parameter (nr__=288344) ! Number of records
character*3243 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*8 Source ! [3000011010100001/3079138010100577]
* Unique source index (SRCID)
character*10 ObsID ! XMM-Newton observation identifier
* (OBS_ID)
integer*4 Nobs ! ? Number of observations involved in
* the stack (N_OBS)
integer*4 Ncont ! ? Number of observations in which the
* source was fitted (N_CONTRIB)
real*8 RAdeg ! (deg) Right ascension (J2000) (RA)
real*8 DEdeg ! (deg) Declination (J2000) (DEC)
real*8 srcML ! (arcsec) Error on position (RADEC_ERR)
real*8 GLON ! (deg) Galactic longitude (LII)
real*8 GLAT ! (deg) Galactic latitude (BII)
real*8 Ximg ! (pix) X coordinate within the re-binned image
* (X_IMA)
real*8 e_Ximg ! (pix) Error on Ximg (X_IMA_ERR)
real*8 Yimg ! (pix) Y coordinate within the re-binned image
* (Y_IMA)
real*8 e_Yimg ! (pix) Error on Yimg (Y_IMA_ERR)
real*8 DistNN ! (arcsec) Distance to the nearest neighbouring
* detection (DIST_NN)
character*17 v3XMM ! Source name, JHHMMSS.s+DDMMSSs, ending on
* an "s" for stacked detection (SRCNAME)
character*21 v3XMMDR7 ! IAU name assigned to the nearest unique
* source in 3XMM-DR7 (IAUNAME_3XMMDR7)
integer*8 SrcidDR7 ! ? Source identifier of DR7 match
* (SRCID_3XMMDR7)
integer*8 DetidDR7 ! ? Detection identifier of DR7 match
* (DETID_3XMMDR7)
real*8 RA7deg ! (deg) ? (SC_)RA of DR7 match (RA_3XMMDR7)
real*8 DE7deg ! (deg) ? (SC_)DEC of DR7 match (DEC_3XMMDR7)
real*8 ePos7 ! (arcsec) ? Position error of DR7 match
* (POSERR_3XMMDR7)
real*8 DistDR7 ! (arcsec) ? Distance to DR7 match (DIST_3XMMDR7)
integer*4 NdetDR7 ! ? Number of DR7 detections
* (NDETECT_3XMMDR7)
real*4 FEP ! (mW/m2) ? All-EPIC flux (EP_FLUX)
real*4 e_FEP ! (mW/m2) ? 1-sigma error on EP_FLUX (EP_FLUX_ERR)
real*4 FEP1 ! (mW/m2) Total flux in energy band 1 (EP_1_FLUX)
real*4 e_FEP1 ! (mW/m2) 1-sigma error on EP_1_FLUX
* (EP_1_FLUX_ERR)
real*4 FEP2 ! (mW/m2) Total flux in energy band 2 (EP_2_FLUX)
real*4 e_FEP2 ! (mW/m2) 1-sigma error on EP_2_FLUX
* (EP_2_FLUX_ERR)
real*4 FEP3 ! (mW/m2) Total flux in energy band 3 (EP_3_FLUX)
real*4 e_FEP3 ! (mW/m2) 1-sigma error on EP_3_FLUX
* (EP_3_FLUX_ERR)
real*4 FEP4 ! (mW/m2) Total flux in energy band 4 (EP_4_FLUX)
real*4 e_FEP4 ! (mW/m2) 1-sigma error on EP_4_FLUX
* (EP_4_FLUX_ERR)
real*4 FEP5 ! (mW/m2) Total flux in energy band 5 (EP_5_FLUX)
real*4 e_FEP5 ! (mW/m2) 1-sigma error on EP_5_FLUX
* (EP_5_FLUX_ERR)
real*4 FPN ! (mW/m2) ? Total EPIC/pn flux (PN_FLUX)
real*4 e_FPN ! (mW/m2) ? 1-sigma error on PN_FLUX (PN_FLUX_ERR)
real*4 FPN1 ! (mW/m2) ? EPIC/pn flux in energy band 1
* (PN_1_FLUX)
real*4 e_FPN1 ! (mW/m2) ? 1-sigma error on PN_1_FLUX
* (PN_1_FLUX_ERR)
real*4 FPN2 ! (mW/m2) ? EPIC/pn flux in energy band 2
* (PN_2_FLUX)
real*4 e_FPN2 ! (mW/m2) ? 1-sigma error on PN_2_FLUX
* (PN_2_FLUX_ERR)
real*4 FPN3 ! (mW/m2) ? EPIC/pn flux in energy band 3
* (PN_3_FLUX)
real*4 e_FPN3 ! (mW/m2) ? 1-sigma error on PN_3_FLUX
* (PN_3_FLUX_ERR)
real*4 FPN4 ! (mW/m2) ? EPIC/pn flux in energy band 4
* (PN_4_FLUX)
real*4 e_FPN4 ! (mW/m2) ? 1-sigma error on PN_4_FLUX
* (PN_4_FLUX_ERR)
real*4 FPN5 ! (mW/m2) ? EPIC/pn flux in energy band 5
* (PN_5_FLUX)
real*4 e_FPN5 ! (mW/m2) ? 1-sigma error on PN_5_FLUX
* (PN_5_FLUX_ERR)
real*4 FM1 ! (mW/m2) ? Total EPIC/MOS1 flux (M1_FLUX)
real*4 e_FM1 ! (mW/m2) ? 1-sigma error on M1_FLUX (M1_FLUX_ERR)
real*4 FM1_1 ! (mW/m2) ? EPIC/MOS1 flux in energy band 1
* (M1_1_FLUX)
real*4 e_FM1_1 ! (mW/m2) ? 1-sigma error on M1_1_FLUX
* (M1_1_FLUX_ERR)
real*4 FM1_2 ! (mW/m2) ? EPIC/MOS1 flux in energy band 2
* (M1_2_FLUX)
real*4 e_FM1_2 ! (mW/m2) ? 1-sigma error on M1_2_FLUX
* (M1_2_FLUX_ERR)
real*4 FM1_3 ! (mW/m2) ? EPIC/MOS1 flux in energy band 3
* (M1_3_FLUX)
real*4 e_FM1_3 ! (mW/m2) ? 1-sigma error on M1_3_FLUX
* (M1_3_FLUX_ERR)
real*4 FM1_4 ! (mW/m2) ? EPIC/MOS1 flux in energy band 4
* (M1_4_FLUX)
real*4 e_FM1_4 ! (mW/m2) ? 1-sigma error on M1_4_FLUX
* (M1_4_FLUX_ERR)
real*4 FM1_5 ! (mW/m2) ? EPIC/MOS1 flux in energy band 5
* (M1_5_FLUX)
real*4 e_FM1_5 ! (mW/m2) ? 1-sigma error on M1_5_FLUX
* (M1_5_FLUX_ERR)
real*4 FM2 ! (mW/m2) ? Total EPIC/MOS2 flux 1 (M2_FLUX)
real*4 e_FM2 ! (mW/m2) ? 1-sigma error on M2_FLUX (M2_FLUX_ERR)
real*4 FM2_1 ! (mW/m2) ? EPIC/MOS2 flux in energy band 1
* (M2_1_FLUX)
real*4 e_FM2_1 ! (mW/m2) ? 1-sigma error on MOS2_1_FLUX
* (M2_1_FLUX_ERR)
real*4 FM2_2 ! (mW/m2) ? EPIC/MOS2 flux in energy band 2
* (M2_2_FLUX)
real*4 e_FM2_2 ! (mW/m2) ? 1-sigma error on MOS2_2_FLUX
* (M2_2_FLUX_ERR)
real*4 FM2_3 ! (mW/m2) ? EPIC/MOS2 flux in energy band 3
* (M2_3_FLUX)
real*4 e_FM2_3 ! (mW/m2) ? 1-sigma error on MOS2_3_FLUX
* (M2_3_FLUX_ERR)
real*4 FM2_4 ! (mW/m2) ? EPIC/MOS2 flux in energy band 4
* (M2_4_FLUX)
real*4 e_FM2_4 ! (mW/m2) ? 1-sigma error on MOS2_4_FLUX
* (M2_4_FLUX_ERR)
real*4 FM2_5 ! (mW/m2) ? EPIC/MOS2 flux in energy band 5
* (M2_5_FLUX)
real*4 e_FM2_5 ! (mW/m2) ? 1-sigma error on MOS2_5_FLUX
* (M2_5_FLUX_ERR)
real*8 EPRate ! (ct/s) All-EPIC count rate (EP_RATE)
real*8 e_EPRate ! (ct/s) 1-sigma error on EP_RATE (EP_RATE_ERR)
real*8 PNRate ! (ct/s) ? Total EPIC/pn count rate (PN_RATE)
real*8 e_PNRate ! (ct/s) ? 1-sigma error on PN_RATE (PN_RATE_ERR)
real*8 PN1Rate ! (ct/s) ? EPIC/pn band 1 count rate (PN_1_RATE)
real*8 e_PN1Rate ! (ct/s) ? 1-sigma error on PN_1_RATE
* (PN_1_RATE_ERR)
real*8 PN2Rate ! (ct/s) ? EPIC/pn band 2 count rate (PN_2_RATE)
real*8 e_PN2Rate ! (ct/s) ? 1-sigma error on PN_2_RATE
* (PN_2_RATE_ERR)
real*8 PN3Rate ! (ct/s) ? EPIC/pn band 3 count rate (PN_3_RATE)
real*8 e_PN3Rate ! (ct/s) ? 1-sigma error on PN_3_RATE
* (PN_3_RATE_ERR)
real*8 PN4Rate ! (ct/s) ? EPIC/pn band 4 count rate (PN_4_RATE)
real*8 e_PN4Rate ! (ct/s) ? 1-sigma error on PN_4_RATE
* (PN_4_RATE_ERR)
real*8 PN5Rate ! (ct/s) ? EPIC/pn band 5 count rate (PN_5_RATE)
real*8 e_PN5Rate ! (ct/s) ? 1-sigma error on PN_5_RATE
* (PN_5_RATE_ERR)
real*8 M1Rate ! (ct/s) ? Total EPIC/MOS1 count rate (M1_RATE)
real*8 e_M1Rate ! (ct/s) ? 1-sigma error on M1_RATE (M1_RATE_ERR)
real*8 M1_1Rate ! (ct/s) ? EPIC/MOS1 band 1 count rate (M1_1_RATE)
real*8 e_M1_1Rate ! (ct/s) ? 1-sigma error on M1_1_RATE
* (M1_1_RATE_ERR)
real*8 M1_2Rate ! (ct/s) ? EPIC/MOS1 band 2 count rate (M1_2_RATE)
real*8 e_M1_2Rate ! (ct/s) ? 1-sigma error on M1_2_RATE
* (M1_2_RATE_ERR)
real*8 M1_3Rate ! (ct/s) ? EPIC/MOS1 band 3 count rate (M1_3_RATE)
real*8 e_M1_3Rate ! (ct/s) ? 1-sigma error on M1_3_RATE
* (M1_3_RATE_ERR)
real*8 M1_4Rate ! (ct/s) ? EPIC/MOS1 band 4 count rate (M1_4_RATE)
real*8 e_M1_4Rate ! (ct/s) ? 1-sigma error on M1_4_RATE
* (M1_4_RATE_ERR)
real*8 M1_5Rate ! (ct/s) ? EPIC/MOS1 band 5 count rate (M1_5_RATE)
real*8 e_M1_5Rate ! (ct/s) ? 1-sigma error on M1_5_RATE
* (M1_5_RATE_ERR)
real*8 M2Rate ! (ct/s) ? Total EPIC/MOS2 count rate (M2_RATE)
real*8 e_M2Rate ! (ct/s) ? 1-sigma error on MOS2_RATE
* (M2_RATE_ERR)
real*8 M2_1Rate ! (ct/s) ? EPIC/MOS2 band 1 count rate (M2_1_RATE)
real*8 e_M2_1Rate ! (ct/s) ? 1-sigma error on MOS2_1_RATE
* (M2_1_RATE_ERR)
real*8 M2_2Rate ! (ct/s) ? EPIC/MOS2 band 2 count rate (M2_2_RATE)
real*8 e_M2_2Rate ! (ct/s) ? 1-sigma error on MOS2_2_RATE
* (M2_2_RATE_ERR)
real*8 M2_3Rate ! (ct/s) ? EPIC/MOS2 band 3 count rate (M2_3_RATE)
real*8 e_M2_3Rate ! (ct/s) ? 1-sigma error on MOS2_3_RATE
* (M2_3_RATE_ERR)
real*8 M2_4Rate ! (ct/s) ? EPIC/MOS2 band 4 count rate (M2_4_RATE)
real*8 e_M2_4Rate ! (ct/s) ? 1-sigma error on MOS2_4_RATE
* (M2_4_RATE_ERR)
real*8 M2_5Rate ! (ct/s) ? EPIC/MOS2 band 5 count rate (M2_5_RATE)
real*8 e_M2_5Rate ! (ct/s) ? 1-sigma error on MOS2_5_RATE
* (M2_5_RATE_ERR)
real*4 EPcts ! (ct) All-EPIC number of counts (EP_CTS)
real*4 e_EPcts ! (ct) 1-sigma error on EP_CTS (EP_CTS_ERR)
real*4 PNcts ! (ct) ? EPIC/pn number of counts (PN_CTS)
real*4 e_PNcts ! (ct) ? 1-sigma error on PN_CTS (PN_CTS_ERR)
real*4 M1cts ! (ct) ? EPIC/MOS1 number of counts (M1_CTS)
real*4 e_M1cts ! (ct) ? 1-sigma error on M1_CTS (M1_CTS_ERR)
real*4 M2cts ! (ct) ? EPIC/MOS2 number of counts (M2_CTS)
real*4 e_M2cts ! (ct) ? 1-sigma error on M2_CTS (M2_CTS_ERR)
real*4 EPdETml ! All-EPIC maximum det. likelihood
* (EP_DET_ML)
real*4 PNdETml ! ? EPIC/pn maximum det. likelihood
* (PN_DET_ML)
real*4 PN_1dETml ! ? EPIC/pn band 1 max.det.likelihood
* (PN_1_DET_ML)
real*4 PN_2dETml ! ? EPIC/pn band 2 max.det.likelihood
* (PN_2_DET_ML)
real*4 PN_3dETml ! ? EPIC/pn band 3 max.det.likelihood
* (PN_3_DET_ML)
real*4 PN_4dETml ! ? EPIC/pn band 4 max.det.likelihood
* (PN_4_DET_ML)
real*4 PN_5dETml ! ? EPIC/pn band 5 max.det.likelihood
* (PN_5_DET_ML)
real*4 M1dETml ! ? EPIC/MOS1 maximum det. likelihood
* (M1_DET_ML)
real*4 M1_1dETml ! ? EPIC/MOS1 band 1 max.det.likelihood
* (M1_1_DET_ML)
real*4 M1_2dETml ! ? EPIC/MOS1 band 2 max.det.likelihood
* (M1_2_DET_ML)
real*4 M1_3dETml ! ? EPIC/MOS1 band 3 max.det.likelihood
* (M1_3_DET_ML)
real*4 M1_4dETml ! ? EPIC/MOS1 band 4 max.det.likelihood
* (M1_4_DET_ML)
real*4 M1_5dETml ! ? EPIC/MOS1 band 5 max.det.likelihood
* (M1_5_DET_ML)
real*4 M2dETml ! ? EPIC/MOS2 maximum det. likelihood
* (M2_DET_ML)
real*4 M2_1dETml ! ? EPIC/MOS2 band 1 max.det.likelihood
* (M2_1_DET_ML)
real*4 M2_2dETml ! ? EPIC/MOS2 band 2 max.det.likelihood
* (M2_2_DET_ML)
real*4 M2_3dETml ! ? EPIC/MOS2 band 3 max.det.likelihood
* (M2_3_DET_ML)
real*4 M2_4dETml ! ? EPIC/MOS2 band 4 max.det.likelihood
* (M2_4_DET_ML)
real*4 M2_5dETml ! ? EPIC/MOS2 band 5 max.det.likelihood
* (M2_5_DET_ML)
real*8 Ext ! (arcsec) Extent radius (EXTENT)
real*8 e_Ext ! (arcsec) 1-sigma error on EXTENT (EXTENT_ERR)
real*4 EtxML ! ? Extent likelihood (EXTENT_ML)
real*8 EPHR1 ! ? All-EPIC hardness ratio 1 (EP_HR1)
real*4 e_EPHR1 ! ? 1-sigma error on EP_HR1 (EP_HR1_ERR)
real*8 EPHR2 ! ? All-EPIC hardness ratio 2 (EP_HR2)
real*4 e_EPHR2 ! ? 1-sigma error on EP_HR2 (EP_HR2_ERR)
real*8 EPHR3 ! ? All-EPIC hardness ratio 3 (EP_HR3)
real*4 e_EPHR3 ! ? 1-sigma error on EP_HR3 (EP_HR3_ERR)
real*8 EPHR4 ! ? All-EPIC hardness ratio 4 (EP_HR4)
real*4 e_EPHR4 ! ? 1-sigma error on EP_HR4 (EP_HR4_ERR)
real*8 PNHR1 ! ? EPIC/pn hardness ratio 1 (PN_HR1)
real*4 e_PNHR1 ! ? 1-sigma error on PN_HR1 (PN_HR1_ERR)
real*8 PNHR2 ! ? EPIC/pn hardness ratio 2 (PN_HR2)
real*4 e_PNHR2 ! ? 1-sigma error on PN_HR2 (PN_HR2_ERR)
real*8 PNHR3 ! ? EPIC/pn hardness ratio 3 (PN_HR3)
real*4 e_PNHR3 ! ? 1-sigma error on PN_HR3 (PN_HR3_ERR)
real*8 PNHR4 ! ? EPIC/pn hardness ratio 4 (PN_HR4)
real*4 e_PNHR4 ! ? 1-sigma error on PN_HR4 (PN_HR4_ERR)
real*8 M1HR1 ! ? EPIC/MOS1 hardness ratio 1 (M1_HR1)
real*4 e_M1HR1 ! ? 1-sigma error on M1_HR1 (M1_HR1_ERR)
real*8 M1HR2 ! ? EPIC/MOS1 hardness ratio 2 (M1_HR2)
real*4 e_M1HR2 ! ? 1-sigma error on M1_HR2 (M1_HR2_ERR)
real*8 M1HR3 ! ? EPIC/MOS1 hardness ratio 3 (M1_HR3)
real*4 e_M1HR3 ! ? 1-sigma error on M1_HR3 (M1_HR3_ERR)
real*8 M1HR4 ! ? EPIC/MOS1 hardness ratio 4 (M1_HR4)
real*4 e_M1HR4 ! ? 1-sigma error on M1_HR4 (M1_HR4_ERR)
real*8 M2HR1 ! ? EPIC/MOS2 hardness ratio 1 (M2_HR1)
real*4 e_M2HR1 ! ? 1-sigma error on M2_HR1 (M2_HR1_ERR)
real*8 M2HR2 ! ? EPIC/MOS2 hardness ratio 2 (M2_HR2)
real*4 e_M2HR2 ! ? 1-sigma error on M2_HR2 (M2_HR2_ERR)
real*8 M2HR3 ! ? EPIC/MOS2 hardness ratio 3 (M2_HR3)
real*4 e_M2HR3 ! ? 1-sigma error on M2_HR3 (M2_HR3_ERR)
real*8 M2HR4 ! ? EPIC/MOS2 hardness ratio 4 (M2_HR4)
real*4 e_M2HR4 ! ? 1-sigma error on M2_HR4 (M2_HR4_ERR)
real*4 PNexp ! (s) ? PSF-weighted pn exposure, all (PN_EXP)
real*4 PN1exp ! (s) ? PSF-weighted pn exposure band 1
* (PN_1_EXP)
real*4 PN2exp ! (s) ? PSF-weighted pn exposure band 2
* (PN_2_EXP)
real*4 PN3exp ! (s) ? PSF-weighted pn exposure band 3
* (PN_3_EXP)
real*4 PN4exp ! (s) ? PSF-weighted pn exposure band 4
* (PN_4_EXP)
real*4 PN5exp ! (s) ? PSF-weighted pn exposure band 5
* (PN_5_EXP)
real*4 M1exp ! (s) ? PSF-weighted MOS1 exposure, all
* (M1_EXP)
real*4 M1_1exp ! (s) ? PSF-weighted MOS1 exposure band 1
* (M1_1_EXP)
real*4 M1_2exp ! (s) ? PSF-weighted MOS1 exposure band 2
* (M1_2_EXP)
real*4 M1_3exp ! (s) ? PSF-weighted MOS1 exposure band 3
* (M1_3_EXP)
real*4 M1_4exp ! (s) ? PSF-weighted MOS1 exposure band 4
* (M1_4_EXP)
real*4 M1_5exp ! (s) ? PSF-weighted MOS1 exposure band 5
* (M1_5_EXP)
real*4 M2exp ! (s) ? PSF-weighted MOS2 exposure, all
* (M2_EXP)
real*4 M2_1exp ! (s) ? PSF-weighted MOS2 exposure band 1
* (M2_1_EXP)
real*4 M2_2exp ! (s) ? PSF-weighted MOS2 exposure band 2
* (M2_2_EXP)
real*4 M2_3exp ! (s) ? PSF-weighted MOS2 exposure band 3
* (M2_3_EXP)
real*4 M2_4exp ! (s) ? PSF-weighted MOS2 exposure band 4
* (M2_4_EXP)
real*4 M2_5exp ! (s) ? PSF-weighted MOS2 exposure band 5
* (M2_5_EXP)
real*8 PNbg ! (ct/pix) ? EPIC/pn background map, all (PN_BG)
real*8 PN1bg ! (ct/pix) ? EPIC/pn background map band 1 (PN_1_BG)
real*8 PN2bg ! (ct/pix) ? EPIC/pn background map band 2 (PN_2_BG)
real*8 PN3bg ! (ct/pix) ? EPIC/pn background map band 3 (PN_3_BG)
real*8 PN4bg ! (ct/pix) ? EPIC/pn background map band 4 (PN_4_BG)
real*8 PN5bg ! (ct/pix) ? EPIC/pn background map band 5 (PN_5_BG)
real*8 M1bg ! (ct/pix) ? EPIC/MOS1 background map, all (M1_BG)
real*8 M1_1bg ! (ct/pix) ? EPIC/MOS1 background map band 1
* (M1_1_BG)
real*8 M1_2bg ! (ct/pix) ? EPIC/MOS1 background map band 2
* (M1_2_BG)
real*8 M1_3bg ! (ct/pix) ? EPIC/MOS1 background map band 3
* (M1_3_BG)
real*8 M1_4bg ! (ct/pix) ? EPIC/MOS1 background map band 4
* (M1_4_BG)
real*8 M1_5bg ! (ct/pix) ? EPIC/MOS1 background map band 5
* (M1_5_BG)
real*8 M2bg ! (ct/pix) ? EPIC/MOS2 background map, all (M2_BG)
real*8 M2_1bg ! (ct/pix) ? EPIC/MOS2 background map band 1
* (M2_1_BG)
real*8 M2_2bg ! (ct/pix) ? EPIC/MOS2 background map band 2
* (M2_2_BG)
real*8 M2_3bg ! (ct/pix) ? EPIC/MOS2 background map band 3
* (M2_3_BG)
real*8 M2_4bg ! (ct/pix) ? EPIC/MOS2 background map band 4
* (M2_4_BG)
real*8 M2_5bg ! (ct/pix) ? EPIC/MOS2 background map band 5
* (M2_5_BG)
real*4 EPontime ! (s) Total good exposure time all-EIC
* (EP_ONTIME)
real*4 PNontime ! (s) Total good pn exposure time (PN_ONTIME)
real*4 M1ontime ! (s) Total good MOS1 exposure time (M1_ONTIME)
real*4 M2ontime ! (s) Total good MOS2 exposure time (M2_ONTIME)
real*8 PNmaskfrac ! ? PSF-weighted CCD coverage pn
* (PN_MASKFRAC)
real*8 M1maskfrac ! ? PSF-weighted CCD coverage MOS1
* (M1_MASKFRAC)
real*8 M2maskfrac ! ? PSF-weighted CCD coverage MOS2
* (M2_MASKFRAC)
real*8 EPoffax ! (arcmin) ? Minimum EPIC offaxis angle (EP_OFFAX)
real*8 PNoffax ! (arcmin) ? pn offaxis angle (PN_OFFAX)
real*8 M1offax ! (arcmin) ? MOS1 offaxis angle (M1_OFFAX)
real*8 M2offax ! (arcmin) ? MOS2 offaxis angle (M2_OFFAX)
real*8 PN1vig ! ? pn vignetting factor band 1 (PN_1_VIG)
real*8 PN2vig ! ? pn vignetting factor band 2 (PN_2_VIG)
real*8 PN3vig ! ? pn vignetting factor band 3 (PN_3_VIG)
real*8 PN4vig ! ? pn vignetting factor band 4 (PN_4_VIG)
real*8 PN5vig ! ? pn vignetting factor band 5 (PN_5_VIG)
real*8 M1_1vig ! ? MOS1 vignetting factor band 1
* (M1_1_VIG)
real*8 M1_2vig ! ? MOS1 vignetting factor band 2
* (M1_2_VIG)
real*8 M1_3vig ! ? MOS1 vignetting factor band 3
* (M1_3_VIG)
real*8 M1_4vig ! ? MOS1 vignetting factor band 4
* (M1_4_VIG)
real*8 M1_5vig ! ? MOS1 vignetting factor band 5
* (M1_5_VIG)
real*8 M2_1vig ! ? MOS2 vignetting factor band 1
* (M2_1_VIG)
real*8 M2_2vig ! ? MOS2 vignetting factor band 2
* (M2_2_VIG)
real*8 M2_3vig ! ? MOS2 vignetting factor band 3
* (M2_3_VIG)
real*8 M2_4vig ! ? MOS2 vignetting factor band 4
* (M2_4_VIG)
real*8 M2_5vig ! ? MOS2 vignetting factor band 5
* (M2_5_VIG)
integer*4 fStack ! Integer quality flag (STACK_FLAG)
character*9 flagEP ! All-EPIC detection flags (EP_FLAG)
character*9 flagPN ! EPIC/pn detection flags (PN_FLAG)
character*9 flagM1 ! EPIC/MOS1 detection flags (M1_FLAG)
character*9 flagM2 ! EPIC/MOS2 detection flags (M2_FLAG)
real*4 Varchi2 ! ? Red. chi2 of all-EPIC variability
* (VAR_CHI2)
real*4 Varchi2_1 ! ? Red. chi2 of band 1 variability
* (VAR_CHI2_1)
real*4 Varchi2_2 ! ? Red. chi2 of band 2 variability
* (VAR_CHI2_2)
real*4 Varchi2_3 ! ? Red. chi2 of band 3 variability
* (VAR_CHI2_3)
real*4 Varchi2_4 ! ? Red. chi2 of band 4 variability
* (VAR_CHI2_4)
real*4 Varchi2_5 ! ? Red. chi2 of band 5 variability
* (VAR_CHI2_5)
real*8 VarProb ! ? Probability to be constant (VAR_PROB)
real*8 VarProb_1 ! ? Band 1 probability to be constant
* (VAR_PROB_1)
real*8 VarProb_2 ! ? Band 2 probability to be constant
* (VAR_PROB_2)
real*8 VarProb_3 ! ? Band 3 probability to be constant
* (VAR_PROB_3)
real*8 VarProb_4 ! ? Band 4 probability to be constant
* (VAR_PROB_4)
real*8 VarProb_5 ! ? Band 5 probability to be constant
* (VAR_PROB_5)
real*4 Fratio ! ? EPIC flux ratio (FRATIO)
real*4 e_Fratio ! ? 1-sigma error on FRATIO (Fratio-ERR)
real*4 Fratio_1 ! ? EPIC band 1 flux ratio (FRATIO_1)
real*4 e_Fratio_1 ! ? 1-sigma error on FRATIO_1
* (FRATIO_1_ERR)
real*4 Fratio_2 ! ? EPIC band 2 flux ratio (FRATIO_2)
real*4 e_Fratio_2 ! ? 1-sigma error on FRATIO_2
* (FRATIO_2_ERR)
real*4 Fratio_3 ! ? EPIC band 3 flux ratio (FRATIO_3)
real*4 e_Fratio_3 ! ? 1-sigma error on FRATIO_3
* (FRATIO_3_ERR)
real*4 Fratio_4 ! ? EPIC band 4 flux ratio (FRATIO_4)
real*4 e_Fratio_4 ! ? 1-sigma error on FRATIO_4
* (FRATIO_4_ERR)
real*4 Fratio_5 ! ? EPIC band 5 flux ratio (FRATIO_5)
real*4 e_Fratio_5 ! ? 1-sigma error on FRATIO_5
* (FRATIO_5_ERR)
real*8 Fluxvar ! ? Largest all-EPIC flux difference
* (FLUXVAR)
real*8 Fluxvar_1 ! ? Largest band 1 flux difference
* (FLUXVAR_1)
real*8 Fluxvar_2 ! ? Largest band 2 flux difference
* (FLUXVAR_2)
real*8 Fluxvar_3 ! ? Largest band 3 flux difference
* (FLUXVAR_3)
real*8 Fluxvar_4 ! ? Largest band 4 flux difference
* (FLUXVAR_4)
real*8 Fluxvar_5 ! ? Largest band 5 flux difference
* (FLUXVAR_5)
real*8 chi2probDR7 ! ? Probability of a DR7 match
* (CHI2PROB_3XMMDR7)
real*8 FvarDR7 ! ? Intra-obs. excess of DR7 match
* (FVAR_3XMMDR7)
real*8 e_FvarDR7 ! ? Intra-obs. excess error (DR7)
* (FVARERR_3XMMDR7)
integer*4 VDR7 ! [0,1] Intra-obs. variab. flag (DR7)
* (VAR_FLAG_3XMMDR7)
integer*4 SDR7 ! ? Quality flag of DR7 match
* (SUM_FLAG_3XMMDR7)
real*8 MJD0 ! (d) MJD of the (first) obs. start (MJD_FIRST)
real*8 MJD1 ! (d) MJD of the (last) obs. end (MJD_LAST)
integer*4 Rev ! ? XMM-Newton revolution number (REVOLUT)
real*8 PApnt ! (deg) ? Mean position angle (PA_PNT)
character*23 PNsubmode ! EPIC/pn submode (PN_SUBMODE)
character*15 M1submode ! EPIC/MOS1 submode (M1_SUBMODE)
character*15 M2submode ! EPIC/MOS2 submode (M2_SUBMODE)
character*6 PNfilter ! EPIC/pn filter (PN_FILTER)
character*6 M1filter ! EPIC/MOS1 filter (M1_FILTER)
character*6 M2filter ! EPIC/MOS2 filter (M2_FILTER)
character*54 URLDR7 ! Web page of the nearest unique source
* in 3XMM-DR7 (URL_3XMMDR7)
C=============================================================================
C Loading file '3xmmdr7s.dat' ! Source catalogue
C Format for file interpretation
1 format(
+ I16,1X,A10,1X,I2,1X,I2,1X,F10.6,1X,F10.6,1X,F9.6,1X,F10.6,1X,
+ F10.6,1X,F11.6,1X,F8.6,1X,F11.6,1X,F8.6,1X,F9.5,6X,A17,1X,A21,
+ 1X,I15,1X,I15,1X,F10.6,1X,F10.6,1X,F9.6,1X,F11.8,1X,I2,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,F9.6,1X,F8.6,1X,F9.6,1X,F8.6,1X,F8.6,1X,
+ F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F9.6,
+ 1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
+ F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,
+ 1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
+ F8.6,1X,F8.6,1X,F8.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,F9.6,1X,F9.6,1X,E12.6,1X,F9.6,
+ 1X,E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,E12.6,
+ 1X,F9.6,1X,E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,
+ E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,E12.6,1X,
+ F9.6,1X,E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,E12.6,1X,F9.6,1X,
+ E12.6,1X,F9.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,F10.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
+ F9.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F9.6,1X,F8.6,
+ 1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,
+ 1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
+ F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,
+ 1X,I1,1X,A9,1X,A9,1X,A9,1X,A9,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,
+ F8.6,1X,F8.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,F10.6,1X,F9.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F8.6,
+ 1X,F10.6,1X,F10.6,1X,I1,1X,I1,1X,F10.4,1X,F10.4,1X,I4,1X,
+ F10.6,1X,A23,1X,A15,1X,A15,1X,A6,1X,A6,1X,A6,1X,A54)
C Effective file loading
open(unit=1,status='old',file=
+'3xmmdr7s.dat')
write(6,*) '....Loading file: 3xmmdr7s.dat'
do i__=1,288344
read(1,'(A3243)')ar__
read(ar__,1)
+ Source,ObsID,Nobs,Ncont,RAdeg,DEdeg,srcML,GLON,GLAT,Ximg,
+ e_Ximg,Yimg,e_Yimg,DistNN,v3XMM,v3XMMDR7,SrcidDR7,DetidDR7,
+ RA7deg,DE7deg,ePos7,DistDR7,NdetDR7,FEP,e_FEP,FEP1,e_FEP1,
+ FEP2,e_FEP2,FEP3,e_FEP3,FEP4,e_FEP4,FEP5,e_FEP5,FPN,e_FPN,
+ FPN1,e_FPN1,FPN2,e_FPN2,FPN3,e_FPN3,FPN4,e_FPN4,FPN5,e_FPN5,
+ FM1,e_FM1,FM1_1,e_FM1_1,FM1_2,e_FM1_2,FM1_3,e_FM1_3,FM1_4,
+ e_FM1_4,FM1_5,e_FM1_5,FM2,e_FM2,FM2_1,e_FM2_1,FM2_2,e_FM2_2,
+ FM2_3,e_FM2_3,FM2_4,e_FM2_4,FM2_5,e_FM2_5,EPRate,e_EPRate,
+ PNRate,e_PNRate,PN1Rate,e_PN1Rate,PN2Rate,e_PN2Rate,PN3Rate,
+ e_PN3Rate,PN4Rate,e_PN4Rate,PN5Rate,e_PN5Rate,M1Rate,e_M1Rate,
+ M1_1Rate,e_M1_1Rate,M1_2Rate,e_M1_2Rate,M1_3Rate,e_M1_3Rate,
+ M1_4Rate,e_M1_4Rate,M1_5Rate,e_M1_5Rate,M2Rate,e_M2Rate,
+ M2_1Rate,e_M2_1Rate,M2_2Rate,e_M2_2Rate,M2_3Rate,e_M2_3Rate,
+ M2_4Rate,e_M2_4Rate,M2_5Rate,e_M2_5Rate,EPcts,e_EPcts,PNcts,
+ e_PNcts,M1cts,e_M1cts,M2cts,e_M2cts,EPdETml,PNdETml,PN_1dETml,
+ PN_2dETml,PN_3dETml,PN_4dETml,PN_5dETml,M1dETml,M1_1dETml,
+ M1_2dETml,M1_3dETml,M1_4dETml,M1_5dETml,M2dETml,M2_1dETml,
+ M2_2dETml,M2_3dETml,M2_4dETml,M2_5dETml,Ext,e_Ext,EtxML,EPHR1,
+ e_EPHR1,EPHR2,e_EPHR2,EPHR3,e_EPHR3,EPHR4,e_EPHR4,PNHR1,
+ e_PNHR1,PNHR2,e_PNHR2,PNHR3,e_PNHR3,PNHR4,e_PNHR4,M1HR1,
+ e_M1HR1,M1HR2,e_M1HR2,M1HR3,e_M1HR3,M1HR4,e_M1HR4,M2HR1,
+ e_M2HR1,M2HR2,e_M2HR2,M2HR3,e_M2HR3,M2HR4,e_M2HR4,PNexp,
+ PN1exp,PN2exp,PN3exp,PN4exp,PN5exp,M1exp,M1_1exp,M1_2exp,
+ M1_3exp,M1_4exp,M1_5exp,M2exp,M2_1exp,M2_2exp,M2_3exp,M2_4exp,
+ M2_5exp,PNbg,PN1bg,PN2bg,PN3bg,PN4bg,PN5bg,M1bg,M1_1bg,M1_2bg,
+ M1_3bg,M1_4bg,M1_5bg,M2bg,M2_1bg,M2_2bg,M2_3bg,M2_4bg,M2_5bg,
+ EPontime,PNontime,M1ontime,M2ontime,PNmaskfrac,M1maskfrac,
+ M2maskfrac,EPoffax,PNoffax,M1offax,M2offax,PN1vig,PN2vig,
+ PN3vig,PN4vig,PN5vig,M1_1vig,M1_2vig,M1_3vig,M1_4vig,M1_5vig,
+ M2_1vig,M2_2vig,M2_3vig,M2_4vig,M2_5vig,fStack,flagEP,flagPN,
+ flagM1,flagM2,Varchi2,Varchi2_1,Varchi2_2,Varchi2_3,Varchi2_4,
+ Varchi2_5,VarProb,VarProb_1,VarProb_2,VarProb_3,VarProb_4,
+ VarProb_5,Fratio,e_Fratio,Fratio_1,e_Fratio_1,Fratio_2,
+ e_Fratio_2,Fratio_3,e_Fratio_3,Fratio_4,e_Fratio_4,Fratio_5,
+ e_Fratio_5,Fluxvar,Fluxvar_1,Fluxvar_2,Fluxvar_3,Fluxvar_4,
+ Fluxvar_5,chi2probDR7,FvarDR7,e_FvarDR7,VDR7,SDR7,MJD0,MJD1,
+ Rev,PApnt,PNsubmode,M1submode,M2submode,PNfilter,M1filter,
+ M2filter,URLDR7
if(ar__(29:30) .EQ. '') Nobs = iNULL__
if(ar__(32:33) .EQ. '') Ncont = iNULL__
if(ar__(186:200) .EQ. '') SrcidDR7 = iNULL__
if(ar__(202:216) .EQ. '') DetidDR7 = iNULL__
if(ar__(218:227) .EQ. '') RA7deg = rNULL__
if(ar__(229:238) .EQ. '') DE7deg = rNULL__
if(ar__(240:248) .EQ. '') ePos7 = rNULL__
if(ar__(250:260) .EQ. '') DistDR7 = rNULL__
if(ar__(262:263) .EQ. '') NdetDR7 = iNULL__
if(ar__(265:276) .EQ. '') FEP = rNULL__
if(ar__(278:289) .EQ. '') e_FEP = rNULL__
if(ar__(421:432) .EQ. '') FPN = rNULL__
if(ar__(434:445) .EQ. '') e_FPN = rNULL__
if(ar__(447:458) .EQ. '') FPN1 = rNULL__
if(ar__(460:471) .EQ. '') e_FPN1 = rNULL__
if(ar__(473:484) .EQ. '') FPN2 = rNULL__
if(ar__(486:497) .EQ. '') e_FPN2 = rNULL__
if(ar__(499:510) .EQ. '') FPN3 = rNULL__
if(ar__(512:523) .EQ. '') e_FPN3 = rNULL__
if(ar__(525:536) .EQ. '') FPN4 = rNULL__
if(ar__(538:549) .EQ. '') e_FPN4 = rNULL__
if(ar__(551:562) .EQ. '') FPN5 = rNULL__
if(ar__(564:575) .EQ. '') e_FPN5 = rNULL__
if(ar__(577:588) .EQ. '') FM1 = rNULL__
if(ar__(590:601) .EQ. '') e_FM1 = rNULL__
if(ar__(603:614) .EQ. '') FM1_1 = rNULL__
if(ar__(616:627) .EQ. '') e_FM1_1 = rNULL__
if(ar__(629:640) .EQ. '') FM1_2 = rNULL__
if(ar__(642:653) .EQ. '') e_FM1_2 = rNULL__
if(ar__(655:666) .EQ. '') FM1_3 = rNULL__
if(ar__(668:679) .EQ. '') e_FM1_3 = rNULL__
if(ar__(681:692) .EQ. '') FM1_4 = rNULL__
if(ar__(694:705) .EQ. '') e_FM1_4 = rNULL__
if(ar__(707:718) .EQ. '') FM1_5 = rNULL__
if(ar__(720:731) .EQ. '') e_FM1_5 = rNULL__
if(ar__(733:744) .EQ. '') FM2 = rNULL__
if(ar__(746:757) .EQ. '') e_FM2 = rNULL__
if(ar__(759:770) .EQ. '') FM2_1 = rNULL__
if(ar__(772:783) .EQ. '') e_FM2_1 = rNULL__
if(ar__(785:796) .EQ. '') FM2_2 = rNULL__
if(ar__(798:809) .EQ. '') e_FM2_2 = rNULL__
if(ar__(811:822) .EQ. '') FM2_3 = rNULL__
if(ar__(824:835) .EQ. '') e_FM2_3 = rNULL__
if(ar__(837:848) .EQ. '') FM2_4 = rNULL__
if(ar__(850:861) .EQ. '') e_FM2_4 = rNULL__
if(ar__(863:874) .EQ. '') FM2_5 = rNULL__
if(ar__(876:887) .EQ. '') e_FM2_5 = rNULL__
if(ar__(908:916) .EQ. '') PNRate = rNULL__
if(ar__(918:925) .EQ. '') e_PNRate = rNULL__
if(ar__(927:934) .EQ. '') PN1Rate = rNULL__
if(ar__(936:943) .EQ. '') e_PN1Rate = rNULL__
if(ar__(945:952) .EQ. '') PN2Rate = rNULL__
if(ar__(954:961) .EQ. '') e_PN2Rate = rNULL__
if(ar__(963:970) .EQ. '') PN3Rate = rNULL__
if(ar__(972:979) .EQ. '') e_PN3Rate = rNULL__
if(ar__(981:988) .EQ. '') PN4Rate = rNULL__
if(ar__(990:997) .EQ. '') e_PN4Rate = rNULL__
if(ar__(999:1007) .EQ. '') PN5Rate = rNULL__
if(ar__(1009:1016) .EQ. '') e_PN5Rate = rNULL__
if(ar__(1018:1025) .EQ. '') M1Rate = rNULL__
if(ar__(1027:1034) .EQ. '') e_M1Rate = rNULL__
if(ar__(1036:1043) .EQ. '') M1_1Rate = rNULL__
if(ar__(1045:1052) .EQ. '') e_M1_1Rate = rNULL__
if(ar__(1054:1061) .EQ. '') M1_2Rate = rNULL__
if(ar__(1063:1070) .EQ. '') e_M1_2Rate = rNULL__
if(ar__(1072:1079) .EQ. '') M1_3Rate = rNULL__
if(ar__(1081:1088) .EQ. '') e_M1_3Rate = rNULL__
if(ar__(1090:1097) .EQ. '') M1_4Rate = rNULL__
if(ar__(1099:1106) .EQ. '') e_M1_4Rate = rNULL__
if(ar__(1108:1115) .EQ. '') M1_5Rate = rNULL__
if(ar__(1117:1124) .EQ. '') e_M1_5Rate = rNULL__
if(ar__(1126:1133) .EQ. '') M2Rate = rNULL__
if(ar__(1135:1142) .EQ. '') e_M2Rate = rNULL__
if(ar__(1144:1151) .EQ. '') M2_1Rate = rNULL__
if(ar__(1153:1160) .EQ. '') e_M2_1Rate = rNULL__
if(ar__(1162:1169) .EQ. '') M2_2Rate = rNULL__
if(ar__(1171:1178) .EQ. '') e_M2_2Rate = rNULL__
if(ar__(1180:1187) .EQ. '') M2_3Rate = rNULL__
if(ar__(1189:1196) .EQ. '') e_M2_3Rate = rNULL__
if(ar__(1198:1205) .EQ. '') M2_4Rate = rNULL__
if(ar__(1207:1214) .EQ. '') e_M2_4Rate = rNULL__
if(ar__(1216:1223) .EQ. '') M2_5Rate = rNULL__
if(ar__(1225:1232) .EQ. '') e_M2_5Rate = rNULL__
if(ar__(1260:1271) .EQ. '') PNcts = rNULL__
if(ar__(1273:1284) .EQ. '') e_PNcts = rNULL__
if(ar__(1286:1297) .EQ. '') M1cts = rNULL__
if(ar__(1299:1310) .EQ. '') e_M1cts = rNULL__
if(ar__(1312:1323) .EQ. '') M2cts = rNULL__
if(ar__(1325:1336) .EQ. '') e_M2cts = rNULL__
if(ar__(1351:1362) .EQ. '') PNdETml = rNULL__
if(ar__(1364:1375) .EQ. '') PN_1dETml = rNULL__
if(ar__(1377:1388) .EQ. '') PN_2dETml = rNULL__
if(ar__(1390:1401) .EQ. '') PN_3dETml = rNULL__
if(ar__(1403:1414) .EQ. '') PN_4dETml = rNULL__
if(ar__(1416:1427) .EQ. '') PN_5dETml = rNULL__
if(ar__(1429:1440) .EQ. '') M1dETml = rNULL__
if(ar__(1442:1453) .EQ. '') M1_1dETml = rNULL__
if(ar__(1455:1466) .EQ. '') M1_2dETml = rNULL__
if(ar__(1468:1479) .EQ. '') M1_3dETml = rNULL__
if(ar__(1481:1492) .EQ. '') M1_4dETml = rNULL__
if(ar__(1494:1505) .EQ. '') M1_5dETml = rNULL__
if(ar__(1507:1518) .EQ. '') M2dETml = rNULL__
if(ar__(1520:1531) .EQ. '') M2_1dETml = rNULL__
if(ar__(1533:1544) .EQ. '') M2_2dETml = rNULL__
if(ar__(1546:1557) .EQ. '') M2_3dETml = rNULL__
if(ar__(1559:1570) .EQ. '') M2_4dETml = rNULL__
if(ar__(1572:1583) .EQ. '') M2_5dETml = rNULL__
if(ar__(1605:1616) .EQ. '') EtxML = rNULL__
if(ar__(1618:1626) .EQ. '') EPHR1 = rNULL__
if(ar__(1628:1639) .EQ. '') e_EPHR1 = rNULL__
if(ar__(1641:1649) .EQ. '') EPHR2 = rNULL__
if(ar__(1651:1662) .EQ. '') e_EPHR2 = rNULL__
if(ar__(1664:1672) .EQ. '') EPHR3 = rNULL__
if(ar__(1674:1685) .EQ. '') e_EPHR3 = rNULL__
if(ar__(1687:1695) .EQ. '') EPHR4 = rNULL__
if(ar__(1697:1708) .EQ. '') e_EPHR4 = rNULL__
if(ar__(1710:1718) .EQ. '') PNHR1 = rNULL__
if(ar__(1720:1731) .EQ. '') e_PNHR1 = rNULL__
if(ar__(1733:1741) .EQ. '') PNHR2 = rNULL__
if(ar__(1743:1754) .EQ. '') e_PNHR2 = rNULL__
if(ar__(1756:1764) .EQ. '') PNHR3 = rNULL__
if(ar__(1766:1777) .EQ. '') e_PNHR3 = rNULL__
if(ar__(1779:1787) .EQ. '') PNHR4 = rNULL__
if(ar__(1789:1800) .EQ. '') e_PNHR4 = rNULL__
if(ar__(1802:1810) .EQ. '') M1HR1 = rNULL__
if(ar__(1812:1823) .EQ. '') e_M1HR1 = rNULL__
if(ar__(1825:1833) .EQ. '') M1HR2 = rNULL__
if(ar__(1835:1846) .EQ. '') e_M1HR2 = rNULL__
if(ar__(1848:1856) .EQ. '') M1HR3 = rNULL__
if(ar__(1858:1869) .EQ. '') e_M1HR3 = rNULL__
if(ar__(1871:1879) .EQ. '') M1HR4 = rNULL__
if(ar__(1881:1892) .EQ. '') e_M1HR4 = rNULL__
if(ar__(1894:1902) .EQ. '') M2HR1 = rNULL__
if(ar__(1904:1915) .EQ. '') e_M2HR1 = rNULL__
if(ar__(1917:1925) .EQ. '') M2HR2 = rNULL__
if(ar__(1927:1938) .EQ. '') e_M2HR2 = rNULL__
if(ar__(1940:1948) .EQ. '') M2HR3 = rNULL__
if(ar__(1950:1961) .EQ. '') e_M2HR3 = rNULL__
if(ar__(1963:1971) .EQ. '') M2HR4 = rNULL__
if(ar__(1973:1984) .EQ. '') e_M2HR4 = rNULL__
if(ar__(1986:1997) .EQ. '') PNexp = rNULL__
if(ar__(1999:2010) .EQ. '') PN1exp = rNULL__
if(ar__(2012:2023) .EQ. '') PN2exp = rNULL__
if(ar__(2025:2036) .EQ. '') PN3exp = rNULL__
if(ar__(2038:2049) .EQ. '') PN4exp = rNULL__
if(ar__(2051:2062) .EQ. '') PN5exp = rNULL__
if(ar__(2064:2075) .EQ. '') M1exp = rNULL__
if(ar__(2077:2088) .EQ. '') M1_1exp = rNULL__
if(ar__(2090:2101) .EQ. '') M1_2exp = rNULL__
if(ar__(2103:2114) .EQ. '') M1_3exp = rNULL__
if(ar__(2116:2127) .EQ. '') M1_4exp = rNULL__
if(ar__(2129:2140) .EQ. '') M1_5exp = rNULL__
if(ar__(2142:2153) .EQ. '') M2exp = rNULL__
if(ar__(2155:2166) .EQ. '') M2_1exp = rNULL__
if(ar__(2168:2179) .EQ. '') M2_2exp = rNULL__
if(ar__(2181:2192) .EQ. '') M2_3exp = rNULL__
if(ar__(2194:2205) .EQ. '') M2_4exp = rNULL__
if(ar__(2207:2218) .EQ. '') M2_5exp = rNULL__
if(ar__(2220:2229) .EQ. '') PNbg = rNULL__
if(ar__(2231:2239) .EQ. '') PN1bg = rNULL__
if(ar__(2241:2249) .EQ. '') PN2bg = rNULL__
if(ar__(2251:2259) .EQ. '') PN3bg = rNULL__
if(ar__(2261:2269) .EQ. '') PN4bg = rNULL__
if(ar__(2271:2279) .EQ. '') PN5bg = rNULL__
if(ar__(2281:2289) .EQ. '') M1bg = rNULL__
if(ar__(2291:2298) .EQ. '') M1_1bg = rNULL__
if(ar__(2300:2307) .EQ. '') M1_2bg = rNULL__
if(ar__(2309:2316) .EQ. '') M1_3bg = rNULL__
if(ar__(2318:2325) .EQ. '') M1_4bg = rNULL__
if(ar__(2327:2334) .EQ. '') M1_5bg = rNULL__
if(ar__(2336:2344) .EQ. '') M2bg = rNULL__
if(ar__(2346:2353) .EQ. '') M2_1bg = rNULL__
if(ar__(2355:2362) .EQ. '') M2_2bg = rNULL__
if(ar__(2364:2371) .EQ. '') M2_3bg = rNULL__
if(ar__(2373:2380) .EQ. '') M2_4bg = rNULL__
if(ar__(2382:2389) .EQ. '') M2_5bg = rNULL__
if(ar__(2443:2450) .EQ. '') PNmaskfrac = rNULL__
if(ar__(2452:2459) .EQ. '') M1maskfrac = rNULL__
if(ar__(2461:2468) .EQ. '') M2maskfrac = rNULL__
if(ar__(2470:2478) .EQ. '') EPoffax = rNULL__
if(ar__(2480:2488) .EQ. '') PNoffax = rNULL__
if(ar__(2490:2498) .EQ. '') M1offax = rNULL__
if(ar__(2500:2508) .EQ. '') M2offax = rNULL__
if(ar__(2510:2517) .EQ. '') PN1vig = rNULL__
if(ar__(2519:2526) .EQ. '') PN2vig = rNULL__
if(ar__(2528:2535) .EQ. '') PN3vig = rNULL__
if(ar__(2537:2544) .EQ. '') PN4vig = rNULL__
if(ar__(2546:2553) .EQ. '') PN5vig = rNULL__
if(ar__(2555:2562) .EQ. '') M1_1vig = rNULL__
if(ar__(2564:2571) .EQ. '') M1_2vig = rNULL__
if(ar__(2573:2580) .EQ. '') M1_3vig = rNULL__
if(ar__(2582:2589) .EQ. '') M1_4vig = rNULL__
if(ar__(2591:2598) .EQ. '') M1_5vig = rNULL__
if(ar__(2600:2607) .EQ. '') M2_1vig = rNULL__
if(ar__(2609:2616) .EQ. '') M2_2vig = rNULL__
if(ar__(2618:2625) .EQ. '') M2_3vig = rNULL__
if(ar__(2627:2634) .EQ. '') M2_4vig = rNULL__
if(ar__(2636:2643) .EQ. '') M2_5vig = rNULL__
if(ar__(2687:2698) .EQ. '') Varchi2 = rNULL__
if(ar__(2700:2711) .EQ. '') Varchi2_1 = rNULL__
if(ar__(2713:2724) .EQ. '') Varchi2_2 = rNULL__
if(ar__(2726:2737) .EQ. '') Varchi2_3 = rNULL__
if(ar__(2739:2750) .EQ. '') Varchi2_4 = rNULL__
if(ar__(2752:2763) .EQ. '') Varchi2_5 = rNULL__
if(ar__(2765:2772) .EQ. '') VarProb = rNULL__
if(ar__(2774:2781) .EQ. '') VarProb_1 = rNULL__
if(ar__(2783:2790) .EQ. '') VarProb_2 = rNULL__
if(ar__(2792:2799) .EQ. '') VarProb_3 = rNULL__
if(ar__(2801:2808) .EQ. '') VarProb_4 = rNULL__
if(ar__(2810:2817) .EQ. '') VarProb_5 = rNULL__
if(ar__(2819:2830) .EQ. '') Fratio = rNULL__
if(ar__(2832:2843) .EQ. '') e_Fratio = rNULL__
if(ar__(2845:2856) .EQ. '') Fratio_1 = rNULL__
if(ar__(2858:2869) .EQ. '') e_Fratio_1 = rNULL__
if(ar__(2871:2882) .EQ. '') Fratio_2 = rNULL__
if(ar__(2884:2895) .EQ. '') e_Fratio_2 = rNULL__
if(ar__(2897:2908) .EQ. '') Fratio_3 = rNULL__
if(ar__(2910:2921) .EQ. '') e_Fratio_3 = rNULL__
if(ar__(2923:2934) .EQ. '') Fratio_4 = rNULL__
if(ar__(2936:2947) .EQ. '') e_Fratio_4 = rNULL__
if(ar__(2949:2960) .EQ. '') Fratio_5 = rNULL__
if(ar__(2962:2973) .EQ. '') e_Fratio_5 = rNULL__
if(ar__(2975:2984) .EQ. '') Fluxvar = rNULL__
if(ar__(2986:2994) .EQ. '') Fluxvar_1 = rNULL__
if(ar__(2996:3005) .EQ. '') Fluxvar_2 = rNULL__
if(ar__(3007:3016) .EQ. '') Fluxvar_3 = rNULL__
if(ar__(3018:3027) .EQ. '') Fluxvar_4 = rNULL__
if(ar__(3029:3038) .EQ. '') Fluxvar_5 = rNULL__
if(ar__(3040:3047) .EQ. '') chi2probDR7 = rNULL__
if(ar__(3049:3058) .EQ. '') FvarDR7 = rNULL__
if(ar__(3060:3069) .EQ. '') e_FvarDR7 = rNULL__
if(ar__(3073:3073) .EQ. '') SDR7 = iNULL__
if(ar__(3097:3100) .EQ. '') Rev = iNULL__
if(ar__(3102:3111) .EQ. '') PApnt = rNULL__
c ..............Just test output...........
write(6,1)
+ Source,ObsID,Nobs,Ncont,RAdeg,DEdeg,srcML,GLON,GLAT,Ximg,
+ e_Ximg,Yimg,e_Yimg,DistNN,v3XMM,v3XMMDR7,SrcidDR7,DetidDR7,
+ RA7deg,DE7deg,ePos7,DistDR7,NdetDR7,FEP,e_FEP,FEP1,e_FEP1,
+ FEP2,e_FEP2,FEP3,e_FEP3,FEP4,e_FEP4,FEP5,e_FEP5,FPN,e_FPN,
+ FPN1,e_FPN1,FPN2,e_FPN2,FPN3,e_FPN3,FPN4,e_FPN4,FPN5,e_FPN5,
+ FM1,e_FM1,FM1_1,e_FM1_1,FM1_2,e_FM1_2,FM1_3,e_FM1_3,FM1_4,
+ e_FM1_4,FM1_5,e_FM1_5,FM2,e_FM2,FM2_1,e_FM2_1,FM2_2,e_FM2_2,
+ FM2_3,e_FM2_3,FM2_4,e_FM2_4,FM2_5,e_FM2_5,EPRate,e_EPRate,
+ PNRate,e_PNRate,PN1Rate,e_PN1Rate,PN2Rate,e_PN2Rate,PN3Rate,
+ e_PN3Rate,PN4Rate,e_PN4Rate,PN5Rate,e_PN5Rate,M1Rate,e_M1Rate,
+ M1_1Rate,e_M1_1Rate,M1_2Rate,e_M1_2Rate,M1_3Rate,e_M1_3Rate,
+ M1_4Rate,e_M1_4Rate,M1_5Rate,e_M1_5Rate,M2Rate,e_M2Rate,
+ M2_1Rate,e_M2_1Rate,M2_2Rate,e_M2_2Rate,M2_3Rate,e_M2_3Rate,
+ M2_4Rate,e_M2_4Rate,M2_5Rate,e_M2_5Rate,EPcts,e_EPcts,PNcts,
+ e_PNcts,M1cts,e_M1cts,M2cts,e_M2cts,EPdETml,PNdETml,PN_1dETml,
+ PN_2dETml,PN_3dETml,PN_4dETml,PN_5dETml,M1dETml,M1_1dETml,
+ M1_2dETml,M1_3dETml,M1_4dETml,M1_5dETml,M2dETml,M2_1dETml,
+ M2_2dETml,M2_3dETml,M2_4dETml,M2_5dETml,Ext,e_Ext,EtxML,EPHR1,
+ e_EPHR1,EPHR2,e_EPHR2,EPHR3,e_EPHR3,EPHR4,e_EPHR4,PNHR1,
+ e_PNHR1,PNHR2,e_PNHR2,PNHR3,e_PNHR3,PNHR4,e_PNHR4,M1HR1,
+ e_M1HR1,M1HR2,e_M1HR2,M1HR3,e_M1HR3,M1HR4,e_M1HR4,M2HR1,
+ e_M2HR1,M2HR2,e_M2HR2,M2HR3,e_M2HR3,M2HR4,e_M2HR4,PNexp,
+ PN1exp,PN2exp,PN3exp,PN4exp,PN5exp,M1exp,M1_1exp,M1_2exp,
+ M1_3exp,M1_4exp,M1_5exp,M2exp,M2_1exp,M2_2exp,M2_3exp,M2_4exp,
+ M2_5exp,PNbg,PN1bg,PN2bg,PN3bg,PN4bg,PN5bg,M1bg,M1_1bg,M1_2bg,
+ M1_3bg,M1_4bg,M1_5bg,M2bg,M2_1bg,M2_2bg,M2_3bg,M2_4bg,M2_5bg,
+ EPontime,PNontime,M1ontime,M2ontime,PNmaskfrac,M1maskfrac,
+ M2maskfrac,EPoffax,PNoffax,M1offax,M2offax,PN1vig,PN2vig,
+ PN3vig,PN4vig,PN5vig,M1_1vig,M1_2vig,M1_3vig,M1_4vig,M1_5vig,
+ M2_1vig,M2_2vig,M2_3vig,M2_4vig,M2_5vig,fStack,flagEP,flagPN,
+ flagM1,flagM2,Varchi2,Varchi2_1,Varchi2_2,Varchi2_3,Varchi2_4,
+ Varchi2_5,VarProb,VarProb_1,VarProb_2,VarProb_3,VarProb_4,
+ VarProb_5,Fratio,e_Fratio,Fratio_1,e_Fratio_1,Fratio_2,
+ e_Fratio_2,Fratio_3,e_Fratio_3,Fratio_4,e_Fratio_4,Fratio_5,
+ e_Fratio_5,Fluxvar,Fluxvar_1,Fluxvar_2,Fluxvar_3,Fluxvar_4,
+ Fluxvar_5,chi2probDR7,FvarDR7,e_FvarDR7,VDR7,SDR7,MJD0,MJD1,
+ Rev,PApnt,PNsubmode,M1submode,M2submode,PNfilter,M1filter,
+ M2filter,URLDR7
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end