FORTRAN Generation
(/./ftp/cats/J/MNRAS/404/1801)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/404/1801 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/MNRAS/404/1801    Horizon MareNostrum cosmological run     (Gay+, 2010)
*================================================================================
*On the filamentary environment of galaxies.
*    Gay C., Pichon C., Borgne D.L., Teyssier R., Sousbie T., Devriendt J.
*   <Mon. Not. R. Astron. Soc., 404, 1801-1816 (2010)>
*   =2010MNRAS.404.1801G
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'mn1d57.dat'	! Catalogue at z=1.56506

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

      real*4        Xpos       (nr__) ! (Mpc) X position (in Mpc/h)
      real*4        Ypos       (nr__) ! (Mpc) Y position (in Mpc/h)
      real*4        Zpos       (nr__) ! (Mpc) Z position (in Mpc/h)
      real*8        z          (nr__) ! Redshift used for the model (1)
      real*8        Zl         (nr__) ! Luminosity-averaged stellar metallicity
      real*8        Zm         (nr__) ! Mass-averaged stellar metallicity
      real*8        Agel       (nr__) ! (Myr) Luminosity-averaged stellar age
      real*8        Agem       (nr__) ! (Myr) Mass-averaged stellar age
      real*8        FUV        (nr__) ! (mag) Observed AB magnitude in FUV band (2)
      real*8        Gmag       (nr__) ! (mag) Observed AB magnitude in G band (2)
      real*8        Rmag       (nr__) ! (mag) Observed AB magnitude in R band (2)
      real*8        Imag       (nr__) ! (mag) Observed AB magnitude in I band (2)
      real*8        zmag       (nr__) ! (mag) Observed AB magnitude in z band (2)
      real*8        Kmag       (nr__) ! (mag) Observed AB magnitude in K band (2)
      real*8        v_3_6_     (nr__) ! (mag) Observed AB magnitude in IRAC 3.6um band (2)
      real*8        v_8_0_     (nr__) ! (mag) Observed AB magnitude in IRAC 8um band (2)
      real*8        FUVM       (nr__) ! (mag) Absolute AB magnitude in FUV band
      real*8        GMAG_1     (nr__) ! (mag) Absolute AB magnitude in G band
      real*8        RMAG_1     (nr__) ! (mag) Absolute AB magnitude in R band
      real*8        IMAG_1     (nr__) ! (mag) Absolute AB magnitude in I band
      real*8        zMAG_1     (nr__) ! (mag) Absolute AB magnitude in z band
      real*8        KMAG_1     (nr__) ! (mag) Absolute AB magnitude in K band
      real*8        v_3_6_M    (nr__) ! (mag) Absolute AB magnitude in IRAC 3.6um band
      real*8        v_8_0_M    (nr__) ! (mag) Absolute AB magnitude in IRAC 8um band
      real*8        Dists      (nr__) ! (Mpc) Distance to the skeleton (in Mpc/h)
      real*8        Distn      (nr__) ! (Mpc) Distance to the nodes of the skeleton
*                                   (in Mpc/h)
      real*8        SFR        (nr__) ! (Msun/yr) Stellar formation rate
      real*4        Mass       (nr__) ! (Msun) Total stellar mass
*Note (1): Redshift values:
*   * for mn1.57.dat, z=1.56506
*   * for mn1.80.dat, z=1.80229
*   * for mn2.10.dat, z=2.09892
*   * for mn2.51.dat, z=2.50993
*   * for mn3.01.dat, z=3.00882
*   * for mn3.53.dat, z=3.53147
*   * for mn3.95.dat, z=3.95021
*Note (2): Available filters:
*  FUV   = FUV Galex, 1520{AA}
*  Gmag  = G, 4810{AA}
*  Rmag  = R, 6980{AA}
*  Imag  = I CFHT 12K, 8130{AA}
*  zmag  = SDSS-z, 8960{AA}
*  Kmag  = Johnson K, 21950{AA}
*  [3.6] = Spitzer IRAC channel 1 IRAC-3.6um, 35610{AA}
*  [8.0] = Spitzer IRAC channel 4 IRAC-8um, 79580{AA}

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

C  Declarations for 'mn1d80.dat'	! Catalogue at z=1.80229

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

      real*4        Xpos_1     (nr__1) ! (Mpc) X position (in Mpc/h)
      real*4        Ypos_1     (nr__1) ! (Mpc) Y position (in Mpc/h)
      real*4        Zpos_1     (nr__1) ! (Mpc) Z position (in Mpc/h)
      real*8        z_1        (nr__1) ! Redshift used for the model (1)
      real*8        Zl_1       (nr__1) ! Luminosity-averaged stellar metallicity
      real*8        Zm_1       (nr__1) ! Mass-averaged stellar metallicity
      real*8        Agel_1     (nr__1) ! (Myr) Luminosity-averaged stellar age
      real*8        Agem_1     (nr__1) ! (Myr) Mass-averaged stellar age
      real*8        FUV_1      (nr__1) ! (mag) Observed AB magnitude in FUV band (2)
      real*8        Gmag_2     (nr__1) ! (mag) Observed AB magnitude in G band (2)
      real*8        Rmag_2     (nr__1) ! (mag) Observed AB magnitude in R band (2)
      real*8        Imag_2     (nr__1) ! (mag) Observed AB magnitude in I band (2)
      real*8        zmag_2     (nr__1) ! (mag) Observed AB magnitude in z band (2)
      real*8        Kmag_2     (nr__1) ! (mag) Observed AB magnitude in K band (2)
      real*8        v_3_6__1   (nr__1) ! (mag) Observed AB magnitude in IRAC 3.6um band (2)
      real*8        v_8_0__1   (nr__1) ! (mag) Observed AB magnitude in IRAC 8um band (2)
      real*8        FUVM_1     (nr__1) ! (mag) Absolute AB magnitude in FUV band
      real*8        GMAG_3     (nr__1) ! (mag) Absolute AB magnitude in G band
      real*8        RMAG_3     (nr__1) ! (mag) Absolute AB magnitude in R band
      real*8        IMAG_3     (nr__1) ! (mag) Absolute AB magnitude in I band
      real*8        zMAG_3     (nr__1) ! (mag) Absolute AB magnitude in z band
      real*8        KMAG_3     (nr__1) ! (mag) Absolute AB magnitude in K band
      real*8        v_3_6_M_1  (nr__1) ! (mag) Absolute AB magnitude in IRAC 3.6um band
      real*8        v_8_0_M_1  (nr__1) ! (mag) Absolute AB magnitude in IRAC 8um band
      real*8        Dists_1    (nr__1) ! (Mpc) Distance to the skeleton (in Mpc/h)
      real*8        Distn_1    (nr__1) ! (Mpc) Distance to the nodes of the skeleton
*                                   (in Mpc/h)
      real*8        SFR_1      (nr__1) ! (Msun/yr) Stellar formation rate
      real*4        Mass_1     (nr__1) ! (Msun) Total stellar mass
*Note (1): Redshift values:
*   * for mn1.57.dat, z=1.56506
*   * for mn1.80.dat, z=1.80229
*   * for mn2.10.dat, z=2.09892
*   * for mn2.51.dat, z=2.50993
*   * for mn3.01.dat, z=3.00882
*   * for mn3.53.dat, z=3.53147
*   * for mn3.95.dat, z=3.95021
*Note (2): Available filters:
*  FUV   = FUV Galex, 1520{AA}
*  Gmag  = G, 4810{AA}
*  Rmag  = R, 6980{AA}
*  Imag  = I CFHT 12K, 8130{AA}
*  zmag  = SDSS-z, 8960{AA}
*  Kmag  = Johnson K, 21950{AA}
*  [3.6] = Spitzer IRAC channel 1 IRAC-3.6um, 35610{AA}
*  [8.0] = Spitzer IRAC channel 4 IRAC-8um, 79580{AA}

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

C  Declarations for 'mn2d10.dat'	! Catalogue at z=2.09892

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

      real*4        Xpos_2     (nr__2) ! (Mpc) X position (in Mpc/h)
      real*4        Ypos_2     (nr__2) ! (Mpc) Y position (in Mpc/h)
      real*4        Zpos_2     (nr__2) ! (Mpc) Z position (in Mpc/h)
      real*8        z_2        (nr__2) ! Redshift used for the model (1)
      real*8        Zl_2       (nr__2) ! Luminosity-averaged stellar metallicity
      real*8        Zm_2       (nr__2) ! Mass-averaged stellar metallicity
      real*8        Agel_2     (nr__2) ! (Myr) Luminosity-averaged stellar age
      real*8        Agem_2     (nr__2) ! (Myr) Mass-averaged stellar age
      real*8        FUV_2      (nr__2) ! (mag) Observed AB magnitude in FUV band (2)
      real*8        Gmag_4     (nr__2) ! (mag) Observed AB magnitude in G band (2)
      real*8        Rmag_4     (nr__2) ! (mag) Observed AB magnitude in R band (2)
      real*8        Imag_4     (nr__2) ! (mag) Observed AB magnitude in I band (2)
      real*8        zmag_4     (nr__2) ! (mag) Observed AB magnitude in z band (2)
      real*8        Kmag_4     (nr__2) ! (mag) Observed AB magnitude in K band (2)
      real*8        v_3_6__2   (nr__2) ! (mag) Observed AB magnitude in IRAC 3.6um band (2)
      real*8        v_8_0__2   (nr__2) ! (mag) Observed AB magnitude in IRAC 8um band (2)
      real*8        FUVM_2     (nr__2) ! (mag) Absolute AB magnitude in FUV band
      real*8        GMAG_5     (nr__2) ! (mag) Absolute AB magnitude in G band
      real*8        RMAG_5     (nr__2) ! (mag) Absolute AB magnitude in R band
      real*8        IMAG_5     (nr__2) ! (mag) Absolute AB magnitude in I band
      real*8        zMAG_5     (nr__2) ! (mag) Absolute AB magnitude in z band
      real*8        KMAG_5     (nr__2) ! (mag) Absolute AB magnitude in K band
      real*8        v_3_6_M_2  (nr__2) ! (mag) Absolute AB magnitude in IRAC 3.6um band
      real*8        v_8_0_M_2  (nr__2) ! (mag) Absolute AB magnitude in IRAC 8um band
      real*8        Dists_2    (nr__2) ! (Mpc) Distance to the skeleton (in Mpc/h)
      real*8        Distn_2    (nr__2) ! (Mpc) Distance to the nodes of the skeleton
*                                   (in Mpc/h)
      real*8        SFR_2      (nr__2) ! (Msun/yr) Stellar formation rate
      real*4        Mass_2     (nr__2) ! (Msun) Total stellar mass
*Note (1): Redshift values:
*   * for mn1.57.dat, z=1.56506
*   * for mn1.80.dat, z=1.80229
*   * for mn2.10.dat, z=2.09892
*   * for mn2.51.dat, z=2.50993
*   * for mn3.01.dat, z=3.00882
*   * for mn3.53.dat, z=3.53147
*   * for mn3.95.dat, z=3.95021
*Note (2): Available filters:
*  FUV   = FUV Galex, 1520{AA}
*  Gmag  = G, 4810{AA}
*  Rmag  = R, 6980{AA}
*  Imag  = I CFHT 12K, 8130{AA}
*  zmag  = SDSS-z, 8960{AA}
*  Kmag  = Johnson K, 21950{AA}
*  [3.6] = Spitzer IRAC channel 1 IRAC-3.6um, 35610{AA}
*  [8.0] = Spitzer IRAC channel 4 IRAC-8um, 79580{AA}

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

C  Declarations for 'mn2d51.dat'	! Catalogue at z=2.50993

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

      real*4        Xpos_3     (nr__3) ! (Mpc) X position (in Mpc/h)
      real*4        Ypos_3     (nr__3) ! (Mpc) Y position (in Mpc/h)
      real*4        Zpos_3     (nr__3) ! (Mpc) Z position (in Mpc/h)
      real*8        z_3        (nr__3) ! Redshift used for the model (1)
      real*8        Zl_3       (nr__3) ! Luminosity-averaged stellar metallicity
      real*8        Zm_3       (nr__3) ! Mass-averaged stellar metallicity
      real*8        Agel_3     (nr__3) ! (Myr) Luminosity-averaged stellar age
      real*8        Agem_3     (nr__3) ! (Myr) Mass-averaged stellar age
      real*8        FUV_3      (nr__3) ! (mag) Observed AB magnitude in FUV band (2)
      real*8        Gmag_6     (nr__3) ! (mag) Observed AB magnitude in G band (2)
      real*8        Rmag_6     (nr__3) ! (mag) Observed AB magnitude in R band (2)
      real*8        Imag_6     (nr__3) ! (mag) Observed AB magnitude in I band (2)
      real*8        zmag_6     (nr__3) ! (mag) Observed AB magnitude in z band (2)
      real*8        Kmag_6     (nr__3) ! (mag) Observed AB magnitude in K band (2)
      real*8        v_3_6__3   (nr__3) ! (mag) Observed AB magnitude in IRAC 3.6um band (2)
      real*8        v_8_0__3   (nr__3) ! (mag) Observed AB magnitude in IRAC 8um band (2)
      real*8        FUVM_3     (nr__3) ! (mag) Absolute AB magnitude in FUV band
      real*8        GMAG_7     (nr__3) ! (mag) Absolute AB magnitude in G band
      real*8        RMAG_7     (nr__3) ! (mag) Absolute AB magnitude in R band
      real*8        IMAG_7     (nr__3) ! (mag) Absolute AB magnitude in I band
      real*8        zMAG_7     (nr__3) ! (mag) Absolute AB magnitude in z band
      real*8        KMAG_7     (nr__3) ! (mag) Absolute AB magnitude in K band
      real*8        v_3_6_M_3  (nr__3) ! (mag) Absolute AB magnitude in IRAC 3.6um band
      real*8        v_8_0_M_3  (nr__3) ! (mag) Absolute AB magnitude in IRAC 8um band
      real*8        Dists_3    (nr__3) ! (Mpc) Distance to the skeleton (in Mpc/h)
      real*8        Distn_3    (nr__3) ! (Mpc) Distance to the nodes of the skeleton
*                                   (in Mpc/h)
      real*8        SFR_3      (nr__3) ! (Msun/yr) Stellar formation rate
      real*4        Mass_3     (nr__3) ! (Msun) Total stellar mass
*Note (1): Redshift values:
*   * for mn1.57.dat, z=1.56506
*   * for mn1.80.dat, z=1.80229
*   * for mn2.10.dat, z=2.09892
*   * for mn2.51.dat, z=2.50993
*   * for mn3.01.dat, z=3.00882
*   * for mn3.53.dat, z=3.53147
*   * for mn3.95.dat, z=3.95021
*Note (2): Available filters:
*  FUV   = FUV Galex, 1520{AA}
*  Gmag  = G, 4810{AA}
*  Rmag  = R, 6980{AA}
*  Imag  = I CFHT 12K, 8130{AA}
*  zmag  = SDSS-z, 8960{AA}
*  Kmag  = Johnson K, 21950{AA}
*  [3.6] = Spitzer IRAC channel 1 IRAC-3.6um, 35610{AA}
*  [8.0] = Spitzer IRAC channel 4 IRAC-8um, 79580{AA}

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

C  Declarations for 'mn3d01.dat'	! Catalogue at z=3.00882

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

      real*4        Xpos_4     (nr__4) ! (Mpc) X position (in Mpc/h)
      real*4        Ypos_4     (nr__4) ! (Mpc) Y position (in Mpc/h)
      real*4        Zpos_4     (nr__4) ! (Mpc) Z position (in Mpc/h)
      real*8        z_4        (nr__4) ! Redshift used for the model (1)
      real*8        Zl_4       (nr__4) ! Luminosity-averaged stellar metallicity
      real*8        Zm_4       (nr__4) ! Mass-averaged stellar metallicity
      real*8        Agel_4     (nr__4) ! (Myr) Luminosity-averaged stellar age
      real*8        Agem_4     (nr__4) ! (Myr) Mass-averaged stellar age
      real*8        FUV_4      (nr__4) ! (mag) Observed AB magnitude in FUV band (2)
      real*8        Gmag_8     (nr__4) ! (mag) Observed AB magnitude in G band (2)
      real*8        Rmag_8     (nr__4) ! (mag) Observed AB magnitude in R band (2)
      real*8        Imag_8     (nr__4) ! (mag) Observed AB magnitude in I band (2)
      real*8        zmag_8     (nr__4) ! (mag) Observed AB magnitude in z band (2)
      real*8        Kmag_8     (nr__4) ! (mag) Observed AB magnitude in K band (2)
      real*8        v_3_6__4   (nr__4) ! (mag) Observed AB magnitude in IRAC 3.6um band (2)
      real*8        v_8_0__4   (nr__4) ! (mag) Observed AB magnitude in IRAC 8um band (2)
      real*8        FUVM_4     (nr__4) ! (mag) Absolute AB magnitude in FUV band
      real*8        GMAG_9     (nr__4) ! (mag) Absolute AB magnitude in G band
      real*8        RMAG_9     (nr__4) ! (mag) Absolute AB magnitude in R band
      real*8        IMAG_9     (nr__4) ! (mag) Absolute AB magnitude in I band
      real*8        zMAG_9     (nr__4) ! (mag) Absolute AB magnitude in z band
      real*8        KMAG_9     (nr__4) ! (mag) Absolute AB magnitude in K band
      real*8        v_3_6_M_4  (nr__4) ! (mag) Absolute AB magnitude in IRAC 3.6um band
      real*8        v_8_0_M_4  (nr__4) ! (mag) Absolute AB magnitude in IRAC 8um band
      real*8        Dists_4    (nr__4) ! (Mpc) Distance to the skeleton (in Mpc/h)
      real*8        Distn_4    (nr__4) ! (Mpc) Distance to the nodes of the skeleton
*                                   (in Mpc/h)
      real*8        SFR_4      (nr__4) ! (Msun/yr) Stellar formation rate
      real*4        Mass_4     (nr__4) ! (Msun) Total stellar mass
*Note (1): Redshift values:
*   * for mn1.57.dat, z=1.56506
*   * for mn1.80.dat, z=1.80229
*   * for mn2.10.dat, z=2.09892
*   * for mn2.51.dat, z=2.50993
*   * for mn3.01.dat, z=3.00882
*   * for mn3.53.dat, z=3.53147
*   * for mn3.95.dat, z=3.95021
*Note (2): Available filters:
*  FUV   = FUV Galex, 1520{AA}
*  Gmag  = G, 4810{AA}
*  Rmag  = R, 6980{AA}
*  Imag  = I CFHT 12K, 8130{AA}
*  zmag  = SDSS-z, 8960{AA}
*  Kmag  = Johnson K, 21950{AA}
*  [3.6] = Spitzer IRAC channel 1 IRAC-3.6um, 35610{AA}
*  [8.0] = Spitzer IRAC channel 4 IRAC-8um, 79580{AA}

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

C  Declarations for 'mn3d53.dat'	! Catalogue at z=3.53147

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

      real*4        Xpos_5     (nr__5) ! (Mpc) X position (in Mpc/h)
      real*4        Ypos_5     (nr__5) ! (Mpc) Y position (in Mpc/h)
      real*4        Zpos_5     (nr__5) ! (Mpc) Z position (in Mpc/h)
      real*8        z_5        (nr__5) ! Redshift used for the model (1)
      real*8        Zl_5       (nr__5) ! Luminosity-averaged stellar metallicity
      real*8        Zm_5       (nr__5) ! Mass-averaged stellar metallicity
      real*8        Agel_5     (nr__5) ! (Myr) Luminosity-averaged stellar age
      real*8        Agem_5     (nr__5) ! (Myr) Mass-averaged stellar age
      real*8        FUV_5      (nr__5) ! (mag) Observed AB magnitude in FUV band (2)
      real*8        Gmag_10    (nr__5) ! (mag) Observed AB magnitude in G band (2)
      real*8        Rmag_10    (nr__5) ! (mag) Observed AB magnitude in R band (2)
      real*8        Imag_10    (nr__5) ! (mag) Observed AB magnitude in I band (2)
      real*8        zmag_10    (nr__5) ! (mag) Observed AB magnitude in z band (2)
      real*8        Kmag_10    (nr__5) ! (mag) Observed AB magnitude in K band (2)
      real*8        v_3_6__5   (nr__5) ! (mag) Observed AB magnitude in IRAC 3.6um band (2)
      real*8        v_8_0__5   (nr__5) ! (mag) Observed AB magnitude in IRAC 8um band (2)
      real*8        FUVM_5     (nr__5) ! (mag) Absolute AB magnitude in FUV band
      real*8        GMAG_11    (nr__5) ! (mag) Absolute AB magnitude in G band
      real*8        RMAG_11    (nr__5) ! (mag) Absolute AB magnitude in R band
      real*8        IMAG_11    (nr__5) ! (mag) Absolute AB magnitude in I band
      real*8        zMAG_11    (nr__5) ! (mag) Absolute AB magnitude in z band
      real*8        KMAG_11    (nr__5) ! (mag) Absolute AB magnitude in K band
      real*8        v_3_6_M_5  (nr__5) ! (mag) Absolute AB magnitude in IRAC 3.6um band
      real*8        v_8_0_M_5  (nr__5) ! (mag) Absolute AB magnitude in IRAC 8um band
      real*8        Dists_5    (nr__5) ! (Mpc) Distance to the skeleton (in Mpc/h)
      real*8        Distn_5    (nr__5) ! (Mpc) Distance to the nodes of the skeleton
*                                   (in Mpc/h)
      real*8        SFR_5      (nr__5) ! (Msun/yr) Stellar formation rate
      real*4        Mass_5     (nr__5) ! (Msun) Total stellar mass
*Note (1): Redshift values:
*   * for mn1.57.dat, z=1.56506
*   * for mn1.80.dat, z=1.80229
*   * for mn2.10.dat, z=2.09892
*   * for mn2.51.dat, z=2.50993
*   * for mn3.01.dat, z=3.00882
*   * for mn3.53.dat, z=3.53147
*   * for mn3.95.dat, z=3.95021
*Note (2): Available filters:
*  FUV   = FUV Galex, 1520{AA}
*  Gmag  = G, 4810{AA}
*  Rmag  = R, 6980{AA}
*  Imag  = I CFHT 12K, 8130{AA}
*  zmag  = SDSS-z, 8960{AA}
*  Kmag  = Johnson K, 21950{AA}
*  [3.6] = Spitzer IRAC channel 1 IRAC-3.6um, 35610{AA}
*  [8.0] = Spitzer IRAC channel 4 IRAC-8um, 79580{AA}

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

C  Declarations for 'mn3d95.dat'	! Catalogue at z=3.95021

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

      real*4        Xpos_6     (nr__6) ! (Mpc) X position (in Mpc/h)
      real*4        Ypos_6     (nr__6) ! (Mpc) Y position (in Mpc/h)
      real*4        Zpos_6     (nr__6) ! (Mpc) Z position (in Mpc/h)
      real*8        z_6        (nr__6) ! Redshift used for the model (1)
      real*8        Zl_6       (nr__6) ! Luminosity-averaged stellar metallicity
      real*8        Zm_6       (nr__6) ! Mass-averaged stellar metallicity
      real*8        Agel_6     (nr__6) ! (Myr) Luminosity-averaged stellar age
      real*8        Agem_6     (nr__6) ! (Myr) Mass-averaged stellar age
      real*8        FUV_6      (nr__6) ! (mag) Observed AB magnitude in FUV band (2)
      real*8        Gmag_12    (nr__6) ! (mag) Observed AB magnitude in G band (2)
      real*8        Rmag_12    (nr__6) ! (mag) Observed AB magnitude in R band (2)
      real*8        Imag_12    (nr__6) ! (mag) Observed AB magnitude in I band (2)
      real*8        zmag_12    (nr__6) ! (mag) Observed AB magnitude in z band (2)
      real*8        Kmag_12    (nr__6) ! (mag) Observed AB magnitude in K band (2)
      real*8        v_3_6__6   (nr__6) ! (mag) Observed AB magnitude in IRAC 3.6um band (2)
      real*8        v_8_0__6   (nr__6) ! (mag) Observed AB magnitude in IRAC 8um band (2)
      real*8        FUVM_6     (nr__6) ! (mag) Absolute AB magnitude in FUV band
      real*8        GMAG_13    (nr__6) ! (mag) Absolute AB magnitude in G band
      real*8        RMAG_13    (nr__6) ! (mag) Absolute AB magnitude in R band
      real*8        IMAG_13    (nr__6) ! (mag) Absolute AB magnitude in I band
      real*8        zMAG_13    (nr__6) ! (mag) Absolute AB magnitude in z band
      real*8        KMAG_13    (nr__6) ! (mag) Absolute AB magnitude in K band
      real*8        v_3_6_M_6  (nr__6) ! (mag) Absolute AB magnitude in IRAC 3.6um band
      real*8        v_8_0_M_6  (nr__6) ! (mag) Absolute AB magnitude in IRAC 8um band
      real*8        Dists_6    (nr__6) ! (Mpc) Distance to the skeleton (in Mpc/h)
      real*8        Distn_6    (nr__6) ! (Mpc) Distance to the nodes of the skeleton
*                                   (in Mpc/h)
      real*8        SFR_6      (nr__6) ! (Msun/yr) Stellar formation rate
      real*4        Mass_6     (nr__6) ! (Msun) Total stellar mass
*Note (1): Redshift values:
*   * for mn1.57.dat, z=1.56506
*   * for mn1.80.dat, z=1.80229
*   * for mn2.10.dat, z=2.09892
*   * for mn2.51.dat, z=2.50993
*   * for mn3.01.dat, z=3.00882
*   * for mn3.53.dat, z=3.53147
*   * for mn3.95.dat, z=3.95021
*Note (2): Available filters:
*  FUV   = FUV Galex, 1520{AA}
*  Gmag  = G, 4810{AA}
*  Rmag  = R, 6980{AA}
*  Imag  = I CFHT 12K, 8130{AA}
*  zmag  = SDSS-z, 8960{AA}
*  Kmag  = Johnson K, 21950{AA}
*  [3.6] = Spitzer IRAC channel 1 IRAC-3.6um, 35610{AA}
*  [8.0] = Spitzer IRAC channel 4 IRAC-8um, 79580{AA}

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

C  Loading file 'mn1d57.dat'	! Catalogue at z=1.56506

C  Format for file interpretation

    1 format(
     +  E11.8,1X,E12.8,1X,E12.8,1X,F7.5,1X,F8.5,1X,F8.5,1X,F11.6,1X,
     +  F11.6,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F9.5,1X,F9.5,1X,F12.8,1X,F15.7,1X,F15.7,1X,E12.9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'mn1d57.dat')
      write(6,*) '....Loading file: mn1d57.dat'
      do i__=1,97563
        read(1,'(A285)')ar__
        read(ar__,1)
     +  Xpos(i__),Ypos(i__),Zpos(i__),z(i__),Zl(i__),Zm(i__),
     +  Agel(i__),Agem(i__),FUV(i__),Gmag(i__),Rmag(i__),Imag(i__),
     +  zmag(i__),Kmag(i__),v_3_6_(i__),v_8_0_(i__),FUVM(i__),
     +  GMAG_1(i__),RMAG_1(i__),IMAG_1(i__),zMAG_1(i__),KMAG_1(i__),
     +  v_3_6_M(i__),v_8_0_M(i__),Dists(i__),Distn(i__),SFR(i__),
     +  Mass(i__)
c    ..............Just test output...........
        write(6,1)
     +  Xpos(i__),Ypos(i__),Zpos(i__),z(i__),Zl(i__),Zm(i__),
     +  Agel(i__),Agem(i__),FUV(i__),Gmag(i__),Rmag(i__),Imag(i__),
     +  zmag(i__),Kmag(i__),v_3_6_(i__),v_8_0_(i__),FUVM(i__),
     +  GMAG_1(i__),RMAG_1(i__),IMAG_1(i__),zMAG_1(i__),KMAG_1(i__),
     +  v_3_6_M(i__),v_8_0_M(i__),Dists(i__),Distn(i__),SFR(i__),
     +  Mass(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'mn1d80.dat'	! Catalogue at z=1.80229

C  Format for file interpretation

    2 format(
     +  E11.8,1X,E12.8,1X,E12.8,1X,F7.5,1X,F8.5,1X,F8.5,1X,F11.6,1X,
     +  F11.6,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F9.5,1X,F9.5,1X,F12.8,1X,F15.7,1X,F15.7,1X,E12.9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'mn1d80.dat')
      write(6,*) '....Loading file: mn1d80.dat'
      do i__=1,103589
        read(1,'(A285)')ar__1
        read(ar__1,2)
     +  Xpos_1(i__),Ypos_1(i__),Zpos_1(i__),z_1(i__),Zl_1(i__),
     +  Zm_1(i__),Agel_1(i__),Agem_1(i__),FUV_1(i__),Gmag_2(i__),
     +  Rmag_2(i__),Imag_2(i__),zmag_2(i__),Kmag_2(i__),v_3_6__1(i__),
     +  v_8_0__1(i__),FUVM_1(i__),GMAG_3(i__),RMAG_3(i__),IMAG_3(i__),
     +  zMAG_3(i__),KMAG_3(i__),v_3_6_M_1(i__),v_8_0_M_1(i__),
     +  Dists_1(i__),Distn_1(i__),SFR_1(i__),Mass_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Xpos_1(i__),Ypos_1(i__),Zpos_1(i__),z_1(i__),Zl_1(i__),
     +  Zm_1(i__),Agel_1(i__),Agem_1(i__),FUV_1(i__),Gmag_2(i__),
     +  Rmag_2(i__),Imag_2(i__),zmag_2(i__),Kmag_2(i__),v_3_6__1(i__),
     +  v_8_0__1(i__),FUVM_1(i__),GMAG_3(i__),RMAG_3(i__),IMAG_3(i__),
     +  zMAG_3(i__),KMAG_3(i__),v_3_6_M_1(i__),v_8_0_M_1(i__),
     +  Dists_1(i__),Distn_1(i__),SFR_1(i__),Mass_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'mn2d10.dat'	! Catalogue at z=2.09892

C  Format for file interpretation

    3 format(
     +  E11.8,1X,E12.8,1X,E12.8,1X,F7.5,1X,F8.5,1X,F8.5,1X,F11.6,1X,
     +  F11.6,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F9.5,1X,F9.5,1X,F12.8,1X,F15.7,1X,F15.7,1X,E12.9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'mn2d10.dat')
      write(6,*) '....Loading file: mn2d10.dat'
      do i__=1,111184
        read(1,'(A285)')ar__2
        read(ar__2,3)
     +  Xpos_2(i__),Ypos_2(i__),Zpos_2(i__),z_2(i__),Zl_2(i__),
     +  Zm_2(i__),Agel_2(i__),Agem_2(i__),FUV_2(i__),Gmag_4(i__),
     +  Rmag_4(i__),Imag_4(i__),zmag_4(i__),Kmag_4(i__),v_3_6__2(i__),
     +  v_8_0__2(i__),FUVM_2(i__),GMAG_5(i__),RMAG_5(i__),IMAG_5(i__),
     +  zMAG_5(i__),KMAG_5(i__),v_3_6_M_2(i__),v_8_0_M_2(i__),
     +  Dists_2(i__),Distn_2(i__),SFR_2(i__),Mass_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  Xpos_2(i__),Ypos_2(i__),Zpos_2(i__),z_2(i__),Zl_2(i__),
     +  Zm_2(i__),Agel_2(i__),Agem_2(i__),FUV_2(i__),Gmag_4(i__),
     +  Rmag_4(i__),Imag_4(i__),zmag_4(i__),Kmag_4(i__),v_3_6__2(i__),
     +  v_8_0__2(i__),FUVM_2(i__),GMAG_5(i__),RMAG_5(i__),IMAG_5(i__),
     +  zMAG_5(i__),KMAG_5(i__),v_3_6_M_2(i__),v_8_0_M_2(i__),
     +  Dists_2(i__),Distn_2(i__),SFR_2(i__),Mass_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'mn2d51.dat'	! Catalogue at z=2.50993

C  Format for file interpretation

    4 format(
     +  E11.8,1X,E12.8,1X,E12.8,1X,F7.5,1X,F8.5,1X,F8.5,1X,F11.6,1X,
     +  F11.6,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F9.5,1X,F9.5,1X,F12.8,1X,F15.7,1X,F15.7,1X,E12.9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'mn2d51.dat')
      write(6,*) '....Loading file: mn2d51.dat'
      do i__=1,119978
        read(1,'(A285)')ar__3
        read(ar__3,4)
     +  Xpos_3(i__),Ypos_3(i__),Zpos_3(i__),z_3(i__),Zl_3(i__),
     +  Zm_3(i__),Agel_3(i__),Agem_3(i__),FUV_3(i__),Gmag_6(i__),
     +  Rmag_6(i__),Imag_6(i__),zmag_6(i__),Kmag_6(i__),v_3_6__3(i__),
     +  v_8_0__3(i__),FUVM_3(i__),GMAG_7(i__),RMAG_7(i__),IMAG_7(i__),
     +  zMAG_7(i__),KMAG_7(i__),v_3_6_M_3(i__),v_8_0_M_3(i__),
     +  Dists_3(i__),Distn_3(i__),SFR_3(i__),Mass_3(i__)
c    ..............Just test output...........
        write(6,4)
     +  Xpos_3(i__),Ypos_3(i__),Zpos_3(i__),z_3(i__),Zl_3(i__),
     +  Zm_3(i__),Agel_3(i__),Agem_3(i__),FUV_3(i__),Gmag_6(i__),
     +  Rmag_6(i__),Imag_6(i__),zmag_6(i__),Kmag_6(i__),v_3_6__3(i__),
     +  v_8_0__3(i__),FUVM_3(i__),GMAG_7(i__),RMAG_7(i__),IMAG_7(i__),
     +  zMAG_7(i__),KMAG_7(i__),v_3_6_M_3(i__),v_8_0_M_3(i__),
     +  Dists_3(i__),Distn_3(i__),SFR_3(i__),Mass_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'mn3d01.dat'	! Catalogue at z=3.00882

C  Format for file interpretation

    5 format(
     +  E11.8,1X,E12.8,1X,E12.8,1X,F7.5,1X,F8.5,1X,F8.5,1X,F11.6,1X,
     +  F11.6,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F9.5,1X,F9.5,1X,F12.8,1X,F15.7,1X,F15.7,1X,E12.9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'mn3d01.dat')
      write(6,*) '....Loading file: mn3d01.dat'
      do i__=1,124642
        read(1,'(A285)')ar__4
        read(ar__4,5)
     +  Xpos_4(i__),Ypos_4(i__),Zpos_4(i__),z_4(i__),Zl_4(i__),
     +  Zm_4(i__),Agel_4(i__),Agem_4(i__),FUV_4(i__),Gmag_8(i__),
     +  Rmag_8(i__),Imag_8(i__),zmag_8(i__),Kmag_8(i__),v_3_6__4(i__),
     +  v_8_0__4(i__),FUVM_4(i__),GMAG_9(i__),RMAG_9(i__),IMAG_9(i__),
     +  zMAG_9(i__),KMAG_9(i__),v_3_6_M_4(i__),v_8_0_M_4(i__),
     +  Dists_4(i__),Distn_4(i__),SFR_4(i__),Mass_4(i__)
c    ..............Just test output...........
        write(6,5)
     +  Xpos_4(i__),Ypos_4(i__),Zpos_4(i__),z_4(i__),Zl_4(i__),
     +  Zm_4(i__),Agel_4(i__),Agem_4(i__),FUV_4(i__),Gmag_8(i__),
     +  Rmag_8(i__),Imag_8(i__),zmag_8(i__),Kmag_8(i__),v_3_6__4(i__),
     +  v_8_0__4(i__),FUVM_4(i__),GMAG_9(i__),RMAG_9(i__),IMAG_9(i__),
     +  zMAG_9(i__),KMAG_9(i__),v_3_6_M_4(i__),v_8_0_M_4(i__),
     +  Dists_4(i__),Distn_4(i__),SFR_4(i__),Mass_4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'mn3d53.dat'	! Catalogue at z=3.53147

C  Format for file interpretation

    6 format(
     +  E11.8,1X,E12.8,1X,E12.8,1X,F7.5,1X,F8.5,1X,F8.5,1X,F11.6,1X,
     +  F11.6,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F9.5,1X,F9.5,1X,F12.8,1X,F15.7,1X,F15.7,1X,E12.9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'mn3d53.dat')
      write(6,*) '....Loading file: mn3d53.dat'
      do i__=1,119612
        read(1,'(A285)')ar__5
        read(ar__5,6)
     +  Xpos_5(i__),Ypos_5(i__),Zpos_5(i__),z_5(i__),Zl_5(i__),
     +  Zm_5(i__),Agel_5(i__),Agem_5(i__),FUV_5(i__),Gmag_10(i__),
     +  Rmag_10(i__),Imag_10(i__),zmag_10(i__),Kmag_10(i__),
     +  v_3_6__5(i__),v_8_0__5(i__),FUVM_5(i__),GMAG_11(i__),
     +  RMAG_11(i__),IMAG_11(i__),zMAG_11(i__),KMAG_11(i__),
     +  v_3_6_M_5(i__),v_8_0_M_5(i__),Dists_5(i__),Distn_5(i__),
     +  SFR_5(i__),Mass_5(i__)
c    ..............Just test output...........
        write(6,6)
     +  Xpos_5(i__),Ypos_5(i__),Zpos_5(i__),z_5(i__),Zl_5(i__),
     +  Zm_5(i__),Agel_5(i__),Agem_5(i__),FUV_5(i__),Gmag_10(i__),
     +  Rmag_10(i__),Imag_10(i__),zmag_10(i__),Kmag_10(i__),
     +  v_3_6__5(i__),v_8_0__5(i__),FUVM_5(i__),GMAG_11(i__),
     +  RMAG_11(i__),IMAG_11(i__),zMAG_11(i__),KMAG_11(i__),
     +  v_3_6_M_5(i__),v_8_0_M_5(i__),Dists_5(i__),Distn_5(i__),
     +  SFR_5(i__),Mass_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'mn3d95.dat'	! Catalogue at z=3.95021

C  Format for file interpretation

    7 format(
     +  E11.8,1X,E12.8,1X,E12.8,1X,F7.5,1X,F8.5,1X,F8.5,1X,F11.6,1X,
     +  F11.6,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F9.5,1X,F9.5,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F9.5,1X,F9.5,1X,F12.8,1X,F15.7,1X,F15.7,1X,E12.9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'mn3d95.dat')
      write(6,*) '....Loading file: mn3d95.dat'
      do i__=1,103187
        read(1,'(A285)')ar__6
        read(ar__6,7)
     +  Xpos_6(i__),Ypos_6(i__),Zpos_6(i__),z_6(i__),Zl_6(i__),
     +  Zm_6(i__),Agel_6(i__),Agem_6(i__),FUV_6(i__),Gmag_12(i__),
     +  Rmag_12(i__),Imag_12(i__),zmag_12(i__),Kmag_12(i__),
     +  v_3_6__6(i__),v_8_0__6(i__),FUVM_6(i__),GMAG_13(i__),
     +  RMAG_13(i__),IMAG_13(i__),zMAG_13(i__),KMAG_13(i__),
     +  v_3_6_M_6(i__),v_8_0_M_6(i__),Dists_6(i__),Distn_6(i__),
     +  SFR_6(i__),Mass_6(i__)
c    ..............Just test output...........
        write(6,7)
     +  Xpos_6(i__),Ypos_6(i__),Zpos_6(i__),z_6(i__),Zl_6(i__),
     +  Zm_6(i__),Agel_6(i__),Agem_6(i__),FUV_6(i__),Gmag_12(i__),
     +  Rmag_12(i__),Imag_12(i__),zmag_12(i__),Kmag_12(i__),
     +  v_3_6__6(i__),v_8_0__6(i__),FUVM_6(i__),GMAG_13(i__),
     +  RMAG_13(i__),IMAG_13(i__),zMAG_13(i__),KMAG_13(i__),
     +  v_3_6_M_6(i__),v_8_0_M_6(i__),Dists_6(i__),Distn_6(i__),
     +  SFR_6(i__),Mass_6(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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