Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/185/1 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-17 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/ApJS/185/1 A catalog of star formation and metallicity (Tojeiro+, 2009) *================================================================================ *A public catalog of stellar masses, star formation and metallicity histories, *and dust content from the Sloan Digital Sky Survey using VESPA. * Tojeiro R., Wilkins S., Heavens A.F., Panter B., Jimenez R. * <Astrophys. J. Suppl. Ser., 185, 1-19 (2009)> * =2009ApJS..185....1T (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal1.dat' ! RunID1 results relative to the galaxy as a whole integer*4 nr__ parameter (nr__=683113) ! Number of records character*121 ar__ ! Full-size record character*12 index (nr__) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID (nr__) ! [1,8] Gives detail of the run (G3) real*4 Mass (nr__) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass (nr__) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb (nr__) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2 (nr__) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N (nr__) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb (nr__) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np (nr__) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z (nr__) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal2.dat' ! RunID2 results relative to the galaxy as a whole integer*4 nr__1 parameter (nr__1=681311) ! Number of records character*121 ar__1 ! Full-size record character*12 index_1 (nr__1) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_1 (nr__1) ! [1,8] Gives detail of the run (G3) real*4 Mass_1 (nr__1) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass_1 (nr__1) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb_1 (nr__1) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2_1 (nr__1) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N_1 (nr__1) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb_1 (nr__1) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np_1 (nr__1) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z_1 (nr__1) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal3.dat' ! RunID3 results relative to the galaxy as a whole integer*4 nr__2 parameter (nr__2=687387) ! Number of records character*121 ar__2 ! Full-size record character*12 index_2 (nr__2) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_2 (nr__2) ! [1,8] Gives detail of the run (G3) real*4 Mass_2 (nr__2) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass_2 (nr__2) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb_2 (nr__2) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2_2 (nr__2) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N_2 (nr__2) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb_2 (nr__2) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np_2 (nr__2) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z_2 (nr__2) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal4.dat' ! RunID4 results relative to the galaxy as a whole integer*4 nr__3 parameter (nr__3=685586) ! Number of records character*121 ar__3 ! Full-size record character*12 index_3 (nr__3) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_3 (nr__3) ! [1,8] Gives detail of the run (G3) real*4 Mass_3 (nr__3) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass_3 (nr__3) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb_3 (nr__3) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2_3 (nr__3) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N_3 (nr__3) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb_3 (nr__3) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np_3 (nr__3) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z_3 (nr__3) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal5.dat' ! RunID5 results relative to the galaxy as a whole integer*4 nr__4 parameter (nr__4=89055) ! Number of records character*121 ar__4 ! Full-size record character*12 index_4 (nr__4) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_4 (nr__4) ! [1,8] Gives detail of the run (G3) real*4 Mass_4 (nr__4) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass_4 (nr__4) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb_4 (nr__4) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2_4 (nr__4) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N_4 (nr__4) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb_4 (nr__4) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np_4 (nr__4) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z_4 (nr__4) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal6.dat' ! RunID6 results relative to the galaxy as a whole integer*4 nr__5 parameter (nr__5=87207) ! Number of records character*121 ar__5 ! Full-size record character*12 index_5 (nr__5) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_5 (nr__5) ! [1,8] Gives detail of the run (G3) real*4 Mass_5 (nr__5) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass_5 (nr__5) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb_5 (nr__5) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2_5 (nr__5) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N_5 (nr__5) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb_5 (nr__5) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np_5 (nr__5) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z_5 (nr__5) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal7.dat' ! RunID7 results relative to the galaxy as a whole integer*4 nr__6 parameter (nr__6=89060) ! Number of records character*121 ar__6 ! Full-size record character*12 index_6 (nr__6) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_6 (nr__6) ! [1,8] Gives detail of the run (G3) real*4 Mass_6 (nr__6) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass_6 (nr__6) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb_6 (nr__6) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2_6 (nr__6) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N_6 (nr__6) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb_6 (nr__6) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np_6 (nr__6) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z_6 (nr__6) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'gal8.dat' ! RunID8 results relative to the galaxy as a whole integer*4 nr__7 parameter (nr__7=89061) ! Number of records character*121 ar__7 ! Full-size record character*12 index_7 (nr__7) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_7 (nr__7) ! [1,8] Gives detail of the run (G3) real*4 Mass_7 (nr__7) ! (Msun) Stellar mass (M_*_) (1) real*4 e_Mass_7 (nr__7) ! (Msun) Statistical error on Mass ({sigma}(M_*_)) (2) real*4 tlb_7 (nr__7) ! (yr) Look-back time of galaxy, assuming a WMAP5 * cosmology real*8 chi2_7 (nr__7) ! {chi}^2^ of the unmasked regions used for * spectral fit real*8 S_N_7 (nr__7) ! Signal-to-noise ratio of the used (unmasked) * spectrum, at the models' resolution integer*4 Nb_7 (nr__7) ! [0,16] Number of recovered bins in galaxy (G2) integer*4 Np_7 (nr__7) ! [0,16] Number of recovered bins with * non-zero mass (G2) real*4 z_7 (nr__7) ! Redshift *Note (1): We use the observed fiber and petrosian magnitudes in the z band to * scale up the stellar mass as: * M_*_=M_*,fiber_*10^0.4(z_p_-fz_p_)^ - Equation (22) * where z_p_ and fz_p_ are the petrosian and fiber magnitudes in the z * band, respectively. *Note (2): We use the full covariance matrix to estimate the statistical * error on the total stellar mass, M_*_: * {sigma}^2^(M_*_)={sigma}^2^(M_*,fiber_){gamma}^2^+ * M^2^_*,fiber_{sigma}^2^({gamma}) - Equation (26) * where {gamma} is the conversion factor between fiber and galaxy mass * of Equation (22) (see note (1)), and {sigma}({gamma}) is the error * associated with this factor, calculated using the errors in the * petrosian and fiber z-band magnitudes. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin1.dat' ! RunID1 results relative to each age bin integer*4 nr__8 parameter (nr__8=5236152) ! Number of records character*88 ar__8 ! Full-size record character*12 index_8 (nr__8) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID (nr__8) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_8 (nr__8) ! [1,8] Gives detail of the run (G3) real*4 Mass_8 (nr__8) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_8 (nr__8) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR (nr__8) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_8 (nr__8) ! ?=0 Metallicity in bin real*4 e_Z (nr__8) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin2.dat' ! RunID2 results relative to each age bin integer*4 nr__9 parameter (nr__9=5140580) ! Number of records character*88 ar__9 ! Full-size record character*12 index_9 (nr__9) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_1 (nr__9) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_9 (nr__9) ! [1,8] Gives detail of the run (G3) real*4 Mass_9 (nr__9) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_9 (nr__9) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR_1 (nr__9) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_9 (nr__9) ! ?=0 Metallicity in bin real*4 e_Z_1 (nr__9) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin3.dat' ! RunID3 results relative to each age bin integer*4 nr__10 parameter (nr__10=5611841) ! Number of records character*88 ar__10 ! Full-size record character*12 index_10 (nr__10) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_2 (nr__10) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_10 (nr__10) ! [1,8] Gives detail of the run (G3) real*4 Mass_10 (nr__10) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_10 (nr__10) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR_2 (nr__10) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_10 (nr__10) ! ?=0 Metallicity in bin real*4 e_Z_2 (nr__10) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin4.dat' ! RunID4 results relative to each age bin integer*4 nr__11 parameter (nr__11=5535239) ! Number of records character*88 ar__11 ! Full-size record character*12 index_11 (nr__11) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_3 (nr__11) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_11 (nr__11) ! [1,8] Gives detail of the run (G3) real*4 Mass_11 (nr__11) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_11 (nr__11) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR_3 (nr__11) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_11 (nr__11) ! ?=0 Metallicity in bin real*4 e_Z_3 (nr__11) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin5.dat' ! RunID5 results relative to each age bin integer*4 nr__12 parameter (nr__12=390329) ! Number of records character*88 ar__12 ! Full-size record character*12 index_12 (nr__12) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_4 (nr__12) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_12 (nr__12) ! [1,8] Gives detail of the run (G3) real*4 Mass_12 (nr__12) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_12 (nr__12) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR_4 (nr__12) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_12 (nr__12) ! ?=0 Metallicity in bin real*4 e_Z_4 (nr__12) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin6.dat' ! RunID6 results relative to each age bin integer*4 nr__13 parameter (nr__13=388292) ! Number of records character*88 ar__13 ! Full-size record character*12 index_13 (nr__13) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_5 (nr__13) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_13 (nr__13) ! [1,8] Gives detail of the run (G3) real*4 Mass_13 (nr__13) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_13 (nr__13) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR_5 (nr__13) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_13 (nr__13) ! ?=0 Metallicity in bin real*4 e_Z_5 (nr__13) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin7.dat' ! RunID7 results relative to each age bin integer*4 nr__14 parameter (nr__14=633730) ! Number of records character*88 ar__14 ! Full-size record character*12 index_14 (nr__14) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_6 (nr__14) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_14 (nr__14) ! [1,8] Gives detail of the run (G3) real*4 Mass_14 (nr__14) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_14 (nr__14) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR_6 (nr__14) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_14 (nr__14) ! ?=0 Metallicity in bin real*4 e_Z_6 (nr__14) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'bin8.dat' ! RunID8 results relative to each age bin integer*4 nr__15 parameter (nr__15=632347) ! Number of records character*88 ar__15 ! Full-size record character*12 index_15 (nr__15) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_7 (nr__15) ! [0,29] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_15 (nr__15) ! [1,8] Gives detail of the run (G3) real*4 Mass_15 (nr__15) ! (Msun) ?=0 Mass formed in bin; corrected for fiber * aperture (1) real*4 e_Mass_15 (nr__15) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); * corrected for fiber aperture (2) real*4 SFR_7 (nr__15) ! (Msun/Gyr) ?=0 Star formation rate in the bin real*8 Z_15 (nr__15) ! ?=0 Metallicity in bin real*4 e_Z_7 (nr__15) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin1.dat' ! RunID1 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__16 parameter (nr__16=10929808) ! Number of records character*74 ar__16 ! Full-size record character*12 index_16 (nr__16) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_8 (nr__16) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_16 (nr__16) ! [1,8] Gives detail of the run (G3) real*4 Mass_16 (nr__16) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_16 (nr__16) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_16 (nr__16) ! ?=0 Metallicity in bin real*4 e_Z_8 (nr__16) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin2.dat' ! RunID2 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__17 parameter (nr__17=10929808) ! Number of records character*74 ar__17 ! Full-size record character*12 index_17 (nr__17) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_9 (nr__17) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_17 (nr__17) ! [1,8] Gives detail of the run (G3) real*4 Mass_17 (nr__17) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_17 (nr__17) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_17 (nr__17) ! ?=0 Metallicity in bin real*4 e_Z_9 (nr__17) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin3.dat' ! RunID3 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__18 parameter (nr__18=10998208) ! Number of records character*74 ar__18 ! Full-size record character*12 index_18 (nr__18) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_10 (nr__18) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_18 (nr__18) ! [1,8] Gives detail of the run (G3) real*4 Mass_18 (nr__18) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_18 (nr__18) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_18 (nr__18) ! ?=0 Metallicity in bin real*4 e_Z_10 (nr__18) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin4.dat' ! RunID4 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__19 parameter (nr__19=10998208) ! Number of records character*74 ar__19 ! Full-size record character*12 index_19 (nr__19) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_11 (nr__19) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_19 (nr__19) ! [1,8] Gives detail of the run (G3) real*4 Mass_19 (nr__19) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_19 (nr__19) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_19 (nr__19) ! ?=0 Metallicity in bin real*4 e_Z_11 (nr__19) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin5.dat' ! RunID5 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__20 parameter (nr__20=1424880) ! Number of records character*74 ar__20 ! Full-size record character*12 index_20 (nr__20) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_12 (nr__20) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_20 (nr__20) ! [1,8] Gives detail of the run (G3) real*4 Mass_20 (nr__20) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_20 (nr__20) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_20 (nr__20) ! ?=0 Metallicity in bin real*4 e_Z_12 (nr__20) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin6.dat' ! RunID6 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__21 parameter (nr__21=1424128) ! Number of records character*74 ar__21 ! Full-size record character*12 index_21 (nr__21) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_13 (nr__21) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_21 (nr__21) ! [1,8] Gives detail of the run (G3) real*4 Mass_21 (nr__21) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_21 (nr__21) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_21 (nr__21) ! ?=0 Metallicity in bin real*4 e_Z_13 (nr__21) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin7.dat' ! RunID7 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__22 parameter (nr__22=1424976) ! Number of records character*74 ar__22 ! Full-size record character*12 index_22 (nr__22) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_14 (nr__22) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_22 (nr__22) ! [1,8] Gives detail of the run (G3) real*4 Mass_22 (nr__22) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_22 (nr__22) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_22 (nr__22) ! ?=0 Metallicity in bin real*4 e_Z_14 (nr__22) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hrbin8.dat' ! RunID8 results relative to each bin into the highest possible resolution (16 bins) integer*4 nr__23 parameter (nr__23=1424976) ! Number of records character*74 ar__23 ! Full-size record character*12 index_23 (nr__23) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 binID_15 (nr__23) ! [0,15] Bin identifier ({alpha}) as given in * Figure 1 (G2) integer*4 runID_23 (nr__23) ! [1,8] Gives detail of the run (G3) real*4 Mass_23 (nr__23) ! (Msun) ?=0 Mass formed in bin; with weights from * Equation (7) and corrected for fiber aperture (1) real*4 e_Mass_23 (nr__23) ! (Msun) ?=0 Mass uncertainty ({sigma}_u({alpha})_); with * weights from Equation (7) and corrected for * fiber aperture (2) real*8 Z_23 (nr__23) ! ?=0 Metallicity in bin real*4 e_Z_15 (nr__23) ! ?=0 {sigma}_Z({alpha})_ as derived from * Equation (24) (3) *Note (1): For any given age bin {alpha}, the mass formed in each bin {alpha} * in units of solar masses is given by the equation (17): * u_{alpha}_=x_{alpha}_4{pi}D_L_^2^(1+z). * See section 2.1.2 for high-resolution bins explanations. *Note (2): For convenience, we also define a covariance matrix of the * unrecycled mass per bin as (using Equation (17)): * C_{alpha}{beta}_(u)=C_{alpha}{beta}_(x)[4{pi}D_L_^2^(1+z)]^2^ * - Equation (25) - from which we can estimate the error in the unscaled * mass formed in bin {alpha} as * {sigma}_u({alpha})_=C_{alpha}{alpha}_(u)^0.5^. * See section 2.1.2 for high-resolution bins explanations. *Note (3): To estimate how much noise affects our recovered solutions we take * a rather empirical approach. For each recovered solution we create * n_error_ random noisy realizations and we apply VESPA to each of these * spectra. In the current runs we have used n_error_=20. We re-bin each * recovered solution in the parametrization of the solution we want to * analyze and estimate the covariance matrices: * C(Z)_{alpha}{beta}=<(Z_{alpha}-Z{bar}_{alpha}_)(Z_{beta}-Z{bar}_{beta}_)> * - equation (24). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust1.dat' ! RunID1 dust information integer*4 nr__24 parameter (nr__24=1366226) ! Number of records character*27 ar__24 ! Full-size record character*12 index_24 (nr__24) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_24 (nr__24) ! [1,8] Gives detail of the run (G3) integer*4 dustID (nr__24) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal (nr__24) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust2.dat' ! RunID2 dust information integer*4 nr__25 parameter (nr__25=1366226) ! Number of records character*27 ar__25 ! Full-size record character*12 index_25 (nr__25) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_25 (nr__25) ! [1,8] Gives detail of the run (G3) integer*4 dustID_1 (nr__25) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal_1 (nr__25) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust3.dat' ! RunID3 dust information integer*4 nr__26 parameter (nr__26=1374776) ! Number of records character*27 ar__26 ! Full-size record character*12 index_26 (nr__26) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_26 (nr__26) ! [1,8] Gives detail of the run (G3) integer*4 dustID_2 (nr__26) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal_2 (nr__26) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust4.dat' ! RunID4 dust information integer*4 nr__27 parameter (nr__27=1374776) ! Number of records character*27 ar__27 ! Full-size record character*12 index_27 (nr__27) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_27 (nr__27) ! [1,8] Gives detail of the run (G3) integer*4 dustID_3 (nr__27) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal_3 (nr__27) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust5.dat' ! RunID5 dust information integer*4 nr__28 parameter (nr__28=178110) ! Number of records character*27 ar__28 ! Full-size record character*12 index_28 (nr__28) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_28 (nr__28) ! [1,8] Gives detail of the run (G3) integer*4 dustID_4 (nr__28) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal_4 (nr__28) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust6.dat' ! RunID6 dust information integer*4 nr__29 parameter (nr__29=178016) ! Number of records character*27 ar__29 ! Full-size record character*12 index_29 (nr__29) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_29 (nr__29) ! [1,8] Gives detail of the run (G3) integer*4 dustID_5 (nr__29) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal_5 (nr__29) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust7.dat' ! RunID7 dust information integer*4 nr__30 parameter (nr__30=178122) ! Number of records character*27 ar__30 ! Full-size record character*12 index_30 (nr__30) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_30 (nr__30) ! [1,8] Gives detail of the run (G3) integer*4 dustID_6 (nr__30) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal_6 (nr__30) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dust8.dat' ! RunID8 dust information integer*4 nr__31 parameter (nr__31=178122) ! Number of records character*27 ar__31 ! Full-size record character*12 index_31 (nr__31) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 runID_31 (nr__31) ! [1,8] Gives detail of the run (G3) integer*4 dustID_7 (nr__31) ! [1/2] Dust identifier: 1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_ real*4 dustVal_7 (nr__31) ! [0,4] Value of dust extinction; either * {tau}^BC^_V_ or {tau}^ISM^_V_, according to * dustID c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'lookup.dat' ! Links between VESPA's unique identifier with SDSS's own identifiers integer*4 nr__32 parameter (nr__32=781788) ! Number of records character*46 ar__32 ! Full-size record character*12 index_32 (nr__32) ! Unique identifier, constructed from SDSS's * plate, MJD and fiberID info (G1) integer*4 MJD (nr__32) ! SDSS-DR7 MJD (=bytes 1-5) integer*4 Plate (nr__32) ! SDSS-DR7 Plate (=bytes 6-9) integer*4 Fiber (nr__32) ! SDSS-DR7 Fiber (=bytes 10-12) character*18 SpecObj (nr__32) ! SDSS-DR7 specObjID * (=plate*2^48^ + MJD*2^32^ + fiber*2^22^) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'runprop.dat' ! Run-specific details (table7 of paper) integer*4 nr__33 parameter (nr__33=8) ! Number of records character*20 ar__33 ! Full-size record integer*4 runID_32 (nr__33) ! [1,8] Unique run ID (G3) character*4 SSP (nr__33) ! Single stellar population (SSP) model used: * BC03 (Bruzual & Charlot 2003MNRAS.344.1000B) or * M05 (Maraston, 2005MNRAS.362..799M) real*4 V1_0 (nr__33) ! [1.0] VESPA code version integer*4 dustID_8 (nr__33) ! [1/2] Dust model used (1 for {tau}^BC^_V_ and * 2 for {tau}^ISM^_V_) character*3 SDSS (nr__33) ! [DR7] SDSS's data release character*3 Sample (nr__33) ! [MGS/LRG] MGS for Main Galaxies Sample or * LRG for Luminous Red Galaxies Sample (LRGS) C============================================================================= C Loading file 'gal1.dat' ! RunID1 results relative to the galaxy as a whole C Format for file interpretation 1 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal1.dat') write(6,*) '....Loading file: gal1.dat' do i__=1,683113 read(1,'(A121)')ar__ read(ar__,1) + index(i__),runID(i__),Mass(i__),e_Mass(i__),tlb(i__), + chi2(i__),S_N(i__),Nb(i__),Np(i__),z(i__) c ..............Just test output........... write(6,1) + index(i__),runID(i__),Mass(i__),e_Mass(i__),tlb(i__), + chi2(i__),S_N(i__),Nb(i__),Np(i__),z(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'gal2.dat' ! RunID2 results relative to the galaxy as a whole C Format for file interpretation 2 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal2.dat') write(6,*) '....Loading file: gal2.dat' do i__=1,681311 read(1,'(A121)')ar__1 read(ar__1,2) + index_1(i__),runID_1(i__),Mass_1(i__),e_Mass_1(i__), + tlb_1(i__),chi2_1(i__),S_N_1(i__),Nb_1(i__),Np_1(i__),z_1(i__) c ..............Just test output........... write(6,2) + index_1(i__),runID_1(i__),Mass_1(i__),e_Mass_1(i__), + tlb_1(i__),chi2_1(i__),S_N_1(i__),Nb_1(i__),Np_1(i__),z_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'gal3.dat' ! RunID3 results relative to the galaxy as a whole C Format for file interpretation 3 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal3.dat') write(6,*) '....Loading file: gal3.dat' do i__=1,687387 read(1,'(A121)')ar__2 read(ar__2,3) + index_2(i__),runID_2(i__),Mass_2(i__),e_Mass_2(i__), + tlb_2(i__),chi2_2(i__),S_N_2(i__),Nb_2(i__),Np_2(i__),z_2(i__) c ..............Just test output........... write(6,3) + index_2(i__),runID_2(i__),Mass_2(i__),e_Mass_2(i__), + tlb_2(i__),chi2_2(i__),S_N_2(i__),Nb_2(i__),Np_2(i__),z_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'gal4.dat' ! RunID4 results relative to the galaxy as a whole C Format for file interpretation 4 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal4.dat') write(6,*) '....Loading file: gal4.dat' do i__=1,685586 read(1,'(A121)')ar__3 read(ar__3,4) + index_3(i__),runID_3(i__),Mass_3(i__),e_Mass_3(i__), + tlb_3(i__),chi2_3(i__),S_N_3(i__),Nb_3(i__),Np_3(i__),z_3(i__) c ..............Just test output........... write(6,4) + index_3(i__),runID_3(i__),Mass_3(i__),e_Mass_3(i__), + tlb_3(i__),chi2_3(i__),S_N_3(i__),Nb_3(i__),Np_3(i__),z_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'gal5.dat' ! RunID5 results relative to the galaxy as a whole C Format for file interpretation 5 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal5.dat') write(6,*) '....Loading file: gal5.dat' do i__=1,89055 read(1,'(A121)')ar__4 read(ar__4,5) + index_4(i__),runID_4(i__),Mass_4(i__),e_Mass_4(i__), + tlb_4(i__),chi2_4(i__),S_N_4(i__),Nb_4(i__),Np_4(i__),z_4(i__) c ..............Just test output........... write(6,5) + index_4(i__),runID_4(i__),Mass_4(i__),e_Mass_4(i__), + tlb_4(i__),chi2_4(i__),S_N_4(i__),Nb_4(i__),Np_4(i__),z_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'gal6.dat' ! RunID6 results relative to the galaxy as a whole C Format for file interpretation 6 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal6.dat') write(6,*) '....Loading file: gal6.dat' do i__=1,87207 read(1,'(A121)')ar__5 read(ar__5,6) + index_5(i__),runID_5(i__),Mass_5(i__),e_Mass_5(i__), + tlb_5(i__),chi2_5(i__),S_N_5(i__),Nb_5(i__),Np_5(i__),z_5(i__) c ..............Just test output........... write(6,6) + index_5(i__),runID_5(i__),Mass_5(i__),e_Mass_5(i__), + tlb_5(i__),chi2_5(i__),S_N_5(i__),Nb_5(i__),Np_5(i__),z_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'gal7.dat' ! RunID7 results relative to the galaxy as a whole C Format for file interpretation 7 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal7.dat') write(6,*) '....Loading file: gal7.dat' do i__=1,89060 read(1,'(A121)')ar__6 read(ar__6,7) + index_6(i__),runID_6(i__),Mass_6(i__),e_Mass_6(i__), + tlb_6(i__),chi2_6(i__),S_N_6(i__),Nb_6(i__),Np_6(i__),z_6(i__) c ..............Just test output........... write(6,7) + index_6(i__),runID_6(i__),Mass_6(i__),e_Mass_6(i__), + tlb_6(i__),chi2_6(i__),S_N_6(i__),Nb_6(i__),Np_6(i__),z_6(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'gal8.dat' ! RunID8 results relative to the galaxy as a whole C Format for file interpretation 8 format( + A12,1X,I1,1X,E14.8,1X,E13.11,1X,E13.8,1X,F13.8,1X,F12.8,1X,I2, + 9X,I2,9X,E14.11) C Effective file loading open(unit=1,status='old',file= +'gal8.dat') write(6,*) '....Loading file: gal8.dat' do i__=1,89061 read(1,'(A121)')ar__7 read(ar__7,8) + index_7(i__),runID_7(i__),Mass_7(i__),e_Mass_7(i__), + tlb_7(i__),chi2_7(i__),S_N_7(i__),Nb_7(i__),Np_7(i__),z_7(i__) c ..............Just test output........... write(6,8) + index_7(i__),runID_7(i__),Mass_7(i__),e_Mass_7(i__), + tlb_7(i__),chi2_7(i__),S_N_7(i__),Nb_7(i__),Np_7(i__),z_7(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin1.dat' ! RunID1 results relative to each age bin C Format for file interpretation 9 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin1.dat') write(6,*) '....Loading file: bin1.dat' do i__=1,5236152 read(1,'(A88)')ar__8 read(ar__8,9) + index_8(i__),binID(i__),runID_8(i__),Mass_8(i__), + e_Mass_8(i__),SFR(i__),Z_8(i__),e_Z(i__) c ..............Just test output........... write(6,9) + index_8(i__),binID(i__),runID_8(i__),Mass_8(i__), + e_Mass_8(i__),SFR(i__),Z_8(i__),e_Z(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin2.dat' ! RunID2 results relative to each age bin C Format for file interpretation 10 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin2.dat') write(6,*) '....Loading file: bin2.dat' do i__=1,5140580 read(1,'(A88)')ar__9 read(ar__9,10) + index_9(i__),binID_1(i__),runID_9(i__),Mass_9(i__), + e_Mass_9(i__),SFR_1(i__),Z_9(i__),e_Z_1(i__) c ..............Just test output........... write(6,10) + index_9(i__),binID_1(i__),runID_9(i__),Mass_9(i__), + e_Mass_9(i__),SFR_1(i__),Z_9(i__),e_Z_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin3.dat' ! RunID3 results relative to each age bin C Format for file interpretation 11 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin3.dat') write(6,*) '....Loading file: bin3.dat' do i__=1,5611841 read(1,'(A88)')ar__10 read(ar__10,11) + index_10(i__),binID_2(i__),runID_10(i__),Mass_10(i__), + e_Mass_10(i__),SFR_2(i__),Z_10(i__),e_Z_2(i__) c ..............Just test output........... write(6,11) + index_10(i__),binID_2(i__),runID_10(i__),Mass_10(i__), + e_Mass_10(i__),SFR_2(i__),Z_10(i__),e_Z_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin4.dat' ! RunID4 results relative to each age bin C Format for file interpretation 12 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin4.dat') write(6,*) '....Loading file: bin4.dat' do i__=1,5535239 read(1,'(A88)')ar__11 read(ar__11,12) + index_11(i__),binID_3(i__),runID_11(i__),Mass_11(i__), + e_Mass_11(i__),SFR_3(i__),Z_11(i__),e_Z_3(i__) c ..............Just test output........... write(6,12) + index_11(i__),binID_3(i__),runID_11(i__),Mass_11(i__), + e_Mass_11(i__),SFR_3(i__),Z_11(i__),e_Z_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin5.dat' ! RunID5 results relative to each age bin C Format for file interpretation 13 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin5.dat') write(6,*) '....Loading file: bin5.dat' do i__=1,390329 read(1,'(A88)')ar__12 read(ar__12,13) + index_12(i__),binID_4(i__),runID_12(i__),Mass_12(i__), + e_Mass_12(i__),SFR_4(i__),Z_12(i__),e_Z_4(i__) c ..............Just test output........... write(6,13) + index_12(i__),binID_4(i__),runID_12(i__),Mass_12(i__), + e_Mass_12(i__),SFR_4(i__),Z_12(i__),e_Z_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin6.dat' ! RunID6 results relative to each age bin C Format for file interpretation 14 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin6.dat') write(6,*) '....Loading file: bin6.dat' do i__=1,388292 read(1,'(A88)')ar__13 read(ar__13,14) + index_13(i__),binID_5(i__),runID_13(i__),Mass_13(i__), + e_Mass_13(i__),SFR_5(i__),Z_13(i__),e_Z_5(i__) c ..............Just test output........... write(6,14) + index_13(i__),binID_5(i__),runID_13(i__),Mass_13(i__), + e_Mass_13(i__),SFR_5(i__),Z_13(i__),e_Z_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin7.dat' ! RunID7 results relative to each age bin C Format for file interpretation 15 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin7.dat') write(6,*) '....Loading file: bin7.dat' do i__=1,633730 read(1,'(A88)')ar__14 read(ar__14,15) + index_14(i__),binID_6(i__),runID_14(i__),Mass_14(i__), + e_Mass_14(i__),SFR_6(i__),Z_14(i__),e_Z_6(i__) c ..............Just test output........... write(6,15) + index_14(i__),binID_6(i__),runID_14(i__),Mass_14(i__), + e_Mass_14(i__),SFR_6(i__),Z_14(i__),e_Z_6(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'bin8.dat' ! RunID8 results relative to each age bin C Format for file interpretation 16 format( + A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,E13.11,1X,F13.11,1X, + E13.11) C Effective file loading open(unit=1,status='old',file= +'bin8.dat') write(6,*) '....Loading file: bin8.dat' do i__=1,632347 read(1,'(A88)')ar__15 read(ar__15,16) + index_15(i__),binID_7(i__),runID_15(i__),Mass_15(i__), + e_Mass_15(i__),SFR_7(i__),Z_15(i__),e_Z_7(i__) c ..............Just test output........... write(6,16) + index_15(i__),binID_7(i__),runID_15(i__),Mass_15(i__), + e_Mass_15(i__),SFR_7(i__),Z_15(i__),e_Z_7(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin1.dat' ! RunID1 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 17 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin1.dat') write(6,*) '....Loading file: hrbin1.dat' do i__=1,10929808 read(1,'(A74)')ar__16 read(ar__16,17) + index_16(i__),binID_8(i__),runID_16(i__),Mass_16(i__), + e_Mass_16(i__),Z_16(i__),e_Z_8(i__) c ..............Just test output........... write(6,17) + index_16(i__),binID_8(i__),runID_16(i__),Mass_16(i__), + e_Mass_16(i__),Z_16(i__),e_Z_8(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin2.dat' ! RunID2 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 18 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin2.dat') write(6,*) '....Loading file: hrbin2.dat' do i__=1,10929808 read(1,'(A74)')ar__17 read(ar__17,18) + index_17(i__),binID_9(i__),runID_17(i__),Mass_17(i__), + e_Mass_17(i__),Z_17(i__),e_Z_9(i__) c ..............Just test output........... write(6,18) + index_17(i__),binID_9(i__),runID_17(i__),Mass_17(i__), + e_Mass_17(i__),Z_17(i__),e_Z_9(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin3.dat' ! RunID3 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 19 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin3.dat') write(6,*) '....Loading file: hrbin3.dat' do i__=1,10998208 read(1,'(A74)')ar__18 read(ar__18,19) + index_18(i__),binID_10(i__),runID_18(i__),Mass_18(i__), + e_Mass_18(i__),Z_18(i__),e_Z_10(i__) c ..............Just test output........... write(6,19) + index_18(i__),binID_10(i__),runID_18(i__),Mass_18(i__), + e_Mass_18(i__),Z_18(i__),e_Z_10(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin4.dat' ! RunID4 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 20 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin4.dat') write(6,*) '....Loading file: hrbin4.dat' do i__=1,10998208 read(1,'(A74)')ar__19 read(ar__19,20) + index_19(i__),binID_11(i__),runID_19(i__),Mass_19(i__), + e_Mass_19(i__),Z_19(i__),e_Z_11(i__) c ..............Just test output........... write(6,20) + index_19(i__),binID_11(i__),runID_19(i__),Mass_19(i__), + e_Mass_19(i__),Z_19(i__),e_Z_11(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin5.dat' ! RunID5 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 21 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin5.dat') write(6,*) '....Loading file: hrbin5.dat' do i__=1,1424880 read(1,'(A74)')ar__20 read(ar__20,21) + index_20(i__),binID_12(i__),runID_20(i__),Mass_20(i__), + e_Mass_20(i__),Z_20(i__),e_Z_12(i__) c ..............Just test output........... write(6,21) + index_20(i__),binID_12(i__),runID_20(i__),Mass_20(i__), + e_Mass_20(i__),Z_20(i__),e_Z_12(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin6.dat' ! RunID6 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 22 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin6.dat') write(6,*) '....Loading file: hrbin6.dat' do i__=1,1424128 read(1,'(A74)')ar__21 read(ar__21,22) + index_21(i__),binID_13(i__),runID_21(i__),Mass_21(i__), + e_Mass_21(i__),Z_21(i__),e_Z_13(i__) c ..............Just test output........... write(6,22) + index_21(i__),binID_13(i__),runID_21(i__),Mass_21(i__), + e_Mass_21(i__),Z_21(i__),e_Z_13(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin7.dat' ! RunID7 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 23 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin7.dat') write(6,*) '....Loading file: hrbin7.dat' do i__=1,1424976 read(1,'(A74)')ar__22 read(ar__22,23) + index_22(i__),binID_14(i__),runID_22(i__),Mass_22(i__), + e_Mass_22(i__),Z_22(i__),e_Z_14(i__) c ..............Just test output........... write(6,23) + index_22(i__),binID_14(i__),runID_22(i__),Mass_22(i__), + e_Mass_22(i__),Z_22(i__),e_Z_14(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hrbin8.dat' ! RunID8 results relative to each bin into the * highest possible resolution (16 bins) C Format for file interpretation 24 format(A12,1X,I2,1X,I1,1X,E14.10,1X,E13.11,1X,F13.11,1X,E13.11) C Effective file loading open(unit=1,status='old',file= +'hrbin8.dat') write(6,*) '....Loading file: hrbin8.dat' do i__=1,1424976 read(1,'(A74)')ar__23 read(ar__23,24) + index_23(i__),binID_15(i__),runID_23(i__),Mass_23(i__), + e_Mass_23(i__),Z_23(i__),e_Z_15(i__) c ..............Just test output........... write(6,24) + index_23(i__),binID_15(i__),runID_23(i__),Mass_23(i__), + e_Mass_23(i__),Z_23(i__),e_Z_15(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust1.dat' ! RunID1 dust information C Format for file interpretation 25 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust1.dat') write(6,*) '....Loading file: dust1.dat' do i__=1,1366226 read(1,'(A27)')ar__24 read(ar__24,25) + index_24(i__),runID_24(i__),dustID(i__),dustVal(i__) c ..............Just test output........... write(6,25) + index_24(i__),runID_24(i__),dustID(i__),dustVal(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust2.dat' ! RunID2 dust information C Format for file interpretation 26 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust2.dat') write(6,*) '....Loading file: dust2.dat' do i__=1,1366226 read(1,'(A27)')ar__25 read(ar__25,26) + index_25(i__),runID_25(i__),dustID_1(i__),dustVal_1(i__) c ..............Just test output........... write(6,26) + index_25(i__),runID_25(i__),dustID_1(i__),dustVal_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust3.dat' ! RunID3 dust information C Format for file interpretation 27 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust3.dat') write(6,*) '....Loading file: dust3.dat' do i__=1,1374776 read(1,'(A27)')ar__26 read(ar__26,27) + index_26(i__),runID_26(i__),dustID_2(i__),dustVal_2(i__) c ..............Just test output........... write(6,27) + index_26(i__),runID_26(i__),dustID_2(i__),dustVal_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust4.dat' ! RunID4 dust information C Format for file interpretation 28 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust4.dat') write(6,*) '....Loading file: dust4.dat' do i__=1,1374776 read(1,'(A27)')ar__27 read(ar__27,28) + index_27(i__),runID_27(i__),dustID_3(i__),dustVal_3(i__) c ..............Just test output........... write(6,28) + index_27(i__),runID_27(i__),dustID_3(i__),dustVal_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust5.dat' ! RunID5 dust information C Format for file interpretation 29 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust5.dat') write(6,*) '....Loading file: dust5.dat' do i__=1,178110 read(1,'(A27)')ar__28 read(ar__28,29) + index_28(i__),runID_28(i__),dustID_4(i__),dustVal_4(i__) c ..............Just test output........... write(6,29) + index_28(i__),runID_28(i__),dustID_4(i__),dustVal_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust6.dat' ! RunID6 dust information C Format for file interpretation 30 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust6.dat') write(6,*) '....Loading file: dust6.dat' do i__=1,178016 read(1,'(A27)')ar__29 read(ar__29,30) + index_29(i__),runID_29(i__),dustID_5(i__),dustVal_5(i__) c ..............Just test output........... write(6,30) + index_29(i__),runID_29(i__),dustID_5(i__),dustVal_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust7.dat' ! RunID7 dust information C Format for file interpretation 31 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust7.dat') write(6,*) '....Loading file: dust7.dat' do i__=1,178122 read(1,'(A27)')ar__30 read(ar__30,31) + index_30(i__),runID_30(i__),dustID_6(i__),dustVal_6(i__) c ..............Just test output........... write(6,31) + index_30(i__),runID_30(i__),dustID_6(i__),dustVal_6(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dust8.dat' ! RunID8 dust information C Format for file interpretation 32 format(A12,1X,I1,1X,I1,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'dust8.dat') write(6,*) '....Loading file: dust8.dat' do i__=1,178122 read(1,'(A27)')ar__31 read(ar__31,32) + index_31(i__),runID_31(i__),dustID_7(i__),dustVal_7(i__) c ..............Just test output........... write(6,32) + index_31(i__),runID_31(i__),dustID_7(i__),dustVal_7(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'lookup.dat' ! Links between VESPA's unique identifier with * SDSS's own identifiers C Format for file interpretation 33 format(A12,1X,I5,1X,I4,1X,I3,1X,A18) C Effective file loading open(unit=1,status='old',file= +'lookup.dat') write(6,*) '....Loading file: lookup.dat' do i__=1,781788 read(1,'(A46)')ar__32 read(ar__32,33) + index_32(i__),MJD(i__),Plate(i__),Fiber(i__),SpecObj(i__) c ..............Just test output........... write(6,33) + index_32(i__),MJD(i__),Plate(i__),Fiber(i__),SpecObj(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'runprop.dat' ! Run-specific details (table7 of paper) C Format for file interpretation 34 format(I1,1X,A4,1X,F3.1,1X,I1,1X,A3,1X,A3) C Effective file loading open(unit=1,status='old',file= +'runprop.dat') write(6,*) '....Loading file: runprop.dat' do i__=1,8 read(1,'(A20)')ar__33 read(ar__33,34) + runID_32(i__),SSP(i__),V1_0(i__),dustID_8(i__),SDSS(i__), + Sample(i__) c ..............Just test output........... write(6,34) + runID_32(i__),SSP(i__),V1_0(i__),dustID_8(i__),SDSS(i__), + Sample(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end