Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/648/A4 into FORTRAN code for reading data files line by line.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Oct-13
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/A+A/648/A4 LoTSS Deep Fields DR1 photometric redshifts (Duncan+, 2021)
*================================================================================
*The LOFAR Two-metre Sky Survey: Deep Fields Data Release 1.
*IV. Photometric redshifts and stellar masses.
* Duncan K.J., Kondapally R., Brown M.J.I., Bonato M., Best P.N.,
* Roettgering H.J.A., Bondi M., Bowler R.A.A., Cochrane R.K., Guerkan G.,
* Hardcastle M.J., Jarvis M.J., Kunert-Bajraszewska M., Leslie S.K., Malek K.,
* Morabito L.K., O'Sullivan S.P., Prandoni I., Sabater J., Shimwell T.W.,
* Smith D.J.B., Wang L., Wolowska A.
* <Astron. Astrophys. 648, A4 (2021)>
* =2021A&A...648A...4D (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'bootes.dat' ! Bootes Field Photometric redshift catalogue
integer*4 nr__
parameter (nr__=2214329) ! Number of records
character*460 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID ! Unique object identifier for multi-wavelength
* catalogue in this field
integer*4 IDOpt ! Identifier for the optical-NIR detected chi2
* catalogue (ID_OPTICAL)
integer*4 IDSpitzer ! ?=-99 Identifier for the Spitzer detected chi2
* catalogue (ID_SPITZER)
real*8 RAdeg ! (deg) Right ascension (J2000) (RA)
real*8 DEdeg ! (deg) Declination (J2000) (DEC)
real*4 E_B_V ! E(B-V) galactic extinction for the source
* based on its position and
* Schlegel+(1998ApJ...500..525S)
real*8 Class ! Stellarity Parameter (CLASSSTAR)
integer*4 overlap ! [0/1] Bit flag indicating the multi-wavelength
* coverage (FLAGOVERLAP)
integer*4 clean ! [0/3] 3-valued flag indicating bright star
* masking (FLAGCLEAN)
real*8 zbest ! Best available redshift estimate (Z_BEST)
integer*4 f_zbest ! [0/1] Source of z_best, 0 = photometric,
* 1 = spectroscopic (Z_BEST_SOURCE)
real*8 zspec ! ?=-99 Literature Spectroscopic Redshift
* (Z_SPEC)
integer*4 r_zspec ! ?=-99 Source of the spectroscopic redshift
* (Z_SOURCE)
integer*4 q_zspec ! ?=-99 Spectroscopic redshift quality: flag Q=3
* means probable, Q>=4 means reliable (Z_QUAL)
real*8 z1med ! Median of the primary redshift peak above
* 80% HPD CI (z1_median)
real*8 z1min ! Lower bound of the primary 80% HPD CI peak
* (z1_min)
real*8 z1max ! Upper bound of the primary 80% HPD CI peak
* (z1_max)
real*8 z1area ! Integrated area of the primary 80% HPD CI peak
* (z1_area)
real*8 z2med ! Median of the secondary redshift peak
* (if present) above 80% HPD CI (z2_median)
real*8 z2min ! Lower bound of the secondary 80% HPD CI peak
* (z2_min)
real*8 z2max ! Upper bound of the secondary 80% HPD CI peak
* (z2_max)
real*8 z2area ! Integrated area of the secondary 80% HPD CI
* peak (z2_area)
integer*4 nfiltEazy ! ?=-99 Number of filters included in EAZY
* template fit (nfilt_eazy) (G1)
integer*4 nfiltAtl ! ?=-99 Number of filters included in Atlas+AGN
* template fit (nfilt_atlas) (G2)
integer*4 nfiltAna ! ?=-99 Number of filters included in Ananna+
* template fit (nfilt_ananna) (G3)
real*8 chirBest ! ?=-99 Chi^2 / nfilt for best-fit galaxy/AGN
* template (any library) (chi_r_best)
real*8 chirStelR ! ?=-99 Chi^2 / nfilt for best-fit stellar
* template (chi_r_stellar)
character*6 StelType ! Stellar type of best-fit stellar template
* (stellar_type)
integer*4 AGN ! [0/1] Sources flagged by any one of
* optAGN/IRAGN/XrayAGN
integer*4 optAGN ! [0/1] Source is included in Million Quasar
* Catalog (Flesch 2015PASA...32...10F)
* compilation or spectroscopically identified
* AGN
integer*4 IRAGN ! [0/1] Source satisfies Donley et al.
* (2012ApJ...748..142D) IR AGN criteria
integer*4 XrayAGN ! [0/1] Source has X-ray counterpart
real*4 XrayFlux ! (W/m2) ?=-9.90E-02 XBootes Soft X-ray Flux
* (XrayFlux_0.5-2)
real*8 HRX ! ?=-99 XBootes X-ray Hardness Ratio
real*8 aptmz ! ?=-99 Estimated z-band total flux correction
* from Legacy Surveys model flux (z)
* (ap_to_model_z)
real*8 e_aptmz ! ?=999 Statistical uncertainty on tot_z
* (ap_to_model_err_z)
real*8 aptmzS ! ?=-99 Estimated z-band total flux correction
* from Legacy Surveys model flux (z_Subaru)
* (ap_to_model_z_Subaru)
real*8 e_aptmzS ! ?=999 Statistical uncertainty on tot_z_Subaru
* (ap_to_model_err_z_Subaru)
real*8 zmodel ! ?=-99 Model grid redshift used in stellar mass
* fit
real*8 chibest ! ?=-99 Minimum Chi^2 in stellar mass fit
* (for Z_BEST) (chi_best)
real*8 Massmed ! ([Msun]) ?=-99 50th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_median)
real*8 Massl68 ! ([Msun]) ?=-99 16th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_l68)
real*8 Massu68 ! ([Msun]) ?=-99 84th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_u68)
integer*4 Nfilts ! ?=-99 Number of bands included in stellar
* mass fit
real*8 umag ! (mag) ?=-99 Rest-frame u-band magnitude for best-fit
* SED (for Z_BEST) (u_rest)
real*8 Bwmag ! (mag) ?=-99 Rest-frame Bw-band magnitude for
* best-fit SED (for Z_BEST) (Bw_rest)
real*8 Rmag ! (mag) ?=-99 Rest-frame R-band magnitude for
* best-fit SED (for Z_BEST) (R_rest)
real*8 imag ! (mag) ?=-99 Rest-frame I-band magnitude for
* best-fit SED (for Z_BEST) (I_rest)
real*8 zmag ! (mag) ?=-99 Rest-frame z-band magnitude for
* best-fit SED (for Z_BEST) (z_rest)
real*8 zSmag ! (mag) ?=-99 Rest-frame z_Subaru-band magnitude for
* best-fit SED (for Z_BEST) (z_Subaru_rest)
real*8 ymag ! (mag) ?=-99 Rest-frame y-band magnitude for
* best-fit SED (for Z_BEST) (y_rest)
real*8 Jmag ! (mag) ?=-99 Rest-frame J-band magnitude for
* best-fit SED (for Z_BEST) (J_rest)
real*8 Hmag ! (mag) ?=-99 Rest-frame H-band magnitude for
* best-fit SED (for Z_BEST) (H_rest)
real*8 Kmag ! (mag) ?=-99 Rest-frame K-band magnitude for
* best-fit SED (for Z_BEST) (K_rest)
real*8 Ksmag ! (mag) ?=-99 Rest-frame Ks-band magnitude for
* best-fit SED (for Z_BEST) (Ks_rest)
real*8 IRAC1 ! (mag) ?=-99 Rest-frame IRAC 3.6micron magnitude for
* best-fit SED (for Z_BEST) (z_rest)
real*8 IRAC2 ! (mag) ?=-99 Rest-frame IRAC 4.5micron magnitude for
* best-fit SED (for Z_BEST) (z_rest)
real*8 IRAC3 ! (mag) ?=-99 Rest-frame IRAC 5.8micron magnitude for
* best-fit SED (for Z_BEST) (z_rest)
real*8 IRAC4 ! (mag) ?=-99 Rest-frame IRAC 8.0micron magnitude for
* best-fit SED (for Z_BEST) (z_rest)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'en1.dat' ! ELAIS-N1 Photometric redshift catalogue
integer*4 nr__1
parameter (nr__1=2105993) ! Number of records
character*505 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID_1 ! Unique object identifier for multi-wavelength
* catalogue in this field
integer*4 IDOpt_1 ! Identifier for the optical-NIR detected chi2
* catalogue (ID_OPTICAL)
integer*4 IDSpitzer_1 ! ?=-99 Identifier for the Spitzer detected chi2
* catalogue (ID_SPITZER)
real*8 RAdeg_1 ! (deg) Right ascension (J2000) (RA)
real*8 DEdeg_1 ! (deg) Declination (J2000) (DEC)
real*4 E_B_V_1 ! (mag) E(B-V) galactic extinction for the source
* based on its position and
* Schlegel et al. (1998ApJ...500..525S)
real*4 Class_1 ! Stellarity Parameter
integer*4 overlap_1 ! [0/7] Bit flag indicating the multi-wavelength
* coverage
integer*4 clean_1 ! [1/3] 3-valued flag indicating bright star
* masking
real*8 zbest_1 ! Best available redshift estimate (Z_BEST)
integer*4 f_zbest_1 ! [0/1] Source of z_best, 0 = photometric,
* 1 = spectroscopic (Z_BEST_SOURCE)
real*8 zspec_1 ! ?=-99 Literature Spectroscopic Redshift
* (Z_SPEC)
integer*4 e_zspec ! ?=-99 Source of the spectroscopic redshift
* (Z_SOURCE)
integer*4 q_zspec_1 ! ?=-99 Spectroscopic redshift quality: flag Q=3
* means probable, Q>=4 means reliable (Z_QUAL)
real*8 z1med_1 ! Median of the primary redshift peak above
* 80% HPD CI (z1_median)
real*8 z1min_1 ! Lower bound of the primary 80% HPD CI peak
* (z1_min)
real*8 z1max_1 ! Upper bound of the primary 80% HPD CI peak
* (z1_max)
real*8 z1area_1 ! Integrated area of the primary 80% HPD CI peak
* (z1_area)
real*8 z2med_1 ! Median of the secondary redshift peak
* (if present) above 80% HPD CI (z2_median)
real*8 z2min_1 ! Lower bound of the secondary 80% HPD CI peak
* (z2_min)
real*8 z2max_1 ! Upper bound of the secondary 80% HPD CI peak
* (z2_max)
real*8 z2area_1 ! Integrated area of the secondary 80% HPD CI
* peak (z2_area)
integer*4 nfiltEazy_1 ! ?=-99 Number of filters included in EAZY
* template fit (nfilt_eazy) (G1)
integer*4 nfiltAtl_1 ! ?=-99 Number of filters included in Atlas+AGN
* template fit (nfilt_atlas) (G2)
integer*4 nfiltAna_1 ! ?=-99 Number of filters included in Ananna+
* template fit (nfilt_ananna) (G3)
real*8 chirBest_1 ! ?=-99 Chi^2 / nfilt for best-fit galaxy/AGN
* template (any library) (chi_r_best)
real*8 chirStel ! ?=-99 Chi^2 / nfilt for best-fit stellar
* template (chi_r_stellar)
character*6 StelType_1 ! Stellar type of best-fit stellar template
* (stellar_type)
integer*4 AGN_1 ! [0/1] Sources flagged by any one of
* optAGN/IRAGN/XrayAGN
integer*4 optAGN_1 ! [0/1] Source is included in Million Quasar
* Catalog (Flesch 2015PASA...32...10F)
* compilation or spectroscopically identified
* AGN
integer*4 IRAGN_1 ! [0/1] Source satisfies Donley et al.
* (2012ApJ...748..142D) IR AGN criteria
integer*4 XrayAGN_1 ! [0/1] Source has X-ray counterpart
character*21 v2RXSID ! ? ID for 2RXS X-ray Catalog counterpart (G4)
character*23 XMMSL2ID ! ? ID for XMMSL2 X-ray Catalog counterpart (G4)
real*8 aptmg ! ?=-99 Estimated g-band total flux correction
* from Legacy Surveys model flux (g)
* (ap_to_model_g)
real*8 e_aptmg ! ?=999 Statistical uncertainty on atpmg
* (ap_to_model_err_g)
real*8 aptmr ! ?=-99 Estimated r-band total flux correction
* from Legacy Surveys model flux (r)
* (ap_to_model_r)
real*8 e_aptmr ! ?=999 Statistical uncertainty on atpmg
* (ap_to_model_err_r)
real*8 aptmz_1 ! ?=-99 Estimated z-band total flux correction
* from Legacy Surveys model flux (z)
* (ap_to_model_z)
real*8 e_aptmz_1 ! ?=999 Statistical uncertainty on atpmg
* (ap_to_model_err_z)
real*8 zmodel_1 ! ?=-99 Model grid redshift used in stellar
* mass fit
real*8 chibest_1 ! ?=-99 Minimum Chi^2 in stellar mass fit
* (for Z_BEST)
* (chi_best)
real*8 Massmed_1 ! ([Msun]) ?=-99 50th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_median)
real*8 Massl68_1 ! ([Msun]) ?=-99 16th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_l68)
real*8 Massu68_1 ! ([Msun]) ?=-99 84th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_u68)
integer*4 Nfilts_1 ! ?=-99 Number of bands included in stellar
* mass fit
real*8 umag_1 ! (mag) ?=-99 Rest-frame u-band magnitude for
* best-fit SED (for Z_BEST) (u_rest)
real*8 gmag ! (mag) ?=-99 Rest-frame g-band magnitude for
* best-fit SED (for Z_BEST) (g_rest)
real*8 rmag_1 ! (mag) ?=-99 Rest-frame r-band magnitude for
* best-fit SED (for Z_BEST) (r_rest)
real*8 imag_1 ! (mag) ?=-99 Rest-frame i-band magnitude for
* best-fit SED (for Z_BEST) (i_rest)
real*8 zmag_1 ! (mag) ?=-99 Rest-frame z-band magnitude for
* best-fit SED (for Z_BEST) (z_rest)
real*8 ymag_1 ! (mag) ?=-99 Rest-frame y-band magnitude for
* best-fit SED (for Z_BEST) (y_rest)
real*8 Jmag_1 ! (mag) ?=-99 Rest-frame J-band magnitude for
* best-fit SED (for Z_BEST) (J_rest)
real*8 Kmag_1 ! (mag) ?=-99 Rest-frame K-band magnitude for
* best-fit SED (for Z_BEST) (K_rest)
real*8 IRAC1vs ! (mag) ?=-99 Rest-frame IRAC 3.6micron magnitude for
* best-fit SED (for Z_BEST) (ch1_servs_rest)
real*8 IRAC2vs ! (mag) ?=-99 Rest-frame IRAC 4.5micron magnitude for
* best-fit SED (for Z_BEST) (ch2_servs_rest)
real*8 IRAC1sw ! (mag) ?=-99 Rest-frame IRAC 3.6micron magnitude for
* best-fit SED (for Z_BEST) (ch1_swire_rest)
real*8 IRAC2sw ! (mag) ?=-99 Rest-frame IRAC 4.5micron magnitude for
* best-fit SED (for Z_BEST) (ch2_swire_rest)
real*8 IRAC3sw ! (mag) ?=-99 Rest-frame IRAC 5.8micron magnitude for
* best-fit SED (for Z_BEST) (ch3_swire_rest)
real*8 IRAC4sw ! (mag) ?=-99 Rest-frame IRAC 8.0micron magnitude for
* best-fit SED (for Z_BEST) (ch4_swire_rest)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lockman.dat' ! Lockman Hole photometric redshift catalogue
integer*4 nr__2
parameter (nr__2=3041794) ! Number of records
character*502 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID_2 ! Unique object identifier for multi-wavelength
* catalogue in this field
integer*4 IDOpt_2 ! Identifier for the optical-NIR detected chi2
* catalogue (ID_OPTICAL)
integer*4 IDSpitzer_2 ! ?=-99 Identifier for the Spitzer detected chi2
* catalogue (ID_SPITZER)
real*8 RAdeg_2 ! (deg) Right ascension (J2000) (RA)
real*8 DEdeg_2 ! (deg) Declination (J2000) (DEC)
real*4 E_B_V_1 ! (mag) E(B-V) galactic extinction for the source
* based on its position and
* Schlegel et al. (1998ApJ...500..525S)
real*4 Class_2 ! Stellarity Parameter (CLASSSTAR)
integer*4 clean_2 ! [1/3] 3-valued flag indicating bright star
* masking (FLAGCLEAN)
real*8 zbest_2 ! ?=-99 Best available redshift estimate
* (Z_BEST)
integer*4 f_zbest_2 ! [0/1] Source of z_best, 0 = photometric,
* 1 = spectroscopic (Z_BEST_SOURCE)
real*8 zspec_2 ! ?=-99 Literature Spectroscopic Redshift
* (Z_SPEC)
real*8 z1med_2 ! ?=-99 Median of the primary redshift peak
* above 80% HPD CI (z1_median)
real*8 z1min_2 ! ?=-99 Lower bound of the primary 80% HPD CI
* peak (z1_min)
real*8 z1max_2 ! ?=-99 Upper bound of the primary 80% HPD CI
* peak (z1_max)
real*8 z1area_2 ! ?=-99 Integrated area of the primary 80%
* HPD CI peak (z1_area)
real*8 z2med_2 ! ?=-99 Median of the secondary redshift peak
* (if present) above 80% HPD CI (z2_median)
real*8 z2min_2 ! ?=-99 Lower bound of the secondary 80% HPD
* CI peak (z2_min)
real*8 z2max_2 ! ?=-99 Upper bound of the secondary 80% HPD
* CI peak (z2_max)
real*8 z2area_2 ! ?=-99 Integrated area of the secondary 80%
* HPD CI peak (z2_area)
integer*4 nfiltEazy_2 ! ?=-99 Number of filters included in EAZY
* template fit (nfilt_eazy) (G1)
integer*4 nfiltAtl_2 ! ?=-99 Number of filters included in Atlas+AGN
* template fit (nfilt_atlas) (G2)
integer*4 nfiltAna_2 ! ?=-99 Number of filters included in Ananna+
* template fit (nfilt_ananna) (G3)
real*8 chirBest_2 ! ?=-99 Chi^2 / nfilt for best-fit galaxy/AGN
* template (any library) (chi_r_best)
real*8 chirStel_1 ! ?=-99 Chi^2 / nfilt for best-fit stellar
* template (chi_r_stellar)
character*6 StelType_2 ! Stellar type of best-fit stellar template
* (stellar_type)
integer*4 AGN_2 ! [0/1] Sources flagged by any one of
* optAGN/IRAGN/XrayAGN
integer*4 optAGN_2 ! [0/1] Source is included in Million Quasar
* Catalog (Flesch 2015PASA...32...10F)
* compilation or spectroscopically identified
* AGN
integer*4 IRAGN_2 ! [0/1] Source satisfies Donley et al.
* (2012ApJ...748..142D) IR AGN criteria
integer*4 XrayAGN_2 ! ?=-99 Source has X-ray counterpart
character*21 v2RXSID_1 ! ? ID for 2RXS X-ray Catalog counterpart
character*23 XMMSL2ID_1 ! ? ID for XMMSL2 X-ray Catalog counterpart
real*8 aptmg_1 ! ?=-99 Estimated g-band total flux correction
* from Legacy Surveys model flux (g)
* (ap_to_model_g)
real*8 e_aptmg_1 ! ?=999 Statistical uncertainty on atpmg
* (ap_to_model_err_g)
real*8 aptmr_1 ! ?=-99 Estimated r-band total flux correction
* from Legacy Surveys model flux (r)
* (ap_to_model_r)
real*8 e_aptmr_1 ! ?=999 Statistical uncertainty on atpmr
* (ap_to_model_err_r)
real*8 aptmz_2 ! ?=-99 Estimated z-band total flux correction
* from Legacy Surveys model flux (z)
* (ap_to_model_z)
real*8 e_aptmz_2 ! ?=999 Statistical uncertainty on atpmz
* (ap_to_model_err_z)
real*8 zmodel_2 ! ?=-99 Model grid redshift used in stellar
* mass fit
real*8 chibest_2 ! ?=-99 Minimum Chi^2 in stellar mass fit
* (for Z_BEST) (chi_best)
real*8 Massmed_2 ! ([Msun]) ?=-99 50th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_median)
real*8 Massl68_2 ! ([Msun]) ?=-99 16th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_l68)
real*8 Massu68_2 ! ([Msun]) ?=-99 84th percentile of the marginalised
* stellar mass posterior (for Z_BEST)
* (Mass_u68)
integer*4 Nfilts_2 ! ?=-99 Number of bands included in stellar
* mass fit
real*8 umag_2 ! (mag) ?=-99 Rest-frame u-band magnitude for
* best-fit SED (for Z_BEST) (u_rest)
real*8 gmag_1 ! (mag) ?=-99 Rest-frame g-band magnitude for
* best-fit SED (for Z_BEST) (g_rest)
real*8 rmag_2 ! (mag) ?=-99 Rest-frame r-band magnitude for
* best-fit SED (for Z_BEST) (r_rest)
real*8 zmag_2 ! (mag) ?=-99 Rest-frame z-band magnitude for
* best-fit SED (for Z_BEST) (z_rest)
real*8 grcsmag ! (mag) ?=-99 Rest-frame g_rcs-band magnitude for
* best-fit SED (for Z_BEST) (g_rcs_rest)
real*8 rrcsmag ! (mag) ?=-99 Rest-frame r_rcs-band magnitude for
* best-fit SED (for Z_BEST) (r_rcs_rest)
real*8 ircsmag ! (mag) ?=-99 Rest-frame i_rcs-band magnitude for
* best-fit SED (for Z_BEST) (i_rcs_rest)
real*8 zrcsmag ! (mag) ?=-99 Rest-frame z_rcs-band magnitude for
* best-fit SED (for Z_BEST) (z_rcs_rest)
real*8 Jmag_2 ! (mag) ?=-99 Rest-frame J-band magnitude for
* best-fit SED (for Z_BEST) (J_rest)
real*8 Kmag_2 ! (mag) ?=-99 Rest-frame K-band magnitude for
* best-fit SED (for Z_BEST) (K_rest)
real*8 IRAC1sv ! (mag) ?=-99 Rest-frame IRAC 3.6micron magnitude for
* best-fit SED (for Z_BEST) (ch1_servs_rest)
real*8 IRAC2sv ! (mag) ?=-99 Rest-frame IRAC 4.5micron magnitude for
* best-fit SED (for Z_BEST) (ch2_servs_rest)
real*8 IRAC1sw_1 ! (mag) ?=-99 Rest-frame IRAC 3.6micron magnitude for
* best-fit SED (for Z_BEST) (ch1_swire_rest)
real*8 IRAC2sw_1 ! (mag) ?=-99 Rest-frame IRAC 4.5micron magnitude for
* best-fit SED (for Z_BEST) (ch2_swire_rest)
real*8 IRAC3sw_1 ! (mag) ?=-99 Rest-frame IRAC 5.8micron magnitude for
* best-fit SED (for Z_BEST) (ch3_swire_rest)
real*8 IRAC4sw_1 ! (mag) ?=-99 Rest-frame IRAC 8.0micron magnitude for
* best-fit SED (for Z_BEST) (ch4_swire_rest)
C=============================================================================
C Loading file 'bootes.dat' ! Bootes Field Photometric redshift catalogue
C Format for file interpretation
1 format(
+ I7,1X,I7,1X,I6,1X,F16.12,1X,F15.12,1X,F6.4,1X,F8.4,1X,I1,1X,
+ I1,1X,F8.4,1X,I1,1X,F8.4,1X,I3,1X,I3,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,I3,1X,I3,1X,I3,1X,
+ F8.4,1X,F8.4,1X,A6,1X,I1,1X,I1,1X,I1,1X,I1,1X,E9.2,1X,F8.4,1X,
+ F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F9.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,I3,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'bootes.dat')
write(6,*) '....Loading file: bootes.dat'
do i__=1,2214329
read(1,'(A460)')ar__
read(ar__,1)
+ ID,IDOpt,IDSpitzer,RAdeg,DEdeg,E_B_V,Class,overlap,clean,
+ zbest,f_zbest,zspec,r_zspec,q_zspec,z1med,z1min,z1max,z1area,
+ z2med,z2min,z2max,z2area,nfiltEazy,nfiltAtl,nfiltAna,chirBest,
+ chirStelR,StelType,AGN,optAGN,IRAGN,XrayAGN,XrayFlux,HRX,
+ aptmz,e_aptmz,aptmzS,e_aptmzS,zmodel,chibest,Massmed,Massl68,
+ Massu68,Nfilts,umag,Bwmag,Rmag,imag,zmag,zSmag,ymag,Jmag,Hmag,
+ Kmag,Ksmag,IRAC1,IRAC2,IRAC3,IRAC4
c ..............Just test output...........
write(6,1)
+ ID,IDOpt,IDSpitzer,RAdeg,DEdeg,E_B_V,Class,overlap,clean,
+ zbest,f_zbest,zspec,r_zspec,q_zspec,z1med,z1min,z1max,z1area,
+ z2med,z2min,z2max,z2area,nfiltEazy,nfiltAtl,nfiltAna,chirBest,
+ chirStelR,StelType,AGN,optAGN,IRAGN,XrayAGN,XrayFlux,HRX,
+ aptmz,e_aptmz,aptmzS,e_aptmzS,zmodel,chibest,Massmed,Massl68,
+ Massu68,Nfilts,umag,Bwmag,Rmag,imag,zmag,zSmag,ymag,Jmag,Hmag,
+ Kmag,Ksmag,IRAC1,IRAC2,IRAC3,IRAC4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'en1.dat' ! ELAIS-N1 Photometric redshift catalogue
C Format for file interpretation
2 format(
+ I7,1X,I7,1X,I6,1X,F16.12,1X,F15.12,1X,F6.4,1X,F6.4,1X,I1,1X,
+ I1,1X,F8.4,1X,I1,1X,F8.4,1X,I3,1X,I3,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,I3,1X,I3,1X,I3,1X,
+ F8.4,1X,F8.4,1X,A6,1X,I1,1X,I1,1X,I1,1X,I1,1X,A21,1X,A23,1X,
+ F8.4,1X,F8.4,1X,F15.4,1X,F12.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F9.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,I3,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'en1.dat')
write(6,*) '....Loading file: en1.dat'
do i__=1,2105993
read(1,'(A505)')ar__1
read(ar__1,2)
+ ID_1,IDOpt_1,IDSpitzer_1,RAdeg_1,DEdeg_1,E_B_V_1,Class_1,
+ overlap_1,clean_1,zbest_1,f_zbest_1,zspec_1,e_zspec,q_zspec_1,
+ z1med_1,z1min_1,z1max_1,z1area_1,z2med_1,z2min_1,z2max_1,
+ z2area_1,nfiltEazy_1,nfiltAtl_1,nfiltAna_1,chirBest_1,
+ chirStel,StelType_1,AGN_1,optAGN_1,IRAGN_1,XrayAGN_1,v2RXSID,
+ XMMSL2ID,aptmg,e_aptmg,aptmr,e_aptmr,aptmz_1,e_aptmz_1,
+ zmodel_1,chibest_1,Massmed_1,Massl68_1,Massu68_1,Nfilts_1,
+ umag_1,gmag,rmag_1,imag_1,zmag_1,ymag_1,Jmag_1,Kmag_1,IRAC1vs,
+ IRAC2vs,IRAC1sw,IRAC2sw,IRAC3sw,IRAC4sw
c ..............Just test output...........
write(6,2)
+ ID_1,IDOpt_1,IDSpitzer_1,RAdeg_1,DEdeg_1,E_B_V_1,Class_1,
+ overlap_1,clean_1,zbest_1,f_zbest_1,zspec_1,e_zspec,q_zspec_1,
+ z1med_1,z1min_1,z1max_1,z1area_1,z2med_1,z2min_1,z2max_1,
+ z2area_1,nfiltEazy_1,nfiltAtl_1,nfiltAna_1,chirBest_1,
+ chirStel,StelType_1,AGN_1,optAGN_1,IRAGN_1,XrayAGN_1,v2RXSID,
+ XMMSL2ID,aptmg,e_aptmg,aptmr,e_aptmr,aptmz_1,e_aptmz_1,
+ zmodel_1,chibest_1,Massmed_1,Massl68_1,Massu68_1,Nfilts_1,
+ umag_1,gmag,rmag_1,imag_1,zmag_1,ymag_1,Jmag_1,Kmag_1,IRAC1vs,
+ IRAC2vs,IRAC1sw,IRAC2sw,IRAC3sw,IRAC4sw
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'lockman.dat' ! Lockman Hole photometric redshift catalogue
C Format for file interpretation
3 format(
+ I7,1X,I7,1X,I6,1X,F16.12,1X,F15.12,1X,F6.4,1X,F6.4,1X,I1,1X,
+ F8.4,1X,I1,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F8.4,1X,F8.4,1X,I3,1X,I3,1X,I3,1X,F8.4,1X,F8.4,1X,A6,
+ 1X,I1,1X,I1,1X,I1,1X,I1,1X,A21,1X,A23,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F9.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,I3,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'lockman.dat')
write(6,*) '....Loading file: lockman.dat'
do i__=1,3041794
read(1,'(A502)')ar__2
read(ar__2,3)
+ ID_2,IDOpt_2,IDSpitzer_2,RAdeg_2,DEdeg_2,E_B_V_1,Class_2,
+ clean_2,zbest_2,f_zbest_2,zspec_2,z1med_2,z1min_2,z1max_2,
+ z1area_2,z2med_2,z2min_2,z2max_2,z2area_2,nfiltEazy_2,
+ nfiltAtl_2,nfiltAna_2,chirBest_2,chirStel_1,StelType_2,AGN_2,
+ optAGN_2,IRAGN_2,XrayAGN_2,v2RXSID_1,XMMSL2ID_1,aptmg_1,
+ e_aptmg_1,aptmr_1,e_aptmr_1,aptmz_2,e_aptmz_2,zmodel_2,
+ chibest_2,Massmed_2,Massl68_2,Massu68_2,Nfilts_2,umag_2,
+ gmag_1,rmag_2,zmag_2,grcsmag,rrcsmag,ircsmag,zrcsmag,Jmag_2,
+ Kmag_2,IRAC1sv,IRAC2sv,IRAC1sw_1,IRAC2sw_1,IRAC3sw_1,IRAC4sw_1
c ..............Just test output...........
write(6,3)
+ ID_2,IDOpt_2,IDSpitzer_2,RAdeg_2,DEdeg_2,E_B_V_1,Class_2,
+ clean_2,zbest_2,f_zbest_2,zspec_2,z1med_2,z1min_2,z1max_2,
+ z1area_2,z2med_2,z2min_2,z2max_2,z2area_2,nfiltEazy_2,
+ nfiltAtl_2,nfiltAna_2,chirBest_2,chirStel_1,StelType_2,AGN_2,
+ optAGN_2,IRAGN_2,XrayAGN_2,v2RXSID_1,XMMSL2ID_1,aptmg_1,
+ e_aptmg_1,aptmr_1,e_aptmr_1,aptmz_2,e_aptmz_2,zmodel_2,
+ chibest_2,Massmed_2,Massl68_2,Massu68_2,Nfilts_2,umag_2,
+ gmag_1,rmag_2,zmag_2,grcsmag,rrcsmag,ircsmag,zrcsmag,Jmag_2,
+ Kmag_2,IRAC1sv,IRAC2sv,IRAC1sw_1,IRAC2sw_1,IRAC3sw_1,IRAC4sw_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end