Conversion of standardized ReadMe file for
file /./ftp/cats/I/356 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-Sep-10
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. I/356 Gaia DR3 Part 2. Extra-galactic (Gaia Collaboration, 2022)
*================================================================================
*Gaia Data Release 3 (Gaia DR3) Part 2: Extra-galactic Gaia sources.
* Gaia collaboration
* <Astron. Astrophys., 674, A41 (2023)>
* =2023A&A...674A..41G
* =2022yCat.1356....0G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'galcand.dat' ! Parameters derived from various modules
dedicated to the classification and
characterisation of sources considered as
galaxy candidates (Galaxy_candidates)
integer*4 nr__
parameter (nr__=4842342) ! Number of records
character*1032 ar__ ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=2016.0)
integer*8 Source (nr__) ! Unique source identifier (unique within a
* particular Data Release) (source_id)
integer*8 SolID (nr__) ! Solution Identifier (solution_id)
character*15 Class (nr__) ! Name of best class (in vcclassd.dat table)
* (vari_best_class_name)
real*4 ClassSc (nr__) ! ? Score of the best class
* (vari_best_class_score)
real*4 PGal (nr__) ! ? Probability from DSC-Combmod of being a
* galaxy (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_galaxy)
real*4 PQSO (nr__) ! ? Probability from DSC-Combmod of being a
* quasar (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_quasar)
character*14 ClassDSCC (nr__) ! Class assigned by DSC based on the
* probability from its Combmod classifier
* (classlabel_dsc)
character*12 ClassDCSSA (nr__) ! Class assigned by DSC based on the
* probability from its Specmod and Allosmod
* classifiers (classlabel_dsc_joint)
character*11 ClassOA (nr__) ! Class assigned by OA the neuron that
* represents the source (classlabel_oa)
real*4 z (nr__) ! ? Redshift from UGC (redshift_ugc)
real*4 zlow (nr__) ! ? Redshift prediction lower limit from UGC
* (redshift_ugc_lower)
real*8 zupp (nr__) ! ? Redshift prediction upper limit from UGC
* (redshift_ugc_upper)
integer*4 Ntrans (nr__) ! Number of transits used for the
* morphological analysis (n_transits)
real*8 PAS (nr__) ! (deg) ? Fitted position angle of the source for
* the Sersic Profile (posangle_sersic)
real*4 e_PAS (nr__) ! (deg) ? Error on the fitted position angle of the
* source for the Sersic Profile
* (posangle_sersic_error)
real*8 IntS (nr__) ! (e-/s) ? Fitted intensity of the source for the
* Sersic Profile (intensity_sersic)
real*4 e_IntS (nr__) ! (e-/s) ? Error on the fitted intensity of the
* source at effective radius radiusSersic
* (intensity_sersic_error)
real*8 RadS (nr__) ! (mas) ? Fitted effective radius of the source for
* the Sersic Profile (radius_sersic)
real*4 e_RadS (nr__) ! (mas) ? Error on the fitted effective radius of
* the source for the Sersic Profile
* (radius_sersic_error)
real*4 ellS (nr__) ! ? Fitted ellipticity of source for the
* Sersic Profile (ellipticity_sersic)
real*4 e_ellS (nr__) ! ? Error on the fitted ellipticity of the
* source for the Sersic Profile
* (ellipticity_sersic_error)
real*8 nS (nr__) ! ? Fitted Sersic Index for Sersic Profile
* (n_sersic)
real*4 e_nS (nr__) ! ? Error on the fitted Sersic Index for
* Sersic Profile (n_sersic_error)
real*8 L2S (nr__) ! ? L2 norm for the Sersic Profile
* (l2_sersic)
real*4 mpcvS0 (nr__) ! ? Parameter 0 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[0])
real*4 mpcvS1 (nr__) ! ? Parameter 1 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[1])
real*4 mpcvS2 (nr__) ! ? Parameter 2 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[2])
real*4 mpcvS3 (nr__) ! ? Parameter 3 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[3])
real*4 mpcvS4 (nr__) ! ? Parameter 4 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[4])
real*4 mpcvS5 (nr__) ! ? Parameter 5 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[5])
real*4 mpcvS6 (nr__) ! ? Parameter 6 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[6])
real*4 mpcvS7 (nr__) ! ? Parameter 7 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[7])
real*4 mpcvS8 (nr__) ! ? Parameter 8 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[8])
real*4 mpcvS9 (nr__) ! ? Parameter 9 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the Sersic Profile
* (morph_params_corr_vec_sersic[9])
integer*4 FlagsS (nr__) ! [0/8]? Flag indicative of processing or
* scientific quality for the morphological
* parameters fitting for the Sersic Profile
* (flags_sersic) (4)
real*8 PAdV (nr__) ! (deg) ? Fitted position angle of the source for
* the de Vaucouleurs Profile
* (posangle_de_vaucouleurs)
real*4 e_PAdV (nr__) ! (deg) ? Error on the fitted position angle of the
* source for the de Vaucouleurs Profile
* (posangle_de_vaucouleurs_error)
real*8 IntdV (nr__) ! (e-/s) ? Fitted intensity of the source for the
* de Vaucouleurs Profile
* (intensity_de_vaucouleurs)
real*4 e_IntdV (nr__) ! (e-/s) ? Error on the fitted intensity of the
* bulge for the de Vaucouleurs Profile
* (intensity_de_vaucouleurs_error)
real*8 RaddV (nr__) ! (mas) ? Fitted effective radius of the source for
* de Vaucouleurs Profile
* (radius_de_vaucouleurs)
real*4 e_RaddV (nr__) ! (mas) ? Error on the fitted effective radius of
* the source for the de Vaucouleurs Profile
* (radius_de_vaucouleurs_error)
real*4 ElldV (nr__) ! ? Fitted ellipticity of source for the
* de Vaucouleurs Profile
* (ellipticity_de_vaucouleurs)
real*4 e_ElldV (nr__) ! ? Error on the fitted ellipticity of the
* source for the de Vaucouleurs Profile
* (ellipticity_de_vaucouleurs_error)
real*8 L2dV (nr__) ! ? L2 norm for the de Vaucouleurs Profile
* (l2_de_vaucouleurs)
real*4 mpcvdV0 (nr__) ! ? Parameter 0 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the de Vaucouleurs
* Profile
* (morph_params_corr_vec_de_vaucouleurs[0])
real*4 mpcvdV1 (nr__) ! ? Parameter 1 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the de Vaucouleurs
* Profile
* (morph_params_corr_vec_de_vaucouleurs[1])
real*4 mpcvdV2 (nr__) ! ? Parameter 2 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the de Vaucouleurs
* Profile
* (morph_params_corr_vec_de_vaucouleurs[2])
real*4 mpcvdV3 (nr__) ! ? Parameter 3 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the de Vaucouleurs
* Profile
* (morph_params_corr_vec_de_vaucouleurs[3])
real*4 mpcvdV4 (nr__) ! ? Parameter 4 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the de Vaucouleurs
* Profile
* (morph_params_corr_vec_de_vaucouleurs[4])
real*4 mpcvdV5 (nr__) ! ? Parameter 5 of vector form of the upper
* triangle of the correlation matrix for the
* fitted parameters for the de Vaucouleurs
* Profile
* (morph_params_corr_vec_de_vaucouleurs[5])
integer*4 FlagsdV (nr__) ! [0/8]? Flag indicative of processing or
* scientific quality for the morphological
* parameters fitting for the de Vaucouleurs
* Profile (flags_de_vaucouleurs) (1)
integer*4 FlagSel (nr__) ! Bit indicative of whether the input data
* from a given module met the source list
* eligibility criteria for the source of
* interest (source_selection_flags) (2)
integer*4 CatID (nr__) ! [0/1]? Catalog id (catalogue_id)
real*4 RAdeg (nr__) ! (deg) Right ascension (ICRS) at Ep=2016.0 (ra)
real*4 DEdeg (nr__) ! (deg) Declination (ICRS) at Ep=2016.0 (dec)
*Note (1): flagsS as follows:
* 1 = Elliptical profile fitted, other source exists at less that 2.5arcsec,
* doubtful solution
* 2 = Circular profile fitted, other source exists at less that 2.5arcsec,
* doubtful solution
* 3 = Elliptical profile fitted, the position angle did not converge and one
* parameter or more is at the bound, Poor solution.
* 4 = Elliptical profile fitted, the position angle did not converge
* 5 = Elliptical profile fitted, one parameter or more is at the bound,
* Poor solution
* 6 = Elliptical profile fitted
* 7 = Circular profile fitted, one parameter or more is at the bound,
* Poor solution
* 8 = Circular profile fitted
*Note (2): The bit is coded as:
* bit 0 = The source meets the eligibility criteria for the output of the
* Surface brightness analysis module
* bit 1 = The source belongs to the {GalaxyCatalogueName} table.
* bit 2 = The source meets the eligibility criteria for the output of the
* classification module based on photometric lightcurves
* bit 3 = The source meets the eligibility criteria for the output of the DSC
* module
* bit 4 = The source meets the eligibility criteria for the redshift
* determined by the UGC module
* bit 5 = The source meets the eligibility criteria for the classification
* output of the UGC module (not yet applicable to DR3).
*Note (4): Sersic flags as follows:
* 1 = Elliptical profile fitted, other source exists at less that 2.5",
* doubtful solution
* 2 = Circular profile fitted, other source exists at less that 2.5",
* doubtful solution
* 3 = Elliptical profile fitted, the position angle did not converge and
* one parameter or more is at the bound, poor solution.
* 4 = Elliptical profile fitted, the position angle did not converge
* 5 = Elliptical profile fitted, one parameter or more is at the bound,
* poor solution
* 6 = Elliptical profile fitted
* 7 = Circular profile fitted, one parameter or more is at the bound,
* poor solution
* 8 = Circular profile fitted
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'qsocand.dat' ! Parameters derived from various modules
dedicated to the classification and
characterisation of sources considered as QSO
candidates (QSO_candidates)
integer*4 nr__1
parameter (nr__1=6649162) ! Number of records
character*995 ar__1 ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=2016.0)
integer*8 Source_1 (nr__1) ! Unique source identifier (unique within a
* particular Data Release) (source_id)
integer*8 SolID_1 (nr__1) ! Solution Identifier (solution_id)
integer*4 ASF (nr__1) ! [0/1] Flag indicating if the source is part
* of the astrometric selection
* (astrometric_selection_flag)
integer*4 GCS (nr__1) ! [0/1] Flag indicative of whether the source
* was used define the Gaia-CRF3
* (gaia_crf_source)
character*26 Class_1 (nr__1) ! Name of best class (in vcclassd.dat table)
* (vari_best_class_name)
real*4 ClassSc_1 (nr__1) ! ? Score of the best class
* (vari_best_class_score)
real*4 fvarG (nr__1) ! ? Fractional variability in the G band
* (fractional_variability_g)
real*4 SFIndex (nr__1) ! ? Index of the first-order structure function
* in the G band (structure_function_index)
real*8 s_SFIndex (nr__1) ! ? Standard deviation of the index of the
* structure function
* (structure_function_index_scatter)
real*4 QSOVar (nr__1) ! ? Quasar variability metric in the G band
* (qso_variability)
real*4 NonQSOVar (nr__1) ! ? Non-quasar variability metric in the G band
* (non_qso_variability)
real*4 VAGNMenSc (nr__1) ! ? Membership score (0=lowest,1=highest) of
* source to be of AGN type
* (vari_agn_membership_score)
real*4 PQSO_1 (nr__1) ! ? Probability from DSC-Combmod of being a
* quasar (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_quasar)
real*4 PGal_1 (nr__1) ! ? Probability from DSC-Combmod of being a
* galaxy (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_galaxy)
character*14 ClassDSCC_1(nr__1) ! Class assigned by DSC based on the
* probability from its Combmod classifier
* (classlabel_dsc)
character*12 ClassDCSSA_1(nr__1) ! Class assigned by DSC based on the
* probability from its Specmod and Allosmod
* classifiers (classlabel_dsc_joint)
character*11 ClassOA_1 (nr__1) ! Class assigned by OA the neuron that
* represents the source (classlabel_oa)
real*8 z_1 (nr__1) ! ? Redshift from QSOC (redshift_qsoc)
real*8 zlow_1 (nr__1) ! ? Redshift lower confidence level from QSOC
* (redshift_qsoc_lower)
real*8 zupp_1 (nr__1) ! ? Redshift upper confidence level from QSOC
* (redshift_qsoc_upper)
real*8 CCF (nr__1) ! ? Value of the cross-correlation function
* used to derive the redshift from QSOC
* (ccfratio_qsoc)
real*4 QSOCSc (nr__1) ! ? Redshift zscore from QSOC (zscore_qsoc)
integer*4 flagsQSOC (nr__1) ! [0/31]? Processing flags for the analysis
* based on BP/RP Spectra from QSOC
* (flags_qsoc) (1)
integer*4 Ntrans_1 (nr__1) ! Number of transits used for the morphological
* analysis (n_transits)
real*8 IntQSO (nr__1) ! (e-/s) ? Fitted intensity of the quasar at its
* centre (intensity_quasar)
real*4 e_IntQSO (nr__1) ! (e-/s) ? Error on the fitted intensity of the quasar
* at its centre (intensity_quasar_error)
real*8 IntHost (nr__1) ! (e-/s) ? Fitted intensity of the host galaxy at the
* effective radius (intensity_hostgalaxy)
real*4 e_IntHost (nr__1) ! (e-/s) ? Error on the fitted intensity of the host
* galaxy at effective radius
* (intensity_hostgalaxy_error)
real*8 RadHost (nr__1) ! (mas) ? Fitted effective radius of the host galaxy
* (radius_hostgalaxy)
real*4 e_RadHost (nr__1) ! (mas) ? Error on the fitted effective radius of the
* host galaxy (radius_hostgalaxy_error)
real*8 n (nr__1) ! ? Fitted sersic Index (sersic_index)
real*4 e_n (nr__1) ! ? Error on the fitted sersic Index
* (sersic_index_error)
real*8 EllHost (nr__1) ! ? Fitted ellipticity of the host galaxy
* (ellipticity_hostgalaxy)
real*4 e_EllHost (nr__1) ! ? Error on the fitted ellipticity of the host
* galaxy (ellipticity_hostgalaxy_error)
real*8 PAHost (nr__1) ! (deg) ? Fitted position angle of the host galaxy
* (posangle_hostgalaxy)
real*4 e_PAHost (nr__1) ! (deg) ? Error on the fitted position angle of the
* host galaxy (posangle_hostgalaxy_error)
integer*4 HostDet (nr__1) ! [0/1] Flag indicating whether a host galaxy
* has been detected (host_galaxy_detected)
real*8 L2 (nr__1) ! ? L2 norm for the fitted Sersic profile
* (l2_norm)
real*4 mpcv0 (nr__1) ! ? Parameter 0 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[0])
real*4 mpcv1 (nr__1) ! ? Parameter 1 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[1])
real*4 mpcv2 (nr__1) ! ? Parameter 2 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[2])
real*4 mpcv3 (nr__1) ! ? Parameter 3 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[3])
real*4 mpcv4 (nr__1) ! ? Parameter 4 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[4])
real*4 mpcv5 (nr__1) ! ? Parameter 5 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[5])
real*4 mpcv6 (nr__1) ! ? Parameter 6 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[6])
real*4 mpcv7 (nr__1) ! ? Parameter 7 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[7])
real*4 mpcv8 (nr__1) ! ? Parameter 8 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[8])
real*4 mpcv9 (nr__1) ! ? Parameter 9 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[9])
real*4 mpcv10 (nr__1) ! ? Parameter 10 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[10])
real*4 mpcv11 (nr__1) ! ? Parameter 11 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[11])
real*4 mpcv12 (nr__1) ! ? Parameter 12 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[12])
real*4 mpcv13 (nr__1) ! ? Parameter 13 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[13])
real*4 mpcv14 (nr__1) ! ? Parameter 14 of vector form of the upper
* triangle of the correlation matrix for the
* fitted morphological parameters
* (morph_params_corr_vec[14])
integer*4 FlagHost (nr__1) ! [0/6]? Flag indicative of processing or
* scientific quality for the morphological
* parameters fitting (host_galaxy_flag) (2)
integer*4 FlagSel_1 (nr__1) ! Bit indicative of whether the input data from
* a given module met the source list
* eligibility criteria for the source of
* interest (source_selection_flags) (3)
integer*4 CatID_1 (nr__1) ! ? Catalog id (catalogue_id)
real*4 RAdeg_1 (nr__1) ! (deg) Right ascension (ICRS) at Ep=2016.0 (ra)
real*4 DEdeg_1 (nr__1) ! (deg) Declination (ICRS) at Ep=2016.0 (dec)
*Note (1): FlagsQSOC (binary combination (binary OR)) as:
* 0 = Z_NOWARNING. The processing of this source raised no warning flag.
* 1 = Z_AMBIGUOUS. The cross correlation function has more than one maximum
* with {ccfratioQsoc}>0.85, meaning that at least two redshifts lead to
* a similar chi^2^ and the solution is ambiguous
* 2 = Z_LOWCCFRATIO. The selected redshift leads to a small value of the
* cross correlation function when compared to the maximum value of the
* cross correlation function. Equivalently, {ccfratioQsoc}<0.9.
* 4 = Z_LOWZSCORE. The selected redshift leads to a low {zscoreQsoc}, meaning
* that at least one emission line is either missing or is strongly
* damped. Equivalently, {zscoreQsoc}<0.9.
* 8 = Z_NOTOPTIMAL. We did not choose the redshift having the lowest chi^2^
* (i.e. {ccfratioQsoc}<1).
* 16 = Z_BADSPEC. Raised if one of the following conditions is met:
* the number of BP or RP spectral transits (N_BP_ and N_RP_ hereafter), is
* lower than 10 transits, or G>20.5mag, or G>19+0.03*(N_BP_-10)mag, or
* G>19+0.03*(N_RP_-10)mag.
* This allows the user to filter out uncertain predictions in a simple way
* without having to explicitly deal with the aforementioned formula,
* nor with N_BP_ and N_RP_.
*Note (2): FlagHost as follows:
* 1 = Host Galaxy measured with circular Sersic profile
* 2 = Host Galaxy measured with elliptical Sersic profile
* 3 = No host galaxy detected
* 4 = Poor solution measured with elliptical Sersic profile
* 5 = No convergence of fitting but host galaxy detected
* 6 = No convergence of fitting or doubtful solution due to presence of a
* secondary source at d<5 arcsec
*Note (3): The bit is coded as:
* bit 0 = The source belongs to the Gaia-CRF3 ({AgnCrossId} table)
* bit 1 = The source belongs to the {FrameRotatorSource} table with either
* {usedForReferenceFrameOrientation} or
* {usedForReferenceFrameSpin} set to True.
* bit 2 = The source meets the eligibility criteria for the output of the
* Surface brightness analysis module
* bit 3 = The source belongs to the {QsoCatalogueName} table.
* bit 4 = The source meets the eligibility criteria for the output of the
* classification module based on photometric lightcurves
* bit 5 = The source meets the eligibility criteria for the output of the SOS
* AGN module
* bit 6 = The source meets the eligibility criteria for the output of the DSC
* module
* bit 7 = The source meets the eligibility criteria for the redshift determined
* by the QSOC module
* bit 8 = The source meets the eligibility criteria for the classification
* output of the QSOC module (not yet applicable to DR3)
C=============================================================================
C Loading file 'galcand.dat' ! Parameters derived from various modules
* dedicated to the classification and
* characterisation of sources considered as
* galaxy candidates (Galaxy_candidates)
C Format for file interpretation
1 format(
+ I19,1X,I19,1X,A15,1X,E13.8,1X,E13.8,1X,E14.8,1X,A14,1X,A12,1X,
+ A11,1X,E13.7,1X,E13.7,1X,F11.9,1X,I3,1X,F23.19,1X,E21.16,1X,
+ F21.19,1X,E22.16,1X,F19.14,1X,E22.16,1X,E21.16,1X,E22.16,1X,
+ F10.8,1X,E14.8,1X,F20.17,1X,E22.16,1X,E22.16,1X,E22.16,1X,
+ E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,
+ E22.16,1X,I1,1X,F22.18,1X,E21.16,1X,F20.18,1X,E22.16,1X,
+ F20.14,1X,E22.16,1X,E21.16,1X,E22.16,1X,F20.17,1X,E22.16,1X,
+ E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,I1,1X,I2,1X,
+ I1,1X,E23.16,1X,E23.16)
C Effective file loading
open(unit=1,status='old',file=
+'galcand.dat')
write(6,*) '....Loading file: galcand.dat'
do i__=1,4842342
read(1,'(A1032)')ar__
read(ar__,1)
+ Source(i__),SolID(i__),Class(i__),ClassSc(i__),PGal(i__),
+ PQSO(i__),ClassDSCC(i__),ClassDCSSA(i__),ClassOA(i__),z(i__),
+ zlow(i__),zupp(i__),Ntrans(i__),PAS(i__),e_PAS(i__),IntS(i__),
+ e_IntS(i__),RadS(i__),e_RadS(i__),ellS(i__),e_ellS(i__),
+ nS(i__),e_nS(i__),L2S(i__),mpcvS0(i__),mpcvS1(i__),
+ mpcvS2(i__),mpcvS3(i__),mpcvS4(i__),mpcvS5(i__),mpcvS6(i__),
+ mpcvS7(i__),mpcvS8(i__),mpcvS9(i__),FlagsS(i__),PAdV(i__),
+ e_PAdV(i__),IntdV(i__),e_IntdV(i__),RaddV(i__),e_RaddV(i__),
+ ElldV(i__),e_ElldV(i__),L2dV(i__),mpcvdV0(i__),mpcvdV1(i__),
+ mpcvdV2(i__),mpcvdV3(i__),mpcvdV4(i__),mpcvdV5(i__),
+ FlagsdV(i__),FlagSel(i__),CatID(i__),RAdeg(i__),DEdeg(i__)
if(ar__(57:69) .EQ. '') ClassSc(i__) = rNULL__
if(ar__(71:83) .EQ. '') PGal(i__) = rNULL__
if(ar__(85:98) .EQ. '') PQSO(i__) = rNULL__
if(ar__(140:152) .EQ. '') z(i__) = rNULL__
if(ar__(154:166) .EQ. '') zlow(i__) = rNULL__
if(ar__(168:178) .EQ. '') zupp(i__) = rNULL__
if(ar__(184:206) .EQ. '') PAS(i__) = rNULL__
if(ar__(208:228) .EQ. '') e_PAS(i__) = rNULL__
if(ar__(230:250) .EQ. '') IntS(i__) = rNULL__
if(ar__(252:273) .EQ. '') e_IntS(i__) = rNULL__
if(ar__(275:293) .EQ. '') RadS(i__) = rNULL__
if(ar__(295:316) .EQ. '') e_RadS(i__) = rNULL__
if(ar__(318:338) .EQ. '') ellS(i__) = rNULL__
if(ar__(340:361) .EQ. '') e_ellS(i__) = rNULL__
if(ar__(363:372) .EQ. '') nS(i__) = rNULL__
if(ar__(374:387) .EQ. '') e_nS(i__) = rNULL__
if(ar__(389:408) .EQ. '') L2S(i__) = rNULL__
if(ar__(410:431) .EQ. '') mpcvS0(i__) = rNULL__
if(ar__(433:454) .EQ. '') mpcvS1(i__) = rNULL__
if(ar__(456:477) .EQ. '') mpcvS2(i__) = rNULL__
if(ar__(479:500) .EQ. '') mpcvS3(i__) = rNULL__
if(ar__(502:523) .EQ. '') mpcvS4(i__) = rNULL__
if(ar__(525:546) .EQ. '') mpcvS5(i__) = rNULL__
if(ar__(548:569) .EQ. '') mpcvS6(i__) = rNULL__
if(ar__(571:592) .EQ. '') mpcvS7(i__) = rNULL__
if(ar__(594:615) .EQ. '') mpcvS8(i__) = rNULL__
if(ar__(617:638) .EQ. '') mpcvS9(i__) = rNULL__
if(ar__(640:640) .EQ. '') FlagsS(i__) = iNULL__
if(ar__(642:663) .EQ. '') PAdV(i__) = rNULL__
if(ar__(665:685) .EQ. '') e_PAdV(i__) = rNULL__
if(ar__(687:706) .EQ. '') IntdV(i__) = rNULL__
if(ar__(708:729) .EQ. '') e_IntdV(i__) = rNULL__
if(ar__(731:750) .EQ. '') RaddV(i__) = rNULL__
if(ar__(752:773) .EQ. '') e_RaddV(i__) = rNULL__
if(ar__(775:795) .EQ. '') ElldV(i__) = rNULL__
if(ar__(797:818) .EQ. '') e_ElldV(i__) = rNULL__
if(ar__(820:839) .EQ. '') L2dV(i__) = rNULL__
if(ar__(841:862) .EQ. '') mpcvdV0(i__) = rNULL__
if(ar__(864:885) .EQ. '') mpcvdV1(i__) = rNULL__
if(ar__(887:908) .EQ. '') mpcvdV2(i__) = rNULL__
if(ar__(910:931) .EQ. '') mpcvdV3(i__) = rNULL__
if(ar__(933:954) .EQ. '') mpcvdV4(i__) = rNULL__
if(ar__(956:977) .EQ. '') mpcvdV5(i__) = rNULL__
if(ar__(979:979) .EQ. '') FlagsdV(i__) = iNULL__
if(ar__(984:984) .EQ. '') CatID(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Source(i__),SolID(i__),Class(i__),ClassSc(i__),PGal(i__),
+ PQSO(i__),ClassDSCC(i__),ClassDCSSA(i__),ClassOA(i__),z(i__),
+ zlow(i__),zupp(i__),Ntrans(i__),PAS(i__),e_PAS(i__),IntS(i__),
+ e_IntS(i__),RadS(i__),e_RadS(i__),ellS(i__),e_ellS(i__),
+ nS(i__),e_nS(i__),L2S(i__),mpcvS0(i__),mpcvS1(i__),
+ mpcvS2(i__),mpcvS3(i__),mpcvS4(i__),mpcvS5(i__),mpcvS6(i__),
+ mpcvS7(i__),mpcvS8(i__),mpcvS9(i__),FlagsS(i__),PAdV(i__),
+ e_PAdV(i__),IntdV(i__),e_IntdV(i__),RaddV(i__),e_RaddV(i__),
+ ElldV(i__),e_ElldV(i__),L2dV(i__),mpcvdV0(i__),mpcvdV1(i__),
+ mpcvdV2(i__),mpcvdV3(i__),mpcvdV4(i__),mpcvdV5(i__),
+ FlagsdV(i__),FlagSel(i__),CatID(i__),RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'qsocand.dat' ! Parameters derived from various modules
* dedicated to the classification and
* characterisation of sources considered as QSO
* candidates (QSO_candidates)
C Format for file interpretation
2 format(
+ I19,1X,I19,1X,I1,1X,I1,1X,A26,1X,E13.8,1X,E12.7,1X,E13.7,1X,
+ F20.18,1X,E14.8,1X,E14.8,1X,E22.16,1X,E14.8,1X,E14.8,1X,A14,
+ 1X,A12,1X,A11,1X,F11.9,1X,F12.9,1X,F11.9,1X,F11.9,1X,E13.8,1X,
+ I2,1X,I3,1X,F21.16,1X,E22.16,1X,F21.18,1X,E22.16,1X,F19.14,1X,
+ E21.16,1X,F10.8,1X,E13.8,1X,F20.18,1X,E21.16,1X,F22.18,1X,
+ E21.16,1X,I1,1X,F19.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,
+ E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,
+ E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,E22.16,1X,
+ I1,1X,I3,1X,I3,1X,E23.16,1X,E23.16)
C Effective file loading
open(unit=1,status='old',file=
+'qsocand.dat')
write(6,*) '....Loading file: qsocand.dat'
do i__=1,6649162
read(1,'(A995)')ar__1
read(ar__1,2)
+ Source_1(i__),SolID_1(i__),ASF(i__),GCS(i__),Class_1(i__),
+ ClassSc_1(i__),fvarG(i__),SFIndex(i__),s_SFIndex(i__),
+ QSOVar(i__),NonQSOVar(i__),VAGNMenSc(i__),PQSO_1(i__),
+ PGal_1(i__),ClassDSCC_1(i__),ClassDCSSA_1(i__),ClassOA_1(i__),
+ z_1(i__),zlow_1(i__),zupp_1(i__),CCF(i__),QSOCSc(i__),
+ flagsQSOC(i__),Ntrans_1(i__),IntQSO(i__),e_IntQSO(i__),
+ IntHost(i__),e_IntHost(i__),RadHost(i__),e_RadHost(i__),
+ n(i__),e_n(i__),EllHost(i__),e_EllHost(i__),PAHost(i__),
+ e_PAHost(i__),HostDet(i__),L2(i__),mpcv0(i__),mpcv1(i__),
+ mpcv2(i__),mpcv3(i__),mpcv4(i__),mpcv5(i__),mpcv6(i__),
+ mpcv7(i__),mpcv8(i__),mpcv9(i__),mpcv10(i__),mpcv11(i__),
+ mpcv12(i__),mpcv13(i__),mpcv14(i__),FlagHost(i__),
+ FlagSel_1(i__),CatID_1(i__),RAdeg_1(i__),DEdeg_1(i__)
if(ar__1(72:84) .EQ. '') ClassSc_1(i__) = rNULL__
if(ar__1(86:97) .EQ. '') fvarG(i__) = rNULL__
if(ar__1(99:111) .EQ. '') SFIndex(i__) = rNULL__
if(ar__1(113:132) .EQ. '') s_SFIndex(i__) = rNULL__
if(ar__1(134:147) .EQ. '') QSOVar(i__) = rNULL__
if(ar__1(149:162) .EQ. '') NonQSOVar(i__) = rNULL__
if(ar__1(164:185) .EQ. '') VAGNMenSc(i__) = rNULL__
if(ar__1(187:200) .EQ. '') PQSO_1(i__) = rNULL__
if(ar__1(202:215) .EQ. '') PGal_1(i__) = rNULL__
if(ar__1(257:267) .EQ. '') z_1(i__) = rNULL__
if(ar__1(269:280) .EQ. '') zlow_1(i__) = rNULL__
if(ar__1(282:292) .EQ. '') zupp_1(i__) = rNULL__
if(ar__1(294:304) .EQ. '') CCF(i__) = rNULL__
if(ar__1(306:318) .EQ. '') QSOCSc(i__) = rNULL__
if(ar__1(320:321) .EQ. '') flagsQSOC(i__) = iNULL__
if(ar__1(327:347) .EQ. '') IntQSO(i__) = rNULL__
if(ar__1(349:370) .EQ. '') e_IntQSO(i__) = rNULL__
if(ar__1(372:392) .EQ. '') IntHost(i__) = rNULL__
if(ar__1(394:415) .EQ. '') e_IntHost(i__) = rNULL__
if(ar__1(417:435) .EQ. '') RadHost(i__) = rNULL__
if(ar__1(437:457) .EQ. '') e_RadHost(i__) = rNULL__
if(ar__1(459:468) .EQ. '') n(i__) = rNULL__
if(ar__1(470:482) .EQ. '') e_n(i__) = rNULL__
if(ar__1(484:503) .EQ. '') EllHost(i__) = rNULL__
if(ar__1(505:525) .EQ. '') e_EllHost(i__) = rNULL__
if(ar__1(527:548) .EQ. '') PAHost(i__) = rNULL__
if(ar__1(550:570) .EQ. '') e_PAHost(i__) = rNULL__
if(ar__1(574:592) .EQ. '') L2(i__) = rNULL__
if(ar__1(594:615) .EQ. '') mpcv0(i__) = rNULL__
if(ar__1(617:638) .EQ. '') mpcv1(i__) = rNULL__
if(ar__1(640:661) .EQ. '') mpcv2(i__) = rNULL__
if(ar__1(663:684) .EQ. '') mpcv3(i__) = rNULL__
if(ar__1(686:707) .EQ. '') mpcv4(i__) = rNULL__
if(ar__1(709:730) .EQ. '') mpcv5(i__) = rNULL__
if(ar__1(732:753) .EQ. '') mpcv6(i__) = rNULL__
if(ar__1(755:776) .EQ. '') mpcv7(i__) = rNULL__
if(ar__1(778:799) .EQ. '') mpcv8(i__) = rNULL__
if(ar__1(801:822) .EQ. '') mpcv9(i__) = rNULL__
if(ar__1(824:845) .EQ. '') mpcv10(i__) = rNULL__
if(ar__1(847:868) .EQ. '') mpcv11(i__) = rNULL__
if(ar__1(870:891) .EQ. '') mpcv12(i__) = rNULL__
if(ar__1(893:914) .EQ. '') mpcv13(i__) = rNULL__
if(ar__1(916:937) .EQ. '') mpcv14(i__) = rNULL__
if(ar__1(939:939) .EQ. '') FlagHost(i__) = iNULL__
if(ar__1(945:947) .EQ. '') CatID_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Source_1(i__),SolID_1(i__),ASF(i__),GCS(i__),Class_1(i__),
+ ClassSc_1(i__),fvarG(i__),SFIndex(i__),s_SFIndex(i__),
+ QSOVar(i__),NonQSOVar(i__),VAGNMenSc(i__),PQSO_1(i__),
+ PGal_1(i__),ClassDSCC_1(i__),ClassDCSSA_1(i__),ClassOA_1(i__),
+ z_1(i__),zlow_1(i__),zupp_1(i__),CCF(i__),QSOCSc(i__),
+ flagsQSOC(i__),Ntrans_1(i__),IntQSO(i__),e_IntQSO(i__),
+ IntHost(i__),e_IntHost(i__),RadHost(i__),e_RadHost(i__),
+ n(i__),e_n(i__),EllHost(i__),e_EllHost(i__),PAHost(i__),
+ e_PAHost(i__),HostDet(i__),L2(i__),mpcv0(i__),mpcv1(i__),
+ mpcv2(i__),mpcv3(i__),mpcv4(i__),mpcv5(i__),mpcv6(i__),
+ mpcv7(i__),mpcv8(i__),mpcv9(i__),mpcv10(i__),mpcv11(i__),
+ mpcv12(i__),mpcv13(i__),mpcv14(i__),FlagHost(i__),
+ FlagSel_1(i__),CatID_1(i__),RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end