FORTRAN Generation
(/./ftp/cats/IX/49)

Conversion of standardized ReadMe file for file /./ftp/cats/IX/49 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 2024-Mar-28
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/49               XXL Survey: First results                   (Pierre+, 2016)
*================================================================================
*The XXL Survey: First results.
*   Pierre M., Pacaud F., Adami C., Alis S., Altieri B., Baran B.,
*   Benoist C., Birkinshaw M., Bongiorno A., Bremer M.N., Brusa M.,
*   Butler A., Ciliegi P., Chiappetti L., Clerc N., Corasaniti P.S.,
*   Coupon J., De Breuck C., Democles J., Desai S., Delhaize J., Devriendt J.,
*   Dubois Y., Eckert D., Elyiv A., Ettori S., Evrard A., Faccioli L.,
*   Farahi A., Ferrari C., Finet F., Fotopoulou S., Fourmanoit N., Gandhi P.,
*   Gastaldello F., Gastaud R., Georgantopoulos I., P.Giles, Guennou L.,
*   Guglielmo V., Horellou C., Husband K., Huynh M., Iovino A., Kilbinger M.,
*   Koulouridis E., Lavoie S., Le Brun A.M.C., LeFevre J.P., Lidman C.,
*   Lieu M., Lin C.A., Mantz A., Maughan B.J., Maurogordato S.,
*   McCarthy I.G., McGee S., Melin J.B., Melnyk O., Menanteau F., Novak M.,
*   Paltani S., Plionis M., Poggianti B.M., Pomarede D., Pompei E.,
*   Ponman T.J., Ramos-Ceja M.E., Ranalli P., Rapetti D., Raychaudury S.,
*   Reiprich T.H., Rottgering H., Rozo E., Ryko E., Sadibekova T.,
*   Santos J., Sauvageot J.L., Schimd C., Sereno M., Smith G.P., Smolcic V.,
*   Snowden S., Spergel D., Stanford S., Surdej J., Valageas P., Valotti A.,
*   Valtchanov I., Vignali C., Willis J., Ziparo F.
*   <Astron. Astrophys. 592, A1-A12 (2016)>
*   =2016yCat.9049....0P
*   =2016A&A...592A...1P
*   +2016A&A...592A...2P
*   +2016A&A...592A...3G
*   +2016A&A...592A...4L
*   +2016A&A...592A...5F
*   +2016A&A...592A...6P
*   +2016A&A...592A...7A
*   +2016A&A...592A...8B
*   +2016A&A...592A...9Z
*   +2016A&A...592A..10S
*   +2016A&A...592A..11K
*   +2016A&A...592A..12E
*   +2016PASA...33....1L
*   +2016MNRAS.462.4141L
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'xxlpoint.dat'	! List of XMM pointings in XXL survey

      integer*4 nr__
      parameter (nr__=622)	! Number of records
      character*65 ar__   	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*10  ObsId      (nr__) ! ESA XMM log observation identifier
      character*11  FieldName  (nr__) ! XXL pointing identifier (1)
      real*8        RAdeg      (nr__) ! (deg) Pointing Right ascension (J2000)
      real*8        DEdeg      (nr__) ! (deg) Pointing Declination (J2000)
      real*4        ExpM1      (nr__) ! (ks) MOS1 exposure time (2)
      real*4        ExpM2      (nr__) ! (ks) MOS2 exposure time (2)
      real*4        Exppn      (nr__) ! (ks) pn exposure time (2)
      character*1   Qual       (nr__) ! [0-3] Quality Flag (3)
      character*1   Badfield   (nr__) ! [0-3] Badfield Flag (4)
      character*1   db         (nr__) ! [X] database presence flag (5)
      character*1   cat        (nr__) ! [X] catalogue presence flag (6)
*Note (1): XXLn or XXLs stand for the XXL-N (XXL-S) field;
*  a,b,c... tags indicate that the same a,b,c... tags indicate that the same
*  sky position has been observed several times in different AOs (consult the ESA
*  XMM log using the ESA ObsId) because the quality of earlier pointings was
*  insufficient; the z tag means that a fictitious pointing has been created
*  combining the events of all usable repeated pointings, in order to improve
*  the quality. In total there are 542 and 81 a,b,c and z pointings respectively.
*  In case of repeated fields, as well as of overlaps from adjacent fields, the
*  X-ray catalogue will remove overlapping detections, and consider only the one
*  from the better pointing, or, in case of equal quality, the object with the
*  smallest off-axis angle.
*Note (2): The exposure for the MOS and pn detectors is the remaining
*  exposure (in ks) after selection of the good-time intervals.
*Note (3): Quality flag as follows:
*  0 = good quality
*  1 = low exposure
*  2 = high background
*  3 = low exposure and high background
*Note (4): Badfield flag as follows:
*  0 = best acceptable observation at a given position
*  1 = for deep/good observation from the archives, not part of XXL proper
*  2 = other acceptable XXL observation at same position
*  3 = bad pointings i.e. quality=3
* This flag is used in the overlap removal procedure
*Note (5): This column is ticked if Xamin detected at least one object in this
*  pointing.
*Note (6): This column is ticked if at least one source in this pointing survived
*  the overlap removal procedure and hence entered the X-ray source catalogue.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'xxl100gc.dat'	! the XXL-100-GC catalogue

      integer*4 nr__1
      parameter (nr__1=107)	! Number of records
      character*128 ar__1  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     XLSSC      (nr__1) ! IAU XLSSC cluster identifier
      real*8        RAdeg_1    (nr__1) ! (deg) Cluster Right ascension (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Cluster Declination (J2000)
      real*4        z          (nr__1) ! Cluster redshift
      integer*4     Ngal       (nr__1) ! Number of galaxies used to determine z
      integer*4     C60        (nr__1) ! (ct) Counts within 60" in [0.5-2]keV
      real*4        F60        (nr__1) ! (10-17W/m2) Flux within 60" in [0.5-2]keV
      real*4        e_F60      (nr__1) ! (10-17W/m2) rms uncertainty on F60
      real*4        Ez         (nr__1) ! ? evolution term E(z)
      real*4        r500MT     (nr__1) ! (Mpc) Radius at overdensity 500 from MT
*                                       relation (1)
      real*4        T300kpc    (nr__1) ! (keV) ? Temperature in aperture 300kpc (2)
      real*4        E_T300kpc  (nr__1) ! (keV) ? Plus side error bar on T300kpc
      real*4        e_T300kpc_1(nr__1) ! (keV) ? Minus side error bar on T300kpc
      real*4        LXXL300kpc (nr__1) ! (10+36W) ? 0.5-2keV luminosity in 300kpc (3)
      real*4        e_LXXL300kpc(nr__1) ! (10+36W) ? rms uncertainty on LXXL300kpc
      real*4        Lbol500MT  (nr__1) ! (10+36W) ? bolometric luminosity in r500MT (1)
      real*4        e_Lbol500MT(nr__1) ! (10+36W) ? rms uncertainty on Lbol500MT
      real*4        M500MT     (nr__1) ! (10+13Msun) ? Mass within r500MT (4)
      real*4        e_M500MT   (nr__1) ! (10+13Msun) ? rms uncertainty on M500MT
      real*4        Mgas500    (nr__1) ! (10+13Msun) ? Gas mass within r500MT (5)
      real*4        e_Mgas500  (nr__1) ! (10+13Msun) ? rms uncertainty on Mgas500
      character*1   Flag       (nr__1) ! [0-2] Supplementary object flag (6)
*Note (1): r500MT from T300kpc and the M500-T300kpc relation of paper IV.
*Note (2): T300kpc: not core excised, from Table 1 of paper III.
*Note (3): LXXL300kpc and  Lbol500MT from Table 1 of Paper III.
*Note (4): M500MT from T300kpc and the M500-T300kpc relation of paper IV.
*Note (5): Mgas500 from Table A.1 of paper XIII.
*Note (6): Supplementary flag as follows:
*    0 = cluster in main sample of 100 brightest clusters (Table D.1 of paper II)
*    1 = supplementary cluster detected in bad pointing (Table D.2 of paper II)
*    2 = supplementary cluster not part of pipeline  (Table D.2 of paper II)

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'xxl1000a.dat'	! The XXL_1000_AGN catalogue

      integer*4 nr__2
      parameter (nr__2=1000)	! Number of records
      character*305 ar__2  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*16  v3XLSS     (nr__2) ! IAU-style 3XLSS catalogue name
*                                      (JHHMMSS.s+DDMMSS)
      real*8        RAdeg_2    (nr__2) ! (deg) X-ray Right ascension (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) X-ray Declination (J2000)
      real*4        Bflux      (nr__2) ! (mW/m2) ? Pipeline flux in band B [0.5-2 keV]
      real*4        CDflux     (nr__2) ! (mW/m2) ? Pipeline flux in band CD [2-10 keV]
      real*4        Fmode      (nr__2) ! (mW/m2) Mode of 2-10keV spectral fit flux
      real*4        e_Fmode    (nr__2) ! (mW/m2) 68% confidence lower limit for fit flux
      real*4        E_Fmode_1  (nr__2) ! (mW/m2) 68% confidence upper limit for fit flux
      real*4        Fmed       (nr__2) ! (mW/m2) Median of 2-10keV spectral fit flux
      real*4        Gmode      (nr__2) ! Mode of fitted photon index (Gammamode)
      real*4        e_Gmode    (nr__2) ! 68% confidence lower limit for Gamma
      real*4        E_Gmode_1  (nr__2) ! 68% confidence upper limit for Gamma
      real*4        Gmed       (nr__2) ! Median of fitted photon index (Gammamedian)
      real*4        NHmode     (nr__2) ! (cm-2) Mode of fitted NH (atom/cm^2^)
      real*4        e_NHmode   (nr__2) ! (cm-2) 68% confidence lower limit for NH
*                                      (atom/cm^2^)
      real*4        E_NHmode_1 (nr__2) ! (cm-2) 68% confidence upper limit for NH
*                                      (atom/cm^2^)
      real*4        NHmed      (nr__2) ! (cm-2) Median of fitted NH (atom/cm^2^)
      character*1   Xflag      (nr__2) ! [0-5] X-ray spectrum quality flag (1)
      real*8        RACdeg     (nr__2) ! (deg) ? Counterpart Right ascension (J2000)
      real*8        DECdeg     (nr__2) ! (deg) ? Counterpart Declination (J2000)
      real*4        gmag       (nr__2) ! (mag) ? g band magnitude
      real*4        e_gmag     (nr__2) ! (mag) ? Error on g band magnitude
      real*4        rmag       (nr__2) ! (mag) ? r band magnitude
      real*4        e_rmag     (nr__2) ! (mag) ? Error on r band magnitude
      real*4        imag       (nr__2) ! (mag) ? i band magnitude
      real*4        e_imag     (nr__2) ! (mag) ? Error on i band magnitude
      real*4        zmag       (nr__2) ! (mag) ? z band magnitude
      real*4        e_zmag     (nr__2) ! (mag) ? Error on z band magnitude
      character*4   r_photo    (nr__2) ! Parent survey(s) for griz photometry (2)
      real*4        zsp        (nr__2) ! ? Spectroscopic redshift
      character*19  r_zsp      (nr__2) ! ? Parent survey (spectroscopy)
      real*4        zph        (nr__2) ! ? Photometric redshift
      real*4        e_zph      (nr__2) ! ? 68% confidence zphot lower limit
      real*4        E_zph_1    (nr__2) ! ? 68% confidence zphot upper limit
      integer*4     zphClass   (nr__2) ! ? Classification for best zphot (3)
      real*4        Pstar      (nr__2) ! [0,1]? Probability to be a star
      real*4        Poutlier   (nr__2) ! [0,1]? Probability to be an outlier (4)
*Note (1): The values of the Xflag are specified in Table 2 of the paper.
*   Table 2: X-ray spectral quality classes defined based on the XMM pointing
*    quality and the signal-to-noise ratio (S/N) of the spectrum.
* ---------------------------------------------------------------------------
* Flag Quality class   N     FXamin/Ffit  PN counts     MOS1 counts  MOS2 counts
*      Pointing S/N       med  mean sigma min  max med  min max med  min max med
* ------------------------------------------------------------------------------
*  1    0       >6    547 1.00 1.0   0.3  25 11194 268  13 2626 105  6 2651 107
*  2    0/1/2   >4    281 1.01 1.1   0.5   8  1879 103   3  244  43  6  415  40
*  3    0/1/2    3-4   75 1.07 2     3     6   246 43    2   49  23  4   81  20
*  4    0/1/2   <3     96 0.78 52   191    3   177 16    1   34   6  1   55   5
*  5    NoData   -      1 -     -     -    -     -  -    2   2    -  1    1   -
* ------------------------------------------------------------------------------
*Note (2): 4-letter code where each character gives the origin of the griz
*   magnitude (blank if counterpart missing or not in the optical band).
*    The 1-character survey code is:
*    C = CFHT
*    S = SDSS
*    B = BCS
*    D = DECam
*Note (3): Photometric class as follows:
*    1 = passive
*    2 = starforming
*    3 = starburst
*    4 = AGN
*    5 = QSO
*Note (4): Probability to have |zph-zsp|>0.15(1+zsp).

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'xxl_vla.dat'	! VLA 3GHz sources in XXL-N

      integer*4 nr__3
      parameter (nr__3=155)	! Number of records
      character*141 ar__3  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     ID         (nr__3) ! Source identifier name
      real*8        RAdeg_3    (nr__3) ! (deg) Right ascension (J2000)
      real*8        DEdeg_3    (nr__3) ! (deg) Declination (J2000)
      integer*4     RAh        (nr__3) ! (h) Right ascension (J2000)
      integer*4     RAm        (nr__3) ! (min) Right ascension (J2000)
      real*4        RAs        (nr__3) ! (s) Right ascension (J2000)
      character*1   DE_        (nr__3) ! Declination sign (J2000)
      integer*4     DEd        (nr__3) ! (deg) Declination (J2000)
      integer*4     DEm        (nr__3) ! (arcmin) Declination (J2000)
      real*4        DEs        (nr__3) ! (arcsec) Declination (J2000)
      real*4        e_RAs      (nr__3) ! (arcsec) ?=-99.99 Positional error on RA (1)
      real*4        e_DEs      (nr__3) ! (arcsec) ?=-99.99 Positional error on DE (1)
      real*8        Speak      (nr__3) ! (uJy/beam) ?=-99.99 Peak surface brightness at 3GHz(1)
      real*4        e_Speak    (nr__3) ! (uJy/beam) ?=-99.99 Error on peak surface brightness (1)
      real*8        Stotal     (nr__3) ! (uJy) Total flux density at 3GHz (2)
      real*4        e_Stotal   (nr__3) ! (uJy) Error on total flux density
      character*1   Res        (nr__3) ! [0-1] Resolved flag (3)
      character*1   SL         (nr__3) ! [0-1] Sidelobe flag (4)
      character*1   Mult       (nr__3) ! [0-1] Multicomponent flag (5)
      real*4        maj        (nr__3) ! (arcsec) ?=0 Deconvolved major axis (1)(6)
      real*4        e_maj      (nr__3) ! (arcsec) ?=0 Deconvolved major axis (1)(6)
      real*4        min        (nr__3) ! (arcsec) ?=0 Deconvolved minor axis (1)(6)
      real*4        e_min      (nr__3) ! (arcsec) ?=0 Deconvolved minor axis (1)(6)
      real*4        PA         (nr__3) ! (deg) ?=0 Deconvolved position angle (1)(6)
      real*4        e_PA       (nr__3) ! (deg) ?=0 Deconvolved position angle (1)(6)
*Note (1): undefined (-99.99) for multicomponent sources
*Note (2): Equal to the peak surface brightness for unresolved sources
*Note (3): Resolved flag as follows:
*           0 = Unresolved
*           1 = Resolved
*Note (4): Sidelobe flag as follows:
*           1 = source located on a sidelobe
*Note (5): Multicomponent flag as follows:
*           0 = Single component
*           1 = Multicomponent
*Note (6): minor and major axes and PA set to 0 for unresolved sources.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'list.dat'	! Information on VLA fits image

      integer*4 nr__4
      parameter (nr__4=1)	! Number of records
      character*99 ar__4  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_4    (nr__4) ! (deg) Right Ascension of center (J2000)
      real*8        DEdeg_4    (nr__4) ! (deg) Declination of center (J2000)
      integer*4     Nx         (nr__4) ! Number of pixels along X-axis
      integer*4     Ny         (nr__4) ! Number of pixels along Y-axis
      character*26  Obs_Date   (nr__4) ! Observation date
      real*4        Freq       (nr__4) ! (Hz) Observed frequency
      integer*4     size       (nr__4) ! (Kibyte) Size of FITS file
      character*12  FileName   (nr__4) ! Name of FITS file
      character*13  Title      (nr__4) ! Title of the FITS file

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'xxl_atca.dat'	! ATCA pilot sample in XXL-S

      integer*4 nr__5
      parameter (nr__5=1389)	! Number of records
      character*159 ar__5  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*16  Name       (nr__5) ! Source identifier name (JHHMMSS.s+DDMMSS)
      real*8        RAdeg_5    (nr__5) ! (deg) Right ascension (J2000)
      real*8        DEdeg_5    (nr__5) ! (deg) Declination (J2000)
      integer*4     RAh_1      (nr__5) ! (h) Right ascension (J2000)
      integer*4     RAm_1      (nr__5) ! (min) Right ascension (J2000)
      real*4        RAs_1      (nr__5) ! (s) Right ascension (J2000)
      character*1   DE__1      (nr__5) ! Declination sign (J2000)
      integer*4     DEd_1      (nr__5) ! (deg) Declination (J2000)
      integer*4     DEm_1      (nr__5) ! (arcmin) Declination (J2000)
      real*4        DEs_1      (nr__5) ! (arcsec) Declination (J2000)
      real*4        e_RAs_1    (nr__5) ! (arcsec) ?=-99 Positional error on RA  (1)
      real*4        e_DEs_1    (nr__5) ! (arcsec) ?=-99 Positional error on DE (1)
      real*8        Speak_1    (nr__5) ! (mJy/beam) ?=-99 Peak surface brightness at 2.1GHz (1)
      real*8        e_Speak_1  (nr__5) ! (mJy/beam) ?=-99 Error on Speak (in mJy/beam) (1)
      real*8        Stotal_1   (nr__5) ! (mJy) Total flux density at 2.1GHz (2)
      real*8        e_Stotal_1 (nr__5) ! (mJy) ?=-99 Error on Stotal (1)
      real*4        rms        (nr__5) ! (mJy) Noise level at the source position
      real*4        maj_1      (nr__5) ! (arcsec) ?=0 Deconvolved major axis (3)
      real*4        min_1      (nr__5) ! (arcsec) ?=0 Deconvolved minor axis (3)
      real*4        PA_1       (nr__5) ! (deg) ?=0 Deconvolved position angle (3)
      character*1   Res_1      (nr__5) ! [0/1] Resolved flag  (4)
      character*1   Mult_1     (nr__5) ! [0/1] Multicomponent flag (5)
      real*8        SNR        (nr__5) ! Signal-to-noise ratio
*Note (1): undefined (-99) for multicomponent sources
*Note (2): Equal to the peak surface brightness for unresolved sources
*Note (3): minor and major axes and PA set to 0 for unresolved sources
*Note (4): Resolved flag as follows:
*           0 = Unresolved
*           1 = Resolved
*Note (5): Multicomponent flag as follows:
*           0 = Single component
*           1 = Multicomponent

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'xxlaaoz.dat'	! XXL AAOmega redshift catalogue

      integer*4 nr__6
      parameter (nr__6=3660)	! Number of records
      character*67 ar__6  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      character*19  XXLname    (nr__6) ! XXL Name (JHHMMSS.ss+DDMMSS.s)
      real*8        RAdeg_6    (nr__6) ! (deg) Right ascension (J2000)
      real*8        DEdeg_6    (nr__6) ! (deg) Declination (J2000)
      real*8        z_1        (nr__6) ! Redshift
      integer*4     q_z        (nr__6) ! Flag for redshift quality (1)
      character*6   Com        (nr__6) ! Comment
*Note (1): Redshift quality flag as follows:
*      6 = star
*      4 = the redshift is likely to be correct more than 99% of the time
*      3 = the redshift is likely to be correct about 95% of the time

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'xxl100bc.dat'	! Summary of XXL-100-GC clusters and BCGs properties

      integer*4 nr__7
      parameter (nr__7=85)	! Number of records
      character*84 ar__7  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     XLSSC_1    (nr__7) ! IAU XLSSC cluster identifier
      real*4        zCl        (nr__7) ! Cluster redshift
      integer*4     M500       (nr__7) ! (10+13Msun) ?=- Mass inside of r500
      integer*4     e_M500     (nr__7) ! (10+13Msun) ? rms uncertainty in M500
      real*4        r500       (nr__7) ! (Mpc) Radius at overdensity 500 from MT
*                                      relation (1)
      real*8        RAdeg_7    (nr__7) ! (deg) BCG right ascension (J2000)
      real*8        DEdeg_7    (nr__7) ! (deg) BCG declination (J2000)
      real*4        zBCG       (nr__7) ! ?=- BCG redshift
      real*4        Offset     (nr__7) ! (arcsec) BCG offset from the X-ray centroid
*                                      in arcsec unit
      real*4        Offset500  (nr__7) ! BCG offset from the X-ray centroid
*                                      in r500 unit
      real*4        BGCMass    (nr__7) ! (10+11Msun) BCG stellar mass
      real*4        E_BGCMass  (nr__7) ! (10+11Msun) Error on BGCMass (upper value)
      real*4        e_BGCMass_1(nr__7) ! (10+11Msun) Error on BGCMass (upper value)
      real*4        Dm12       (nr__7) ! (mag) z-band magnitude difference between the
*                                      brightest and second brightest cluster
*                                      members
*Note (1): r500 from T300kpc and the M500-T300kpc relation of paper IV
*  (Lieu et al., 2016A&A...592A...4L).

C=============================================================================

C  Loading file 'xxlpoint.dat'	! List of XMM pointings in XXL survey

C  Format for file interpretation

    1 format(
     +  A10,1X,A11,1X,F9.5,1X,F9.5,1X,F4.1,1X,F4.1,1X,F4.1,1X,A1,1X,
     +  A1,1X,A1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'xxlpoint.dat')
      write(6,*) '....Loading file: xxlpoint.dat'
      do i__=1,622
        read(1,'(A65)')ar__
        read(ar__,1)
     +  ObsId(i__),FieldName(i__),RAdeg(i__),DEdeg(i__),ExpM1(i__),
     +  ExpM2(i__),Exppn(i__),Qual(i__),Badfield(i__),db(i__),cat(i__)
c    ..............Just test output...........
        write(6,1)
     +  ObsId(i__),FieldName(i__),RAdeg(i__),DEdeg(i__),ExpM1(i__),
     +  ExpM2(i__),Exppn(i__),Qual(i__),Badfield(i__),db(i__),cat(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'xxl100gc.dat'	! the XXL-100-GC catalogue

C  Format for file interpretation

    2 format(
     +  6X,I3,1X,F7.3,1X,F7.3,1X,F5.3,1X,I2,1X,I4,1X,F5.2,1X,F4.2,1X,
     +  F4.2,1X,F4.2,1X,F5.2,2X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F6.2,
     +  1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'xxl100gc.dat')
      write(6,*) '....Loading file: xxl100gc.dat'
      do i__=1,107
        read(1,'(A128)')ar__1
        read(ar__1,2)
     +  XLSSC(i__),RAdeg_1(i__),DEdeg_1(i__),z(i__),Ngal(i__),
     +  C60(i__),F60(i__),e_F60(i__),Ez(i__),r500MT(i__),T300kpc(i__),
     +  E_T300kpc(i__),e_T300kpc_1(i__),LXXL300kpc(i__),
     +  e_LXXL300kpc(i__),Lbol500MT(i__),e_Lbol500MT(i__),M500MT(i__),
     +  e_M500MT(i__),Mgas500(i__),e_Mgas500(i__),Flag(i__)
        if(ar__1(52:55) .EQ. '') Ez(i__) = rNULL__
        if(ar__1(62:66) .EQ. '') T300kpc(i__) = rNULL__
        if(ar__1(69:73) .EQ. '') E_T300kpc(i__) = rNULL__
        if(ar__1(75:79) .EQ. '') e_T300kpc_1(i__) = rNULL__
        if(ar__1(81:85) .EQ. '') LXXL300kpc(i__) = rNULL__
        if(ar__1(87:90) .EQ. '') e_LXXL300kpc(i__) = rNULL__
        if(ar__1(92:97) .EQ. '') Lbol500MT(i__) = rNULL__
        if(ar__1(99:102) .EQ. '') e_Lbol500MT(i__) = rNULL__
        if(ar__1(104:108) .EQ. '') M500MT(i__) = rNULL__
        if(ar__1(110:114) .EQ. '') e_M500MT(i__) = rNULL__
        if(ar__1(116:120) .EQ. '') Mgas500(i__) = rNULL__
        if(ar__1(122:126) .EQ. '') e_Mgas500(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  XLSSC(i__),RAdeg_1(i__),DEdeg_1(i__),z(i__),Ngal(i__),
     +  C60(i__),F60(i__),e_F60(i__),Ez(i__),r500MT(i__),T300kpc(i__),
     +  E_T300kpc(i__),e_T300kpc_1(i__),LXXL300kpc(i__),
     +  e_LXXL300kpc(i__),Lbol500MT(i__),e_Lbol500MT(i__),M500MT(i__),
     +  e_M500MT(i__),Mgas500(i__),e_Mgas500(i__),Flag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'xxl1000a.dat'	! The XXL_1000_AGN catalogue

C  Format for file interpretation

    3 format(
     +  6X,A16,1X,F8.4,1X,F8.4,1X,E11.5,1X,E11.5,1X,E11.5,1X,E11.5,1X,
     +  E11.5,1X,E11.5,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,E11.5,1X,
     +  E11.5,1X,E11.5,1X,E11.5,1X,A1,1X,F8.4,1X,F8.4,1X,F5.2,1X,F4.2,
     +  1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A4,1X,F5.3,
     +  1X,A19,1X,F5.3,1X,F5.3,1X,F5.3,1X,I1,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'xxl1000a.dat')
      write(6,*) '....Loading file: xxl1000a.dat'
      do i__=1,1000
        read(1,'(A305)')ar__2
        read(ar__2,3)
     +  v3XLSS(i__),RAdeg_2(i__),DEdeg_2(i__),Bflux(i__),CDflux(i__),
     +  Fmode(i__),e_Fmode(i__),E_Fmode_1(i__),Fmed(i__),Gmode(i__),
     +  e_Gmode(i__),E_Gmode_1(i__),Gmed(i__),NHmode(i__),
     +  e_NHmode(i__),E_NHmode_1(i__),NHmed(i__),Xflag(i__),
     +  RACdeg(i__),DECdeg(i__),gmag(i__),e_gmag(i__),rmag(i__),
     +  e_rmag(i__),imag(i__),e_imag(i__),zmag(i__),e_zmag(i__),
     +  r_photo(i__),zsp(i__),r_zsp(i__),zph(i__),e_zph(i__),
     +  E_zph_1(i__),zphClass(i__),Pstar(i__),Poutlier(i__)
        if(ar__2(42:52) .EQ. '') Bflux(i__) = rNULL__
        if(ar__2(54:64) .EQ. '') CDflux(i__) = rNULL__
        if(ar__2(184:191) .EQ. '') RACdeg(i__) = rNULL__
        if(ar__2(193:200) .EQ. '') DECdeg(i__) = rNULL__
        if(ar__2(202:206) .EQ. '') gmag(i__) = rNULL__
        if(ar__2(208:211) .EQ. '') e_gmag(i__) = rNULL__
        if(ar__2(213:217) .EQ. '') rmag(i__) = rNULL__
        if(ar__2(219:222) .EQ. '') e_rmag(i__) = rNULL__
        if(ar__2(224:228) .EQ. '') imag(i__) = rNULL__
        if(ar__2(230:233) .EQ. '') e_imag(i__) = rNULL__
        if(ar__2(235:239) .EQ. '') zmag(i__) = rNULL__
        if(ar__2(241:244) .EQ. '') e_zmag(i__) = rNULL__
        if(ar__2(251:255) .EQ. '') zsp(i__) = rNULL__
        if(ar__2(277:281) .EQ. '') zph(i__) = rNULL__
        if(ar__2(283:287) .EQ. '') e_zph(i__) = rNULL__
        if(ar__2(289:293) .EQ. '') E_zph_1(i__) = rNULL__
        if(ar__2(295:295) .EQ. '') zphClass(i__) = iNULL__
        if(ar__2(297:300) .EQ. '') Pstar(i__) = rNULL__
        if(ar__2(302:305) .EQ. '') Poutlier(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  v3XLSS(i__),RAdeg_2(i__),DEdeg_2(i__),Bflux(i__),CDflux(i__),
     +  Fmode(i__),e_Fmode(i__),E_Fmode_1(i__),Fmed(i__),Gmode(i__),
     +  e_Gmode(i__),E_Gmode_1(i__),Gmed(i__),NHmode(i__),
     +  e_NHmode(i__),E_NHmode_1(i__),NHmed(i__),Xflag(i__),
     +  RACdeg(i__),DECdeg(i__),gmag(i__),e_gmag(i__),rmag(i__),
     +  e_rmag(i__),imag(i__),e_imag(i__),zmag(i__),e_zmag(i__),
     +  r_photo(i__),zsp(i__),r_zsp(i__),zph(i__),e_zph(i__),
     +  E_zph_1(i__),zphClass(i__),Pstar(i__),Poutlier(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'xxl_vla.dat'	! VLA 3GHz sources in XXL-N

C  Format for file interpretation

    4 format(
     +  I3,1X,F10.7,1X,F10.7,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,
     +  F5.2,1X,F6.2,1X,F6.2,1X,F7.2,1X,F6.2,1X,F7.2,1X,F5.2,1X,A1,1X,
     +  A1,1X,A1,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'xxl_vla.dat')
      write(6,*) '....Loading file: xxl_vla.dat'
      do i__=1,155
        read(1,'(A141)')ar__3
        read(ar__3,4)
     +  ID(i__),RAdeg_3(i__),DEdeg_3(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),e_RAs(i__),e_DEs(i__),
     +  Speak(i__),e_Speak(i__),Stotal(i__),e_Stotal(i__),Res(i__),
     +  SL(i__),Mult(i__),maj(i__),e_maj(i__),min(i__),e_min(i__),
     +  PA(i__),e_PA(i__)
c    ..............Just test output...........
        write(6,4)
     +  ID(i__),RAdeg_3(i__),DEdeg_3(i__),RAh(i__),RAm(i__),RAs(i__),
     +  DE_(i__),DEd(i__),DEm(i__),DEs(i__),e_RAs(i__),e_DEs(i__),
     +  Speak(i__),e_Speak(i__),Stotal(i__),e_Stotal(i__),Res(i__),
     +  SL(i__),Mult(i__),maj(i__),e_maj(i__),min(i__),e_min(i__),
     +  PA(i__),e_PA(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'list.dat'	! Information on VLA fits image

C  Format for file interpretation

    5 format(F9.5,F9.5,1X,I4,1X,I4,1X,A26,1X,E9.4,1X,I6,1X,A12,1X,A13)

C  Effective file loading

      open(unit=1,status='old',file=
     +'list.dat')
      write(6,*) '....Loading file: list.dat'
      do i__=1,1
        read(1,'(A99)')ar__4
        read(ar__4,5)
     +  RAdeg_4(i__),DEdeg_4(i__),Nx(i__),Ny(i__),Obs_Date(i__),
     +  Freq(i__),size(i__),FileName(i__),Title(i__)
c    ..............Just test output...........
        write(6,5)
     +  RAdeg_4(i__),DEdeg_4(i__),Nx(i__),Ny(i__),Obs_Date(i__),
     +  Freq(i__),size(i__),FileName(i__),Title(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'xxl_atca.dat'	! ATCA pilot sample in XXL-S

C  Format for file interpretation

    6 format(
     +  9X,A16,1X,F11.7,1X,F11.7,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,
     +  1X,F6.3,1X,F6.2,1X,F6.2,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,
     +  F5.3,1X,F5.2,1X,F5.2,1X,F6.2,1X,A1,1X,A1,1X,F7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'xxl_atca.dat')
      write(6,*) '....Loading file: xxl_atca.dat'
      do i__=1,1389
        read(1,'(A159)')ar__5
        read(ar__5,6)
     +  Name(i__),RAdeg_5(i__),DEdeg_5(i__),RAh_1(i__),RAm_1(i__),
     +  RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
     +  e_RAs_1(i__),e_DEs_1(i__),Speak_1(i__),e_Speak_1(i__),
     +  Stotal_1(i__),e_Stotal_1(i__),rms(i__),maj_1(i__),min_1(i__),
     +  PA_1(i__),Res_1(i__),Mult_1(i__),SNR(i__)
c    ..............Just test output...........
        write(6,6)
     +  Name(i__),RAdeg_5(i__),DEdeg_5(i__),RAh_1(i__),RAm_1(i__),
     +  RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
     +  e_RAs_1(i__),e_DEs_1(i__),Speak_1(i__),e_Speak_1(i__),
     +  Stotal_1(i__),e_Stotal_1(i__),rms(i__),maj_1(i__),min_1(i__),
     +  PA_1(i__),Res_1(i__),Mult_1(i__),SNR(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'xxlaaoz.dat'	! XXL AAOmega redshift catalogue

C  Format for file interpretation

    7 format(12X,A19,1X,F7.3,1X,F8.4,1X,F9.6,1X,I1,1X,A6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'xxlaaoz.dat')
      write(6,*) '....Loading file: xxlaaoz.dat'
      do i__=1,3660
        read(1,'(A67)')ar__6
        read(ar__6,7)
     +  XXLname(i__),RAdeg_6(i__),DEdeg_6(i__),z_1(i__),q_z(i__),
     +  Com(i__)
c    ..............Just test output...........
        write(6,7)
     +  XXLname(i__),RAdeg_6(i__),DEdeg_6(i__),z_1(i__),q_z(i__),
     +  Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'xxl100bc.dat'	! Summary of XXL-100-GC clusters and BCGs properties

C  Format for file interpretation

    8 format(
     +  6X,I3,1X,F5.3,1X,I2,1X,I2,1X,F5.3,1X,F8.4,1X,F8.4,1X,F5.3,1X,
     +  F4.1,1X,F5.3,1X,F5.2,2X,F4.2,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'xxl100bc.dat')
      write(6,*) '....Loading file: xxl100bc.dat'
      do i__=1,85
        read(1,'(A84)')ar__7
        read(ar__7,8)
     +  XLSSC_1(i__),zCl(i__),M500(i__),e_M500(i__),r500(i__),
     +  RAdeg_7(i__),DEdeg_7(i__),zBCG(i__),Offset(i__),
     +  Offset500(i__),BGCMass(i__),E_BGCMass(i__),e_BGCMass_1(i__),
     +  Dm12(i__)
        if (M500(i__) .EQ. 45) M500(i__) =  iNULL__
        if(ar__7(20:21) .EQ. '') e_M500(i__) = iNULL__
        if (idig(ar__7(47:51)).EQ.0) zBCG(i__) =  rNULL__
c    ..............Just test output...........
        write(6,8)
     +  XLSSC_1(i__),zCl(i__),M500(i__),e_M500(i__),r500(i__),
     +  RAdeg_7(i__),DEdeg_7(i__),zBCG(i__),Offset(i__),
     +  Offset500(i__),BGCMass(i__),E_BGCMass(i__),e_BGCMass_1(i__),
     +  Dm12(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end