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