FORTRAN Generation
(/./ftp/cats/J/A_A/506/999)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/506/999 into FORTRAN code for reading data files line by line.

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



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-18
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/506/999       Photometric Calibration of SNLS            (Regnault+, 2009)
*================================================================================
*Photometric Calibration of the SuperNova Legacy Survey Fields.
*    Regnault N., Conley A., Guy J., Sullivan M., Cuillandre J.-C., Astier P.,
*    Balland C., Basa S., Carlberg R.G., Fouchez D., Hardin D., Hook I.M.,
*    Howell D.A., Pain R., Perrett K., Pritchet C.J.
*   <Astron. Astrophys. 506, 999 (2009)>
*   =2009A&A...506..999R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table19.dat'	! Tertiary rgiz catalog for the SNLS field D1

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg       ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg       ! (deg) Declination, in decimal degrees (J2000)
      real*4        gmag        ! (mag) ?=-1 g-band magnitude (Vega-like system)
      real*4        e_gmag      ! (mag) ?=-1 Uncertainty on the g-band magnitude
      integer*4     o_gmag      ! [1,40]?=-1 Number of observations
      real*4        dkg         ! ?=-1 Grid g-band dk color term for this object
      real*4        rmag        ! (mag) ?=-1 r-band magnitude(Vega-like system)
      real*4        e_rmag      ! (mag) ?=-1 Uncertainty on the r-band magnitude
      integer*4     o_rmag      ! [1,40]?=-1 Number of observations
      real*4        dkr         ! ?=-1 Grid r-band dk color term for this object
      real*4        imag        ! (mag) ?=-1 i-band magnitude(Vega-like system)
      real*4        e_imag      ! (mag) ?=-1 Uncertainty on the i-band magnitude
      integer*4     o_imag      ! [1,40]?=-1 Number of observations
      real*4        dki         ! ?=-1 Grid i-band dk color term for this object
      real*4        zmag        ! (mag) ?=-1 z-band magnitude(Vega-like system)
      real*4        e_zmag      ! (mag) ?=-1 Uncertainty on the z-band magnitude
      integer*4     o_zmag      ! [1,40]?=-1 Number of observations
      real*4        dkz         ! ?=-1 Grid z-band dk color term for this object

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

C  Declarations for 'table20.dat'	! Tertiary rgiz catalog for the SNLS field D2

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_1     ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_1     ! (deg) Declination, in decimal degrees (J2000)
      real*4        gmag_1      ! (mag) ?=-1 g-band magnitude (Vega-like system)
      real*4        e_gmag_1    ! (mag) ?=-1 Uncertainty on the g-band magnitude
      integer*4     o_gmag_1    ! [1,40]?=-1 Number of observations
      real*4        dkg_1       ! ?=-1 Grid g-band dk color term for this object
      real*4        rmag_1      ! (mag) ?=-1 r-band magnitude(Vega-like system)
      real*4        e_rmag_1    ! (mag) ?=-1 Uncertainty on the r-band magnitude
      integer*4     o_rmag_1    ! [1,40]?=-1 Number of observations
      real*4        dkr_1       ! ?=-1 Grid r-band dk color term for this object
      real*4        imag_1      ! (mag) ?=-1 i-band magnitude(Vega-like system)
      real*4        e_imag_1    ! (mag) ?=-1 Uncertainty on the i-band magnitude
      integer*4     o_imag_1    ! [1,40]?=-1 Number of observations
      real*4        dki_1       ! ?=-1 Grid i-band dk color term for this object
      real*4        zmag_1      ! (mag) ?=-1 z-band magnitude(Vega-like system)
      real*4        e_zmag_1    ! (mag) ?=-1 Uncertainty on the z-band magnitude
      integer*4     o_zmag_1    ! [1,40]?=-1 Number of observations
      real*4        dkz_1       ! ?=-1 Grid z-band dk color term for this object

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

C  Declarations for 'table21.dat'	! Tertiary rgiz catalog for the SNLS field D3

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_2     ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_2     ! (deg) Declination, in decimal degrees (J2000)
      real*4        gmag_2      ! (mag) ?=-1 g-band magnitude (Vega-like system)
      real*4        e_gmag_2    ! (mag) ?=-1 Uncertainty on the g-band magnitude
      integer*4     o_gmag_2    ! [1,40]?=-1 Number of observations
      real*4        dkg_2       ! ?=-1 Grid g-band dk color term for this object
      real*4        rmag_2      ! (mag) ?=-1 r-band magnitude(Vega-like system)
      real*4        e_rmag_2    ! (mag) ?=-1 Uncertainty on the r-band magnitude
      integer*4     o_rmag_2    ! [1,40]?=-1 Number of observations
      real*4        dkr_2       ! ?=-1 Grid r-band dk color term for this object
      real*4        imag_2      ! (mag) ?=-1 i-band magnitude(Vega-like system)
      real*4        e_imag_2    ! (mag) ?=-1 Uncertainty on the i-band magnitude
      integer*4     o_imag_2    ! [1,40]?=-1 Number of observations
      real*4        dki_2       ! ?=-1 Grid i-band dk color term for this object
      real*4        zmag_2      ! (mag) ?=-1 z-band magnitude(Vega-like system)
      real*4        e_zmag_2    ! (mag) ?=-1 Uncertainty on the z-band magnitude
      integer*4     o_zmag_2    ! [1,40]?=-1 Number of observations
      real*4        dkz_2       ! ?=-1 Grid z-band dk color term for this object

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

C  Declarations for 'table22.dat'	! Tertiary rgiz catalog for the SNLS field D4

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_3     ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_3     ! (deg) Declination, in decimal degrees (J2000)
      real*4        gmag_3      ! (mag) ?=-1 g-band magnitude (Vega-like system)
      real*4        e_gmag_3    ! (mag) ?=-1 Uncertainty on the g-band magnitude
      integer*4     o_gmag_3    ! [1,40]?=-1 Number of observations
      real*4        dkg_3       ! ?=-1 Grid g-band dk color term for this object
      real*4        rmag_3      ! (mag) ?=-1 r-band magnitude(Vega-like system)
      real*4        e_rmag_3    ! (mag) ?=-1 Uncertainty on the r-band magnitude
      integer*4     o_rmag_3    ! [1,40]?=-1 Number of observations
      real*4        dkr_3       ! ?=-1 Grid r-band dk color term for this object
      real*4        imag_3      ! (mag) ?=-1 i-band magnitude(Vega-like system)
      real*4        e_imag_3    ! (mag) ?=-1 Uncertainty on the i-band magnitude
      integer*4     o_imag_3    ! [1,40]?=-1 Number of observations
      real*4        dki_3       ! ?=-1 Grid i-band dk color term for this object
      real*4        zmag_3      ! (mag) ?=-1 z-band magnitude(Vega-like system)
      real*4        e_zmag_3    ! (mag) ?=-1 Uncertainty on the z-band magnitude
      integer*4     o_zmag_3    ! [1,40]?=-1 Number of observations
      real*4        dkz_3       ! ?=-1 Grid z-band dk color term for this object

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

C  Declarations for 'd1u.dat'	! D1 tertiary stars u-band magnitudes

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_4     ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_4     ! (deg) Declination, in decimal degrees (J2000)
      real*4        umag        ! (mag) u-band magnitude
      real*4        e_umag      ! (mag) uncertainty on the u-band magnitude
      integer*4     o_umag      ! number of observations
      real*4        dku         ! grid u-band dk color term for this object

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

C  Declarations for 'd2u.dat'	! D2 tertiary stars u-band magnitudes

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_5     ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_5     ! (deg) Declination, in decimal degrees (J2000)
      real*4        umag_1      ! (mag) u-band magnitude
      real*4        e_umag_1    ! (mag) uncertainty on the u-band magnitude
      integer*4     o_umag_1    ! number of observations
      real*4        dku_1       ! grid u-band dk color term for this object

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

C  Declarations for 'd3u.dat'	! D3 tertiary stars u-band magnitudes

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_6     ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_6     ! (deg) Declination, in decimal degrees (J2000)
      real*4        umag_2      ! (mag) u-band magnitude
      real*4        e_umag_2    ! (mag) uncertainty on the u-band magnitude
      integer*4     o_umag_2    ! number of observations
      real*4        dku_2       ! grid u-band dk color term for this object

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

C  Declarations for 'd4u.dat'	! D4 tertiary stars u-band magnitudes

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_7     ! (deg) Right ascension, in decimal degrees (J2000)
      real*8        DEdeg_7     ! (deg) Declination, in decimal degrees (J2000)
      real*4        umag_3      ! (mag) u-band magnitude
      real*4        e_umag_3    ! (mag) uncertainty on the u-band magnitude
      integer*4     o_umag_3    ! number of observations
      real*4        dku_3       ! grid u-band dk color term for this object

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

C  Declarations for 'tabled1.dat'	! Full covariance matrix of the BD+17 4708 MegaCam
                               and Landolt magnitude uncertainties

      integer*4 nr__8
      parameter (nr__8=9)	! Number of records
      character*123 ar__8  	! Full-size record

      real*4        C1          ! Covariance matrix coefficient ({Delta}g)
      real*4        C2          ! Covariance matrix coefficient ({Delta}r)
      real*4        C3          ! Covariance matrix coefficient ({Delta}i)
      real*4        C4          ! Covariance matrix coefficient ({Delta}z)
      real*4        C5          ! Covariance matrix coefficient (U_ref_)
      real*4        C6          ! Covariance matrix coefficient (B_ref_)
      real*4        C7          ! Covariance matrix coefficient (V_ref_)
      real*4        C8          ! Covariance matrix coefficient (R_ref_)
      real*4        C9          ! Covariance matrix coefficient (I_ref_)

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

C  Declarations for 'tabled2.dat'	! Covariance matrix of the uncertainties affecting (1)
                             the differences between the tertiary MegaCam
                             magnitudes and the MegaCam magnitudes of BD+17 4708
                             and (2) the Landolt magnitudes of BD+17 4708

      integer*4 nr__9
      parameter (nr__9=9)	! Number of records
      character*123 ar__9  	! Full-size record

      real*4        C1_1        ! Covariance matrix coefficient ({Delta}g)
      real*4        C2_1        ! Covariance matrix coefficient ({Delta}r)
      real*4        C3_1        ! Covariance matrix coefficient ({Delta}i)
      real*4        C4_1        ! Covariance matrix coefficient ({Delta}z)
      real*4        C5_1        ! Covariance matrix coefficient (U_ref_)
      real*4        C6_1        ! Covariance matrix coefficient (B_ref_)
      real*4        C7_1        ! Covariance matrix coefficient (V_ref_)
      real*4        C8_1        ! Covariance matrix coefficient (R_ref_)
      real*4        C9_1        ! Covariance matrix coefficient (I_ref_)

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

C  Declarations for 'tabled3.dat'	! Covariance matrix of the uncertainties affecting the
                             synthetic instrumental magnitudes of BD+17 4708

      integer*4 nr__10
      parameter (nr__10=9)	! Number of records
      character*123 ar__10 	! Full-size record

      real*4        C1_2        ! Covariance matrix coefficient ({Delta}g)
      real*4        C2_2        ! Covariance matrix coefficient ({Delta}r)
      real*4        C3_2        ! Covariance matrix coefficient ({Delta}i)
      real*4        C4_2        ! Covariance matrix coefficient ({Delta}z)
      real*4        C5_2        ! Covariance matrix coefficient (U_ref_)
      real*4        C6_2        ! Covariance matrix coefficient (B_ref_)
      real*4        C7_2        ! Covariance matrix coefficient (V_ref_)
      real*4        C8_2        ! Covariance matrix coefficient (R_ref_)
      real*4        C9_2        ! Covariance matrix coefficient (I_ref_)

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

C  Declarations for 'table13.dat'	! Open transmission of the MegaCam imager (no filters)

      integer*4 nr__11
      parameter (nr__11=7190)	! Number of records
      character*15 ar__11 	! Full-size record

      real*8        WL          ! (0.1nm) Wavelength
      real*8        OT          ! [0,1] Open Transmission (no filters) (1)
*Note (1): product of (a) the transmission of the optical system, including
*     the camera window, (b) the mirror reflectivity (c) the average
*     transmission of the atmosphere above Mauna Kea and (c) the average
*     quantum efficiency of the MegaCam CCDs.

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

C  Declarations for 'table14.dat'	! u_M filter scan (Sagem/REOSC)

      integer*4 nr__12
      parameter (nr__12=1545)	! Number of records
      character*95 ar__12 	! Full-size record

      integer*4     WL_1        ! (0.1nm) Wavelength
      real*4        T000        ! [0,1] Transmission at the center of the filter
      real*4        T023        ! [0,1] Transmission at 23mm from the center
      real*4        T047        ! [0,1] Transmission at 47mm from the center
      real*4        T070        ! [0,1] Transmission at 70mm from the center
      real*4        T093        ! [0,1] Transmission at 93mm from the center
      real*4        T117        ! [0,1] Transmission at 117mm from the center
      real*4        T140        ! [0,1] Transmission at 140mm from the center
      real*4        T163        ! [0,1] Transmission at 163mm from the center
      real*4        T186        ! [0,1] Transmission at 186mm from the center
      real*4        T210        ! [0,1] Transmission at 210mm from the center

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

C  Declarations for 'table15.dat'	! g_M filter scan (Sagem/REOSC)

      integer*4 nr__13
      parameter (nr__13=371)	! Number of records
      character*95 ar__13 	! Full-size record

      integer*4     WL_2        ! (0.1nm) Wavelength
      real*4        T000_1      ! [0,1] Transmission at the center of the filter
      real*4        T023_1      ! [0,1] Transmission at 23mm from the center
      real*4        T047_1      ! [0,1] Transmission at 47mm from the center
      real*4        T070_1      ! [0,1] Transmission at 70mm from the center
      real*4        T093_1      ! [0,1] Transmission at 93mm from the center
      real*4        T117_1      ! [0,1] Transmission at 117mm from the center
      real*4        T140_1      ! [0,1] Transmission at 140mm from the center
      real*4        T163_1      ! [0,1] Transmission at 163mm from the center
      real*4        T186_1      ! [0,1] Transmission at 186mm from the center
      real*4        T210_1      ! [0,1] Transmission at 210mm from the center

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

C  Declarations for 'table16.dat'	! r_M filter scan (Sagem/REOSC)

      integer*4 nr__14
      parameter (nr__14=357)	! Number of records
      character*95 ar__14 	! Full-size record

      integer*4     WL_3        ! (0.1nm) Wavelength
      real*4        T000_2      ! [0,1] Transmission at the center of the filter
      real*4        T023_2      ! [0,1] Transmission at 23mm from the center
      real*4        T047_2      ! [0,1] Transmission at 47mm from the center
      real*4        T070_2      ! [0,1] Transmission at 70mm from the center
      real*4        T093_2      ! [0,1] Transmission at 93mm from the center
      real*4        T117_2      ! [0,1] Transmission at 117mm from the center
      real*4        T140_2      ! [0,1] Transmission at 140mm from the center
      real*4        T163_2      ! [0,1] Transmission at 163mm from the center
      real*4        T186_2      ! [0,1] Transmission at 186mm from the center
      real*4        T210_2      ! [0,1] Transmission at 210mm from the center

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

C  Declarations for 'table17.dat'	! i_M filter scan (Sagem/REOSC)

      integer*4 nr__15
      parameter (nr__15=446)	! Number of records
      character*95 ar__15 	! Full-size record

      integer*4     WL_4        ! (0.1nm) Wavelength
      real*4        T000_3      ! [0,1] Transmission at the center of the filter
      real*4        T023_3      ! [0,1] Transmission at 23mm from the center
      real*4        T047_3      ! [0,1] Transmission at 47mm from the center
      real*4        T070_3      ! [0,1] Transmission at 70mm from the center
      real*4        T093_3      ! [0,1] Transmission at 93mm from the center
      real*4        T117_3      ! [0,1] Transmission at 117mm from the center
      real*4        T140_3      ! [0,1] Transmission at 140mm from the center
      real*4        T163_3      ! [0,1] Transmission at 163mm from the center
      real*4        T186_3      ! [0,1] Transmission at 186mm from the center
      real*4        T210_3      ! [0,1] Transmission at 210mm from the center

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

C  Declarations for 'table18.dat'	! z_M filter scan (Sagem/REOSC)

      integer*4 nr__16
      parameter (nr__16=632)	! Number of records
      character*95 ar__16 	! Full-size record

      integer*4     WL_5        ! (0.1nm) Wavelength
      real*4        T000_4      ! [0,1] Transmission at the center of the filter
      real*4        T023_4      ! [0,1] Transmission at 23mm from the center
      real*4        T047_4      ! [0,1] Transmission at 47mm from the center
      real*4        T070_4      ! [0,1] Transmission at 70mm from the center
      real*4        T093_4      ! [0,1] Transmission at 93mm from the center
      real*4        T117_4      ! [0,1] Transmission at 117mm from the center
      real*4        T140_4      ! [0,1] Transmission at 140mm from the center
      real*4        T163_4      ! [0,1] Transmission at 163mm from the center
      real*4        T186_4      ! [0,1] Transmission at 186mm from the center
      real*4        T210_4      ! [0,1] Transmission at 210mm from the center

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

C  Loading file 'table19.dat'	! Tertiary rgiz catalog for the SNLS field D1

C  Format for file interpretation

    1 format(
     +  F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table19.dat')
      write(6,*) '....Loading file: table19.dat'
      do i__=1,2546
        read(1,'(A153)')ar__
        read(ar__,1)
     +  RAdeg,DEdeg,gmag,e_gmag,o_gmag,dkg,rmag,e_rmag,o_rmag,dkr,
     +  imag,e_imag,o_imag,dki,zmag,e_zmag,o_zmag,dkz
c    ..............Just test output...........
        write(6,1)
     +  RAdeg,DEdeg,gmag,e_gmag,o_gmag,dkg,rmag,e_rmag,o_rmag,dkr,
     +  imag,e_imag,o_imag,dki,zmag,e_zmag,o_zmag,dkz
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table20.dat'	! Tertiary rgiz catalog for the SNLS field D2

C  Format for file interpretation

    2 format(
     +  F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table20.dat')
      write(6,*) '....Loading file: table20.dat'
      do i__=1,3878
        read(1,'(A153)')ar__1
        read(ar__1,2)
     +  RAdeg_1,DEdeg_1,gmag_1,e_gmag_1,o_gmag_1,dkg_1,rmag_1,
     +  e_rmag_1,o_rmag_1,dkr_1,imag_1,e_imag_1,o_imag_1,dki_1,zmag_1,
     +  e_zmag_1,o_zmag_1,dkz_1
c    ..............Just test output...........
        write(6,2)
     +  RAdeg_1,DEdeg_1,gmag_1,e_gmag_1,o_gmag_1,dkg_1,rmag_1,
     +  e_rmag_1,o_rmag_1,dkr_1,imag_1,e_imag_1,o_imag_1,dki_1,zmag_1,
     +  e_zmag_1,o_zmag_1,dkz_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table21.dat'	! Tertiary rgiz catalog for the SNLS field D3

C  Format for file interpretation

    3 format(
     +  F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table21.dat')
      write(6,*) '....Loading file: table21.dat'
      do i__=1,2827
        read(1,'(A153)')ar__2
        read(ar__2,3)
     +  RAdeg_2,DEdeg_2,gmag_2,e_gmag_2,o_gmag_2,dkg_2,rmag_2,
     +  e_rmag_2,o_rmag_2,dkr_2,imag_2,e_imag_2,o_imag_2,dki_2,zmag_2,
     +  e_zmag_2,o_zmag_2,dkz_2
c    ..............Just test output...........
        write(6,3)
     +  RAdeg_2,DEdeg_2,gmag_2,e_gmag_2,o_gmag_2,dkg_2,rmag_2,
     +  e_rmag_2,o_rmag_2,dkr_2,imag_2,e_imag_2,o_imag_2,dki_2,zmag_2,
     +  e_zmag_2,o_zmag_2,dkz_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table22.dat'	! Tertiary rgiz catalog for the SNLS field D4

C  Format for file interpretation

    4 format(
     +  F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3,3X,F6.3,3X,F6.3,3X,
     +  I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table22.dat')
      write(6,*) '....Loading file: table22.dat'
      do i__=1,5294
        read(1,'(A153)')ar__3
        read(ar__3,4)
     +  RAdeg_3,DEdeg_3,gmag_3,e_gmag_3,o_gmag_3,dkg_3,rmag_3,
     +  e_rmag_3,o_rmag_3,dkr_3,imag_3,e_imag_3,o_imag_3,dki_3,zmag_3,
     +  e_zmag_3,o_zmag_3,dkz_3
c    ..............Just test output...........
        write(6,4)
     +  RAdeg_3,DEdeg_3,gmag_3,e_gmag_3,o_gmag_3,dkg_3,rmag_3,
     +  e_rmag_3,o_rmag_3,dkr_3,imag_3,e_imag_3,o_imag_3,dki_3,zmag_3,
     +  e_zmag_3,o_zmag_3,dkz_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'd1u.dat'	! D1 tertiary stars u-band magnitudes

C  Format for file interpretation

    5 format(F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'d1u.dat')
      write(6,*) '....Loading file: d1u.dat'
      do i__=1,1394
        read(1,'(A54)')ar__4
        read(ar__4,5)RAdeg_4,DEdeg_4,umag,e_umag,o_umag,dku
c    ..............Just test output...........
        write(6,5)RAdeg_4,DEdeg_4,umag,e_umag,o_umag,dku
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'd2u.dat'	! D2 tertiary stars u-band magnitudes

C  Format for file interpretation

    6 format(F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'d2u.dat')
      write(6,*) '....Loading file: d2u.dat'
      do i__=1,2349
        read(1,'(A54)')ar__5
        read(ar__5,6)RAdeg_5,DEdeg_5,umag_1,e_umag_1,o_umag_1,dku_1
c    ..............Just test output...........
        write(6,6)RAdeg_5,DEdeg_5,umag_1,e_umag_1,o_umag_1,dku_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'd3u.dat'	! D3 tertiary stars u-band magnitudes

C  Format for file interpretation

    7 format(F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'d3u.dat')
      write(6,*) '....Loading file: d3u.dat'
      do i__=1,714
        read(1,'(A54)')ar__6
        read(ar__6,7)RAdeg_6,DEdeg_6,umag_2,e_umag_2,o_umag_2,dku_2
c    ..............Just test output...........
        write(6,7)RAdeg_6,DEdeg_6,umag_2,e_umag_2,o_umag_2,dku_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'd4u.dat'	! D4 tertiary stars u-band magnitudes

C  Format for file interpretation

    8 format(F9.5,3X,F9.5,3X,F6.3,3X,F6.3,3X,I3,3X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'d4u.dat')
      write(6,*) '....Loading file: d4u.dat'
      do i__=1,3051
        read(1,'(A54)')ar__7
        read(ar__7,8)RAdeg_7,DEdeg_7,umag_3,e_umag_3,o_umag_3,dku_3
c    ..............Just test output...........
        write(6,8)RAdeg_7,DEdeg_7,umag_3,e_umag_3,o_umag_3,dku_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tabled1.dat'	! Full covariance matrix of the BD+17 4708 MegaCam
*                               and Landolt magnitude uncertainties

C  Format for file interpretation

    9 format(
     +  E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,
     +  3X,E11.5,3X,E11.5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tabled1.dat')
      write(6,*) '....Loading file: tabled1.dat'
      do i__=1,9
        read(1,'(A123)')ar__8
        read(ar__8,9)C1,C2,C3,C4,C5,C6,C7,C8,C9
c    ..............Just test output...........
        write(6,9)C1,C2,C3,C4,C5,C6,C7,C8,C9
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tabled2.dat'	! Covariance matrix of the uncertainties affecting (1)
*                             the differences between the tertiary MegaCam
*                             magnitudes and the MegaCam magnitudes of BD+17 4708
*                             and (2) the Landolt magnitudes of BD+17 4708

C  Format for file interpretation

   10 format(
     +  E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,
     +  3X,E11.5,3X,E11.5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tabled2.dat')
      write(6,*) '....Loading file: tabled2.dat'
      do i__=1,9
        read(1,'(A123)')ar__9
        read(ar__9,10)C1_1,C2_1,C3_1,C4_1,C5_1,C6_1,C7_1,C8_1,C9_1
c    ..............Just test output...........
        write(6,10)C1_1,C2_1,C3_1,C4_1,C5_1,C6_1,C7_1,C8_1,C9_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tabled3.dat'	! Covariance matrix of the uncertainties affecting the
*                             synthetic instrumental magnitudes of BD+17 4708

C  Format for file interpretation

   11 format(
     +  E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5,
     +  3X,E11.5,3X,E11.5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tabled3.dat')
      write(6,*) '....Loading file: tabled3.dat'
      do i__=1,9
        read(1,'(A123)')ar__10
        read(ar__10,11)C1_2,C2_2,C3_2,C4_2,C5_2,C6_2,C7_2,C8_2,C9_2
c    ..............Just test output...........
        write(6,11)C1_2,C2_2,C3_2,C4_2,C5_2,C6_2,C7_2,C8_2,C9_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table13.dat'	! Open transmission of the MegaCam imager (no filters)

C  Format for file interpretation

   12 format(F7.1,1X,F7.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table13.dat')
      write(6,*) '....Loading file: table13.dat'
      do i__=1,7190
        read(1,'(A15)')ar__11
        read(ar__11,12)WL,OT
c    ..............Just test output...........
        write(6,12)WL,OT
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table14.dat'	! u_M filter scan (Sagem/REOSC)

C  Format for file interpretation

   13 format(
     +  I5,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,
     +  F6.4,3X,F6.4,3X,F6.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table14.dat')
      write(6,*) '....Loading file: table14.dat'
      do i__=1,1545
        read(1,'(A95)')ar__12
        read(ar__12,13)
     +  WL_1,T000,T023,T047,T070,T093,T117,T140,T163,T186,T210
c    ..............Just test output...........
        write(6,13)
     +  WL_1,T000,T023,T047,T070,T093,T117,T140,T163,T186,T210
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table15.dat'	! g_M filter scan (Sagem/REOSC)

C  Format for file interpretation

   14 format(
     +  I5,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,
     +  F6.4,3X,F6.4,3X,F6.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table15.dat')
      write(6,*) '....Loading file: table15.dat'
      do i__=1,371
        read(1,'(A95)')ar__13
        read(ar__13,14)
     +  WL_2,T000_1,T023_1,T047_1,T070_1,T093_1,T117_1,T140_1,T163_1,
     +  T186_1,T210_1
c    ..............Just test output...........
        write(6,14)
     +  WL_2,T000_1,T023_1,T047_1,T070_1,T093_1,T117_1,T140_1,T163_1,
     +  T186_1,T210_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table16.dat'	! r_M filter scan (Sagem/REOSC)

C  Format for file interpretation

   15 format(
     +  I5,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,
     +  F6.4,3X,F6.4,3X,F6.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table16.dat')
      write(6,*) '....Loading file: table16.dat'
      do i__=1,357
        read(1,'(A95)')ar__14
        read(ar__14,15)
     +  WL_3,T000_2,T023_2,T047_2,T070_2,T093_2,T117_2,T140_2,T163_2,
     +  T186_2,T210_2
c    ..............Just test output...........
        write(6,15)
     +  WL_3,T000_2,T023_2,T047_2,T070_2,T093_2,T117_2,T140_2,T163_2,
     +  T186_2,T210_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table17.dat'	! i_M filter scan (Sagem/REOSC)

C  Format for file interpretation

   16 format(
     +  I5,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,
     +  F6.4,3X,F6.4,3X,F6.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table17.dat')
      write(6,*) '....Loading file: table17.dat'
      do i__=1,446
        read(1,'(A95)')ar__15
        read(ar__15,16)
     +  WL_4,T000_3,T023_3,T047_3,T070_3,T093_3,T117_3,T140_3,T163_3,
     +  T186_3,T210_3
c    ..............Just test output...........
        write(6,16)
     +  WL_4,T000_3,T023_3,T047_3,T070_3,T093_3,T117_3,T140_3,T163_3,
     +  T186_3,T210_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table18.dat'	! z_M filter scan (Sagem/REOSC)

C  Format for file interpretation

   17 format(
     +  I5,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,F6.4,3X,
     +  F6.4,3X,F6.4,3X,F6.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table18.dat')
      write(6,*) '....Loading file: table18.dat'
      do i__=1,632
        read(1,'(A95)')ar__16
        read(ar__16,17)
     +  WL_5,T000_4,T023_4,T047_4,T070_4,T093_4,T117_4,T140_4,T163_4,
     +  T186_4,T210_4
c    ..............Just test output...........
        write(6,17)
     +  WL_5,T000_4,T023_4,T047_4,T070_4,T093_4,T117_4,T140_4,T163_4,
     +  T186_4,T210_4
c    .......End.of.Just test output...........
      end do
      close(1)

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