FORTRAN Generation
(/./ftp/cats/J/A_A/569/A124)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/569/A124 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-May-19
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/569/A124      GUViCS. Ultraviolet Source Catalogs          (Voyer+, 2014)
*================================================================================
*The GALEX Ultraviolet Virgo Cluster Survey (GUViCS).
*III. The Ultraviolet Source Catalogs.
*    Voyer E.N., Boselli A., Boissier S., Heinis S., Cortese L., Ferrarese L.,
*    Cote P., Cuillandre J.-C., Gwyn S.D.J., Peng E.W., Zhang H., Liu C.
*   <Astron. Astrophys. 569, A124 (2014)>
*   =2014A&A...569A.124V
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'uv_ves.dat'	! GUViCS UV Virgo Cluster extended sources (H2)

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VCC        (nr__) ! ? Virgo Cluster Catalog ID (Binggeli et al.
*                                    1985, 1993, J/AJ/90/1681; McDonald et al.
*                                    2011, VII/266)
      integer*4     NGC        (nr__) ! ? NGC identification
      character*14  ID         (nr__) ! Identification from another catalog (2)
      real*8        RAdeg      (nr__) ! (deg) Right ascension (J2000)
      real*8        DEdeg      (nr__) ! (deg) [-1/185] Declination (J2000)
      character*28  Ntile      (nr__) ! Deepest GALEX tile available in NUV filter
      character*28  Ftile      (nr__) ! Deepest GALEX tile available in FUV filter
      real*8        Nexp       (nr__) ! (s) ?=-999 Exposure time of NUV GALEX image
      real*8        Fexp       (nr__) ! (s) ?=-999 Exposure time of FUV GALEX image
      real*4        Bmag       (nr__) ! (mag) ?=-99 BT magnitude from Binggeli et al.
*                                        (J/AJ/90/1681)
      integer*4     HRV        (nr__) ! (km/s) Velocity from Binggeli et al. (J/AJ/90/1681)
      integer*4     a          (nr__) ! (arcsec) Semi-major axis of galaxy aperture (1)
      integer*4     b          (nr__) ! (arcsec) Semi-minor axis of galaxy aperture (1)
      real*4        PA         (nr__) ! (deg) [-90/360] Position angle of galaxy
      integer*4     TT         (nr__) ! ?=-99 Galaxy type from GOLDMine,
*                                        where available (3)
      real*8        Nflx       (nr__) ! (ct/s) NUV integrated flux
      real*4        e_Nflx     (nr__) ! (ct/s) Error on NUV integrated flux
      real*8        Fflx       (nr__) ! (ct/s) FUV integrated flux
      real*8        e_Fflx     (nr__) ! (ct/s) ?=-999 Error on FUV integrated flux
      real*4        NUVmag     (nr__) ! (mag) NUV integrated magnitude
      real*4        e_NUVmag   (nr__) ! (mag) ?=-999 Error on NUV integrated magnitude
      real*8        FUVmag     (nr__) ! (mag) ?=-999 FUV integrated magnitude
      real*8        e_FUVmag   (nr__) ! (mag) ?=-999 Error on FUV integrated magnitude
*Note (1): The galaxy aperture is defined in the NUV galaxy image,
*     but used for both NUV and FUV integrated photometry.
*Note (2): If the galaxy does not have a VCC or NGC ID, it is given the next
*     most relevant ID from the GOLDMine or NED databases, depending where
*     it was selected from.
*Note (3): Galaxy types are listed where available and are from the GOLDMine
*  database Binggeli et al. (1985, 1993, Cat. J/AJ/90/1681).
*  Possible values are -3 through 20:
*    -3 = dS0
*    -2 = dE/dS0
*    -1 = dE (d:E)
*     0 = E-E/S0
*     1 = S0
*     2 = S0a-S0/Sa
*     3 = Sa
*     4 = Sab
*     5 = Sb
*     6 = Sbc
*     7 = Sc (dSc)
*     8 = Scd
*     9 = Sd
*    10 = Sdm-Sd/Sm
*    11 = Sm
*    12 = Im (Im/S)
*    13 = pec
*    14 = S/BCD (dS/BCD dS0/BCD Sd/BCD)
*    15 = Sm/BCD
*    16 = Im/BCD
*    17 = BCD
*    18 = S (dS)
*    19 = dIm/dE
*    20 = ?

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

C  Declarations for 'uv_vps.dat'	! GUViCS UV Virgo Cluster point-like sources (I2)

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VPS        (nr__1) ! [1/1231331] ID number from this paper (VPS_ID)
      real*8        RAdeg_1    (nr__1) ! (deg) Right ascension at UV galaxy center (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination at UV galaxy center (J2000)
      character*35  NUVfld     (nr__1) ! Deepest NUV GALEX field of source
*                                  (NUV_OrgPipFld)
      character*35  FUVfld     (nr__1) ! Deepest FUV GALEX field of source
*                                  (FUV_OrgPipFld)
      real*8        Nexp_1     (nr__1) ! (s) ?=-99 Exposure time of NUVfld (NUV_Expt)
      real*8        Fexp_1     (nr__1) ! (s) ?=-99 Exposure time of FUVfld (FUV_Expt)
      integer*4     No         (nr__1) ! [0/10]?=-99 Number of shallower NUV fields
*                                        (NUV_GUVfldOvlps) (1)
      integer*4     Fo         (nr__1) ! [0/6]?=-99 Number of shallower FUV fields
*                                        (FUV_GUVfldOvlps) (1)
      integer*4     Npf        (nr__1) ! [0/1]?=-99 NUV photometry flag, 0=good, 1=bad
*                                       (NUV_GoodPhotFlag)
      integer*4     Fpf        (nr__1) ! [0/1]?=-99 FUV photometry flag, 0=good, 1=bad
*                                       (FUV_GoodPhotFlag)
      integer*4     Nsf        (nr__1) ! [0/3] NUV source flag (NUV_SrcRemFlag) (3)
      integer*4     Fsf        (nr__1) ! [0/3]?=-1 FUV source flag (FUV_SrcRemFlag) (3)
      real*8        E_B_V      (nr__1) ! (mag) ?=-99 Galactic reddening as E(B-V) (E_BV) (4)
      real*8        Nrad       (nr__1) ! (deg) ?=-99 Source distance from center of FOV
*                                       (NUV_FOV_RADIUS)
      real*8        Frad       (nr__1) ! (deg) ?=-99 Source distance from center of FOV
*                                       (FUV_FOV_RADIUS)
      real*8        Nflx_1     (nr__1) ! (ct/s) ?=-99 Calibrated flux from AUTO aperture
*                                       (NUV_FLUX)
      real*8        e_Nflx_1   (nr__1) ! (ct/s) ?=-99 Error on NUV flux (NUV_FLUXERR)
      real*8        NUVmag_1   (nr__1) ! (mag) ?=-99 Calibrated AB mag from AUTO aperture
*                                       (NUV_MAG)
      real*8        e_NUVmag_1 (nr__1) ! (mag) ?=-99 Error on NUV mag (NUV_MAGERR)
      real*8        Fflx_1     (nr__1) ! (ct/s) ?=-99 Calibrated flux from AUTO aperture
*                                       (FUV_FLUX)
      real*8        e_Fflx_1   (nr__1) ! (ct/s) ?=-99 Error on Fflx (FUV_FLUXERR)
      real*8        FUVmag_1   (nr__1) ! (mag) ?=-99 Calibrated AB mag from AUTO aperture
*                                       (FUV_MAG)
      real*8        e_FUVmag_1 (nr__1) ! (mag) ?=-99 Error on FUVmag (FUV_MAGERR)
      real*8        NS_N       (nr__1) ! ?=-99 signal-to-noise from NUV_FLUX value
*                                       (NUV_S2N)
      real*8        FS_N       (nr__1) ! ?=-99 signal-to-noise from FUV_FLUX value
*                                       (FUV_S2N)
      real*8        Fflx1      (nr__1) ! (ct/s) ?=-99 FUV flux derived using NUV position
*                                       (FUV_NCAT_FLUX)
      real*8        e_Fflx1    (nr__1) ! (ct/s) ?=-99 Error on Fflx1 (FUV_NCAT_FLUXERR)
      real*8        FUVmag1    (nr__1) ! (mag) ?=-99 calibrated AB mag from FFUV1
*                                       (FUV_NCAT_MAG)
      real*8        e_FUVmag1  (nr__1) ! (mag) ?=-99 Error on FUVmag1 (FUV_NCAT_MAGERR)
      real*8        FS_N1      (nr__1) ! ?=-99 Signal to noise for Fflx1 (FUV_NCAT_S2N)
      real*8        Nbg        (nr__1) ! (s-1/arcsec2) ?=-99 Background value in photons/s/arcsec^2^
*                                       (NUV_SKYBG) (5)
      real*8        Fbg        (nr__1) ! (s-1/arcsec2) ?=-99 Background value in photons/s/arcsec^2^
*                                       (FUV_SKYBG) (6)
      integer*4     Naf        (nr__1) ! NUV artifact flags (NUV_ARTIFACT) (7)
      integer*4     Faf        (nr__1) ! ?=-999 FUV artifact flags (FUV_ARTIFACT) (8)
      integer*4     Nnum       (nr__1) ! [0/27029] Running GALEX NUV object number
*                                           (NUV_NUMBER)
      real*8        Nr2        (nr__1) ! (pix) [-209/262] NUV half-light radius from AUTO
*                                             aperture (NUV_FLUX_RADIUS_2)
      real*8        Nr4        (nr__1) ! (pix) [-41748/62879] NUV 90%-light radius from AUTO
*                                                 aperture (NUV_FLUX_RADIUS_4)
      real*8        NrK        (nr__1) ! [3.5/9.8]?=-99 NUV Kron apertures in units
*                                       of Na or Nb (NUV_KRON_RADIUS)
      real*8        Nmu        (nr__1) ! [-2.4/9.2]?=-99 Peak surface brightness above
*                                                  backgrd (NUV_MU_MAX)
      real*8        Na         (nr__1) ! (pix) [0.7/148]?=-99 Profile RMS along NUV
*                                       major axis (NUV_A_IMAGE)
      real*8        Nb         (nr__1) ! (pix) [0.2/86]?=-99 Profile RMS along NUV
*                                       minor axis (NUV_B_IMAGE)
      real*8        Ntheta     (nr__1) ! (deg) [-90/90]?=-99 Position angle (NUV_THETA_IMAGE)
      real*8        e_Na       (nr__1) ! (pix) ?=-99 RMS position error along major axis
*                                       (NUV_ERRA_IMAGE)
      real*8        e_Nb       (nr__1) ! (pix) ?=-99 RMS position error along minor axis
*                                       (NUV_ERRB_IMAGE)
      real*8        e_Ntheta   (nr__1) ! (deg) [-90/90]?=-99 Error ellipse NUV position angle
*                                       (NUV_ERRTHETA_IMAGE)
      real*8        Nr         (nr__1) ! (deg) ?=-99 NUV FWHM assuming a gaussian core
*                                       (NUV_FWHM_WORLD)
      integer*4     NSf_1      (nr__1) ! [0/31] NUV extraction flags (NUV_FLAGS)
      real*8        Ns_g       (nr__1) ! [0/1] SExtractor star/galaxy classifier
*                                        (1=point source) (NUV_CLASS_STAR)
      integer*4     Fnum       (nr__1) ! [0/10720] Running GALEX FUV object number
*                                           (FUV_NUMBER)
      real*8        Fr2        (nr__1) ! (pix) [0.5/52]?=-99 Half-light radius from AUTO
*                                          aperture (FUV_FLUX_RADIUS_2)
      real*8        Fr4        (nr__1) ! (pix) 90% radius from AUTO aperture
*                                       (FUV_FLUX_RADIUS_4)
      real*8        FrK        (nr__1) ! [3.5/9.9]?=-99 FUV Kron apertures in units of
*                                       A or B (FUV_KRON_RADIUS)
      real*8        Fmu        (nr__1) ! [-2.2/10.4]?=-99 Peak surface brightness above
*                                       backgrd (FUV_MU_MAX)
      real*8        Fa         (nr__1) ! (pix) [0.5/65]?=-99 Profile RMS along FUV major axis
*                                       (FUV_A_IMAGE)
      real*8        Fb         (nr__1) ! (pix) [0.3/36]?=-99 Profile RMS along FUV minor axis
*                                       (FUV_B_IMAGE)
      real*8        Ftheta     (nr__1) ! (deg) [-90/90]?=-99 Error ellipse FUV osition angle
*                                         (FUV_THETA_IMAGE)
      real*8        e_Fa       (nr__1) ! (pix) ?=-99 RMS position error along major axis
*                                       (FUV_ERRA_IMAGE)
      real*8        e_Fb       (nr__1) ! (pix) ?=-99 RMS position error along minor axis
*                                       (FUV_ERRB_IMAGE)
      real*8        e_Ftheta   (nr__1) ! (deg) [-90/90]?=-99 Error ellipse FUV position angle
*                                       (FUV_ERRTHETA_IMAGE)
      real*8        Fr         (nr__1) ! (deg) ?=-99 FUV FWHM assuming a gaussian core
*                                       (FUV_FWHM_WORLD)
      integer*4     FSf_1      (nr__1) ! [0/31] FUV extraction flags (FUV_FLAGS)
      real*8        Fs_g       (nr__1) ! [0/1] SExtractor star/galaxy classifier
*                                        (1=point source) (FUV_CLASS_STAR)
      integer*4     m          (nr__1) ! Number of SDSS matches to source
*                                      (SDSS_NumMtchs)
      character*10  mUV        (nr__1) ! Multiple UV match to optical source
*                                       (SDSS_multUV) (9)
      integer*8     ObjID      (nr__1) ! ?=-9999 SDSS object identifier (SDSS_ObjID)
      real*8        Sep        (nr__1) ! (deg) ?=-99 Angular separation from match
*                                       (SDSS_MtchSep)
      integer*4     ph         (nr__1) ! [0/1]?=-9999 SDSS photometry clean=1,
*                                       not clean=0 (SDSS_phot_Clean)
      integer*4     cl         (nr__1) ! [0/6]?=-9999 SDSS Type classification of
*                                       object (SDSS_phot_ObjType) (G1)
      character*6   SDSS       (nr__1) ! SDSS survey name: sdss, boss, segue1 or segue2
*                                      (SDSS_survey)
      real*8        zsp        (nr__1) ! ?=-99 SDSS spectroscopic z (SDSS_specz)
      real*8        e_zsp      (nr__1) ! [-4/14]?=-99 Error on zsp (SDSS_specz_Err)
      real*8        zNN        (nr__1) ! ?=-99 SDSS photometric z (SDSS_photzNN) (11)
      real*8        e_zNN      (nr__1) ! ?=-99 Error on zNN (SDSS_photzNN_Err)
      real*8        zRF        (nr__1) ! ?=-99 SDSS photometric z (SDSS_photzRF) (12)
      real*8        e_zRF      (nr__1) ! ?=-99 Error on zRF (SDSS_photzRF_Err)
      character*2   VM         (nr__1) ! [VM -] Flag for Virgo Cluster members
*                                       (SDSS_VM) (13)
      character*2   Vbg        (nr__1) ! [BK -] Flag for Virgo background sources
*                                       (SDSS_BCKD) (14)
      character*4   Bcl        (nr__1) ! Flag indicates classification method
*                                       (SDSS_BCKDclass) (15)
      character*8   f_Bcl      (nr__1) ! Flag for classification breakdown
*                                       (SDSS_BCKDclassBrkdn) (16)
      character*1   st         (nr__1) ! [*- ] * if source is SDSS determined star
*                                       (SDSS_STAR)
      character*4   n_st       (nr__1) ! Flag indicates classification method
*                                       (SDSS_STARclass) (17)
      character*1   TBD        (nr__1) ! [?-] Flag for objects to be determined
*                                       (SDSS_TBD) (18)
      integer*4     NGVS       (nr__1) ! Number of NGVS matches to source
*                                       (NGVS_NumMtchs) (19)
      character*3   f_NGVS     (nr__1) ! [mUV -] multiple UV match to optical NGVS
*                                       source (NGVS_multUV) (19)
      character*1   n_NGVS     (nr__1) ! [*- ] * if source is NGVS determined star
*                                       (NGVS_STAR) (19)
      integer*4     NED        (nr__1) ! Number of NED matches to source
*                                       (NED_NumMtchs)
      character*3   f_NED      (nr__1) ! [mUV -] multiple UV match to optical
*                                        source (NED_multUV) (20)
      character*30  NED_ID     (nr__1) ! Primary Object identifier from NED
*                                       (NED_ObjID)
      real*8        NEDRAdeg   (nr__1) ! (deg) ?=-99 Right ascension of NED optical
*                                       source (NED_RAdeg)
      real*8        NEDDEdeg   (nr__1) ! (deg) ?=-99 Declination of NED optical source
*                                       (NED_DEdeg)
      real*8        NEDsep     (nr__1) ! (deg) ?=-99 Angular separation from match
*                                       (NED_MtchSep)
      real*8        NEDz       (nr__1) ! ?=-99 Spectroscopic z of the NED source
*                                       (NED_z)
      character*4   q_NEDz     (nr__1) ! The NED quality flag on NED_z (NED_zQF)
      character*2   VM1        (nr__1) ! [VM -] Flag for Virgo Cluster members
*                                       (NED_VM) (21)
      character*2   Vbg1       (nr__1) ! [BK -] Flag for Virgo background sources
*                                       (NED_BCKD) (22)
      character*1   st1        (nr__1) ! [*-] Flag for stellar sources (NED_STAR) (23)
*Note (1): Number of shallower NUV fields that also contain this source.
*Note (3): Flag as follows:
*   -1 = No image taken in this band
*    0 = Source is detected in this band
*    1 = Source never detected in this band. Null values for this band in the
*        catalog are from the original GALEX pipeline and set to -99.
*    2 = Source detected in this band but deeper GALEX image overlaps source
*        location. Null values for this band in the catalog are set to -99.
*    3 = Spurious source detection in this band. Section 3.2 describes the
*        exposure time dependent spurious detection limits applied based on
*        Morrissey et al. (2007ApJS..173..682M). Null values for this band in
*        the catalog are set to -99.
*Note (4): Galactic reddening expressed as E(B-V) from the Galactic extinctions
*    maps Schlegel et al. (1998ApJ...500..525S).
*Note (5): Background value in photons per second per square arcsecond at the
*   source position taken in GALEX pipeline -nd-skybg.fits file image.
*Note (6): Background value in photons per second per square arcsecond at the
*   source position taken in GALEX pipeline -fd-skybg.fits file image.
*Note (7): Logical OR of artifact flags for pixels within a 3x3 pixel box in
*    GALEX pipeline -nd-flags.fits image.
*Note (8): Logical OR of artifact flags for pixels within a 3x3 pixel box in
*    GALEX pipeline -fd-flags.fits image
*Note (9): Flag is set to 'mUV' if this UV sources shares a SDSS match with one
*     or more other UV sources.
*Note (11): Photometric redshift from the SDSS 'neural network' algorithm.
*Note (12): Photometric redshift from the SDSS 'random forest' algorithm.
*Note (13): Flag value is 'VM' if the source is a member of the Virgo Cluster
*    based on SDSS photometric z or SDSS spectroscopic z.
*Note (14): Flag value is 'BK' if the source is located in the background of the
*    Virgo Cluster based on its SDSS photometric or spectroscopic z.
*Note (15): If the sources is a background source, as indicated by BK, then this
*     flag indicates how this was determined:
*    PHOT = background source based on photometric z
*    SPEC = background source based on spectroscopic z
*    BOTH = background source based on photometric and spectroscopic z
*Note (16): Flag value is various combinations of 'PNN','PRF', and 'SPEC' to
*    provide further details on which photometric z and/or spectroscopic z was
*    used to classify the source as a background source.
*Note (17): If the source is a star, this flag value indicates how that was
*    determined in the SDSS data:
*    PHOT = star from photo point spread function
*    SPEC = star from spectra
*    BOTH = star from both PHOT and SPEC methods
*Note (18): Flag value is 'TBD' if UV sources is matched to an SDSS source, but
*    the sources is uncategorizable based on available SDSS optical parameters.
*Note (19): NGVS is the Next Generation Virgo cluster Survey.
*    The flag is set to 'mUV' if this UV sources shares a NGVS match with one
*    or more other UV sources.
*Note (20): Flag is set to 'mUV' if this UV sources shares a NED match with one
*    or more other UV sources.
*Note (21): Flag value is 'VM' if the source is a member of the Virgo Cluster
*    based on NED spectroscopic z.
*Note (22): Flag value is 'BK' if the source is located in the background of the
*    Virgo Cluster based on the NED spectroscopic z.
*Note (23): Flag is set to '*' if the object is a star based on visual
*    investigation of NED data.

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

C  Declarations for 'guv_sdss.dat'	! All GUViCS-SDSS matched point-like sources (J2)

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VPS_1      (nr__2) ! [1/1231331] ID number from this paper (VPS_ID)
      integer*8     ObjID_1    (nr__2) ! SDSS object identifier (SDSS_ObjID)
      integer*8     SpObjID    (nr__2) ! SDSS ID points to object spectrum
*                                  (SDSS_SpecObjID)
      real*8        RAdeg_2    (nr__2) ! (deg) Right ascension at UV center (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination at UV center (J2000)
      real*8        Sep_1      (nr__2) ! (deg) Angular separation to match
      real*8        SRAdeg     (nr__2) ! (deg) Right ascension SDSS source (SDSS_RAdeg) (1)
      real*8        SDEdeg     (nr__2) ! (deg) Declination SDSS source (SDSS_DEdeg) (1)
      integer*4     ph_1       (nr__2) ! [0/1] clean photometry clean=1, not clean=0
*                                       (SDSS_phot_Clean)
      integer*4     cl_1       (nr__2) ! [0/6] Type classification of object
*                                       (SDSS_phot_ObjType) (G1)
      real*8        umag       (nr__2) ! (mag) ?=-99 u' Petrosian magnitude (u_petroMag) (1)
      real*8        e_umag     (nr__2) ! (mag) ?=-99 Error on umag (u_petroMagErr)
      real*8        umag1      (nr__2) ! (mag) ?=-99 u' magnitude (u_modelMag) (2)
      real*8        e_umag1    (nr__2) ! (mag) ?=-99 Error on umag1 (u_modelMagErr)
      real*8        umag2      (nr__2) ! (mag) ?=-99 deVaucouleurs+exponential mag
*                                       (u_cModelMag) (3)
      real*8        e_umag2    (nr__2) ! (mag) ?=-99 Error on umag2 (u_cModelMagErr)
      real*8        ur         (nr__2) ! (arcsec) ?=-99 Petrosian radius in u' (u_petroRad)
      real*8        e_ur       (nr__2) ! (arcsec) ?=-99 Error on ur (u_petroRadErr)
      real*8        ur50       (nr__2) ! (arcsec) ?=-99 Radius with 50% of Petrosian flux
*                                      (u_petroR50)
      real*8        e_ur50     (nr__2) ! (arcsec) ?=-99 Error on ur50 (u_petroR50Err)
      real*8        ur90       (nr__2) ! (arcsec) ?=-99 Radius with 90% of Petrosian flux
*                                      (u_petroR90)
      real*8        e_ur90     (nr__2) ! (arcsec) ?=-99 Error on ur90 (u_petroR90Err)
      real*8        gmag       (nr__2) ! (mag) ?=-99 Petrosian magnitude (g_petroMag) (1)
      real*8        e_gmag     (nr__2) ! (mag) ?=-99 Error on gmag (g_petorMagErr)
      real*8        gmag1      (nr__2) ! (mag) ?=-99 g' magnitude (g_modelMag) (2)
      real*8        e_gmag1    (nr__2) ! (mag) ?=-99 Error on gmag1 (g_modelMagErr)
      real*8        gmag2      (nr__2) ! (mag) ?=-99 DeVaucouleurs+exponential mag
*                                      (g_cModelMag) (3)
      real*8        e_gmag2    (nr__2) ! (mag) ?=-99 Error on gmag2 (g_cModelMagErr)
      real*8        gr         (nr__2) ! (arcsec) ?=-99 Petrosian radius in g' (g_petroRad)
      real*8        e_gr       (nr__2) ! (arcsec) ?=-99 Error on gr (g_petroRadErr)
      real*8        gr50       (nr__2) ! (arcsec) ?=-99 Radius with 50% of Petrosian flux
*                                        (g_petroR50)
      real*8        e_gr50     (nr__2) ! (arcsec) ?=-99 Error on gr50 (g_petroR50Err)
      real*8        gr90       (nr__2) ! (arcsec) ?=-99 Radius with 90% of Petrosian flux
*                                      (g_petroR90)
      real*8        e_gr90     (nr__2) ! (arcsec) ?=-99 Error on gr90 (g_petroR90Err)
      real*8        rmag       (nr__2) ! (mag) ?=-99 Petrosian magnitude for r' (1)
*                                      (r_petroMag)
      real*8        e_rmag     (nr__2) ! (mag) ?=-99 Error on rmag (r_petroMagErr)
      real*8        rmag1      (nr__2) ! (mag) ?=-99 r' magnitude (r_modelMag) (2)
      real*8        e_rmag1    (nr__2) ! (mag) ?=-99 Error on rmag1 (r_modelMagErr)
      real*8        rmag2      (nr__2) ! (mag) ?=-99 DeVaucouleurs+exponential mag
*                                      (r_cModelMag) (3)
      real*8        e_rmag2    (nr__2) ! (mag) ?=-99 Error on rmag2 (r_cModelMagErr)
      real*8        rr         (nr__2) ! (arcsec) ?=-99 Petrosian radius in r' (r_petroRad)
      real*8        e_rr       (nr__2) ! (arcsec) ?=-99 Error on rr (r_petroRadErr)
      real*8        rr50       (nr__2) ! (arcsec) ?=-99 Radius with 50% of Petrosian flux
*                                      (r_petroR50)
      real*8        e_rr50     (nr__2) ! (arcsec) ?=-99 Error on rr50 (r_petroR50Err)
      real*8        rr90       (nr__2) ! (arcsec) ?=-99 Radius with 90% of Petrosian flux
*                                      (r_petroR90)
      real*8        e_rr90     (nr__2) ! (arcsec) ?=-99 Error on rr90 (r_petroR90Err)
      real*8        imag       (nr__2) ! (mag) ?=-99 Petrosian magnitude for i'
      real*8        e_imag     (nr__2) ! (mag) ?=-99 Error on imag (i_petroMagErr)
      real*8        imag1      (nr__2) ! (mag) ?=-99 i' magnitude (i_modelMag) (2)
      real*8        e_imag1    (nr__2) ! (mag) ?=-99 Error on imag1 (i_modelMagErr)
      real*8        imag2      (nr__2) ! (mag) ?=-99 DeVaucouleurs+exponential mag
*                                      (i_cModelMag) (3)
      real*8        e_imag2    (nr__2) ! (mag) ?=-99 Error on imag2 (i_cModelMagErr)
      real*8        ir         (nr__2) ! (arcsec) ?=-99 Petrosian radius in i' (i_petroRad)
      real*8        e_ir       (nr__2) ! (arcsec) ?=-99 Error on ir (i_petroRadErr)
      real*8        ir50       (nr__2) ! (arcsec) ?=-99 Radius with 50% of Petrosian flux
*                                      (i_petroR50)
      real*8        e_ir50     (nr__2) ! (arcsec) ?=-99 Error on ir50 (i_petroR50Err)
      real*8        ir90       (nr__2) ! (arcsec) ?=-99 Radius with 90% of Petrosian flux
*                                      (i_petroR90)
      real*8        e_ir90     (nr__2) ! (arcsec) ?=-99 Error on ir90 (i_petroR90Err)
      real*8        zmag       (nr__2) ! (mag) ?=-99 Petrosian magnitude for z' in
*                                      (z_petroMag) (1)
      real*8        e_zmag     (nr__2) ! (mag) ?=-99 Error on zmag (z_petroMagErr)
      real*8        zmag1      (nr__2) ! (mag) ?=-99 z' magnitude (z_modelMag) (2)
      real*8        e_zmag1    (nr__2) ! (mag) ?=-99 Error on zmag1 (z_modelMagErr)
      real*8        zmag2      (nr__2) ! (mag) ?=-99 DeVaucouleurs+exponential mag
*                                      (z_cModelMag) (3)
      real*8        e_zmag2    (nr__2) ! (mag) ?=-99 Error on zmag2 (z_cModelMagErr)
      real*8        zr         (nr__2) ! (arcsec) ?=-99 Petrosian radius in z' (z_petroRad)
      real*8        e_zr       (nr__2) ! (arcsec) ?=-99 Error on zr (z_petroRadErr)
      real*8        zr50       (nr__2) ! (arcsec) ?=-99 Radius with 50% of Petrosian flux
*                                      (z_petroR50)
      real*8        e_zr50     (nr__2) ! (arcsec) ?=-99 Error on zr50 (z_petroR50Err)
      real*8        zr90       (nr__2) ! (arcsec) ?=-99 Radius with 90% of Petrosian flux
*                                      (z_petroR90)
      real*8        e_zr90     (nr__2) ! (arcsec) ?=-99 Error on zr90 (z_petroR90Err)
      character*6   SDSS_1     (nr__2) ! SDSS survey name: sdss, boss, segue1 or segue2
*                                      (SDSS_survey)
      real*8        zsp_1      (nr__2) ! ?=-99 Best spectroscopic z (SDSS_specz)
      real*8        e_zsp_1    (nr__2) ! [-4/15]?=-99 Error on zsp (SDSS_specz_Err)
      character*6   spcl       (nr__2) ! Spectroscopic classification of SDSS source
*                                      (SDSS_specz_class) (6)
      integer*4     f_spcl     (nr__2) ! ?=-888 Bitmask of warning values, 0=good
*                                      (SDSS_specz_Warning)
      real*8        zsp1       (nr__2) ! ?=-99 Best SDSS spectroscopic z
*                                      (SDSS_specz_noQSO) (4)
      real*8        e_zsp1     (nr__2) ! [-4/14]?=-99 Error on zsp1
*                                      (SDSS_specz_Err_noQSO)
      integer*4     f_zsp1     (nr__2) ! ?=-888 Bitmask of warning values, 0=good
*                                      (SDSS_specz_Warning_noQSO)
      character*6   spcl1      (nr__2) ! Classification of SDSS source (6)
*                                      (SDSS_specz_class_noQSO)
      real*8        zNN_1      (nr__2) ! ?=-99 Photometric z (SDSS_photzNN) (7)
      real*8        e_zNN_1    (nr__2) ! ?=-99 Error on zNN (SDSS_photzNN_Err)
      real*8        zRF_1      (nr__2) ! ?=-99 Photometric z (SDSS_photzRF) (8)
      real*8        e_zRF_1    (nr__2) ! ?=-99 Error on zRF (SDSS_photzRF_Err)
*Note (1): These SDSS parameters or measurements are taken from the r' data
*     or made in the r' aperture.
*Note (2): Better of De Vaucouleurs/exponential magnitudes in the given band.
*Note (3): Linear combination of De Vaucouleurs+exponential magnitudes.
*Note (4): Parameter is calculated when excluding the QSO fit.
*Note (6): Classifications are based on features in the objects spectrum, and
*    can be one of the following classifications codes: GALAXY, QSO, or STAR.
*Note (7): Photometric redshift from the SDSS 'neural network' algorithm.
*Note (8): Photometric redshift from the SDSS 'random forest' algorithm.

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

C  Loading file 'uv_ves.dat'	! GUViCS UV Virgo Cluster extended sources (H2)

C  Format for file interpretation

    1 format(
     +  4X,I4,5X,I4,1X,A14,1X,F9.5,1X,F9.5,1X,A28,1X,A28,1X,F8.2,1X,
     +  F8.2,1X,F6.2,1X,I6,1X,I3,1X,I3,1X,F6.2,1X,I3,1X,F8.3,1X,F5.3,
     +  1X,F8.3,1X,F8.3,1X,F6.3,1X,F5.3,1X,F8.3,1X,F8.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'uv_ves.dat')
      write(6,*) '....Loading file: uv_ves.dat'
      do i__=1,1770
        read(1,'(A225)')ar__
        read(ar__,1)
     +  VCC(i__),NGC(i__),ID(i__),RAdeg(i__),DEdeg(i__),Ntile(i__),
     +  Ftile(i__),Nexp(i__),Fexp(i__),Bmag(i__),HRV(i__),a(i__),
     +  b(i__),PA(i__),TT(i__),Nflx(i__),e_Nflx(i__),Fflx(i__),
     +  e_Fflx(i__),NUVmag(i__),e_NUVmag(i__),FUVmag(i__),
     +  e_FUVmag(i__)
        if(ar__(5:8) .EQ. '') VCC(i__) = iNULL__
        if(ar__(14:17) .EQ. '') NGC(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  VCC(i__),NGC(i__),ID(i__),RAdeg(i__),DEdeg(i__),Ntile(i__),
     +  Ftile(i__),Nexp(i__),Fexp(i__),Bmag(i__),HRV(i__),a(i__),
     +  b(i__),PA(i__),TT(i__),Nflx(i__),e_Nflx(i__),Fflx(i__),
     +  e_Fflx(i__),NUVmag(i__),e_NUVmag(i__),FUVmag(i__),
     +  e_FUVmag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'uv_vps.dat'	! GUViCS UV Virgo Cluster point-like sources (I2)

C  Format for file interpretation

    2 format(
     +  I7,1X,F12.8,1X,F11.8,1X,A35,1X,A35,1X,F8.2,1X,F8.2,1X,I4,1X,
     +  I4,1X,I4,1X,I4,1X,I1,1X,I2,2X,F12.8,2X,F12.8,2X,F12.8,1X,
     +  F10.4,1X,F9.4,1X,F9.4,1X,F9.4,1X,F10.4,1X,F9.4,1X,F9.4,1X,
     +  F9.4,1X,F9.4,1X,F9.4,1X,F11.4,1X,F9.4,1X,F9.4,1X,F9.4,1X,
     +  F10.4,2X,F12.8,2X,F12.8,1X,I3,1X,I4,1X,I5,1X,F10.5,1X,F12.5,
     +  1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,
     +  F9.5,1X,F9.5,1X,I2,1X,F7.5,1X,I5,1X,F9.5,1X,F11.5,1X,F9.5,1X,
     +  F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,
     +  1X,I2,1X,F7.5,1X,I1,1X,A10,1X,I19,3X,F12.8,1X,I5,1X,I5,1X,A6,
     +  27X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,1X,A2,
     +  4X,A2,9X,A4,7X,A8,3X,A1,5X,A4,2X,A1,5X,I2,1X,A3,3X,A1,1X,I1,
     +  1X,A3,5X,A30,3X,F9.5,3X,F9.5,3X,F12.8,3X,F9.5,1X,A4,5X,A2,1X,
     +  A2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'uv_vps.dat')
      write(6,*) '....Loading file: uv_vps.dat'
      do i__=1,1230855
        read(1,'(A987)')ar__1
        read(ar__1,2)
     +  VPS(i__),RAdeg_1(i__),DEdeg_1(i__),NUVfld(i__),FUVfld(i__),
     +  Nexp_1(i__),Fexp_1(i__),No(i__),Fo(i__),Npf(i__),Fpf(i__),
     +  Nsf(i__),Fsf(i__),E_B_V(i__),Nrad(i__),Frad(i__),Nflx_1(i__),
     +  e_Nflx_1(i__),NUVmag_1(i__),e_NUVmag_1(i__),Fflx_1(i__),
     +  e_Fflx_1(i__),FUVmag_1(i__),e_FUVmag_1(i__),NS_N(i__),
     +  FS_N(i__),Fflx1(i__),e_Fflx1(i__),FUVmag1(i__),e_FUVmag1(i__),
     +  FS_N1(i__),Nbg(i__),Fbg(i__),Naf(i__),Faf(i__),Nnum(i__),
     +  Nr2(i__),Nr4(i__),NrK(i__),Nmu(i__),Na(i__),Nb(i__),
     +  Ntheta(i__),e_Na(i__),e_Nb(i__),e_Ntheta(i__),Nr(i__),
     +  NSf_1(i__),Ns_g(i__),Fnum(i__),Fr2(i__),Fr4(i__),FrK(i__),
     +  Fmu(i__),Fa(i__),Fb(i__),Ftheta(i__),e_Fa(i__),e_Fb(i__),
     +  e_Ftheta(i__),Fr(i__),FSf_1(i__),Fs_g(i__),m(i__),mUV(i__),
     +  ObjID(i__),Sep(i__),ph(i__),cl(i__),SDSS(i__),zsp(i__),
     +  e_zsp(i__),zNN(i__),e_zNN(i__),zRF(i__),e_zRF(i__),VM(i__),
     +  Vbg(i__),Bcl(i__),f_Bcl(i__),st(i__),n_st(i__),TBD(i__),
     +  NGVS(i__),f_NGVS(i__),n_NGVS(i__),NED(i__),f_NED(i__),
     +  NED_ID(i__),NEDRAdeg(i__),NEDDEdeg(i__),NEDsep(i__),NEDz(i__),
     +  q_NEDz(i__),VM1(i__),Vbg1(i__),st1(i__)
c    ..............Just test output...........
        write(6,2)
     +  VPS(i__),RAdeg_1(i__),DEdeg_1(i__),NUVfld(i__),FUVfld(i__),
     +  Nexp_1(i__),Fexp_1(i__),No(i__),Fo(i__),Npf(i__),Fpf(i__),
     +  Nsf(i__),Fsf(i__),E_B_V(i__),Nrad(i__),Frad(i__),Nflx_1(i__),
     +  e_Nflx_1(i__),NUVmag_1(i__),e_NUVmag_1(i__),Fflx_1(i__),
     +  e_Fflx_1(i__),FUVmag_1(i__),e_FUVmag_1(i__),NS_N(i__),
     +  FS_N(i__),Fflx1(i__),e_Fflx1(i__),FUVmag1(i__),e_FUVmag1(i__),
     +  FS_N1(i__),Nbg(i__),Fbg(i__),Naf(i__),Faf(i__),Nnum(i__),
     +  Nr2(i__),Nr4(i__),NrK(i__),Nmu(i__),Na(i__),Nb(i__),
     +  Ntheta(i__),e_Na(i__),e_Nb(i__),e_Ntheta(i__),Nr(i__),
     +  NSf_1(i__),Ns_g(i__),Fnum(i__),Fr2(i__),Fr4(i__),FrK(i__),
     +  Fmu(i__),Fa(i__),Fb(i__),Ftheta(i__),e_Fa(i__),e_Fb(i__),
     +  e_Ftheta(i__),Fr(i__),FSf_1(i__),Fs_g(i__),m(i__),mUV(i__),
     +  ObjID(i__),Sep(i__),ph(i__),cl(i__),SDSS(i__),zsp(i__),
     +  e_zsp(i__),zNN(i__),e_zNN(i__),zRF(i__),e_zRF(i__),VM(i__),
     +  Vbg(i__),Bcl(i__),f_Bcl(i__),st(i__),n_st(i__),TBD(i__),
     +  NGVS(i__),f_NGVS(i__),n_NGVS(i__),NED(i__),f_NED(i__),
     +  NED_ID(i__),NEDRAdeg(i__),NEDDEdeg(i__),NEDsep(i__),NEDz(i__),
     +  q_NEDz(i__),VM1(i__),Vbg1(i__),st1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'guv_sdss.dat'	! All GUViCS-SDSS matched point-like sources (J2)

C  Format for file interpretation

    3 format(
     +  I7,1X,I19,1X,I19,1X,F9.5,1X,F8.5,1X,F10.8,1X,F12.8,1X,F11.8,
     +  1X,I1,1X,I1,3X,F9.5,2X,F10.5,3X,F9.5,3X,F9.5,3X,F9.5,1X,F11.5,
     +  1X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F9.5,
     +  1X,F11.5,3X,F9.5,3X,F9.5,3X,F9.5,1X,F11.5,1X,F10.6,3X,F10.6,
     +  3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F9.5,1X,F11.5,3X,F9.5,
     +  3X,F9.5,3X,F9.5,1X,F11.5,1X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,
     +  3X,F10.6,3X,F10.6,3X,F9.5,1X,F11.5,3X,F9.5,3X,F9.5,3X,F9.5,1X,
     +  F11.5,1X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,
     +  3X,F9.5,2X,F10.5,3X,F9.5,3X,F9.5,3X,F9.5,1X,F11.5,1X,F10.6,3X,
     +  F10.6,3X,F10.6,3X,F10.6,3X,F10.6,3X,F10.6,1X,A6,27X,F10.6,3X,
     +  F10.6,1X,A6,26X,I4,3X,F10.6,3X,F10.6,2X,I4,1X,A6,27X,F10.6,3X,
     +  F10.6,3X,F10.6,3X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'guv_sdss.dat')
      write(6,*) '....Loading file: guv_sdss.dat'
      do i__=1,1056793
        read(1,'(A1055)')ar__2
        read(ar__2,3)
     +  VPS_1(i__),ObjID_1(i__),SpObjID(i__),RAdeg_2(i__),
     +  DEdeg_2(i__),Sep_1(i__),SRAdeg(i__),SDEdeg(i__),ph_1(i__),
     +  cl_1(i__),umag(i__),e_umag(i__),umag1(i__),e_umag1(i__),
     +  umag2(i__),e_umag2(i__),ur(i__),e_ur(i__),ur50(i__),
     +  e_ur50(i__),ur90(i__),e_ur90(i__),gmag(i__),e_gmag(i__),
     +  gmag1(i__),e_gmag1(i__),gmag2(i__),e_gmag2(i__),gr(i__),
     +  e_gr(i__),gr50(i__),e_gr50(i__),gr90(i__),e_gr90(i__),
     +  rmag(i__),e_rmag(i__),rmag1(i__),e_rmag1(i__),rmag2(i__),
     +  e_rmag2(i__),rr(i__),e_rr(i__),rr50(i__),e_rr50(i__),
     +  rr90(i__),e_rr90(i__),imag(i__),e_imag(i__),imag1(i__),
     +  e_imag1(i__),imag2(i__),e_imag2(i__),ir(i__),e_ir(i__),
     +  ir50(i__),e_ir50(i__),ir90(i__),e_ir90(i__),zmag(i__),
     +  e_zmag(i__),zmag1(i__),e_zmag1(i__),zmag2(i__),e_zmag2(i__),
     +  zr(i__),e_zr(i__),zr50(i__),e_zr50(i__),zr90(i__),e_zr90(i__),
     +  SDSS_1(i__),zsp_1(i__),e_zsp_1(i__),spcl(i__),f_spcl(i__),
     +  zsp1(i__),e_zsp1(i__),f_zsp1(i__),spcl1(i__),zNN_1(i__),
     +  e_zNN_1(i__),zRF_1(i__),e_zRF_1(i__)
c    ..............Just test output...........
        write(6,3)
     +  VPS_1(i__),ObjID_1(i__),SpObjID(i__),RAdeg_2(i__),
     +  DEdeg_2(i__),Sep_1(i__),SRAdeg(i__),SDEdeg(i__),ph_1(i__),
     +  cl_1(i__),umag(i__),e_umag(i__),umag1(i__),e_umag1(i__),
     +  umag2(i__),e_umag2(i__),ur(i__),e_ur(i__),ur50(i__),
     +  e_ur50(i__),ur90(i__),e_ur90(i__),gmag(i__),e_gmag(i__),
     +  gmag1(i__),e_gmag1(i__),gmag2(i__),e_gmag2(i__),gr(i__),
     +  e_gr(i__),gr50(i__),e_gr50(i__),gr90(i__),e_gr90(i__),
     +  rmag(i__),e_rmag(i__),rmag1(i__),e_rmag1(i__),rmag2(i__),
     +  e_rmag2(i__),rr(i__),e_rr(i__),rr50(i__),e_rr50(i__),
     +  rr90(i__),e_rr90(i__),imag(i__),e_imag(i__),imag1(i__),
     +  e_imag1(i__),imag2(i__),e_imag2(i__),ir(i__),e_ir(i__),
     +  ir50(i__),e_ir50(i__),ir90(i__),e_ir90(i__),zmag(i__),
     +  e_zmag(i__),zmag1(i__),e_zmag1(i__),zmag2(i__),e_zmag2(i__),
     +  zr(i__),e_zr(i__),zr50(i__),e_zr50(i__),zr90(i__),e_zr90(i__),
     +  SDSS_1(i__),zsp_1(i__),e_zsp_1(i__),spcl(i__),f_spcl(i__),
     +  zsp1(i__),e_zsp1(i__),f_zsp1(i__),spcl1(i__),zNN_1(i__),
     +  e_zNN_1(i__),zRF_1(i__),e_zRF_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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