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