FORTRAN Generation
(/./ftp/cats/J/AJ/131/1648)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/131/1648 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-Mar-29
C=============================================================================
*  This code was generated from the ReadMe file documenting a catalogue
*  according to the "Standard for Documentation of Astronomical Catalogues"
*  currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
*  (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
*  Please report problems or questions to   
C=============================================================================

      implicit none
*  Unspecified or NULL values, generally corresponding to blank columns,
*  are assigned one of the following special values:
*     rNULL__    for unknown or NULL floating-point values
*     iNULL__    for unknown or NULL   integer      values
      real*4     rNULL__
      integer*4  iNULL__
      parameter  (rNULL__=--2147483648.)  	! NULL real number
      parameter  (iNULL__=(-2147483647-1))	! NULL int  number
      integer    idig			! testing NULL number

C=============================================================================
Cat. J/AJ/131/1648    Type Ia supernovae at high and low redshifts   (Blondin+, 2006)
*================================================================================
*Using line profiles to test the fraternity of type Ia supernovae at high and low
*redshifts.
*    Blondin S., Dessart L., Leibundgut B., Branch D., Hoflich P., Tonry J.L.,
*    Matheson T., Foley R.J., Chornock R., Filippenko A.V., Sollerman J.,
*    Spyromilio J., Kirshner R.P., Wood-Vasey W.M., Clocchiatti A., Aguilera C.,
*    Barris B., Becker A.C., Challis P., Covarrubias R., Davis T.M.,
*    Garnavich P., Hicken M., Jha S., Krisciunas K., Li W., Miceli A.,
*    Miknaitis G., Pignata G., Prieto J.L., Rest A., Riess A.G., Salvo M.E.,
*    Schmidt B.P., Smith R.C., Stubbs C.W., Suntzeff N.B.
*   <Astron. J., 131, 1648-1666 (2006)>
*   =2006AJ....131.1648B
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Local SN Ia data

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

      character*6   SN         (nr__) ! The IAU Name of the supernova
      character*13  Galaxy     (nr__) ! Host galaxy (1)
      character*3   M          (nr__) ! Messier name for the host galaxy
      integer*4     cz         (nr__) ! (km/s) Heliocentric radial velocity of host galaxy (2)
      character*1   n_cz       (nr__) ! [a] Indicates cz from Riess et al.
*                                  (1999AJ....117..707R), Table 3
      real*4        E_B_V      (nr__) ! (mag) Average host-galaxy reddening (3)
      real*4        e_E_B_V    (nr__) ! (mag) ? The 1{sigma} error in E(B-V)
      character*1   n_E_B_V    (nr__) ! Note on E(B-V) (4)
      real*8        BMAG       (nr__) ! (mag) Absolute B band magnitude at maximum (6)
      real*4        e_BMAG     (nr__) ! (mag) The 1{sigma} error in BMag
      real*4        Delm15     (nr__) ! (mag) Decline in B between maximum and +15 days (7)
      real*4        e_Delm15   (nr__) ! (mag) ? The 1{sigma} error in Delm15
      character*1   n_Delm15   (nr__) ! Flag on Delm15 (4)
      real*8        HJD        (nr__) ! (d) Heliocentric Julian Date of maximum (8)
      real*4        e_HJD      (nr__) ! (d) The 1{sigma} error in HJD
      character*47  Phase      (nr__) ! Rest-frame phase(s) (G2) (9)
      character*5   Refs       (nr__) ! Reference(s) (10)
*Note (1): Anon, or "Anonymous", means the galaxy has no official designation.
*Note (2): As quoted in the NASA/IPAC Extragalactic Database (NED); we assume
*          an error of 300km/s in the galaxy redshift.
*Note (3): As quoted in Table 2, column (7) of Phillips et al.
*          (1999AJ....118.1766P).
*Note (4): Flag on E(B6V), defined as follows:
*    b = From Li et al. (1999AJ....117.2709L); the host galaxy E(B-V) was
*        determined from the total E(B-V) quoted in Li et al.
*        (1999AJ....117.2709L) and a Galactic reddening of 0.11mag, as derived
*        using the dust IR emission maps of Schlegel et al.
*        (1998ApJ...500..525S).
*    c = From Turatto et al (1998MNRAS.283....1T).
*    d = From Saha et al. (2001ApJ...562..314S). Note: this value assumes the
*        supernova experiences the same reddening as the Cepheids in NGC 5584;
*        see Saha et al. 2001ApJ...562..314S for a thorough discussion of this
*        issue.
*    e = From Riess et al. (1999AJ....118.2675R).
*    f = From Krisciunas et al. (2000ApJ...539..658K).
*    g = From Garnavich et al. (2004ApJ...613.1120G); see this same paper for
*        a long discussion on the low/null reddening in the host galaxy.
*    h = From Stritzinger et al. (2002AJ....124.2100S).
*    i = From Li et al. (2001PASP..113.1178L).
*    j = From Benetti et al. (2004MNRAS.348..261B); consult this article and
*        Krisciunas et al. (2004AJ....128.3034K) for an extensive discussion
*        on the abnormally high E(B-V) value obtained for SN 2002bo.
*Note (6): Derived from an MLCS2k2 fit to the rest-frame B-band light curves
*          (Jha 2002, Ph.D. thesis, Harvard Univ.; Jha et al. 2007, ApJ,
*          submitted, astro-ph/0612666).
*Note (7): As quoted in Table 2, column (2) of Phillips et al.
*          (1999AJ....118.1766P).
*Note (8): Derived from an MLCS2k2 fit to the rest-frame B-band light curves.
*Note (9): Adjacent phases are listed in between square brackets; a "(2)"
*          indicates that two spectra correspond to a same rounded phase.
*Note (10): Reference of refereed articles presenting optical spectroscopic
*          data in the range -15<=Phase<=+30 [days from B-band maximum], 
*          as found on the NASA ADS server (see refs.dat); may be blank.

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

C  Declarations for 'table3.dat'	! High-z SN Ia data

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

      character*6   SN_1       (nr__1) ! The IAU name
      real*4        z          (nr__1) ! Redshift (1)
      real*4        E_B_V_1    (nr__1) ! (mag) The (B-V) color excess (2)
      real*4        e_E_B_V_1  (nr__1) ! (mag) The 1{sigma} error in E(B-V)
      real*4        Delm15_1   (nr__1) ! (mag) Decline in B between maximum and +15 days (2)
      real*4        e_Delm15_1 (nr__1) ! (mag) The 1{sigma} error in Delm15
      real*8        MJD        (nr__1) ! (d) Modified Julian Date of maximum (-2400000) (2)
      real*4        e_MJD      (nr__1) ! (d) The 1{sigma} error in MJD
      character*12  Phase_1    (nr__1) ! Rest-frame phase(s) (G2)
*Note (1): For redshifts determined from narrow lines in the host galaxy (quoted
*          with three decimal places), we assume an error of ~0.001; for those
*          determined via cross-correlations with local SN Ia spectra templates
*          Matheson et al. (2005AJ....129.2352M), we assume an error of 0.01.
*Note (2): Light-curve parameters output by the fitting routine of
*          Prieto et al. (2006ApJ...647..501P). E(B-V) values correspond to
*          host-galaxy extinction only.

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

C  Declarations for 'refs.dat'	! References

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

      integer*4     Ref        (nr__2) ! Reference number
      character*19  Bibcode    (nr__2) ! Bibcode
      character*22  Aut        (nr__2) ! Authors
      character*47  Com        (nr__2) ! Comments

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

C  Declarations for 'table5.dat'	! Absorption velocities in Local SNe Ia (10^3^km/s)

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

      character*6   SN_2       (nr__3) ! SN identification
      character*1   f_SN       (nr__3) ! [cD] Flag on SN (G1)
      integer*4     Phase_2    (nr__3) ! Rest-frame phase (G2)
      real*4        CaII_B     (nr__3) ! (Mm/s) ? CaII(3945) Blue absorption velocity (1)
      real*4        E_CaII_B   (nr__3) ! (Mm/s) ? Upper limit error in CaII-B
      real*4        e_CaII_B_1 (nr__3) ! (Mm/s) ? Lower limit error or error in CaII-B
      character*1   f_CaII_B   (nr__3) ! Flag on CaII-B (G3)
      real*4        CaII_R     (nr__3) ! (Mm/s) ? CaII(3945) Red/Single absorption velocity (1)
      real*4        E_CaII_R   (nr__3) ! (Mm/s) ? Upper limit error in CaII-R
      real*4        e_CaII_R_1 (nr__3) ! (Mm/s) ? Lower limit error or error in CaII-R
      character*1   f_CaII_R   (nr__3) ! Flag on CaII-R (G3)
      real*4        SII5454    (nr__3) ! (Mm/s) ? SII(5454) absorption velocity
      real*4        e_SII5454  (nr__3) ! (Mm/s) ? Error in SII5454
      character*1   f_SII5454  (nr__3) ! Flag on SII5454 (G3)
      real*4        SII5640    (nr__3) ! (Mm/s) ? SII(5640) absorption velocity
      real*4        E_SII5640  (nr__3) ! (Mm/s) ? Upper limit error in SII5640
      real*4        e_SII5640_1(nr__3) ! (Mm/s) ? Lower limit error or error in SII5640
      character*1   f_SII5640  (nr__3) ! Flag on SII5640 (G3)
      real*4        SiII6355   (nr__3) ! (Mm/s) ? SiII(6355) absorption velocity
      real*4        E_SiII6355 (nr__3) ! (Mm/s) ? Upper limit error in SiII6355
      real*4        e_SiII6355_1(nr__3) ! (Mm/s) ? Lower limit error or error in SiII6355
      character*1   f_SiII6355 (nr__3) ! Flag on SiII6355 (G3)
*Note (1): v_abs_ measurements for CaII(3945) are separated into "blue" and
*          "red" components (when a double absorption is present). Both "red"
*          and single-absorption velocities are reported in the same column.

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

C  Declarations for 'table6.dat'	! Absorption velocities in ESSENCE high-z SN Ia

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

      character*6   SN_3       (nr__4) ! SN identification
      character*1   f_SN_1     (nr__4) ! [cD] Flag on SN (G1)
      integer*4     Phase_3    (nr__4) ! Rest-frame phase (G2)
      real*4        CaII_B_1   (nr__4) ! (Mm/s) ? CaII(3945) Blue absorption velocity (1)
      real*4        E_CaII_B_2 (nr__4) ! (Mm/s) ? Upper limit error in CaII-B
      real*4        e_CaII_B_3 (nr__4) ! (Mm/s) ? Lower limit error or error in CaII-B
      character*1   f_CaII_B_1 (nr__4) ! Flag on CaII-B (G3)
      real*4        CaII_R_1   (nr__4) ! (Mm/s) ? CaII(3945) Red/Single absorption velocity (1)
      real*4        E_CaII_R_2 (nr__4) ! (Mm/s) ? Upper limit error in CaII-R
      real*4        e_CaII_R_3 (nr__4) ! (Mm/s) ? Lower limit error or error in CaII-R
      character*1   f_CaII_R_1 (nr__4) ! Flag on CaII-R (G3)
      real*4        SII5454_1  (nr__4) ! (Mm/s) ? SII(5454) absorption velocity
      real*4        e_SII5454_1(nr__4) ! (Mm/s) ? Error in SII5454
      character*1   f_SII5454_1(nr__4) ! Flag on SII5454 (G3)
      real*4        SII5640_1  (nr__4) ! (Mm/s) ? SII(5640) absorption velocity
      real*4        E_SII5640_2(nr__4) ! (Mm/s) ? Upper limit error in SII5640
      real*4        e_SII5640_3(nr__4) ! (Mm/s) ? Lower limit error or error in SII5640
      character*1   f_SII5640_1(nr__4) ! Flag on SII5640 (G3)
      real*4        SiII6355_1 (nr__4) ! (Mm/s) ? SiII(6355) absorption velocity
      real*4        E_SiII6355_2(nr__4) ! (Mm/s) ? Upper limit error in SiII6355
      real*4        e_SiII6355_3(nr__4) ! (Mm/s) ? Lower limit error or error in SiII6355
      character*1   f_SiII6355_1(nr__4) ! Flag on SiII6355 (G3)
*Note (1): v_abs_ measurements for CaII(3945) are separated into "blue" and
*          "red" components (when a double absorption is present). Both "red"
*          and single-absorption velocities are reported in the same column.

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

C  Declarations for 'table7.dat'	! Emission-peak velocities in local SN Ia

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

      character*6   SN_4       (nr__5) ! SN identification
      character*1   f_SN_2     (nr__5) ! [D] Flag on SN (G1)
      integer*4     Phase_4    (nr__5) ! Rest-frame phase (G2)
      real*4        CaII3945   (nr__5) ! (Mm/s) ? Ca II (3945) emission-peak velocity
      real*4        e_CaII3945 (nr__5) ! (Mm/s) ? Error in CaII3945
      character*1   f_CaII3945 (nr__5) ! Flag on CaII3945 (G3)
      real*4        SII5454_2  (nr__5) ! (Mm/s) ? S II (5454) emission-peak velocity
      real*4        e_SII5454_2(nr__5) ! (Mm/s) ? Error in SII5454
      character*1   f_SII5454_2(nr__5) ! Flag on SII5454 (G3)
      real*4        SII5640_2  (nr__5) ! (Mm/s) ? S II (5640) emission-peak velocity
      real*4        e_SII5640_4(nr__5) ! (Mm/s) ? Error in SII5640
      character*1   f_SII5640_2(nr__5) ! Flag on SII5640 (G3)
      real*4        SiII6355_2 (nr__5) ! (Mm/s) ? Si II (6355) emission-peak velocity
      real*4        e_SiII6355_4(nr__5) ! (Mm/s) ? Error in SiII6355
      character*1   f_SiII6355_2(nr__5) ! Flag on SiII6355 (G3)

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

C  Declarations for 'table8.dat'	! Emission-peak velocities in ESSENCE high-z SN Ia

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

      character*6   SN_5       (nr__6) ! SN identification
      character*1   f_SN_3     (nr__6) ! [D] Flag on SN (G1)
      integer*4     Phase_5    (nr__6) ! Rest-frame phase (G2)
      real*4        CaII3945_1 (nr__6) ! (Mm/s) ? Ca II (3945) emission-peak velocity
      real*4        e_CaII3945_1(nr__6) ! (Mm/s) ? Error in CaII3945
      character*1   f_CaII3945_1(nr__6) ! Flag on CaII3945 (G3)
      real*4        SII5454_3  (nr__6) ! (Mm/s) ? S II (5454) emission-peak velocity
      real*4        e_SII5454_3(nr__6) ! (Mm/s) ? Error in SII5454
      character*1   f_SII5454_3(nr__6) ! Flag on SII5454 (G3)
      real*4        SII5640_3  (nr__6) ! (Mm/s) ? S II (5640) emission-peak velocity
      real*4        e_SII5640_5(nr__6) ! (Mm/s) ? Error in SII5640
      character*1   f_SII5640_3(nr__6) ! Flag on SII5640 (G3)
      real*4        SiII6355_3 (nr__6) ! (Mm/s) ? Si II (6355) emission-peak velocity
      real*4        e_SiII6355_5(nr__6) ! (Mm/s) ? Error in SiII6355
      character*1   f_SiII6355_3(nr__6) ! Flag on SiII6355 (G3)

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

C  Loading file 'table2.dat'	! Local SN Ia data

C  Format for file interpretation

    1 format(
     +  A6,1X,A13,1X,A3,1X,I5,1X,A1,1X,F4.2,1X,F4.2,1X,A1,1X,F7.3,1X,
     +  F5.3,1X,F4.2,1X,F4.2,1X,A1,1X,F10.2,1X,F4.2,1X,A47,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,30
        read(1,'(A140)')ar__
        read(ar__,1)
     +  SN(i__),Galaxy(i__),M(i__),cz(i__),n_cz(i__),E_B_V(i__),
     +  e_E_B_V(i__),n_E_B_V(i__),BMAG(i__),e_BMAG(i__),Delm15(i__),
     +  e_Delm15(i__),n_Delm15(i__),HJD(i__),e_HJD(i__),Phase(i__),
     +  Refs(i__)
        if(ar__(39:42) .EQ. '') e_E_B_V(i__) = rNULL__
        if(ar__(65:68) .EQ. '') e_Delm15(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  SN(i__),Galaxy(i__),M(i__),cz(i__),n_cz(i__),E_B_V(i__),
     +  e_E_B_V(i__),n_E_B_V(i__),BMAG(i__),e_BMAG(i__),Delm15(i__),
     +  e_Delm15(i__),n_Delm15(i__),HJD(i__),e_HJD(i__),Phase(i__),
     +  Refs(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! High-z SN Ia data

C  Format for file interpretation

    2 format(
     +  A6,1X,F5.3,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F8.2,1X,F4.2,1X,
     +  A12)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,38
        read(1,'(A59)')ar__1
        read(ar__1,2)
     +  SN_1(i__),z(i__),E_B_V_1(i__),e_E_B_V_1(i__),Delm15_1(i__),
     +  e_Delm15_1(i__),MJD(i__),e_MJD(i__),Phase_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  SN_1(i__),z(i__),E_B_V_1(i__),e_E_B_V_1(i__),Delm15_1(i__),
     +  e_Delm15_1(i__),MJD(i__),e_MJD(i__),Phase_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    3 format(I2,1X,A19,1X,A22,1X,A47)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,37
        read(1,'(A93)')ar__2
        read(ar__2,3)Ref(i__),Bibcode(i__),Aut(i__),Com(i__)
c    ..............Just test output...........
        write(6,3)Ref(i__),Bibcode(i__),Aut(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Absorption velocities in Local SNe Ia (10^3^km/s)

C  Format for file interpretation

    4 format(
     +  A6,1X,A1,1X,I3,1X,F5.1,1X,F3.1,1X,F3.1,1X,A1,1X,F5.1,1X,F3.1,
     +  1X,F3.1,1X,A1,1X,F5.1,1X,F3.1,1X,A1,1X,F5.1,1X,F3.1,1X,F3.1,
     +  1X,A1,1X,F5.1,1X,F3.1,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,229
        read(1,'(A88)')ar__3
        read(ar__3,4)
     +  SN_2(i__),f_SN(i__),Phase_2(i__),CaII_B(i__),E_CaII_B(i__),
     +  e_CaII_B_1(i__),f_CaII_B(i__),CaII_R(i__),E_CaII_R(i__),
     +  e_CaII_R_1(i__),f_CaII_R(i__),SII5454(i__),e_SII5454(i__),
     +  f_SII5454(i__),SII5640(i__),E_SII5640(i__),e_SII5640_1(i__),
     +  f_SII5640(i__),SiII6355(i__),E_SiII6355(i__),
     +  e_SiII6355_1(i__),f_SiII6355(i__)
        if(ar__3(14:18) .EQ. '') CaII_B(i__) = rNULL__
        if(ar__3(20:22) .EQ. '') E_CaII_B(i__) = rNULL__
        if(ar__3(24:26) .EQ. '') e_CaII_B_1(i__) = rNULL__
        if(ar__3(30:34) .EQ. '') CaII_R(i__) = rNULL__
        if(ar__3(36:38) .EQ. '') E_CaII_R(i__) = rNULL__
        if(ar__3(40:42) .EQ. '') e_CaII_R_1(i__) = rNULL__
        if(ar__3(46:50) .EQ. '') SII5454(i__) = rNULL__
        if(ar__3(52:54) .EQ. '') e_SII5454(i__) = rNULL__
        if(ar__3(58:62) .EQ. '') SII5640(i__) = rNULL__
        if(ar__3(64:66) .EQ. '') E_SII5640(i__) = rNULL__
        if(ar__3(68:70) .EQ. '') e_SII5640_1(i__) = rNULL__
        if(ar__3(74:78) .EQ. '') SiII6355(i__) = rNULL__
        if(ar__3(80:82) .EQ. '') E_SiII6355(i__) = rNULL__
        if(ar__3(84:86) .EQ. '') e_SiII6355_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  SN_2(i__),f_SN(i__),Phase_2(i__),CaII_B(i__),E_CaII_B(i__),
     +  e_CaII_B_1(i__),f_CaII_B(i__),CaII_R(i__),E_CaII_R(i__),
     +  e_CaII_R_1(i__),f_CaII_R(i__),SII5454(i__),e_SII5454(i__),
     +  f_SII5454(i__),SII5640(i__),E_SII5640(i__),e_SII5640_1(i__),
     +  f_SII5640(i__),SiII6355(i__),E_SiII6355(i__),
     +  e_SiII6355_1(i__),f_SiII6355(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Absorption velocities in ESSENCE high-z SN Ia

C  Format for file interpretation

    5 format(
     +  A6,1X,A1,1X,I3,1X,F5.1,1X,F3.1,1X,F3.1,1X,A1,1X,F5.1,1X,F3.1,
     +  1X,F3.1,1X,A1,1X,F5.1,1X,F3.1,1X,A1,1X,F5.1,1X,F3.1,1X,F3.1,
     +  1X,A1,1X,F5.1,1X,F3.1,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,48
        read(1,'(A88)')ar__4
        read(ar__4,5)
     +  SN_3(i__),f_SN_1(i__),Phase_3(i__),CaII_B_1(i__),
     +  E_CaII_B_2(i__),e_CaII_B_3(i__),f_CaII_B_1(i__),CaII_R_1(i__),
     +  E_CaII_R_2(i__),e_CaII_R_3(i__),f_CaII_R_1(i__),
     +  SII5454_1(i__),e_SII5454_1(i__),f_SII5454_1(i__),
     +  SII5640_1(i__),E_SII5640_2(i__),e_SII5640_3(i__),
     +  f_SII5640_1(i__),SiII6355_1(i__),E_SiII6355_2(i__),
     +  e_SiII6355_3(i__),f_SiII6355_1(i__)
        if(ar__4(14:18) .EQ. '') CaII_B_1(i__) = rNULL__
        if(ar__4(20:22) .EQ. '') E_CaII_B_2(i__) = rNULL__
        if(ar__4(24:26) .EQ. '') e_CaII_B_3(i__) = rNULL__
        if(ar__4(30:34) .EQ. '') CaII_R_1(i__) = rNULL__
        if(ar__4(36:38) .EQ. '') E_CaII_R_2(i__) = rNULL__
        if(ar__4(40:42) .EQ. '') e_CaII_R_3(i__) = rNULL__
        if(ar__4(46:50) .EQ. '') SII5454_1(i__) = rNULL__
        if(ar__4(52:54) .EQ. '') e_SII5454_1(i__) = rNULL__
        if(ar__4(58:62) .EQ. '') SII5640_1(i__) = rNULL__
        if(ar__4(64:66) .EQ. '') E_SII5640_2(i__) = rNULL__
        if(ar__4(68:70) .EQ. '') e_SII5640_3(i__) = rNULL__
        if(ar__4(74:78) .EQ. '') SiII6355_1(i__) = rNULL__
        if(ar__4(80:82) .EQ. '') E_SiII6355_2(i__) = rNULL__
        if(ar__4(84:86) .EQ. '') e_SiII6355_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  SN_3(i__),f_SN_1(i__),Phase_3(i__),CaII_B_1(i__),
     +  E_CaII_B_2(i__),e_CaII_B_3(i__),f_CaII_B_1(i__),CaII_R_1(i__),
     +  E_CaII_R_2(i__),e_CaII_R_3(i__),f_CaII_R_1(i__),
     +  SII5454_1(i__),e_SII5454_1(i__),f_SII5454_1(i__),
     +  SII5640_1(i__),E_SII5640_2(i__),e_SII5640_3(i__),
     +  f_SII5640_1(i__),SiII6355_1(i__),E_SiII6355_2(i__),
     +  e_SiII6355_3(i__),f_SiII6355_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Emission-peak velocities in local SN Ia

C  Format for file interpretation

    6 format(
     +  A6,1X,A1,1X,I3,1X,F4.1,1X,F3.1,1X,A1,1X,F4.1,1X,F3.1,1X,A1,1X,
     +  F4.1,1X,F3.1,1X,A1,1X,F4.1,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,178
        read(1,'(A56)')ar__5
        read(ar__5,6)
     +  SN_4(i__),f_SN_2(i__),Phase_4(i__),CaII3945(i__),
     +  e_CaII3945(i__),f_CaII3945(i__),SII5454_2(i__),
     +  e_SII5454_2(i__),f_SII5454_2(i__),SII5640_2(i__),
     +  e_SII5640_4(i__),f_SII5640_2(i__),SiII6355_2(i__),
     +  e_SiII6355_4(i__),f_SiII6355_2(i__)
        if(ar__5(14:17) .EQ. '') CaII3945(i__) = rNULL__
        if(ar__5(19:21) .EQ. '') e_CaII3945(i__) = rNULL__
        if(ar__5(25:28) .EQ. '') SII5454_2(i__) = rNULL__
        if(ar__5(30:32) .EQ. '') e_SII5454_2(i__) = rNULL__
        if(ar__5(36:39) .EQ. '') SII5640_2(i__) = rNULL__
        if(ar__5(41:43) .EQ. '') e_SII5640_4(i__) = rNULL__
        if(ar__5(47:50) .EQ. '') SiII6355_2(i__) = rNULL__
        if(ar__5(52:54) .EQ. '') e_SiII6355_4(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  SN_4(i__),f_SN_2(i__),Phase_4(i__),CaII3945(i__),
     +  e_CaII3945(i__),f_CaII3945(i__),SII5454_2(i__),
     +  e_SII5454_2(i__),f_SII5454_2(i__),SII5640_2(i__),
     +  e_SII5640_4(i__),f_SII5640_2(i__),SiII6355_2(i__),
     +  e_SiII6355_4(i__),f_SiII6355_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Emission-peak velocities in ESSENCE high-z SN Ia

C  Format for file interpretation

    7 format(
     +  A6,1X,A1,1X,I3,1X,F4.1,1X,F3.1,1X,A1,1X,F4.1,1X,F3.1,1X,A1,1X,
     +  F4.1,1X,F3.1,1X,A1,1X,F4.1,1X,F3.1,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,39
        read(1,'(A56)')ar__6
        read(ar__6,7)
     +  SN_5(i__),f_SN_3(i__),Phase_5(i__),CaII3945_1(i__),
     +  e_CaII3945_1(i__),f_CaII3945_1(i__),SII5454_3(i__),
     +  e_SII5454_3(i__),f_SII5454_3(i__),SII5640_3(i__),
     +  e_SII5640_5(i__),f_SII5640_3(i__),SiII6355_3(i__),
     +  e_SiII6355_5(i__),f_SiII6355_3(i__)
        if(ar__6(14:17) .EQ. '') CaII3945_1(i__) = rNULL__
        if(ar__6(19:21) .EQ. '') e_CaII3945_1(i__) = rNULL__
        if(ar__6(25:28) .EQ. '') SII5454_3(i__) = rNULL__
        if(ar__6(30:32) .EQ. '') e_SII5454_3(i__) = rNULL__
        if(ar__6(36:39) .EQ. '') SII5640_3(i__) = rNULL__
        if(ar__6(41:43) .EQ. '') e_SII5640_5(i__) = rNULL__
        if(ar__6(47:50) .EQ. '') SiII6355_3(i__) = rNULL__
        if(ar__6(52:54) .EQ. '') e_SiII6355_5(i__) = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  SN_5(i__),f_SN_3(i__),Phase_5(i__),CaII3945_1(i__),
     +  e_CaII3945_1(i__),f_CaII3945_1(i__),SII5454_3(i__),
     +  e_SII5454_3(i__),f_SII5454_3(i__),SII5640_3(i__),
     +  e_SII5640_5(i__),f_SII5640_3(i__),SiII6355_3(i__),
     +  e_SiII6355_5(i__),f_SiII6355_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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