FORTRAN Generation
(/./ftp/cats/J/A_A/705/A84)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/705/A84 into FORTRAN code for loading all data files into arrays.

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



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Jun-08
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/705/A84             Galaxy size-mass relation                (Liao+, 2026)
*================================================================================
*Environmental invariance of the Galaxy size-mass relation.
*    Liao L.-W., Cooper A.
*    <Astron. Astrophys. 705, A84 (2026)>
*    =2026A&A...705A..84L        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'fig4a.dat'	! Spectroscopic-redshift sample,
                                 circular R50 sizes, in a given magnitude

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

      character*2   Sample     (nr__) ! [SZ PZ] Sample name (G1)
      character*6   Param      (nr__) ! Parameter name (G2)
      real*4        Mbin       (nr__) ! ([mag]) Stellar-magnitude bin (log10(rMAG))
      real*8        logmuall   (nr__) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall (nr__) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured   (nr__) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured (nr__) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue  (nr__) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue(nr__) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell   (nr__) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell (nr__) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi   (nr__) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi (nr__) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig4b.dat'	! Spectroscopic-redshift sample,
                                 major-axis sizes, in a given magnitude

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

      character*2   Sample_1   (nr__1) ! [SZ PZ] Sample name (G1)
      character*6   Param_1    (nr__1) ! Parameter name (G2)
      real*4        Mbin_1     (nr__1) ! ([mag]) Stellar-magnitude bin (log10(rMAG))
      real*8        logmuall_1 (nr__1) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall_1(nr__1) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured_1 (nr__1) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured_1(nr__1) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue_1(nr__1) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue_1(nr__1) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell_1 (nr__1) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell_1(nr__1) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi_1 (nr__1) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi_1(nr__1) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig4c.dat'	! Photometric-redshift sample,
                                 circular R50 sizes, in a given magnitude

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

      character*2   Sample_2   (nr__2) ! [SZ PZ] Sample name (G1)
      character*6   Param_2    (nr__2) ! Parameter name (G2)
      real*4        Mbin_2     (nr__2) ! ([mag]) Stellar-magnitude bin (log10(rMAG))
      real*8        logmuall_2 (nr__2) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall_2(nr__2) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured_2 (nr__2) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured_2(nr__2) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue_2(nr__2) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue_2(nr__2) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell_2 (nr__2) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell_2(nr__2) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi_2 (nr__2) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi_2(nr__2) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig4d.dat'	! Photometric-redshift sample,
                                 major-axis sizes, in a given magnitude

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

      character*2   Sample_3   (nr__3) ! [SZ PZ] Sample name (G1)
      character*6   Param_3    (nr__3) ! Parameter name (G2)
      real*4        Mbin_3     (nr__3) ! ([mag]) Stellar-magnitude bin (log10(rMAG))
      real*8        logmuall_3 (nr__3) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall_3(nr__3) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured_3 (nr__3) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured_3(nr__3) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue_3(nr__3) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue_3(nr__3) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell_3 (nr__3) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell_3(nr__3) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi_3 (nr__3) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi_3(nr__3) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig5a.dat'	! Spectroscopic-redshift sample,
                                 circularized half-light radii,
                                 in a given stellar mass bin

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

      character*2   Sample_4   (nr__4) ! [SZ PZ] Sample name (G1)
      character*4   Param_4    (nr__4) ! Parameter name (G2)
      real*4        sm         (nr__4) ! ([Msun]) Log10 stellar mass (bin center)
      real*8        logmuall_4 (nr__4) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall_4(nr__4) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured_4 (nr__4) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured_4(nr__4) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue_4(nr__4) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue_4(nr__4) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell_4 (nr__4) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell_4(nr__4) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi_4 (nr__4) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi_4(nr__4) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig5b.dat'	! Spectroscopic-redshift sample,
                                 major-axis half-light radii,
                                 in a given stellar mass bin

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

      character*2   Sample_5   (nr__5) ! [SZ PZ] Sample name (G1)
      character*4   Param_5    (nr__5) ! Parameter name (G2)
      real*4        sm_1       (nr__5) ! ([Msun]) Log10 stellar mass (bin center)
      real*8        logmuall_5 (nr__5) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall_5(nr__5) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured_5 (nr__5) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured_5(nr__5) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue_5(nr__5) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue_5(nr__5) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell_5 (nr__5) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell_5(nr__5) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi_5 (nr__5) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi_5(nr__5) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig5c.dat'	! Photometric-redshift sample,
                                 circularized half-light radii,
                                 in a given stellar mass bin

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

      character*2   Sample_6   (nr__6) ! [SZ PZ] Sample name (G1)
      character*4   Param_6    (nr__6) ! Parameter name (G2)
      real*4        sm_2       (nr__6) ! ([Msun]) Log10 stellar mass (bin center)
      real*8        logmuall_6 (nr__6) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall_6(nr__6) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured_6 (nr__6) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured_6(nr__6) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue_6(nr__6) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue_6(nr__6) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell_6 (nr__6) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell_6(nr__6) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi_6 (nr__6) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi_6(nr__6) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig5d.dat'	! Photometric-redshift sample,
                                 major-axis half-light radii,
                                 in a given stellar mass bin

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

      character*2   Sample_7   (nr__7) ! [SZ PZ] Sample name (G1)
      character*4   Param_7    (nr__7) ! Parameter name (G2)
      real*4        sm_3       (nr__7) ! ([Msun]) Log10 stellar mass (bin center)
      real*8        logmuall_7 (nr__7) ! ([kpc]) ?=-99 Param median size for all galaxies
      real*8        e_logmuall_7(nr__7) ! ([kpc]) ?=-99 Param 1-sigma error for all galaxies
      real*8        logmured_7 (nr__7) ! ([kpc]) ?=-99 Param median size for red subsample
      real*8        e_logmured_7(nr__7) ! ([kpc]) ?=-99 Param 1-sigma error for red subsample
      real*8        logmublue_7(nr__7) ! ([kpc]) ?=-99 Param median size for blue subsample
      real*8        e_logmublue_7(nr__7) ! ([kpc]) ?=-99 Param 1-sigma error for blue subsample
      real*8        logmuell_7 (nr__7) ! ([kpc]) ?=-99 Param median size for
*                                   elliptical subsample
      real*8        e_logmuell_7(nr__7) ! ([kpc]) ?=-99 Param 1-sigma error for
*                                   elliptical subsample
      real*8        logmuspi_7 (nr__7) ! ([kpc]) ?=-99 Param median size for spiral subsample
      real*8        e_logmuspi_7(nr__7) ! ([kpc]) ?=-99 Param 1-sigma error for spiral subsample

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

C  Declarations for 'fig11.dat'	! Galaxy size-mass relation measured within a
                                 fixed 2Mpc aperture at z=0.05

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

      real*4        sm_4       (nr__8) ! ([Msun]) Log10 stellar mass (bin center)
      real*8        logmured_8 (nr__8) ! ([kpc]) ?=-99 Median size of red galaxies
      real*8        e_logmured_8(nr__8) ! ([kpc]) ?=-99 1-sigma error for red galaxies
      real*8        logmublue_8(nr__8) ! ([kpc]) ?=-99 Median size of blue galaxies
      real*8        e_logmublue_8(nr__8) ! ([kpc]) ?=-99 1-sigma error for blue galaxies
      real*8        logmuell_8 (nr__8) ! ([kpc]) ?=-99 Median size of elliptical galaxies
      real*8        e_logmuell_8(nr__8) ! ([kpc]) ?=-99 1-sigma error for elliptical galaxies
      real*8        logmuspi_8 (nr__8) ! ([kpc]) ?=-99 Median size of spiral galaxies
      real*8        e_logmuspi_8(nr__8) ! ([kpc]) ?=-99 1-sigma error for spiral galaxies
      real*8        logmuredell(nr__8) ! ([kpc]) ?=-99 Median size of red ellipticals
      real*8        e_logmuredell(nr__8) ! ([kpc]) ?=-99 1-sigma error for red ellipticals
      real*8        logmuredspi(nr__8) ! ([kpc]) ?=-99 Median size of red spirals
      real*8        e_logmuredspi(nr__8) ! ([kpc]) ?=-99 1-sigma error for red spirals
      integer*4     logmublueell(nr__8) ! ([kpc]) [-99] Median size of blue ellipticals
      integer*4     e_logmublueell(nr__8) ! ([kpc]) [-99] 1-sigma error for blue ellipticals
      real*8        logmubluespi(nr__8) ! ([kpc]) ?=-99 Median size of blue spirals
      real*8        e_logmubluespi(nr__8) ! ([kpc]) ?=-99 1-sigma error for blue spirals

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

C  Loading file 'fig4a.dat'	! Spectroscopic-redshift sample,
*                                 circular R50 sizes, in a given magnitude

C  Format for file interpretation

    1 format(
     +  A2,1X,A6,1X,F5.1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig4a.dat')
      write(6,*) '....Loading file: fig4a.dat'
      do i__=1,11
        read(1,'(A125)')ar__
        read(ar__,1)
     +  Sample(i__),Param(i__),Mbin(i__),logmuall(i__),
     +  e_logmuall(i__),logmured(i__),e_logmured(i__),logmublue(i__),
     +  e_logmublue(i__),logmuell(i__),e_logmuell(i__),logmuspi(i__),
     +  e_logmuspi(i__)
c    ..............Just test output...........
        write(6,1)
     +  Sample(i__),Param(i__),Mbin(i__),logmuall(i__),
     +  e_logmuall(i__),logmured(i__),e_logmured(i__),logmublue(i__),
     +  e_logmublue(i__),logmuell(i__),e_logmuell(i__),logmuspi(i__),
     +  e_logmuspi(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig4b.dat'	! Spectroscopic-redshift sample,
*                                 major-axis sizes, in a given magnitude

C  Format for file interpretation

    2 format(
     +  A2,1X,A6,1X,F5.1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig4b.dat')
      write(6,*) '....Loading file: fig4b.dat'
      do i__=1,11
        read(1,'(A125)')ar__1
        read(ar__1,2)
     +  Sample_1(i__),Param_1(i__),Mbin_1(i__),logmuall_1(i__),
     +  e_logmuall_1(i__),logmured_1(i__),e_logmured_1(i__),
     +  logmublue_1(i__),e_logmublue_1(i__),logmuell_1(i__),
     +  e_logmuell_1(i__),logmuspi_1(i__),e_logmuspi_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Sample_1(i__),Param_1(i__),Mbin_1(i__),logmuall_1(i__),
     +  e_logmuall_1(i__),logmured_1(i__),e_logmured_1(i__),
     +  logmublue_1(i__),e_logmublue_1(i__),logmuell_1(i__),
     +  e_logmuell_1(i__),logmuspi_1(i__),e_logmuspi_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig4c.dat'	! Photometric-redshift sample,
*                                 circular R50 sizes, in a given magnitude

C  Format for file interpretation

    3 format(
     +  A2,1X,A6,1X,F5.1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig4c.dat')
      write(6,*) '....Loading file: fig4c.dat'
      do i__=1,15
        read(1,'(A125)')ar__2
        read(ar__2,3)
     +  Sample_2(i__),Param_2(i__),Mbin_2(i__),logmuall_2(i__),
     +  e_logmuall_2(i__),logmured_2(i__),e_logmured_2(i__),
     +  logmublue_2(i__),e_logmublue_2(i__),logmuell_2(i__),
     +  e_logmuell_2(i__),logmuspi_2(i__),e_logmuspi_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  Sample_2(i__),Param_2(i__),Mbin_2(i__),logmuall_2(i__),
     +  e_logmuall_2(i__),logmured_2(i__),e_logmured_2(i__),
     +  logmublue_2(i__),e_logmublue_2(i__),logmuell_2(i__),
     +  e_logmuell_2(i__),logmuspi_2(i__),e_logmuspi_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig4d.dat'	! Photometric-redshift sample,
*                                 major-axis sizes, in a given magnitude

C  Format for file interpretation

    4 format(
     +  A2,1X,A6,1X,F5.1,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig4d.dat')
      write(6,*) '....Loading file: fig4d.dat'
      do i__=1,15
        read(1,'(A125)')ar__3
        read(ar__3,4)
     +  Sample_3(i__),Param_3(i__),Mbin_3(i__),logmuall_3(i__),
     +  e_logmuall_3(i__),logmured_3(i__),e_logmured_3(i__),
     +  logmublue_3(i__),e_logmublue_3(i__),logmuell_3(i__),
     +  e_logmuell_3(i__),logmuspi_3(i__),e_logmuspi_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  Sample_3(i__),Param_3(i__),Mbin_3(i__),logmuall_3(i__),
     +  e_logmuall_3(i__),logmured_3(i__),e_logmured_3(i__),
     +  logmublue_3(i__),e_logmublue_3(i__),logmuell_3(i__),
     +  e_logmuell_3(i__),logmuspi_3(i__),e_logmuspi_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig5a.dat'	! Spectroscopic-redshift sample,
*                                 circularized half-light radii,
*                                 in a given stellar mass bin

C  Format for file interpretation

    5 format(
     +  A2,1X,A4,1X,F5.2,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig5a.dat')
      write(6,*) '....Loading file: fig5a.dat'
      do i__=1,14
        read(1,'(A123)')ar__4
        read(ar__4,5)
     +  Sample_4(i__),Param_4(i__),sm(i__),logmuall_4(i__),
     +  e_logmuall_4(i__),logmured_4(i__),e_logmured_4(i__),
     +  logmublue_4(i__),e_logmublue_4(i__),logmuell_4(i__),
     +  e_logmuell_4(i__),logmuspi_4(i__),e_logmuspi_4(i__)
c    ..............Just test output...........
        write(6,5)
     +  Sample_4(i__),Param_4(i__),sm(i__),logmuall_4(i__),
     +  e_logmuall_4(i__),logmured_4(i__),e_logmured_4(i__),
     +  logmublue_4(i__),e_logmublue_4(i__),logmuell_4(i__),
     +  e_logmuell_4(i__),logmuspi_4(i__),e_logmuspi_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig5b.dat'	! Spectroscopic-redshift sample,
*                                 major-axis half-light radii,
*                                 in a given stellar mass bin

C  Format for file interpretation

    6 format(
     +  A2,1X,A4,1X,F5.2,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig5b.dat')
      write(6,*) '....Loading file: fig5b.dat'
      do i__=1,14
        read(1,'(A123)')ar__5
        read(ar__5,6)
     +  Sample_5(i__),Param_5(i__),sm_1(i__),logmuall_5(i__),
     +  e_logmuall_5(i__),logmured_5(i__),e_logmured_5(i__),
     +  logmublue_5(i__),e_logmublue_5(i__),logmuell_5(i__),
     +  e_logmuell_5(i__),logmuspi_5(i__),e_logmuspi_5(i__)
c    ..............Just test output...........
        write(6,6)
     +  Sample_5(i__),Param_5(i__),sm_1(i__),logmuall_5(i__),
     +  e_logmuall_5(i__),logmured_5(i__),e_logmured_5(i__),
     +  logmublue_5(i__),e_logmublue_5(i__),logmuell_5(i__),
     +  e_logmuell_5(i__),logmuspi_5(i__),e_logmuspi_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig5c.dat'	! Photometric-redshift sample,
*                                 circularized half-light radii,
*                                 in a given stellar mass bin

C  Format for file interpretation

    7 format(
     +  A2,1X,A4,1X,F5.2,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig5c.dat')
      write(6,*) '....Loading file: fig5c.dat'
      do i__=1,14
        read(1,'(A123)')ar__6
        read(ar__6,7)
     +  Sample_6(i__),Param_6(i__),sm_2(i__),logmuall_6(i__),
     +  e_logmuall_6(i__),logmured_6(i__),e_logmured_6(i__),
     +  logmublue_6(i__),e_logmublue_6(i__),logmuell_6(i__),
     +  e_logmuell_6(i__),logmuspi_6(i__),e_logmuspi_6(i__)
c    ..............Just test output...........
        write(6,7)
     +  Sample_6(i__),Param_6(i__),sm_2(i__),logmuall_6(i__),
     +  e_logmuall_6(i__),logmured_6(i__),e_logmured_6(i__),
     +  logmublue_6(i__),e_logmublue_6(i__),logmuell_6(i__),
     +  e_logmuell_6(i__),logmuspi_6(i__),e_logmuspi_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig5d.dat'	! Photometric-redshift sample,
*                                 major-axis half-light radii,
*                                 in a given stellar mass bin

C  Format for file interpretation

    8 format(
     +  A2,1X,A4,1X,F5.2,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig5d.dat')
      write(6,*) '....Loading file: fig5d.dat'
      do i__=1,14
        read(1,'(A123)')ar__7
        read(ar__7,8)
     +  Sample_7(i__),Param_7(i__),sm_3(i__),logmuall_7(i__),
     +  e_logmuall_7(i__),logmured_7(i__),e_logmured_7(i__),
     +  logmublue_7(i__),e_logmublue_7(i__),logmuell_7(i__),
     +  e_logmuell_7(i__),logmuspi_7(i__),e_logmuspi_7(i__)
c    ..............Just test output...........
        write(6,8)
     +  Sample_7(i__),Param_7(i__),sm_3(i__),logmuall_7(i__),
     +  e_logmuall_7(i__),logmured_7(i__),e_logmured_7(i__),
     +  logmublue_7(i__),e_logmublue_7(i__),logmuell_7(i__),
     +  e_logmuell_7(i__),logmuspi_7(i__),e_logmuspi_7(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'fig11.dat'	! Galaxy size-mass relation measured within a
*                                 fixed 2Mpc aperture at z=0.05

C  Format for file interpretation

    9 format(
     +  F5.2,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,
     +  F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,I3,1X,
     +  I3,1X,F10.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'fig11.dat')
      write(6,*) '....Loading file: fig11.dat'
      do i__=1,12
        read(1,'(A167)')ar__8
        read(ar__8,9)
     +  sm_4(i__),logmured_8(i__),e_logmured_8(i__),logmublue_8(i__),
     +  e_logmublue_8(i__),logmuell_8(i__),e_logmuell_8(i__),
     +  logmuspi_8(i__),e_logmuspi_8(i__),logmuredell(i__),
     +  e_logmuredell(i__),logmuredspi(i__),e_logmuredspi(i__),
     +  logmublueell(i__),e_logmublueell(i__),logmubluespi(i__),
     +  e_logmubluespi(i__)
c    ..............Just test output...........
        write(6,9)
     +  sm_4(i__),logmured_8(i__),e_logmured_8(i__),logmublue_8(i__),
     +  e_logmublue_8(i__),logmuell_8(i__),e_logmuell_8(i__),
     +  logmuspi_8(i__),e_logmuspi_8(i__),logmuredell(i__),
     +  e_logmuredell(i__),logmuredspi(i__),e_logmuredspi(i__),
     +  logmublueell(i__),e_logmublueell(i__),logmubluespi(i__),
     +  e_logmubluespi(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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