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