Conversion of standardized ReadMe file for
file /./ftp/cats/I/356 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-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 ! Unique source identifier (unique within a
* particular Data Release) (source_id)
integer*8 SolID ! Solution Identifier (solution_id)
character*15 Class ! Name of best class (in vcclassd.dat table)
* (vari_best_class_name)
real*4 ClassSc ! ? Score of the best class
* (vari_best_class_score)
real*4 PGal ! ? Probability from DSC-Combmod of being a
* galaxy (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_galaxy)
real*4 PQSO ! ? Probability from DSC-Combmod of being a
* quasar (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_quasar)
character*14 ClassDSCC ! Class assigned by DSC based on the
* probability from its Combmod classifier
* (classlabel_dsc)
character*12 ClassDCSSA ! Class assigned by DSC based on the
* probability from its Specmod and Allosmod
* classifiers (classlabel_dsc_joint)
character*11 ClassOA ! Class assigned by OA the neuron that
* represents the source (classlabel_oa)
real*4 z ! ? Redshift from UGC (redshift_ugc)
real*4 zlow ! ? Redshift prediction lower limit from UGC
* (redshift_ugc_lower)
real*8 zupp ! ? Redshift prediction upper limit from UGC
* (redshift_ugc_upper)
integer*4 Ntrans ! Number of transits used for the
* morphological analysis (n_transits)
real*8 PAS ! (deg) ? Fitted position angle of the source for
* the Sersic Profile (posangle_sersic)
real*4 e_PAS ! (deg) ? Error on the fitted position angle of the
* source for the Sersic Profile
* (posangle_sersic_error)
real*8 IntS ! (e-/s) ? Fitted intensity of the source for the
* Sersic Profile (intensity_sersic)
real*4 e_IntS ! (e-/s) ? Error on the fitted intensity of the
* source at effective radius radiusSersic
* (intensity_sersic_error)
real*8 RadS ! (mas) ? Fitted effective radius of the source for
* the Sersic Profile (radius_sersic)
real*4 e_RadS ! (mas) ? Error on the fitted effective radius of
* the source for the Sersic Profile
* (radius_sersic_error)
real*4 ellS ! ? Fitted ellipticity of source for the
* Sersic Profile (ellipticity_sersic)
real*4 e_ellS ! ? Error on the fitted ellipticity of the
* source for the Sersic Profile
* (ellipticity_sersic_error)
real*8 nS ! ? Fitted Sersic Index for Sersic Profile
* (n_sersic)
real*4 e_nS ! ? Error on the fitted Sersic Index for
* Sersic Profile (n_sersic_error)
real*8 L2S ! ? L2 norm for the Sersic Profile
* (l2_sersic)
real*4 mpcvS0 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! [0/8]? Flag indicative of processing or
* scientific quality for the morphological
* parameters fitting for the Sersic Profile
* (flags_sersic) (4)
real*8 PAdV ! (deg) ? Fitted position angle of the source for
* the de Vaucouleurs Profile
* (posangle_de_vaucouleurs)
real*4 e_PAdV ! (deg) ? Error on the fitted position angle of the
* source for the de Vaucouleurs Profile
* (posangle_de_vaucouleurs_error)
real*8 IntdV ! (e-/s) ? Fitted intensity of the source for the
* de Vaucouleurs Profile
* (intensity_de_vaucouleurs)
real*4 e_IntdV ! (e-/s) ? Error on the fitted intensity of the
* bulge for the de Vaucouleurs Profile
* (intensity_de_vaucouleurs_error)
real*8 RaddV ! (mas) ? Fitted effective radius of the source for
* de Vaucouleurs Profile
* (radius_de_vaucouleurs)
real*4 e_RaddV ! (mas) ? Error on the fitted effective radius of
* the source for the de Vaucouleurs Profile
* (radius_de_vaucouleurs_error)
real*4 ElldV ! ? Fitted ellipticity of source for the
* de Vaucouleurs Profile
* (ellipticity_de_vaucouleurs)
real*4 e_ElldV ! ? Error on the fitted ellipticity of the
* source for the de Vaucouleurs Profile
* (ellipticity_de_vaucouleurs_error)
real*8 L2dV ! ? L2 norm for the de Vaucouleurs Profile
* (l2_de_vaucouleurs)
real*4 mpcvdV0 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! [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 ! 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 ! [0/1]? Catalog id (catalogue_id)
real*4 RAdeg ! (deg) Right ascension (ICRS) at Ep=2016.0 (ra)
real*4 DEdeg ! (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 ! Unique source identifier (unique within a
* particular Data Release) (source_id)
integer*8 SolID_1 ! Solution Identifier (solution_id)
integer*4 ASF ! [0/1] Flag indicating if the source is part
* of the astrometric selection
* (astrometric_selection_flag)
integer*4 GCS ! [0/1] Flag indicative of whether the source
* was used define the Gaia-CRF3
* (gaia_crf_source)
character*26 Class_1 ! Name of best class (in vcclassd.dat table)
* (vari_best_class_name)
real*4 ClassSc_1 ! ? Score of the best class
* (vari_best_class_score)
real*4 fvarG ! ? Fractional variability in the G band
* (fractional_variability_g)
real*4 SFIndex ! ? Index of the first-order structure function
* in the G band (structure_function_index)
real*8 s_SFIndex ! ? Standard deviation of the index of the
* structure function
* (structure_function_index_scatter)
real*4 QSOVar ! ? Quasar variability metric in the G band
* (qso_variability)
real*4 NonQSOVar ! ? Non-quasar variability metric in the G band
* (non_qso_variability)
real*4 VAGNMenSc ! ? Membership score (0=lowest,1=highest) of
* source to be of AGN type
* (vari_agn_membership_score)
real*4 PQSO_1 ! ? Probability from DSC-Combmod of being a
* quasar (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_quasar)
real*4 PGal_1 ! ? Probability from DSC-Combmod of being a
* galaxy (data used: BP/RP spectrum,
* photometry, astrometry)
* (classprob_dsc_combmod_galaxy)
character*14 ClassDSCC_1 ! Class assigned by DSC based on the
* probability from its Combmod classifier
* (classlabel_dsc)
character*12 ClassDCSSA_1 ! Class assigned by DSC based on the
* probability from its Specmod and Allosmod
* classifiers (classlabel_dsc_joint)
character*11 ClassOA_1 ! Class assigned by OA the neuron that
* represents the source (classlabel_oa)
real*8 z_1 ! ? Redshift from QSOC (redshift_qsoc)
real*8 zlow_1 ! ? Redshift lower confidence level from QSOC
* (redshift_qsoc_lower)
real*8 zupp_1 ! ? Redshift upper confidence level from QSOC
* (redshift_qsoc_upper)
real*8 CCF ! ? Value of the cross-correlation function
* used to derive the redshift from QSOC
* (ccfratio_qsoc)
real*4 QSOCSc ! ? Redshift zscore from QSOC (zscore_qsoc)
integer*4 flagsQSOC ! [0/31]? Processing flags for the analysis
* based on BP/RP Spectra from QSOC
* (flags_qsoc) (1)
integer*4 Ntrans_1 ! Number of transits used for the morphological
* analysis (n_transits)
real*8 IntQSO ! (e-/s) ? Fitted intensity of the quasar at its
* centre (intensity_quasar)
real*4 e_IntQSO ! (e-/s) ? Error on the fitted intensity of the quasar
* at its centre (intensity_quasar_error)
real*8 IntHost ! (e-/s) ? Fitted intensity of the host galaxy at the
* effective radius (intensity_hostgalaxy)
real*4 e_IntHost ! (e-/s) ? Error on the fitted intensity of the host
* galaxy at effective radius
* (intensity_hostgalaxy_error)
real*8 RadHost ! (mas) ? Fitted effective radius of the host galaxy
* (radius_hostgalaxy)
real*4 e_RadHost ! (mas) ? Error on the fitted effective radius of the
* host galaxy (radius_hostgalaxy_error)
real*8 n ! ? Fitted sersic Index (sersic_index)
real*4 e_n ! ? Error on the fitted sersic Index
* (sersic_index_error)
real*8 EllHost ! ? Fitted ellipticity of the host galaxy
* (ellipticity_hostgalaxy)
real*4 e_EllHost ! ? Error on the fitted ellipticity of the host
* galaxy (ellipticity_hostgalaxy_error)
real*8 PAHost ! (deg) ? Fitted position angle of the host galaxy
* (posangle_hostgalaxy)
real*4 e_PAHost ! (deg) ? Error on the fitted position angle of the
* host galaxy (posangle_hostgalaxy_error)
integer*4 HostDet ! [0/1] Flag indicating whether a host galaxy
* has been detected (host_galaxy_detected)
real*8 L2 ! ? L2 norm for the fitted Sersic profile
* (l2_norm)
real*4 mpcv0 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! ? 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 ! [0/6]? Flag indicative of processing or
* scientific quality for the morphological
* parameters fitting (host_galaxy_flag) (2)
integer*4 FlagSel_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 ! ? Catalog id (catalogue_id)
real*4 RAdeg_1 ! (deg) Right ascension (ICRS) at Ep=2016.0 (ra)
real*4 DEdeg_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,SolID,Class,ClassSc,PGal,PQSO,ClassDSCC,ClassDCSSA,
+ ClassOA,z,zlow,zupp,Ntrans,PAS,e_PAS,IntS,e_IntS,RadS,e_RadS,
+ ellS,e_ellS,nS,e_nS,L2S,mpcvS0,mpcvS1,mpcvS2,mpcvS3,mpcvS4,
+ mpcvS5,mpcvS6,mpcvS7,mpcvS8,mpcvS9,FlagsS,PAdV,e_PAdV,IntdV,
+ e_IntdV,RaddV,e_RaddV,ElldV,e_ElldV,L2dV,mpcvdV0,mpcvdV1,
+ mpcvdV2,mpcvdV3,mpcvdV4,mpcvdV5,FlagsdV,FlagSel,CatID,RAdeg,
+ DEdeg
if(ar__(57:69) .EQ. '') ClassSc = rNULL__
if(ar__(71:83) .EQ. '') PGal = rNULL__
if(ar__(85:98) .EQ. '') PQSO = rNULL__
if(ar__(140:152) .EQ. '') z = rNULL__
if(ar__(154:166) .EQ. '') zlow = rNULL__
if(ar__(168:178) .EQ. '') zupp = rNULL__
if(ar__(184:206) .EQ. '') PAS = rNULL__
if(ar__(208:228) .EQ. '') e_PAS = rNULL__
if(ar__(230:250) .EQ. '') IntS = rNULL__
if(ar__(252:273) .EQ. '') e_IntS = rNULL__
if(ar__(275:293) .EQ. '') RadS = rNULL__
if(ar__(295:316) .EQ. '') e_RadS = rNULL__
if(ar__(318:338) .EQ. '') ellS = rNULL__
if(ar__(340:361) .EQ. '') e_ellS = rNULL__
if(ar__(363:372) .EQ. '') nS = rNULL__
if(ar__(374:387) .EQ. '') e_nS = rNULL__
if(ar__(389:408) .EQ. '') L2S = rNULL__
if(ar__(410:431) .EQ. '') mpcvS0 = rNULL__
if(ar__(433:454) .EQ. '') mpcvS1 = rNULL__
if(ar__(456:477) .EQ. '') mpcvS2 = rNULL__
if(ar__(479:500) .EQ. '') mpcvS3 = rNULL__
if(ar__(502:523) .EQ. '') mpcvS4 = rNULL__
if(ar__(525:546) .EQ. '') mpcvS5 = rNULL__
if(ar__(548:569) .EQ. '') mpcvS6 = rNULL__
if(ar__(571:592) .EQ. '') mpcvS7 = rNULL__
if(ar__(594:615) .EQ. '') mpcvS8 = rNULL__
if(ar__(617:638) .EQ. '') mpcvS9 = rNULL__
if(ar__(640:640) .EQ. '') FlagsS = iNULL__
if(ar__(642:663) .EQ. '') PAdV = rNULL__
if(ar__(665:685) .EQ. '') e_PAdV = rNULL__
if(ar__(687:706) .EQ. '') IntdV = rNULL__
if(ar__(708:729) .EQ. '') e_IntdV = rNULL__
if(ar__(731:750) .EQ. '') RaddV = rNULL__
if(ar__(752:773) .EQ. '') e_RaddV = rNULL__
if(ar__(775:795) .EQ. '') ElldV = rNULL__
if(ar__(797:818) .EQ. '') e_ElldV = rNULL__
if(ar__(820:839) .EQ. '') L2dV = rNULL__
if(ar__(841:862) .EQ. '') mpcvdV0 = rNULL__
if(ar__(864:885) .EQ. '') mpcvdV1 = rNULL__
if(ar__(887:908) .EQ. '') mpcvdV2 = rNULL__
if(ar__(910:931) .EQ. '') mpcvdV3 = rNULL__
if(ar__(933:954) .EQ. '') mpcvdV4 = rNULL__
if(ar__(956:977) .EQ. '') mpcvdV5 = rNULL__
if(ar__(979:979) .EQ. '') FlagsdV = iNULL__
if(ar__(984:984) .EQ. '') CatID = iNULL__
c ..............Just test output...........
write(6,1)
+ Source,SolID,Class,ClassSc,PGal,PQSO,ClassDSCC,ClassDCSSA,
+ ClassOA,z,zlow,zupp,Ntrans,PAS,e_PAS,IntS,e_IntS,RadS,e_RadS,
+ ellS,e_ellS,nS,e_nS,L2S,mpcvS0,mpcvS1,mpcvS2,mpcvS3,mpcvS4,
+ mpcvS5,mpcvS6,mpcvS7,mpcvS8,mpcvS9,FlagsS,PAdV,e_PAdV,IntdV,
+ e_IntdV,RaddV,e_RaddV,ElldV,e_ElldV,L2dV,mpcvdV0,mpcvdV1,
+ mpcvdV2,mpcvdV3,mpcvdV4,mpcvdV5,FlagsdV,FlagSel,CatID,RAdeg,
+ DEdeg
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,SolID_1,ASF,GCS,Class_1,ClassSc_1,fvarG,SFIndex,
+ s_SFIndex,QSOVar,NonQSOVar,VAGNMenSc,PQSO_1,PGal_1,
+ ClassDSCC_1,ClassDCSSA_1,ClassOA_1,z_1,zlow_1,zupp_1,CCF,
+ QSOCSc,flagsQSOC,Ntrans_1,IntQSO,e_IntQSO,IntHost,e_IntHost,
+ RadHost,e_RadHost,n,e_n,EllHost,e_EllHost,PAHost,e_PAHost,
+ HostDet,L2,mpcv0,mpcv1,mpcv2,mpcv3,mpcv4,mpcv5,mpcv6,mpcv7,
+ mpcv8,mpcv9,mpcv10,mpcv11,mpcv12,mpcv13,mpcv14,FlagHost,
+ FlagSel_1,CatID_1,RAdeg_1,DEdeg_1
if(ar__1(72:84) .EQ. '') ClassSc_1 = rNULL__
if(ar__1(86:97) .EQ. '') fvarG = rNULL__
if(ar__1(99:111) .EQ. '') SFIndex = rNULL__
if(ar__1(113:132) .EQ. '') s_SFIndex = rNULL__
if(ar__1(134:147) .EQ. '') QSOVar = rNULL__
if(ar__1(149:162) .EQ. '') NonQSOVar = rNULL__
if(ar__1(164:185) .EQ. '') VAGNMenSc = rNULL__
if(ar__1(187:200) .EQ. '') PQSO_1 = rNULL__
if(ar__1(202:215) .EQ. '') PGal_1 = rNULL__
if(ar__1(257:267) .EQ. '') z_1 = rNULL__
if(ar__1(269:280) .EQ. '') zlow_1 = rNULL__
if(ar__1(282:292) .EQ. '') zupp_1 = rNULL__
if(ar__1(294:304) .EQ. '') CCF = rNULL__
if(ar__1(306:318) .EQ. '') QSOCSc = rNULL__
if(ar__1(320:321) .EQ. '') flagsQSOC = iNULL__
if(ar__1(327:347) .EQ. '') IntQSO = rNULL__
if(ar__1(349:370) .EQ. '') e_IntQSO = rNULL__
if(ar__1(372:392) .EQ. '') IntHost = rNULL__
if(ar__1(394:415) .EQ. '') e_IntHost = rNULL__
if(ar__1(417:435) .EQ. '') RadHost = rNULL__
if(ar__1(437:457) .EQ. '') e_RadHost = rNULL__
if(ar__1(459:468) .EQ. '') n = rNULL__
if(ar__1(470:482) .EQ. '') e_n = rNULL__
if(ar__1(484:503) .EQ. '') EllHost = rNULL__
if(ar__1(505:525) .EQ. '') e_EllHost = rNULL__
if(ar__1(527:548) .EQ. '') PAHost = rNULL__
if(ar__1(550:570) .EQ. '') e_PAHost = rNULL__
if(ar__1(574:592) .EQ. '') L2 = rNULL__
if(ar__1(594:615) .EQ. '') mpcv0 = rNULL__
if(ar__1(617:638) .EQ. '') mpcv1 = rNULL__
if(ar__1(640:661) .EQ. '') mpcv2 = rNULL__
if(ar__1(663:684) .EQ. '') mpcv3 = rNULL__
if(ar__1(686:707) .EQ. '') mpcv4 = rNULL__
if(ar__1(709:730) .EQ. '') mpcv5 = rNULL__
if(ar__1(732:753) .EQ. '') mpcv6 = rNULL__
if(ar__1(755:776) .EQ. '') mpcv7 = rNULL__
if(ar__1(778:799) .EQ. '') mpcv8 = rNULL__
if(ar__1(801:822) .EQ. '') mpcv9 = rNULL__
if(ar__1(824:845) .EQ. '') mpcv10 = rNULL__
if(ar__1(847:868) .EQ. '') mpcv11 = rNULL__
if(ar__1(870:891) .EQ. '') mpcv12 = rNULL__
if(ar__1(893:914) .EQ. '') mpcv13 = rNULL__
if(ar__1(916:937) .EQ. '') mpcv14 = rNULL__
if(ar__1(939:939) .EQ. '') FlagHost = iNULL__
if(ar__1(945:947) .EQ. '') CatID_1 = iNULL__
c ..............Just test output...........
write(6,2)
+ Source_1,SolID_1,ASF,GCS,Class_1,ClassSc_1,fvarG,SFIndex,
+ s_SFIndex,QSOVar,NonQSOVar,VAGNMenSc,PQSO_1,PGal_1,
+ ClassDSCC_1,ClassDCSSA_1,ClassOA_1,z_1,zlow_1,zupp_1,CCF,
+ QSOCSc,flagsQSOC,Ntrans_1,IntQSO,e_IntQSO,IntHost,e_IntHost,
+ RadHost,e_RadHost,n,e_n,EllHost,e_EllHost,PAHost,e_PAHost,
+ HostDet,L2,mpcv0,mpcv1,mpcv2,mpcv3,mpcv4,mpcv5,mpcv6,mpcv7,
+ mpcv8,mpcv9,mpcv10,mpcv11,mpcv12,mpcv13,mpcv14,FlagHost,
+ FlagSel_1,CatID_1,RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end