FORTRAN Generation
(/./ftp/cats/II/286)

Conversion of standardized ReadMe file for file /./ftp/cats/II/286 into FORTRAN code for loading all data files into arrays.

Note that special values are assigned to unknown or unspecified numbers (also called NULL numbers); when necessary, the coordinate components making up the right ascension and declination are converted into floating-point numbers representing these angles in degrees.



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-20
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. II/286        VIRMOS deep imaging survey. VVDS-F02 catalog  (VIRMOS team+, 2008)
*================================================================================
*The VIRMOS deep imaging survey.
*II. CFH12K BVRI optical data for the 0226-04 deep field.
*III. ESO/WFI deep U-band imaging of the 0226-04 deep field.
*IV. Near-infrared observations.
*V. The K-band follow-up in the 0226-04 field.
*    McCracken H.J.,  Radovich M., Iovino A., Temporin S.,
*    Adami C., Arnaboldi M., Arnouts S., Bardelli S., Bertin E., Bolzonella M.,
*    Bondi M., Bongiorno A., Bottini D., Brinchmann J., Busarello G., Cappi A.,
*    Charlot S., Ciliegi P., Contini T., Cucciati O., Cuillandre J.-C.,
*    Dantel-Fort M., De La Torre S., Foucaud S., Franzetti P., Garilli B.,
*    Gavignaud I., Gregorini L., Guzzo L., Gwyn S., Ilbert O., Lamareille F.,
*    Le Brun V., Le Fevre O., MacCagni D., Marano B., Marinoni C., Massarotti M.,
*    Mathez G., Mazure A., Mellier Y., Meneux B., Merighi R., Merluzzi P.,
*    Paioro L., Paltani S., Pello R., Picat J.P., Pollo A., Pozzetti L.,
*    Ripepi V., Rizzo D., Saracco P., Scaramella R., Scodeggio M., Tresse L.,
*    Vergani D., Vettolani G., Walcher C.J., Zamorani G., Zanichelli A., Zucca E.
*   <Astron. Astrophys., 410, 17-32 (2003); 417, 51-60 (2004);
*    442, 423-436 (2005); 482, 81-95 (2008)>
*   =2003A&A...410...17M
*   =2004A&A...417...51R
*   =2005A&A...442..423I
*   =2008A&A...482...81T
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'morph.dat'	! Morphological catalog (in filter I)

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS       (nr__) ! [20035000/20472316] Sequential number
      real*8        RAdeg      (nr__) ! (deg) Right ascension in decimal degrees, J2000
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees, J2000
      real*8        Xi         (nr__) ! (pix) Image X position (X_IMAGE)
      real*8        Yi         (nr__) ! (pix) Image Y position (Y_IMAGE)
      real*4        Ai         (nr__) ! (pix) Image Semi-major axis of image (A_IMAGE)
      real*4        e_Ai       (nr__) ! (pix) rms uncertainty on Ai (ERR_A_IMAGE)
      real*4        Bi         (nr__) ! (pix) Image Semi-minor axis of image (B_IMAGE)
      real*4        e_Bi       (nr__) ! (pix) rms uncertainty on Bi (ERR_B_IMAGE)
      real*4        PAi        (nr__) ! (deg) Image position angle (THETA_IMAGE)
      real*4        ePAi       (nr__) ! (deg) rms uncertainty on PAi (ERR_THETA_IMAGE)
      real*8        X2i        (nr__) ! (pix2) Image X2 position (X2_IMAGE)
      real*4        e_X2i      (nr__) ! (pix2) rms uncertainty on X2i (ERR_X2_IMAGE)
      real*8        Y2i        (nr__) ! (pix2) Image Y2 position (Y2_IMAGE)
      real*4        e_Y2i      (nr__) ! (pix2) rms uncertainty on Y2i (ERR_Y2_IMAGE)
      real*8        XYi        (nr__) ! (pix2) Image XY position (XY_IMAGE)
      real*4        eXYi       (nr__) ! (pix2) rms uncertainty on XYi (ERR_XY_IMAGE)
      real*4        Aw         (nr__) ! (arcsec) World semi-major axis (A_WORLD)
      real*4        e_Aw       (nr__) ! (arcsec) rms uncertainty on Aw (ERR_A_WORLD)
      real*4        Bw         (nr__) ! (arcsec) World semi-minor axis (B_WORLD)
      real*4        e_Bw       (nr__) ! (arcsec) rms uncertainty on Bw (ERR_B_WORLD)
      real*4        PAw        (nr__) ! (deg) World position angle (THETA_J2000)
      real*4        ePAw       (nr__) ! (deg) rms uncertainty on PAw (ERR_THETA_WORLD)
      real*8        X2w        (nr__) ! (arcsec+2) World X2 position (X2_WORLD)
      real*4        e_X2w      (nr__) ! (arcsec+2) rms uncertainty on X2w (ERR_X2_WORLD)
      real*8        Y2w        (nr__) ! (arcsec+2) World X2 position (Y2_WORLD)
      real*4        e_Y2w      (nr__) ! (arcsec+2) rms uncertainty on Y2w (ERR_Y2_WORLD)
      real*8        XYw        (nr__) ! (arcsec+2) World XY position (XY_WORLD)
      real*4        eXYw       (nr__) ! (arcsec+2) rms uncertainty on XYw (ERR_XY_WORLD)
      real*4        S_G        (nr__) ! [0/1] Star/Galaxy class (1=star) (CLASS_STAR)
      integer*4     Flag       (nr__) ! Undetailed flags (18, 22 or 26)

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

C  Declarations for 'ucat.dat'	! *U-band catalog (Ueso and Uloiano bands)

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS_1     (nr__1) ! [20035000,20472316] Sequential number
      real*8        RAdeg_1    (nr__1) ! (deg) Rights ascension J2000 (ALPHA_J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination J2000 (DELTA_J2000)
      real*4        mauto      (nr__1) ! (mag) ?=99.999 Auto AB magnitude (MAG_AUTO)
      real*4        e_mauto    (nr__1) ! (mag) ?=99.999 rms uncertainty on mauto
*                                     (ERR_MAG_AUTO)
      real*4        mautoc     (nr__1) ! (mag) ?=99.999 Corrected auto AB magnitude
*                                     (MAG_AUTO_CORR)
      real*4        mag1ap     (nr__1) ! (mag) ?=99.999 Aperture AB magnitude 1 (MAG1_APER)
      real*4        e_mag1ap   (nr__1) ! (mag) ?=99.999 rms uncertainty on mag1ap
*                                    (ERR_MAG1_APER)
      real*4        mag2ap     (nr__1) ! (mag) ?=99.999 Aperture AB magnitude 2 (MAG2_APER)
      real*4        e_mag2ap   (nr__1) ! (mag) ?=99.999 rms uncertainty on mag2ap
*                                    (ERR_MAG2_APER)
      real*4        mautol     (nr__1) ! (mag) ?=99.999 Autoloiano magnitude (MAG_AUTOLOIANO)
      real*4        e_mautol   (nr__1) ! (mag) ?=99.999 rms uncertainty on mautol
*                                    (ERR_MAG_AUTOLOIANO)
      real*4        m1autol    (nr__1) ! (mag) ?=99.999 Auto Uloiano AB magnitude 1
*                                    (MAG1_APERLOIANO)
      real*4        e_m1autol  (nr__1) ! (mag) ?=99.999 rms uncertainty on m1autol
*                                    (ERR_MAG1_APERLOIANO)
      real*4        m2autol    (nr__1) ! (mag) ?=99.999 Autoloiano AB magnitude 2
*                                    (MAG2_APERLOIANO)
      real*4        e_m2autol  (nr__1) ! (mag) ?=99.999 rms uncertainty on m2autol
*                                    (ERR_MAG2_APERLOIANO)
      integer*4     isoA       (nr__1) ! (pix) ?=99999 Image isophotal area (ISOAREA_IMAGE)
      real*4        thr        (nr__1) ! ?=999.9 Threshold (THRESHOLD)
      real*4        muthr      (nr__1) ! (mag/arcsec2) ?=99.999 Surface brightness threshold
*                                     (MU_THRESHOLD)
      real*4        backg      (nr__1) ! ?=99.999 Backgroung flux (BACKGROUND)
      real*8        Fmax       (nr__1) ! ?=9999.999 Maximum flux (FLUX_MAX)
      real*4        muMax      (nr__1) ! (mag/arcsec2) ?=99.999 Maximum Surface brightness (MU_MAX)
      real*4        Krad       (nr__1) ! (arcsec) ?=99.99 Kron radius (KRON_RADIUS)
      real*8        F1rad      (nr__1) ! (arcsec) ?=9999.99 Flux1 radius (FLUX1_RADIUS)
      real*8        F2rad      (nr__1) ! (arcsec) ?=9999.99 Flux2 radius (FLUX2_RADIUS)
      real*8        F3rad      (nr__1) ! (arcsec) ?=9999.99 Flux3 radius (FLUX3_RADIUS)
      real*4        miso       (nr__1) ! (mag) ?=99.999 Isophotal AB magnitude (MAG_ISO)
      real*4        e_miso     (nr__1) ! (mag) ?=99.999 rms uncertainty on miso (ERR_MAG_ISO)
      real*8        Fiso       (nr__1) ! ?=9999999.9 Isophotal flux (FLUX_ISO)
      real*4        e_Fiso     (nr__1) ! ?=9999.9 rms uncertainty on Fiso
*                                         (ERR_FLUX_ISO)
      real*4        misoc      (nr__1) ! (mag) ?=99.999 Corrected isophotal AB magnitude
*                                        (MAG_ISOCOR)
      real*4        e_misoc    (nr__1) ! (mag) ?=99.999 rms uncertainty on misoc
*                                        (ERR_MAG_ISOCOR)
      real*8        Fisoc      (nr__1) ! ?=9999999.9 Corrected isophotal flux
*                                        (FLUX_ISOCOR)
      real*4        e_Fisoc    (nr__1) ! ?=9999.9 rms uncertainty on Fisoc
*                                        (ERR_FLUX_ISOCOR)
      real*8        Fauto      (nr__1) ! ?=9999999.9 Flux Auto (FLUX_AUTO)
      real*4        e_Fauto    (nr__1) ! ?=9999.9 rms uncertainty on Fauto
*                                        (ERR_FLUX_AUTO)
      real*4        mbest      (nr__1) ! (mag) ?=99.999 Best value of AB magnitude (MAG_BEST)
      real*4        e_mbest    (nr__1) ! (mag) ?=99.999 rms uncertainty on mbest
*                                        (ERR_MAG_BEST)
      real*8        Fbest      (nr__1) ! ?=9999999.9 Best value for flux (FLUX_BEST)
      real*4        e_Fbest    (nr__1) ! ?=9999.9 rms uncertainty on Fbest
*                                        (ERR_FLUX_BEST)
      real*8        F1ap       (nr__1) ! ?=9999999.9 Aperture flux 1 (FLUX1_APER)
      real*4        e_F1ap     (nr__1) ! ?=9999.9 rms uncertainty on F1ap
*                                        (ERR_FLUX1_APER)
      real*8        F2ap       (nr__1) ! ?=9999999.9 Aperture flux 2 (FLUX2_APER)
      real*4        e_F2ap     (nr__1) ! ?=9999.9 rms uncertainty on F2ap
*                                        (ERR_FLUX2_APER)
      integer*4     Sflg       (nr__1) ! ?=99 Sextrator flags (FLAGS)

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

C  Declarations for 'bcat.dat'	! B-band catalog (Johnson)

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS_2     (nr__2) ! [20035000,20472316] Sequential number
      real*8        RAdeg_2    (nr__2) ! (deg) Rights ascension J2000 (ALPHA_J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination J2000 (DELTA_J2000)
      real*4        mauto_1    (nr__2) ! (mag) ?=99.999 Auto AB magnitude (MAG_AUTO)
      real*4        e_mauto_1  (nr__2) ! (mag) ?=99.999 rms uncertainty on mauto
*                                     (ERR_MAG_AUTO)
      real*4        mautoc_1   (nr__2) ! (mag) ?=99.999 Corrected auto AB magnitude
*                                     (MAG_AUTO_CORR)
      real*4        mag1ap_1   (nr__2) ! (mag) ?=99.999 Aperture AB magnitude 1 (MAG1_APER)
      real*4        e_mag1ap_1 (nr__2) ! (mag) ?=99.999 rms uncertainty on mag1ap
*                                    (ERR_MAG1_APER)
      real*4        mag2ap_1   (nr__2) ! (mag) ?=99.999 Aperture AB magnitude 2 (MAG2_APER)
      real*4        e_mag2ap_1 (nr__2) ! (mag) ?=99.999 rms uncertainty on mag2ap
*                                    (ERR_MAG2_APER)
      integer*4     isoA_1     (nr__2) ! (pix) ?=99999 Image isophotal area (ISOAREA_IMAGE)
      real*4        thr_1      (nr__2) ! ?=999.9 Threshold (THRESHOLD)
      real*4        muthr_1    (nr__2) ! (mag/arcsec2) ?=99.999 Surface brightness threshold
*                                     (MU_THRESHOLD)
      real*4        backg_1    (nr__2) ! ?=99.999 Backgroung flux (BACKGROUND)
      real*8        Fmax_1     (nr__2) ! ?=9999.999 Maximum flux (FLUX_MAX)
      real*4        muMax_1    (nr__2) ! (mag/arcsec2) ?=99.999 Maximum Surface brightness (MU_MAX)
      real*4        Krad_1     (nr__2) ! (arcsec) ?=99.99 Kron radius (KRON_RADIUS)
      real*8        F1rad_1    (nr__2) ! (arcsec) ?=9999.99 Flux1 radius (FLUX1_RADIUS)
      real*8        F2rad_1    (nr__2) ! (arcsec) ?=9999.99 Flux2 radius (FLUX2_RADIUS)
      real*8        F3rad_1    (nr__2) ! (arcsec) ?=9999.99 Flux3 radius (FLUX3_RADIUS)
      real*4        miso_1     (nr__2) ! (mag) ?=99.999 Isophotal AB magnitude (MAG_ISO)
      real*4        e_miso_1   (nr__2) ! (mag) ?=99.999 rms uncertainty on miso (ERR_MAG_ISO)
      real*8        Fiso_1     (nr__2) ! ?=9999999.9 Isophotal flux (FLUX_ISO)
      real*4        e_Fiso_1   (nr__2) ! ?=9999.9 rms uncertainty on Fiso
*                                         (ERR_FLUX_ISO)
      real*4        misoc_1    (nr__2) ! (mag) ?=99.999 Corrected isophotal AB magnitude
*                                        (MAG_ISOCOR)
      real*4        e_misoc_1  (nr__2) ! (mag) ?=99.999 rms uncertainty on misoc
*                                        (ERR_MAG_ISOCOR)
      real*8        Fisoc_1    (nr__2) ! ?=9999999.9 Corrected isophotal flux
*                                        (FLUX_ISOCOR)
      real*4        e_Fisoc_1  (nr__2) ! ?=9999.9 rms uncertainty on Fisoc
*                                        (ERR_FLUX_ISOCOR)
      real*8        Fauto_1    (nr__2) ! ?=9999999.9 Flux Auto (FLUX_AUTO)
      real*4        e_Fauto_1  (nr__2) ! ?=9999.9 rms uncertainty on Fauto
*                                        (ERR_FLUX_AUTO)
      real*4        mbest_1    (nr__2) ! (mag) ?=99.999 Best value of AB magnitude (MAG_BEST)
      real*4        e_mbest_1  (nr__2) ! (mag) ?=99.999 rms uncertainty on mbest
*                                        (ERR_MAG_BEST)
      real*8        Fbest_1    (nr__2) ! ?=9999999.9 Best value for flux (FLUX_BEST)
      real*4        e_Fbest_1  (nr__2) ! ?=9999.9 rms uncertainty on Fbest
*                                        (ERR_FLUX_BEST)
      real*8        F1ap_1     (nr__2) ! ?=9999999.9 Aperture flux 1 (FLUX1_APER)
      real*4        e_F1ap_1   (nr__2) ! ?=9999.9 rms uncertainty on F1ap
*                                        (ERR_FLUX1_APER)
      real*8        F2ap_1     (nr__2) ! ?=9999999.9 Aperture flux 2 (FLUX2_APER)
      real*4        e_F2ap_1   (nr__2) ! ?=9999.9 rms uncertainty on F2ap
*                                        (ERR_FLUX2_APER)
      integer*4     Sflg_1     (nr__2) ! ?=99 Sextrator flags (FLAGS)

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

C  Declarations for 'vcat.dat'	! V-band catalog (Johnson)

      integer*4 nr__3
      parameter (nr__3=253588)	! Number of records
      character*297 ar__3  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS_3     (nr__3) ! [20035000,20472316] Sequential number
      real*8        RAdeg_3    (nr__3) ! (deg) Rights ascension J2000 (ALPHA_J2000)
      real*8        DEdeg_3    (nr__3) ! (deg) Declination J2000 (DELTA_J2000)
      real*4        mauto_2    (nr__3) ! (mag) ?=99.999 Auto AB magnitude (MAG_AUTO)
      real*4        e_mauto_2  (nr__3) ! (mag) ?=99.999 rms uncertainty on mauto
*                                     (ERR_MAG_AUTO)
      real*4        mautoc_2   (nr__3) ! (mag) ?=99.999 Corrected auto AB magnitude
*                                     (MAG_AUTO_CORR)
      real*4        mag1ap_2   (nr__3) ! (mag) ?=99.999 Aperture AB magnitude 1 (MAG1_APER)
      real*4        e_mag1ap_2 (nr__3) ! (mag) ?=99.999 rms uncertainty on mag1ap
*                                    (ERR_MAG1_APER)
      real*4        mag2ap_2   (nr__3) ! (mag) ?=99.999 Aperture AB magnitude 2 (MAG2_APER)
      real*4        e_mag2ap_2 (nr__3) ! (mag) ?=99.999 rms uncertainty on mag2ap
*                                    (ERR_MAG2_APER)
      integer*4     isoA_2     (nr__3) ! (pix) ?=99999 Image isophotal area (ISOAREA_IMAGE)
      real*4        thr_2      (nr__3) ! ?=999.9 Threshold (THRESHOLD)
      real*4        muthr_2    (nr__3) ! (mag/arcsec2) ?=99.999 Surface brightness threshold
*                                     (MU_THRESHOLD)
      real*4        backg_2    (nr__3) ! ?=99.999 Backgroung flux (BACKGROUND)
      real*8        Fmax_2     (nr__3) ! ?=9999.999 Maximum flux (FLUX_MAX)
      real*4        muMax_2    (nr__3) ! (mag/arcsec2) ?=99.999 Maximum Surface brightness (MU_MAX)
      real*4        Krad_2     (nr__3) ! (arcsec) ?=99.99 Kron radius (KRON_RADIUS)
      real*8        F1rad_2    (nr__3) ! (arcsec) ?=9999.99 Flux1 radius (FLUX1_RADIUS)
      real*8        F2rad_2    (nr__3) ! (arcsec) ?=9999.99 Flux2 radius (FLUX2_RADIUS)
      real*8        F3rad_2    (nr__3) ! (arcsec) ?=9999.99 Flux3 radius (FLUX3_RADIUS)
      real*4        miso_2     (nr__3) ! (mag) ?=99.999 Isophotal AB magnitude (MAG_ISO)
      real*4        e_miso_2   (nr__3) ! (mag) ?=99.999 rms uncertainty on miso (ERR_MAG_ISO)
      real*8        Fiso_2     (nr__3) ! ?=9999999.9 Isophotal flux (FLUX_ISO)
      real*4        e_Fiso_2   (nr__3) ! ?=9999.9 rms uncertainty on Fiso
*                                         (ERR_FLUX_ISO)
      real*4        misoc_2    (nr__3) ! (mag) ?=99.999 Corrected isophotal AB magnitude
*                                        (MAG_ISOCOR)
      real*4        e_misoc_2  (nr__3) ! (mag) ?=99.999 rms uncertainty on misoc
*                                        (ERR_MAG_ISOCOR)
      real*8        Fisoc_2    (nr__3) ! ?=9999999.9 Corrected isophotal flux
*                                        (FLUX_ISOCOR)
      real*4        e_Fisoc_2  (nr__3) ! ?=9999.9 rms uncertainty on Fisoc
*                                        (ERR_FLUX_ISOCOR)
      real*8        Fauto_2    (nr__3) ! ?=9999999.9 Flux Auto (FLUX_AUTO)
      real*4        e_Fauto_2  (nr__3) ! ?=9999.9 rms uncertainty on Fauto
*                                        (ERR_FLUX_AUTO)
      real*4        mbest_2    (nr__3) ! (mag) ?=99.999 Best value of AB magnitude (MAG_BEST)
      real*4        e_mbest_2  (nr__3) ! (mag) ?=99.999 rms uncertainty on mbest
*                                        (ERR_MAG_BEST)
      real*8        Fbest_2    (nr__3) ! ?=9999999.9 Best value for flux (FLUX_BEST)
      real*4        e_Fbest_2  (nr__3) ! ?=9999.9 rms uncertainty on Fbest
*                                        (ERR_FLUX_BEST)
      real*8        F1ap_2     (nr__3) ! ?=9999999.9 Aperture flux 1 (FLUX1_APER)
      real*4        e_F1ap_2   (nr__3) ! ?=9999.9 rms uncertainty on F1ap
*                                        (ERR_FLUX1_APER)
      real*8        F2ap_2     (nr__3) ! ?=9999999.9 Aperture flux 2 (FLUX2_APER)
      real*4        e_F2ap_2   (nr__3) ! ?=9999.9 rms uncertainty on F2ap
*                                        (ERR_FLUX2_APER)
      integer*4     Sflg_2     (nr__3) ! ?=99 Sextrator flags (FLAGS)

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

C  Declarations for 'rcat.dat'	! R-band catalog (Cousins)

      integer*4 nr__4
      parameter (nr__4=253588)	! Number of records
      character*297 ar__4  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS_4     (nr__4) ! [20035000,20472316] Sequential number
      real*8        RAdeg_4    (nr__4) ! (deg) Rights ascension J2000 (ALPHA_J2000)
      real*8        DEdeg_4    (nr__4) ! (deg) Declination J2000 (DELTA_J2000)
      real*4        mauto_3    (nr__4) ! (mag) ?=99.999 Auto AB magnitude (MAG_AUTO)
      real*4        e_mauto_3  (nr__4) ! (mag) ?=99.999 rms uncertainty on mauto
*                                     (ERR_MAG_AUTO)
      real*4        mautoc_3   (nr__4) ! (mag) ?=99.999 Corrected auto AB magnitude
*                                     (MAG_AUTO_CORR)
      real*4        mag1ap_3   (nr__4) ! (mag) ?=99.999 Aperture AB magnitude 1 (MAG1_APER)
      real*4        e_mag1ap_3 (nr__4) ! (mag) ?=99.999 rms uncertainty on mag1ap
*                                    (ERR_MAG1_APER)
      real*4        mag2ap_3   (nr__4) ! (mag) ?=99.999 Aperture AB magnitude 2 (MAG2_APER)
      real*4        e_mag2ap_3 (nr__4) ! (mag) ?=99.999 rms uncertainty on mag2ap
*                                    (ERR_MAG2_APER)
      integer*4     isoA_3     (nr__4) ! (pix) ?=99999 Image isophotal area (ISOAREA_IMAGE)
      real*4        thr_3      (nr__4) ! ?=999.9 Threshold (THRESHOLD)
      real*4        muthr_3    (nr__4) ! (mag/arcsec2) ?=99.999 Surface brightness threshold
*                                     (MU_THRESHOLD)
      real*4        backg_3    (nr__4) ! ?=99.999 Backgroung flux (BACKGROUND)
      real*8        Fmax_3     (nr__4) ! ?=9999.999 Maximum flux (FLUX_MAX)
      real*4        muMax_3    (nr__4) ! (mag/arcsec2) ?=99.999 Maximum Surface brightness (MU_MAX)
      real*4        Krad_3     (nr__4) ! (arcsec) ?=99.99 Kron radius (KRON_RADIUS)
      real*8        F1rad_3    (nr__4) ! (arcsec) ?=9999.99 Flux1 radius (FLUX1_RADIUS)
      real*8        F2rad_3    (nr__4) ! (arcsec) ?=9999.99 Flux2 radius (FLUX2_RADIUS)
      real*8        F3rad_3    (nr__4) ! (arcsec) ?=9999.99 Flux3 radius (FLUX3_RADIUS)
      real*4        miso_3     (nr__4) ! (mag) ?=99.999 Isophotal AB magnitude (MAG_ISO)
      real*4        e_miso_3   (nr__4) ! (mag) ?=99.999 rms uncertainty on miso (ERR_MAG_ISO)
      real*8        Fiso_3     (nr__4) ! ?=9999999.9 Isophotal flux (FLUX_ISO)
      real*4        e_Fiso_3   (nr__4) ! ?=9999.9 rms uncertainty on Fiso
*                                         (ERR_FLUX_ISO)
      real*4        misoc_3    (nr__4) ! (mag) ?=99.999 Corrected isophotal AB magnitude
*                                        (MAG_ISOCOR)
      real*4        e_misoc_3  (nr__4) ! (mag) ?=99.999 rms uncertainty on misoc
*                                        (ERR_MAG_ISOCOR)
      real*8        Fisoc_3    (nr__4) ! ?=9999999.9 Corrected isophotal flux
*                                        (FLUX_ISOCOR)
      real*4        e_Fisoc_3  (nr__4) ! ?=9999.9 rms uncertainty on Fisoc
*                                        (ERR_FLUX_ISOCOR)
      real*8        Fauto_3    (nr__4) ! ?=9999999.9 Flux Auto (FLUX_AUTO)
      real*4        e_Fauto_3  (nr__4) ! ?=9999.9 rms uncertainty on Fauto
*                                        (ERR_FLUX_AUTO)
      real*4        mbest_3    (nr__4) ! (mag) ?=99.999 Best value of AB magnitude (MAG_BEST)
      real*4        e_mbest_3  (nr__4) ! (mag) ?=99.999 rms uncertainty on mbest
*                                        (ERR_MAG_BEST)
      real*8        Fbest_3    (nr__4) ! ?=9999999.9 Best value for flux (FLUX_BEST)
      real*4        e_Fbest_3  (nr__4) ! ?=9999.9 rms uncertainty on Fbest
*                                        (ERR_FLUX_BEST)
      real*8        F1ap_3     (nr__4) ! ?=9999999.9 Aperture flux 1 (FLUX1_APER)
      real*4        e_F1ap_3   (nr__4) ! ?=9999.9 rms uncertainty on F1ap
*                                        (ERR_FLUX1_APER)
      real*8        F2ap_3     (nr__4) ! ?=9999999.9 Aperture flux 2 (FLUX2_APER)
      real*4        e_F2ap_3   (nr__4) ! ?=9999.9 rms uncertainty on F2ap
*                                        (ERR_FLUX2_APER)
      integer*4     Sflg_3     (nr__4) ! ?=99 Sextrator flags (FLAGS)

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

C  Declarations for 'icat.dat'	! I-band catalog (Cousins)

      integer*4 nr__5
      parameter (nr__5=253588)	! Number of records
      character*297 ar__5  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS_5     (nr__5) ! [20035000,20472316] Sequential number
      real*8        RAdeg_5    (nr__5) ! (deg) Rights ascension J2000 (ALPHA_J2000)
      real*8        DEdeg_5    (nr__5) ! (deg) Declination J2000 (DELTA_J2000)
      real*4        mauto_4    (nr__5) ! (mag) ?=99.999 Auto AB magnitude (MAG_AUTO)
      real*4        e_mauto_4  (nr__5) ! (mag) ?=99.999 rms uncertainty on mauto
*                                     (ERR_MAG_AUTO)
      real*4        mautoc_4   (nr__5) ! (mag) ?=99.999 Corrected auto AB magnitude
*                                     (MAG_AUTO_CORR)
      real*4        mag1ap_4   (nr__5) ! (mag) ?=99.999 Aperture AB magnitude 1 (MAG1_APER)
      real*4        e_mag1ap_4 (nr__5) ! (mag) ?=99.999 rms uncertainty on mag1ap
*                                    (ERR_MAG1_APER)
      real*4        mag2ap_4   (nr__5) ! (mag) ?=99.999 Aperture AB magnitude 2 (MAG2_APER)
      real*4        e_mag2ap_4 (nr__5) ! (mag) ?=99.999 rms uncertainty on mag2ap
*                                    (ERR_MAG2_APER)
      integer*4     isoA_4     (nr__5) ! (pix) ?=99999 Image isophotal area (ISOAREA_IMAGE)
      real*4        thr_4      (nr__5) ! ?=999.9 Threshold (THRESHOLD)
      real*4        muthr_4    (nr__5) ! (mag/arcsec2) ?=99.999 Surface brightness threshold
*                                     (MU_THRESHOLD)
      real*4        backg_4    (nr__5) ! ?=99.999 Backgroung flux (BACKGROUND)
      real*8        Fmax_4     (nr__5) ! ?=9999.999 Maximum flux (FLUX_MAX)
      real*4        muMax_4    (nr__5) ! (mag/arcsec2) ?=99.999 Maximum Surface brightness (MU_MAX)
      real*4        Krad_4     (nr__5) ! (arcsec) ?=99.99 Kron radius (KRON_RADIUS)
      real*8        F1rad_4    (nr__5) ! (arcsec) ?=9999.99 Flux1 radius (FLUX1_RADIUS)
      real*8        F2rad_4    (nr__5) ! (arcsec) ?=9999.99 Flux2 radius (FLUX2_RADIUS)
      real*8        F3rad_4    (nr__5) ! (arcsec) ?=9999.99 Flux3 radius (FLUX3_RADIUS)
      real*4        miso_4     (nr__5) ! (mag) ?=99.999 Isophotal AB magnitude (MAG_ISO)
      real*4        e_miso_4   (nr__5) ! (mag) ?=99.999 rms uncertainty on miso (ERR_MAG_ISO)
      real*8        Fiso_4     (nr__5) ! ?=9999999.9 Isophotal flux (FLUX_ISO)
      real*4        e_Fiso_4   (nr__5) ! ?=9999.9 rms uncertainty on Fiso
*                                         (ERR_FLUX_ISO)
      real*4        misoc_4    (nr__5) ! (mag) ?=99.999 Corrected isophotal AB magnitude
*                                        (MAG_ISOCOR)
      real*4        e_misoc_4  (nr__5) ! (mag) ?=99.999 rms uncertainty on misoc
*                                        (ERR_MAG_ISOCOR)
      real*8        Fisoc_4    (nr__5) ! ?=9999999.9 Corrected isophotal flux
*                                        (FLUX_ISOCOR)
      real*4        e_Fisoc_4  (nr__5) ! ?=9999.9 rms uncertainty on Fisoc
*                                        (ERR_FLUX_ISOCOR)
      real*8        Fauto_4    (nr__5) ! ?=9999999.9 Flux Auto (FLUX_AUTO)
      real*4        e_Fauto_4  (nr__5) ! ?=9999.9 rms uncertainty on Fauto
*                                        (ERR_FLUX_AUTO)
      real*4        mbest_4    (nr__5) ! (mag) ?=99.999 Best value of AB magnitude (MAG_BEST)
      real*4        e_mbest_4  (nr__5) ! (mag) ?=99.999 rms uncertainty on mbest
*                                        (ERR_MAG_BEST)
      real*8        Fbest_4    (nr__5) ! ?=9999999.9 Best value for flux (FLUX_BEST)
      real*4        e_Fbest_4  (nr__5) ! ?=9999.9 rms uncertainty on Fbest
*                                        (ERR_FLUX_BEST)
      real*8        F1ap_4     (nr__5) ! ?=9999999.9 Aperture flux 1 (FLUX1_APER)
      real*4        e_F1ap_4   (nr__5) ! ?=9999.9 rms uncertainty on F1ap
*                                        (ERR_FLUX1_APER)
      real*8        F2ap_4     (nr__5) ! ?=9999999.9 Aperture flux 2 (FLUX2_APER)
      real*4        e_F2ap_4   (nr__5) ! ?=9999.9 rms uncertainty on F2ap
*                                        (ERR_FLUX2_APER)
      integer*4     Sflg_4     (nr__5) ! ?=99 Sextrator flags (FLAGS)

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

C  Declarations for 'jcat.dat'	! J-band catalog (Johnson)

      integer*4 nr__6
      parameter (nr__6=253588)	! Number of records
      character*297 ar__6  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS_6     (nr__6) ! [20035000,20472316] Sequential number
      real*8        RAdeg_6    (nr__6) ! (deg) Rights ascension J2000 (ALPHA_J2000)
      real*8        DEdeg_6    (nr__6) ! (deg) Declination J2000 (DELTA_J2000)
      real*4        mauto_5    (nr__6) ! (mag) ?=99.999 Auto AB magnitude (MAG_AUTO)
      real*4        e_mauto_5  (nr__6) ! (mag) ?=99.999 rms uncertainty on mauto
*                                     (ERR_MAG_AUTO)
      real*4        mautoc_5   (nr__6) ! (mag) ?=99.999 Corrected auto AB magnitude
*                                     (MAG_AUTO_CORR)
      real*4        mag1ap_5   (nr__6) ! (mag) ?=99.999 Aperture AB magnitude 1 (MAG1_APER)
      real*4        e_mag1ap_5 (nr__6) ! (mag) ?=99.999 rms uncertainty on mag1ap
*                                    (ERR_MAG1_APER)
      real*4        mag2ap_5   (nr__6) ! (mag) ?=99.999 Aperture AB magnitude 2 (MAG2_APER)
      real*4        e_mag2ap_5 (nr__6) ! (mag) ?=99.999 rms uncertainty on mag2ap
*                                    (ERR_MAG2_APER)
      integer*4     isoA_5     (nr__6) ! (pix) ?=99999 Image isophotal area (ISOAREA_IMAGE)
      real*4        thr_5      (nr__6) ! ?=999.9 Threshold (THRESHOLD)
      real*4        muthr_5    (nr__6) ! (mag/arcsec2) ?=99.999 Surface brightness threshold
*                                     (MU_THRESHOLD)
      real*4        backg_5    (nr__6) ! ?=99.999 Backgroung flux (BACKGROUND)
      real*8        Fmax_5     (nr__6) ! ?=9999.999 Maximum flux (FLUX_MAX)
      real*4        muMax_5    (nr__6) ! (mag/arcsec2) ?=99.999 Maximum Surface brightness (MU_MAX)
      real*4        Krad_5     (nr__6) ! (arcsec) ?=99.99 Kron radius (KRON_RADIUS)
      real*8        F1rad_5    (nr__6) ! (arcsec) ?=9999.99 Flux1 radius (FLUX1_RADIUS)
      real*8        F2rad_5    (nr__6) ! (arcsec) ?=9999.99 Flux2 radius (FLUX2_RADIUS)
      real*8        F3rad_5    (nr__6) ! (arcsec) ?=9999.99 Flux3 radius (FLUX3_RADIUS)
      real*4        miso_5     (nr__6) ! (mag) ?=99.999 Isophotal AB magnitude (MAG_ISO)
      real*4        e_miso_5   (nr__6) ! (mag) ?=99.999 rms uncertainty on miso (ERR_MAG_ISO)
      real*8        Fiso_5     (nr__6) ! ?=9999999.9 Isophotal flux (FLUX_ISO)
      real*4        e_Fiso_5   (nr__6) ! ?=9999.9 rms uncertainty on Fiso
*                                         (ERR_FLUX_ISO)
      real*4        misoc_5    (nr__6) ! (mag) ?=99.999 Corrected isophotal AB magnitude
*                                        (MAG_ISOCOR)
      real*4        e_misoc_5  (nr__6) ! (mag) ?=99.999 rms uncertainty on misoc
*                                        (ERR_MAG_ISOCOR)
      real*8        Fisoc_5    (nr__6) ! ?=9999999.9 Corrected isophotal flux
*                                        (FLUX_ISOCOR)
      real*4        e_Fisoc_5  (nr__6) ! ?=9999.9 rms uncertainty on Fisoc
*                                        (ERR_FLUX_ISOCOR)
      real*8        Fauto_5    (nr__6) ! ?=9999999.9 Flux Auto (FLUX_AUTO)
      real*4        e_Fauto_5  (nr__6) ! ?=9999.9 rms uncertainty on Fauto
*                                        (ERR_FLUX_AUTO)
      real*4        mbest_5    (nr__6) ! (mag) ?=99.999 Best value of AB magnitude (MAG_BEST)
      real*4        e_mbest_5  (nr__6) ! (mag) ?=99.999 rms uncertainty on mbest
*                                        (ERR_MAG_BEST)
      real*8        Fbest_5    (nr__6) ! ?=9999999.9 Best value for flux (FLUX_BEST)
      real*4        e_Fbest_5  (nr__6) ! ?=9999.9 rms uncertainty on Fbest
*                                        (ERR_FLUX_BEST)
      real*8        F1ap_5     (nr__6) ! ?=9999999.9 Aperture flux 1 (FLUX1_APER)
      real*4        e_F1ap_5   (nr__6) ! ?=9999.9 rms uncertainty on F1ap
*                                        (ERR_FLUX1_APER)
      real*8        F2ap_5     (nr__6) ! ?=9999999.9 Aperture flux 2 (FLUX2_APER)
      real*4        e_F2ap_5   (nr__6) ! ?=9999.9 rms uncertainty on F2ap
*                                        (ERR_FLUX2_APER)
      integer*4     Sflg_5     (nr__6) ! ?=99 Sextrator flags (FLAGS)

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

C  Declarations for 'kcat.dat'	! K-band catalog (Johnson)

      integer*4 nr__7
      parameter (nr__7=253588)	! Number of records
      character*297 ar__7  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     VVDS_7     (nr__7) ! [20035000,20472316] Sequential number
      real*8        RAdeg_7    (nr__7) ! (deg) Rights ascension J2000 (ALPHA_J2000)
      real*8        DEdeg_7    (nr__7) ! (deg) Declination J2000 (DELTA_J2000)
      real*4        mauto_6    (nr__7) ! (mag) ?=99.999 Auto AB magnitude (MAG_AUTO)
      real*4        e_mauto_6  (nr__7) ! (mag) ?=99.999 rms uncertainty on mauto
*                                     (ERR_MAG_AUTO)
      real*4        mautoc_6   (nr__7) ! (mag) ?=99.999 Corrected auto AB magnitude
*                                     (MAG_AUTO_CORR)
      real*4        mag1ap_6   (nr__7) ! (mag) ?=99.999 Aperture AB magnitude 1 (MAG1_APER)
      real*4        e_mag1ap_6 (nr__7) ! (mag) ?=99.999 rms uncertainty on mag1ap
*                                    (ERR_MAG1_APER)
      real*4        mag2ap_6   (nr__7) ! (mag) ?=99.999 Aperture AB magnitude 2 (MAG2_APER)
      real*4        e_mag2ap_6 (nr__7) ! (mag) ?=99.999 rms uncertainty on mag2ap
*                                    (ERR_MAG2_APER)
      integer*4     isoA_6     (nr__7) ! (pix) ?=99999 Image isophotal area (ISOAREA_IMAGE)
      real*4        thr_6      (nr__7) ! ?=999.9 Threshold (THRESHOLD)
      real*4        muthr_6    (nr__7) ! (mag/arcsec2) ?=99.999 Surface brightness threshold
*                                     (MU_THRESHOLD)
      real*4        backg_6    (nr__7) ! ?=99.999 Backgroung flux (BACKGROUND)
      real*8        Fmax_6     (nr__7) ! ?=9999.999 Maximum flux (FLUX_MAX)
      real*4        muMax_6    (nr__7) ! (mag/arcsec2) ?=99.999 Maximum Surface brightness (MU_MAX)
      real*4        Krad_6     (nr__7) ! (arcsec) ?=99.99 Kron radius (KRON_RADIUS)
      real*8        F1rad_6    (nr__7) ! (arcsec) ?=9999.99 Flux1 radius (FLUX1_RADIUS)
      real*8        F2rad_6    (nr__7) ! (arcsec) ?=9999.99 Flux2 radius (FLUX2_RADIUS)
      real*8        F3rad_6    (nr__7) ! (arcsec) ?=9999.99 Flux3 radius (FLUX3_RADIUS)
      real*4        miso_6     (nr__7) ! (mag) ?=99.999 Isophotal AB magnitude (MAG_ISO)
      real*4        e_miso_6   (nr__7) ! (mag) ?=99.999 rms uncertainty on miso (ERR_MAG_ISO)
      real*8        Fiso_6     (nr__7) ! ?=9999999.9 Isophotal flux (FLUX_ISO)
      real*4        e_Fiso_6   (nr__7) ! ?=9999.9 rms uncertainty on Fiso
*                                         (ERR_FLUX_ISO)
      real*4        misoc_6    (nr__7) ! (mag) ?=99.999 Corrected isophotal AB magnitude
*                                        (MAG_ISOCOR)
      real*4        e_misoc_6  (nr__7) ! (mag) ?=99.999 rms uncertainty on misoc
*                                        (ERR_MAG_ISOCOR)
      real*8        Fisoc_6    (nr__7) ! ?=9999999.9 Corrected isophotal flux
*                                        (FLUX_ISOCOR)
      real*4        e_Fisoc_6  (nr__7) ! ?=9999.9 rms uncertainty on Fisoc
*                                        (ERR_FLUX_ISOCOR)
      real*8        Fauto_6    (nr__7) ! ?=9999999.9 Flux Auto (FLUX_AUTO)
      real*4        e_Fauto_6  (nr__7) ! ?=9999.9 rms uncertainty on Fauto
*                                        (ERR_FLUX_AUTO)
      real*4        mbest_6    (nr__7) ! (mag) ?=99.999 Best value of AB magnitude (MAG_BEST)
      real*4        e_mbest_6  (nr__7) ! (mag) ?=99.999 rms uncertainty on mbest
*                                        (ERR_MAG_BEST)
      real*8        Fbest_6    (nr__7) ! ?=9999999.9 Best value for flux (FLUX_BEST)
      real*4        e_Fbest_6  (nr__7) ! ?=9999.9 rms uncertainty on Fbest
*                                        (ERR_FLUX_BEST)
      real*8        F1ap_6     (nr__7) ! ?=9999999.9 Aperture flux 1 (FLUX1_APER)
      real*4        e_F1ap_6   (nr__7) ! ?=9999.9 rms uncertainty on F1ap
*                                        (ERR_FLUX1_APER)
      real*8        F2ap_6     (nr__7) ! ?=9999999.9 Aperture flux 2 (FLUX2_APER)
      real*4        e_F2ap_6   (nr__7) ! ?=9999.9 rms uncertainty on F2ap
*                                        (ERR_FLUX2_APER)
      integer*4     Sflg_6     (nr__7) ! ?=99 Sextrator flags (FLAGS)

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

C  Loading file 'morph.dat'	! Morphological catalog (in filter I)

C  Format for file interpretation

    1 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F8.2,1X,F8.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,
     +  F5.2,1X,F5.1,1X,F5.1,1X,F8.2,1X,F5.2,1X,F8.2,1X,F5.2,1X,F8.2,
     +  1X,F5.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,F5.1,1X,F5.1,1X,
     +  F8.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,F4.2,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'morph.dat')
      write(6,*) '....Loading file: morph.dat'
      do i__=1,253588
        read(1,'(A224)')ar__
        read(ar__,1)
     +  VVDS(i__),RAdeg(i__),DEdeg(i__),Xi(i__),Yi(i__),Ai(i__),
     +  e_Ai(i__),Bi(i__),e_Bi(i__),PAi(i__),ePAi(i__),X2i(i__),
     +  e_X2i(i__),Y2i(i__),e_Y2i(i__),XYi(i__),eXYi(i__),Aw(i__),
     +  e_Aw(i__),Bw(i__),e_Bw(i__),PAw(i__),ePAw(i__),X2w(i__),
     +  e_X2w(i__),Y2w(i__),e_Y2w(i__),XYw(i__),eXYw(i__),S_G(i__),
     +  Flag(i__)
c    ..............Just test output...........
        write(6,1)
     +  VVDS(i__),RAdeg(i__),DEdeg(i__),Xi(i__),Yi(i__),Ai(i__),
     +  e_Ai(i__),Bi(i__),e_Bi(i__),PAi(i__),ePAi(i__),X2i(i__),
     +  e_X2i(i__),Y2i(i__),e_Y2i(i__),XYi(i__),eXYi(i__),Aw(i__),
     +  e_Aw(i__),Bw(i__),e_Bw(i__),PAw(i__),ePAw(i__),X2w(i__),
     +  e_X2w(i__),Y2w(i__),e_Y2w(i__),XYw(i__),eXYw(i__),S_G(i__),
     +  Flag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'ucat.dat'	! *U-band catalog (Ueso and Uloiano bands)

C  Format for file interpretation

    2 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,
     +  1X,I5,1X,F5.1,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,F5.2,1X,F7.2,
     +  1X,F7.2,1X,F7.2,1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,1X,F6.3,1X,
     +  F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F6.3,1X,F6.3,1X,F9.1,
     +  1X,F6.1,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'ucat.dat')
      write(6,*) '....Loading file: ucat.dat'
      do i__=1,253588
        read(1,'(A339)')ar__1
        read(ar__1,2)
     +  VVDS_1(i__),RAdeg_1(i__),DEdeg_1(i__),mauto(i__),e_mauto(i__),
     +  mautoc(i__),mag1ap(i__),e_mag1ap(i__),mag2ap(i__),
     +  e_mag2ap(i__),mautol(i__),e_mautol(i__),m1autol(i__),
     +  e_m1autol(i__),m2autol(i__),e_m2autol(i__),isoA(i__),thr(i__),
     +  muthr(i__),backg(i__),Fmax(i__),muMax(i__),Krad(i__),
     +  F1rad(i__),F2rad(i__),F3rad(i__),miso(i__),e_miso(i__),
     +  Fiso(i__),e_Fiso(i__),misoc(i__),e_misoc(i__),Fisoc(i__),
     +  e_Fisoc(i__),Fauto(i__),e_Fauto(i__),mbest(i__),e_mbest(i__),
     +  Fbest(i__),e_Fbest(i__),F1ap(i__),e_F1ap(i__),F2ap(i__),
     +  e_F2ap(i__),Sflg(i__)
c    ..............Just test output...........
        write(6,2)
     +  VVDS_1(i__),RAdeg_1(i__),DEdeg_1(i__),mauto(i__),e_mauto(i__),
     +  mautoc(i__),mag1ap(i__),e_mag1ap(i__),mag2ap(i__),
     +  e_mag2ap(i__),mautol(i__),e_mautol(i__),m1autol(i__),
     +  e_m1autol(i__),m2autol(i__),e_m2autol(i__),isoA(i__),thr(i__),
     +  muthr(i__),backg(i__),Fmax(i__),muMax(i__),Krad(i__),
     +  F1rad(i__),F2rad(i__),F3rad(i__),miso(i__),e_miso(i__),
     +  Fiso(i__),e_Fiso(i__),misoc(i__),e_misoc(i__),Fisoc(i__),
     +  e_Fisoc(i__),Fauto(i__),e_Fauto(i__),mbest(i__),e_mbest(i__),
     +  Fbest(i__),e_Fbest(i__),F1ap(i__),e_F1ap(i__),F2ap(i__),
     +  e_F2ap(i__),Sflg(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'bcat.dat'	! B-band catalog (Johnson)

C  Format for file interpretation

    3 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,I5,1X,F5.1,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,
     +  F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,
     +  1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F6.3,1X,
     +  F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'bcat.dat')
      write(6,*) '....Loading file: bcat.dat'
      do i__=1,253588
        read(1,'(A297)')ar__2
        read(ar__2,3)
     +  VVDS_2(i__),RAdeg_2(i__),DEdeg_2(i__),mauto_1(i__),
     +  e_mauto_1(i__),mautoc_1(i__),mag1ap_1(i__),e_mag1ap_1(i__),
     +  mag2ap_1(i__),e_mag2ap_1(i__),isoA_1(i__),thr_1(i__),
     +  muthr_1(i__),backg_1(i__),Fmax_1(i__),muMax_1(i__),
     +  Krad_1(i__),F1rad_1(i__),F2rad_1(i__),F3rad_1(i__),
     +  miso_1(i__),e_miso_1(i__),Fiso_1(i__),e_Fiso_1(i__),
     +  misoc_1(i__),e_misoc_1(i__),Fisoc_1(i__),e_Fisoc_1(i__),
     +  Fauto_1(i__),e_Fauto_1(i__),mbest_1(i__),e_mbest_1(i__),
     +  Fbest_1(i__),e_Fbest_1(i__),F1ap_1(i__),e_F1ap_1(i__),
     +  F2ap_1(i__),e_F2ap_1(i__),Sflg_1(i__)
c    ..............Just test output...........
        write(6,3)
     +  VVDS_2(i__),RAdeg_2(i__),DEdeg_2(i__),mauto_1(i__),
     +  e_mauto_1(i__),mautoc_1(i__),mag1ap_1(i__),e_mag1ap_1(i__),
     +  mag2ap_1(i__),e_mag2ap_1(i__),isoA_1(i__),thr_1(i__),
     +  muthr_1(i__),backg_1(i__),Fmax_1(i__),muMax_1(i__),
     +  Krad_1(i__),F1rad_1(i__),F2rad_1(i__),F3rad_1(i__),
     +  miso_1(i__),e_miso_1(i__),Fiso_1(i__),e_Fiso_1(i__),
     +  misoc_1(i__),e_misoc_1(i__),Fisoc_1(i__),e_Fisoc_1(i__),
     +  Fauto_1(i__),e_Fauto_1(i__),mbest_1(i__),e_mbest_1(i__),
     +  Fbest_1(i__),e_Fbest_1(i__),F1ap_1(i__),e_F1ap_1(i__),
     +  F2ap_1(i__),e_F2ap_1(i__),Sflg_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'vcat.dat'	! V-band catalog (Johnson)

C  Format for file interpretation

    4 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,I5,1X,F5.1,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,
     +  F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,
     +  1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F6.3,1X,
     +  F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'vcat.dat')
      write(6,*) '....Loading file: vcat.dat'
      do i__=1,253588
        read(1,'(A297)')ar__3
        read(ar__3,4)
     +  VVDS_3(i__),RAdeg_3(i__),DEdeg_3(i__),mauto_2(i__),
     +  e_mauto_2(i__),mautoc_2(i__),mag1ap_2(i__),e_mag1ap_2(i__),
     +  mag2ap_2(i__),e_mag2ap_2(i__),isoA_2(i__),thr_2(i__),
     +  muthr_2(i__),backg_2(i__),Fmax_2(i__),muMax_2(i__),
     +  Krad_2(i__),F1rad_2(i__),F2rad_2(i__),F3rad_2(i__),
     +  miso_2(i__),e_miso_2(i__),Fiso_2(i__),e_Fiso_2(i__),
     +  misoc_2(i__),e_misoc_2(i__),Fisoc_2(i__),e_Fisoc_2(i__),
     +  Fauto_2(i__),e_Fauto_2(i__),mbest_2(i__),e_mbest_2(i__),
     +  Fbest_2(i__),e_Fbest_2(i__),F1ap_2(i__),e_F1ap_2(i__),
     +  F2ap_2(i__),e_F2ap_2(i__),Sflg_2(i__)
c    ..............Just test output...........
        write(6,4)
     +  VVDS_3(i__),RAdeg_3(i__),DEdeg_3(i__),mauto_2(i__),
     +  e_mauto_2(i__),mautoc_2(i__),mag1ap_2(i__),e_mag1ap_2(i__),
     +  mag2ap_2(i__),e_mag2ap_2(i__),isoA_2(i__),thr_2(i__),
     +  muthr_2(i__),backg_2(i__),Fmax_2(i__),muMax_2(i__),
     +  Krad_2(i__),F1rad_2(i__),F2rad_2(i__),F3rad_2(i__),
     +  miso_2(i__),e_miso_2(i__),Fiso_2(i__),e_Fiso_2(i__),
     +  misoc_2(i__),e_misoc_2(i__),Fisoc_2(i__),e_Fisoc_2(i__),
     +  Fauto_2(i__),e_Fauto_2(i__),mbest_2(i__),e_mbest_2(i__),
     +  Fbest_2(i__),e_Fbest_2(i__),F1ap_2(i__),e_F1ap_2(i__),
     +  F2ap_2(i__),e_F2ap_2(i__),Sflg_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'rcat.dat'	! R-band catalog (Cousins)

C  Format for file interpretation

    5 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,I5,1X,F5.1,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,
     +  F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,
     +  1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F6.3,1X,
     +  F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'rcat.dat')
      write(6,*) '....Loading file: rcat.dat'
      do i__=1,253588
        read(1,'(A297)')ar__4
        read(ar__4,5)
     +  VVDS_4(i__),RAdeg_4(i__),DEdeg_4(i__),mauto_3(i__),
     +  e_mauto_3(i__),mautoc_3(i__),mag1ap_3(i__),e_mag1ap_3(i__),
     +  mag2ap_3(i__),e_mag2ap_3(i__),isoA_3(i__),thr_3(i__),
     +  muthr_3(i__),backg_3(i__),Fmax_3(i__),muMax_3(i__),
     +  Krad_3(i__),F1rad_3(i__),F2rad_3(i__),F3rad_3(i__),
     +  miso_3(i__),e_miso_3(i__),Fiso_3(i__),e_Fiso_3(i__),
     +  misoc_3(i__),e_misoc_3(i__),Fisoc_3(i__),e_Fisoc_3(i__),
     +  Fauto_3(i__),e_Fauto_3(i__),mbest_3(i__),e_mbest_3(i__),
     +  Fbest_3(i__),e_Fbest_3(i__),F1ap_3(i__),e_F1ap_3(i__),
     +  F2ap_3(i__),e_F2ap_3(i__),Sflg_3(i__)
c    ..............Just test output...........
        write(6,5)
     +  VVDS_4(i__),RAdeg_4(i__),DEdeg_4(i__),mauto_3(i__),
     +  e_mauto_3(i__),mautoc_3(i__),mag1ap_3(i__),e_mag1ap_3(i__),
     +  mag2ap_3(i__),e_mag2ap_3(i__),isoA_3(i__),thr_3(i__),
     +  muthr_3(i__),backg_3(i__),Fmax_3(i__),muMax_3(i__),
     +  Krad_3(i__),F1rad_3(i__),F2rad_3(i__),F3rad_3(i__),
     +  miso_3(i__),e_miso_3(i__),Fiso_3(i__),e_Fiso_3(i__),
     +  misoc_3(i__),e_misoc_3(i__),Fisoc_3(i__),e_Fisoc_3(i__),
     +  Fauto_3(i__),e_Fauto_3(i__),mbest_3(i__),e_mbest_3(i__),
     +  Fbest_3(i__),e_Fbest_3(i__),F1ap_3(i__),e_F1ap_3(i__),
     +  F2ap_3(i__),e_F2ap_3(i__),Sflg_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'icat.dat'	! I-band catalog (Cousins)

C  Format for file interpretation

    6 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,I5,1X,F5.1,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,
     +  F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,
     +  1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F6.3,1X,
     +  F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'icat.dat')
      write(6,*) '....Loading file: icat.dat'
      do i__=1,253588
        read(1,'(A297)')ar__5
        read(ar__5,6)
     +  VVDS_5(i__),RAdeg_5(i__),DEdeg_5(i__),mauto_4(i__),
     +  e_mauto_4(i__),mautoc_4(i__),mag1ap_4(i__),e_mag1ap_4(i__),
     +  mag2ap_4(i__),e_mag2ap_4(i__),isoA_4(i__),thr_4(i__),
     +  muthr_4(i__),backg_4(i__),Fmax_4(i__),muMax_4(i__),
     +  Krad_4(i__),F1rad_4(i__),F2rad_4(i__),F3rad_4(i__),
     +  miso_4(i__),e_miso_4(i__),Fiso_4(i__),e_Fiso_4(i__),
     +  misoc_4(i__),e_misoc_4(i__),Fisoc_4(i__),e_Fisoc_4(i__),
     +  Fauto_4(i__),e_Fauto_4(i__),mbest_4(i__),e_mbest_4(i__),
     +  Fbest_4(i__),e_Fbest_4(i__),F1ap_4(i__),e_F1ap_4(i__),
     +  F2ap_4(i__),e_F2ap_4(i__),Sflg_4(i__)
c    ..............Just test output...........
        write(6,6)
     +  VVDS_5(i__),RAdeg_5(i__),DEdeg_5(i__),mauto_4(i__),
     +  e_mauto_4(i__),mautoc_4(i__),mag1ap_4(i__),e_mag1ap_4(i__),
     +  mag2ap_4(i__),e_mag2ap_4(i__),isoA_4(i__),thr_4(i__),
     +  muthr_4(i__),backg_4(i__),Fmax_4(i__),muMax_4(i__),
     +  Krad_4(i__),F1rad_4(i__),F2rad_4(i__),F3rad_4(i__),
     +  miso_4(i__),e_miso_4(i__),Fiso_4(i__),e_Fiso_4(i__),
     +  misoc_4(i__),e_misoc_4(i__),Fisoc_4(i__),e_Fisoc_4(i__),
     +  Fauto_4(i__),e_Fauto_4(i__),mbest_4(i__),e_mbest_4(i__),
     +  Fbest_4(i__),e_Fbest_4(i__),F1ap_4(i__),e_F1ap_4(i__),
     +  F2ap_4(i__),e_F2ap_4(i__),Sflg_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'jcat.dat'	! J-band catalog (Johnson)

C  Format for file interpretation

    7 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,I5,1X,F5.1,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,
     +  F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,
     +  1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F6.3,1X,
     +  F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'jcat.dat')
      write(6,*) '....Loading file: jcat.dat'
      do i__=1,253588
        read(1,'(A297)')ar__6
        read(ar__6,7)
     +  VVDS_6(i__),RAdeg_6(i__),DEdeg_6(i__),mauto_5(i__),
     +  e_mauto_5(i__),mautoc_5(i__),mag1ap_5(i__),e_mag1ap_5(i__),
     +  mag2ap_5(i__),e_mag2ap_5(i__),isoA_5(i__),thr_5(i__),
     +  muthr_5(i__),backg_5(i__),Fmax_5(i__),muMax_5(i__),
     +  Krad_5(i__),F1rad_5(i__),F2rad_5(i__),F3rad_5(i__),
     +  miso_5(i__),e_miso_5(i__),Fiso_5(i__),e_Fiso_5(i__),
     +  misoc_5(i__),e_misoc_5(i__),Fisoc_5(i__),e_Fisoc_5(i__),
     +  Fauto_5(i__),e_Fauto_5(i__),mbest_5(i__),e_mbest_5(i__),
     +  Fbest_5(i__),e_Fbest_5(i__),F1ap_5(i__),e_F1ap_5(i__),
     +  F2ap_5(i__),e_F2ap_5(i__),Sflg_5(i__)
c    ..............Just test output...........
        write(6,7)
     +  VVDS_6(i__),RAdeg_6(i__),DEdeg_6(i__),mauto_5(i__),
     +  e_mauto_5(i__),mautoc_5(i__),mag1ap_5(i__),e_mag1ap_5(i__),
     +  mag2ap_5(i__),e_mag2ap_5(i__),isoA_5(i__),thr_5(i__),
     +  muthr_5(i__),backg_5(i__),Fmax_5(i__),muMax_5(i__),
     +  Krad_5(i__),F1rad_5(i__),F2rad_5(i__),F3rad_5(i__),
     +  miso_5(i__),e_miso_5(i__),Fiso_5(i__),e_Fiso_5(i__),
     +  misoc_5(i__),e_misoc_5(i__),Fisoc_5(i__),e_Fisoc_5(i__),
     +  Fauto_5(i__),e_Fauto_5(i__),mbest_5(i__),e_mbest_5(i__),
     +  Fbest_5(i__),e_Fbest_5(i__),F1ap_5(i__),e_F1ap_5(i__),
     +  F2ap_5(i__),e_F2ap_5(i__),Sflg_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'kcat.dat'	! K-band catalog (Johnson)

C  Format for file interpretation

    8 format(
     +  I8,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,I5,1X,F5.1,1X,F6.3,1X,F6.3,1X,F8.3,1X,F6.3,1X,
     +  F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,
     +  1X,F6.3,1X,F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F6.3,1X,
     +  F6.3,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,F9.1,1X,F6.1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'kcat.dat')
      write(6,*) '....Loading file: kcat.dat'
      do i__=1,253588
        read(1,'(A297)')ar__7
        read(ar__7,8)
     +  VVDS_7(i__),RAdeg_7(i__),DEdeg_7(i__),mauto_6(i__),
     +  e_mauto_6(i__),mautoc_6(i__),mag1ap_6(i__),e_mag1ap_6(i__),
     +  mag2ap_6(i__),e_mag2ap_6(i__),isoA_6(i__),thr_6(i__),
     +  muthr_6(i__),backg_6(i__),Fmax_6(i__),muMax_6(i__),
     +  Krad_6(i__),F1rad_6(i__),F2rad_6(i__),F3rad_6(i__),
     +  miso_6(i__),e_miso_6(i__),Fiso_6(i__),e_Fiso_6(i__),
     +  misoc_6(i__),e_misoc_6(i__),Fisoc_6(i__),e_Fisoc_6(i__),
     +  Fauto_6(i__),e_Fauto_6(i__),mbest_6(i__),e_mbest_6(i__),
     +  Fbest_6(i__),e_Fbest_6(i__),F1ap_6(i__),e_F1ap_6(i__),
     +  F2ap_6(i__),e_F2ap_6(i__),Sflg_6(i__)
c    ..............Just test output...........
        write(6,8)
     +  VVDS_7(i__),RAdeg_7(i__),DEdeg_7(i__),mauto_6(i__),
     +  e_mauto_6(i__),mautoc_6(i__),mag1ap_6(i__),e_mag1ap_6(i__),
     +  mag2ap_6(i__),e_mag2ap_6(i__),isoA_6(i__),thr_6(i__),
     +  muthr_6(i__),backg_6(i__),Fmax_6(i__),muMax_6(i__),
     +  Krad_6(i__),F1rad_6(i__),F2rad_6(i__),F3rad_6(i__),
     +  miso_6(i__),e_miso_6(i__),Fiso_6(i__),e_Fiso_6(i__),
     +  misoc_6(i__),e_misoc_6(i__),Fisoc_6(i__),e_Fisoc_6(i__),
     +  Fauto_6(i__),e_Fauto_6(i__),mbest_6(i__),e_mbest_6(i__),
     +  Fbest_6(i__),e_Fbest_6(i__),F1ap_6(i__),e_F1ap_6(i__),
     +  F2ap_6(i__),e_F2ap_6(i__),Sflg_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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