FORTRAN Generation
(/./ftp/cats/J/AJ/125/1817)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/125/1817 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-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/AJ/125/1817       SDSS Early-Type Galaxies Catalog         (Bernardi+, 2003)
*================================================================================
*Early-type galaxies in the Sloan Digital Sky Survey. I. The sample.
*    Bernardi M., Sheth R.K., Annis J., Burles S., Eisenstein D.J.,
*    Finkbeiner D.P., Hogg D.W., Lupton R.H., Schlegel D.J., Subbarao M.,
*    Bahcall N.A., Blakeslee J.P., Brinkmann J., Castander F.J., Connolly A.J.,
*    Csabai I., Doi M., Fukugita M., Frieman J., Heckman T., Hennessy G.S.,
*    Ivezic Z., Knapp G.R., Lamb D.Q., McKay T., Munn J.A., Nichol R.,
*    Okamura S., Schneider D.P., Thakar A.R., York D.G.
*   <Astron. J. 125, 1817 (2003)>
*   =2003AJ....125.1817B
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2a.dat'	! SDSS Early-Type Galaxies (SDSS ETG) Catalog:
                                                  Observed Parameters (r* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg      (nr__) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees (J2000)
      real*8        z          (nr__) ! Redshift
      real*8        e_z        (nr__) ! Error estimate in z
      integer*4     S_N        (nr__) ! Signal-to-noise ratio of the spectrum
      character*2   Band       (nr__) ! [r* g* i* z*] Band
      real*4        rDev       (nr__) ! (arcsec) The Band apparent half-light radius (1)
      real*4        e_rDev     (nr__) ! (arcsec) Error estimate in rDev
      real*4        magDev     (nr__) ! (mag) de Vaucouleurs Band apparent magnitude (1)
      real*4        e_magDev   (nr__) ! (mag) Error estimate in magDev
      real*4        magMod     (nr__) ! (mag) Model Band magnitude
      real*4        e_magMod   (nr__) ! (mag) Error estimate in magMod
      real*4        b_a        (nr__) ! The Band axis ratio
      real*4        e_b_a      (nr__) ! Error estimate in b/a
*Note (1): From the r^1/4^ fit.

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

C  Declarations for 'table2b.dat'	! SDSS ETG Catalog: Observed Parameters (g* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_1    (nr__1) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination in decimal degrees (J2000)
      real*8        z_1        (nr__1) ! Redshift
      real*8        e_z_1      (nr__1) ! Error estimate in z
      integer*4     S_N_1      (nr__1) ! Signal-to-noise ratio of the spectrum
      character*2   Band_1     (nr__1) ! [r* g* i* z*] Band
      real*4        rDev_1     (nr__1) ! (arcsec) The Band apparent half-light radius (1)
      real*4        e_rDev_1   (nr__1) ! (arcsec) Error estimate in rDev
      real*4        magDev_1   (nr__1) ! (mag) de Vaucouleurs Band apparent magnitude (1)
      real*4        e_magDev_1 (nr__1) ! (mag) Error estimate in magDev
      real*4        magMod_1   (nr__1) ! (mag) Model Band magnitude
      real*4        e_magMod_1 (nr__1) ! (mag) Error estimate in magMod
      real*4        b_a_1      (nr__1) ! The Band axis ratio
      real*4        e_b_a_1    (nr__1) ! Error estimate in b/a
*Note (1): From the r^1/4^ fit.

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

C  Declarations for 'table2c.dat'	! SDSS ETG Catalog: Observed Parameters (i* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_2    (nr__2) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination in decimal degrees (J2000)
      real*8        z_2        (nr__2) ! Redshift
      real*8        e_z_2      (nr__2) ! Error estimate in z
      integer*4     S_N_2      (nr__2) ! Signal-to-noise ratio of the spectrum
      character*2   Band_2     (nr__2) ! [r* g* i* z*] Band
      real*4        rDev_2     (nr__2) ! (arcsec) The Band apparent half-light radius (1)
      real*4        e_rDev_2   (nr__2) ! (arcsec) Error estimate in rDev
      real*4        magDev_2   (nr__2) ! (mag) de Vaucouleurs Band apparent magnitude (1)
      real*4        e_magDev_2 (nr__2) ! (mag) Error estimate in magDev
      real*4        magMod_2   (nr__2) ! (mag) Model Band magnitude
      real*4        e_magMod_2 (nr__2) ! (mag) Error estimate in magMod
      real*4        b_a_2      (nr__2) ! The Band axis ratio
      real*4        e_b_a_2    (nr__2) ! Error estimate in b/a
*Note (1): From the r^1/4^ fit.

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

C  Declarations for 'table2d.dat'	! SDSS ETG Catalog: Observed Parameters (z* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_3    (nr__3) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_3    (nr__3) ! (deg) Declination in decimal degrees (J2000)
      real*8        z_3        (nr__3) ! Redshift
      real*8        e_z_3      (nr__3) ! Error estimate in z
      integer*4     S_N_3      (nr__3) ! Signal-to-noise ratio of the spectrum
      character*2   Band_3     (nr__3) ! [r* g* i* z*] Band
      real*4        rDev_3     (nr__3) ! (arcsec) The Band apparent half-light radius (1)
      real*4        e_rDev_3   (nr__3) ! (arcsec) Error estimate in rDev
      real*4        magDev_3   (nr__3) ! (mag) de Vaucouleurs Band apparent magnitude (1)
      real*4        e_magDev_3 (nr__3) ! (mag) Error estimate in magDev
      real*4        magMod_3   (nr__3) ! (mag) Model Band magnitude
      real*4        e_magMod_3 (nr__3) ! (mag) Error estimate in magMod
      real*4        b_a_3      (nr__3) ! The Band axis ratio
      real*4        e_b_a_3    (nr__3) ! Error estimate in b/a
*Note (1): From the r^1/4^ fit.

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

C  Declarations for 'table3a.dat'	! SDSS ETG Catalog: Physical Parameters (r* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_4    (nr__4) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_4    (nr__4) ! (deg) Declination in decimal degrees (J2000)
      real*4        logVDis    (nr__4) ! ([km/s]) Log of the aperture-corrected velocity
*                                       dispersion
      real*4        e_logVDis  (nr__4) ! ([km/s]) Error estimate in logVDis
      character*2   Band_4     (nr__4) ! [r* g* i* z*] Band
      real*4        logRad     (nr__4) ! ([kpc]) Log of Band physical half-light radius (1)
      real*4        e_logRad   (nr__4) ! ([kpc]) Error estimate in logRad  (1)
      real*4        SB         (nr__4) ! (mag/arcsec2) Band surface brightness (2)
      real*4        e_SB       (nr__4) ! (mag/arcsec2) Error estimate in SB
      real*8        MagDev_4   (nr__4) ! (mag) de Vaucouleurs Band absolute magnitude
      real*4        e_MagDev_4 (nr__4) ! (mag) Error estimate in MagDev
      real*8        MagMod_4   (nr__4) ! (mag) Model Band absolute magnitude
      real*4        e_MagMod_4 (nr__4) ! (mag) Error estimate in MagMod
      real*4        Kcorr      (nr__4) ! (mag) The band K-correction
      integer*4     Nden       (nr__4) ! Number of near neighbors used as an
*                                       estimate of the local density (3)
*Note (1): In units of kpc*h_70_^-1^.
*Note (2): Corrected for the cosmological (1+z)^4^ dimming.
*Note (3): 1001 = An object for which the number of neighbors was not determined;
*           102 = An object with >=100 neighbors.

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

C  Declarations for 'table3b.dat'	! SDSS ETG Catalog: Physical Parameters (g* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_5    (nr__5) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_5    (nr__5) ! (deg) Declination in decimal degrees (J2000)
      real*4        logVDis_1  (nr__5) ! ([km/s]) Log of the aperture-corrected velocity
*                                       dispersion
      real*4        e_logVDis_1(nr__5) ! ([km/s]) Error estimate in logVDis
      character*2   Band_5     (nr__5) ! [r* g* i* z*] Band
      real*4        logRad_1   (nr__5) ! ([kpc]) Log of Band physical half-light radius (1)
      real*4        e_logRad_1 (nr__5) ! ([kpc]) Error estimate in logRad  (1)
      real*4        SB_1       (nr__5) ! (mag/arcsec2) Band surface brightness (2)
      real*4        e_SB_1     (nr__5) ! (mag/arcsec2) Error estimate in SB
      real*8        MagDev_5   (nr__5) ! (mag) de Vaucouleurs Band absolute magnitude
      real*4        e_MagDev_5 (nr__5) ! (mag) Error estimate in MagDev
      real*8        MagMod_5   (nr__5) ! (mag) Model Band absolute magnitude
      real*4        e_MagMod_5 (nr__5) ! (mag) Error estimate in MagMod
      real*4        Kcorr_1    (nr__5) ! (mag) The band K-correction
      integer*4     Nden_1     (nr__5) ! Number of near neighbors used as an
*                                       estimate of the local density (3)
*Note (1): In units of kpc*h_70_^-1^.
*Note (2): Corrected for the cosmological (1+z)^4^ dimming.
*Note (3): 1001 = An object for which the number of neighbors was not determined;
*           102 = An object with >=100 neighbors.

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

C  Declarations for 'table3c.dat'	! SDSS ETG Catalog: Physical Parameters (i* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_6    (nr__6) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_6    (nr__6) ! (deg) Declination in decimal degrees (J2000)
      real*4        logVDis_2  (nr__6) ! ([km/s]) Log of the aperture-corrected velocity
*                                       dispersion
      real*4        e_logVDis_2(nr__6) ! ([km/s]) Error estimate in logVDis
      character*2   Band_6     (nr__6) ! [r* g* i* z*] Band
      real*4        logRad_2   (nr__6) ! ([kpc]) Log of Band physical half-light radius (1)
      real*4        e_logRad_2 (nr__6) ! ([kpc]) Error estimate in logRad  (1)
      real*4        SB_2       (nr__6) ! (mag/arcsec2) Band surface brightness (2)
      real*4        e_SB_2     (nr__6) ! (mag/arcsec2) Error estimate in SB
      real*8        MagDev_6   (nr__6) ! (mag) de Vaucouleurs Band absolute magnitude
      real*4        e_MagDev_6 (nr__6) ! (mag) Error estimate in MagDev
      real*8        MagMod_6   (nr__6) ! (mag) Model Band absolute magnitude
      real*4        e_MagMod_6 (nr__6) ! (mag) Error estimate in MagMod
      real*4        Kcorr_2    (nr__6) ! (mag) The band K-correction
      integer*4     Nden_2     (nr__6) ! Number of near neighbors used as an
*                                       estimate of the local density (3)
*Note (1): In units of kpc*h_70_^-1^.
*Note (2): Corrected for the cosmological (1+z)^4^ dimming.
*Note (3): 1001 = An object for which the number of neighbors was not determined;
*           102 = An object with >=100 neighbors.

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

C  Declarations for 'table3d.dat'	! SDSS ETG Catalog: Physical Parameters (z* Band)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_7    (nr__7) ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg_7    (nr__7) ! (deg) Declination in decimal degrees (J2000)
      real*4        logVDis_3  (nr__7) ! ([km/s]) Log of the aperture-corrected velocity
*                                       dispersion
      real*4        e_logVDis_3(nr__7) ! ([km/s]) Error estimate in logVDis
      character*2   Band_7     (nr__7) ! [r* g* i* z*] Band
      real*4        logRad_3   (nr__7) ! ([kpc]) Log of Band physical half-light radius (1)
      real*4        e_logRad_3 (nr__7) ! ([kpc]) Error estimate in logRad  (1)
      real*4        SB_3       (nr__7) ! (mag/arcsec2) Band surface brightness (2)
      real*4        e_SB_3     (nr__7) ! (mag/arcsec2) Error estimate in SB
      real*8        MagDev_7   (nr__7) ! (mag) de Vaucouleurs Band absolute magnitude
      real*4        e_MagDev_7 (nr__7) ! (mag) Error estimate in MagDev
      real*8        MagMod_7   (nr__7) ! (mag) Model Band absolute magnitude
      real*4        e_MagMod_7 (nr__7) ! (mag) Error estimate in MagMod
      real*4        Kcorr_3    (nr__7) ! (mag) The band K-correction
      integer*4     Nden_3     (nr__7) ! Number of near neighbors used as an
*                                       estimate of the local density (3)
*Note (1): In units of kpc*h_70_^-1^.
*Note (2): Corrected for the cosmological (1+z)^4^ dimming.
*Note (3): 1001 = An object for which the number of neighbors was not determined;
*           102 = An object with >=100 neighbors.

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

C  Loading file 'table2a.dat'	! SDSS Early-Type Galaxies (SDSS ETG) Catalog:
*                                                  Observed Parameters (r* Band)

C  Format for file interpretation

    1 format(
     +  F10.6,2X,F9.6,1X,F7.5,1X,F7.5,2X,I2,1X,A2,1X,F5.2,2X,F4.2,2X,
     +  F6.3,3X,F5.3,2X,F6.3,3X,F5.3,1X,F5.3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2a.dat')
      write(6,*) '....Loading file: table2a.dat'
      do i__=1,8666
        read(1,'(A100)')ar__
        read(ar__,1)
     +  RAdeg(i__),DEdeg(i__),z(i__),e_z(i__),S_N(i__),Band(i__),
     +  rDev(i__),e_rDev(i__),magDev(i__),e_magDev(i__),magMod(i__),
     +  e_magMod(i__),b_a(i__),e_b_a(i__)
c    ..............Just test output...........
        write(6,1)
     +  RAdeg(i__),DEdeg(i__),z(i__),e_z(i__),S_N(i__),Band(i__),
     +  rDev(i__),e_rDev(i__),magDev(i__),e_magDev(i__),magMod(i__),
     +  e_magMod(i__),b_a(i__),e_b_a(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2b.dat'	! SDSS ETG Catalog: Observed Parameters (g* Band)

C  Format for file interpretation

    2 format(
     +  F10.6,2X,F9.6,1X,F7.5,1X,F7.5,2X,I2,1X,A2,1X,F5.2,2X,F4.2,2X,
     +  F6.3,3X,F5.3,2X,F6.3,3X,F5.3,1X,F5.3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2b.dat')
      write(6,*) '....Loading file: table2b.dat'
      do i__=1,8666
        read(1,'(A100)')ar__1
        read(ar__1,2)
     +  RAdeg_1(i__),DEdeg_1(i__),z_1(i__),e_z_1(i__),S_N_1(i__),
     +  Band_1(i__),rDev_1(i__),e_rDev_1(i__),magDev_1(i__),
     +  e_magDev_1(i__),magMod_1(i__),e_magMod_1(i__),b_a_1(i__),
     +  e_b_a_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  RAdeg_1(i__),DEdeg_1(i__),z_1(i__),e_z_1(i__),S_N_1(i__),
     +  Band_1(i__),rDev_1(i__),e_rDev_1(i__),magDev_1(i__),
     +  e_magDev_1(i__),magMod_1(i__),e_magMod_1(i__),b_a_1(i__),
     +  e_b_a_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2c.dat'	! SDSS ETG Catalog: Observed Parameters (i* Band)

C  Format for file interpretation

    3 format(
     +  F10.6,2X,F9.6,1X,F7.5,1X,F7.5,2X,I2,1X,A2,1X,F5.2,2X,F4.2,2X,
     +  F6.3,3X,F5.3,2X,F6.3,3X,F5.3,1X,F5.3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2c.dat')
      write(6,*) '....Loading file: table2c.dat'
      do i__=1,8666
        read(1,'(A100)')ar__2
        read(ar__2,3)
     +  RAdeg_2(i__),DEdeg_2(i__),z_2(i__),e_z_2(i__),S_N_2(i__),
     +  Band_2(i__),rDev_2(i__),e_rDev_2(i__),magDev_2(i__),
     +  e_magDev_2(i__),magMod_2(i__),e_magMod_2(i__),b_a_2(i__),
     +  e_b_a_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  RAdeg_2(i__),DEdeg_2(i__),z_2(i__),e_z_2(i__),S_N_2(i__),
     +  Band_2(i__),rDev_2(i__),e_rDev_2(i__),magDev_2(i__),
     +  e_magDev_2(i__),magMod_2(i__),e_magMod_2(i__),b_a_2(i__),
     +  e_b_a_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2d.dat'	! SDSS ETG Catalog: Observed Parameters (z* Band)

C  Format for file interpretation

    4 format(
     +  F10.6,2X,F9.6,1X,F7.5,1X,F7.5,2X,I2,1X,A2,1X,F5.2,2X,F4.2,2X,
     +  F6.3,3X,F5.3,2X,F6.3,3X,F5.3,1X,F5.3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2d.dat')
      write(6,*) '....Loading file: table2d.dat'
      do i__=1,8666
        read(1,'(A100)')ar__3
        read(ar__3,4)
     +  RAdeg_3(i__),DEdeg_3(i__),z_3(i__),e_z_3(i__),S_N_3(i__),
     +  Band_3(i__),rDev_3(i__),e_rDev_3(i__),magDev_3(i__),
     +  e_magDev_3(i__),magMod_3(i__),e_magMod_3(i__),b_a_3(i__),
     +  e_b_a_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  RAdeg_3(i__),DEdeg_3(i__),z_3(i__),e_z_3(i__),S_N_3(i__),
     +  Band_3(i__),rDev_3(i__),e_rDev_3(i__),magDev_3(i__),
     +  e_magDev_3(i__),magMod_3(i__),e_magMod_3(i__),b_a_3(i__),
     +  e_b_a_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3a.dat'	! SDSS ETG Catalog: Physical Parameters (r* Band)

C  Format for file interpretation

    5 format(
     +  F10.6,2X,F9.6,2X,F5.3,2X,F5.3,1X,A2,1X,F6.3,2X,F5.3,1X,F6.3,
     +  2X,F5.3,1X,F7.3,2X,F6.3,1X,F7.3,2X,F6.3,3X,F5.3,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3a.dat')
      write(6,*) '....Loading file: table3a.dat'
      do i__=1,8666
        read(1,'(A110)')ar__4
        read(ar__4,5)
     +  RAdeg_4(i__),DEdeg_4(i__),logVDis(i__),e_logVDis(i__),
     +  Band_4(i__),logRad(i__),e_logRad(i__),SB(i__),e_SB(i__),
     +  MagDev_4(i__),e_MagDev_4(i__),MagMod_4(i__),e_MagMod_4(i__),
     +  Kcorr(i__),Nden(i__)
c    ..............Just test output...........
        write(6,5)
     +  RAdeg_4(i__),DEdeg_4(i__),logVDis(i__),e_logVDis(i__),
     +  Band_4(i__),logRad(i__),e_logRad(i__),SB(i__),e_SB(i__),
     +  MagDev_4(i__),e_MagDev_4(i__),MagMod_4(i__),e_MagMod_4(i__),
     +  Kcorr(i__),Nden(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3b.dat'	! SDSS ETG Catalog: Physical Parameters (g* Band)

C  Format for file interpretation

    6 format(
     +  F10.6,2X,F9.6,2X,F5.3,2X,F5.3,1X,A2,1X,F6.3,2X,F5.3,1X,F6.3,
     +  2X,F5.3,1X,F7.3,2X,F6.3,1X,F7.3,2X,F6.3,3X,F5.3,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3b.dat')
      write(6,*) '....Loading file: table3b.dat'
      do i__=1,8666
        read(1,'(A110)')ar__5
        read(ar__5,6)
     +  RAdeg_5(i__),DEdeg_5(i__),logVDis_1(i__),e_logVDis_1(i__),
     +  Band_5(i__),logRad_1(i__),e_logRad_1(i__),SB_1(i__),
     +  e_SB_1(i__),MagDev_5(i__),e_MagDev_5(i__),MagMod_5(i__),
     +  e_MagMod_5(i__),Kcorr_1(i__),Nden_1(i__)
c    ..............Just test output...........
        write(6,6)
     +  RAdeg_5(i__),DEdeg_5(i__),logVDis_1(i__),e_logVDis_1(i__),
     +  Band_5(i__),logRad_1(i__),e_logRad_1(i__),SB_1(i__),
     +  e_SB_1(i__),MagDev_5(i__),e_MagDev_5(i__),MagMod_5(i__),
     +  e_MagMod_5(i__),Kcorr_1(i__),Nden_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3c.dat'	! SDSS ETG Catalog: Physical Parameters (i* Band)

C  Format for file interpretation

    7 format(
     +  F10.6,2X,F9.6,2X,F5.3,2X,F5.3,1X,A2,1X,F6.3,2X,F5.3,1X,F6.3,
     +  2X,F5.3,1X,F7.3,2X,F6.3,1X,F7.3,2X,F6.3,3X,F5.3,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3c.dat')
      write(6,*) '....Loading file: table3c.dat'
      do i__=1,8666
        read(1,'(A110)')ar__6
        read(ar__6,7)
     +  RAdeg_6(i__),DEdeg_6(i__),logVDis_2(i__),e_logVDis_2(i__),
     +  Band_6(i__),logRad_2(i__),e_logRad_2(i__),SB_2(i__),
     +  e_SB_2(i__),MagDev_6(i__),e_MagDev_6(i__),MagMod_6(i__),
     +  e_MagMod_6(i__),Kcorr_2(i__),Nden_2(i__)
c    ..............Just test output...........
        write(6,7)
     +  RAdeg_6(i__),DEdeg_6(i__),logVDis_2(i__),e_logVDis_2(i__),
     +  Band_6(i__),logRad_2(i__),e_logRad_2(i__),SB_2(i__),
     +  e_SB_2(i__),MagDev_6(i__),e_MagDev_6(i__),MagMod_6(i__),
     +  e_MagMod_6(i__),Kcorr_2(i__),Nden_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3d.dat'	! SDSS ETG Catalog: Physical Parameters (z* Band)

C  Format for file interpretation

    8 format(
     +  F10.6,2X,F9.6,2X,F5.3,2X,F5.3,1X,A2,1X,F6.3,2X,F5.3,1X,F6.3,
     +  2X,F5.3,1X,F7.3,2X,F6.3,1X,F7.3,2X,F6.3,3X,F5.3,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3d.dat')
      write(6,*) '....Loading file: table3d.dat'
      do i__=1,8666
        read(1,'(A110)')ar__7
        read(ar__7,8)
     +  RAdeg_7(i__),DEdeg_7(i__),logVDis_3(i__),e_logVDis_3(i__),
     +  Band_7(i__),logRad_3(i__),e_logRad_3(i__),SB_3(i__),
     +  e_SB_3(i__),MagDev_7(i__),e_MagDev_7(i__),MagMod_7(i__),
     +  e_MagMod_7(i__),Kcorr_3(i__),Nden_3(i__)
c    ..............Just test output...........
        write(6,8)
     +  RAdeg_7(i__),DEdeg_7(i__),logVDis_3(i__),e_logVDis_3(i__),
     +  Band_7(i__),logRad_3(i__),e_logRad_3(i__),SB_3(i__),
     +  e_SB_3(i__),MagDev_7(i__),e_MagDev_7(i__),MagMod_7(i__),
     +  e_MagMod_7(i__),Kcorr_3(i__),Nden_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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