Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/609/A134 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/A+A/609/A134 CSP-I photometry for 34 supernovae (Stritzinger+, 2018) *================================================================================ *The Carnegie Supernova Project I: *photometry data release of low-redshift stripped-envelope supernovae. * Stritzinger M.D., Anderson J.P., Contreras C., Heinrich-Josties E., * Morrell N., Phillips M.M., Anais J., Boldt L., Busta L., Burns C.R., * Campillay A., Corco C., Castellon S., Folatelli G., Gonzalez C., Holmbo S., * Hsiao E.Y., Krzeminski W., Salgado F., Seron J., Torres-Robledo S., * Freedman W.L., Hamuy M., Krisciunas K., Madore B.F., Persson S.E., Roth M., * Suntzeff N.B., Taddia F., Li W., Filippenko A.V. * <Astron. Astrophys. 609, A134 (2018)> * =2018A&A...609A.134S (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! General properties of 34 Type Ibc and Type IIb supernovae integer*4 nr__ parameter (nr__=34) ! Number of records character*149 ar__ ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*6 SN (nr__) ! SN name integer*4 RAh (nr__) ! (h) Right ascension (J2000) integer*4 RAm (nr__) ! (min) Right ascension (J2000) real*4 RAs (nr__) ! (s) Right ascension (J2000) character*1 DE_ (nr__) ! Declination sign (J2000) integer*4 DEd (nr__) ! (deg) Declination (J2000) integer*4 DEm (nr__) ! (arcmin) Declination (J2000) real*4 DEs (nr__) ! (arcsec) Declination (J2000) character*23 Host (nr__) ! Host galaxy name character*19 Morph (nr__) ! Morphology (1) real*4 z (nr__) ! Heliocentric redshift character*1 n_z (nr__) ! [y] Note on z (2) character*9 RDisc (nr__) ! Discovery reference character*18 Disc (nr__) ! Discovery group (3) character*6 Type (nr__) ! Spectral type integer*4 Phase1 (nr__) ! Phase range integer*4 Phase2 (nr__) ! Phase range real*8 T_B_max (nr__) ! (d) Time of maximum brightness real*4 e_T_B_max (nr__) ! (d) rms uncertainty on T(B)max character*1 n_T_B_max (nr__) ! [d] Note on T(B)max (4) integer*4 T5 (nr__) ! ? indicates optical photometry in du Pont * natural system (in table5) integer*4 T6 (nr__) ! ? indicates NIR photometry in Swope natural * system (in table6) integer*4 T7 (nr__) ! ? indicates NIR photometry in Swope natural * system (in table7) integer*4 T8 (nr__) ! ? indicates NIR photometry in du Pont natural * system (in table8) *Note (1): Morphologies are as listed in the NASA/IPAC Extragalactic * Database (NED). *Note (2): Heliocentric redshifts are as given in NED or if not listed in NED, * as determined through host galaxy emission lines in visual-wavelength * spectra of the SN. Measurements made by the authors are indicated with * the y symbol. *Note (3): References/URLs as: * LOSS = see Filippenko 2005, in The Fate of the Most Massive Stars, * ed. R. Humphreys & K. Stanek (San Francisco: ASP), p. 33, * and references therein * Puckett = http://www.cometwatch.com * SDSS-II = SDSS-II SN survey (Frieman et al., 2008AJ....135..338F) * SNFactory = http://snfactory.lbl.gov/ * CHASE = http://www.das.uchile.cl/proyectoCHASE *Note (4): d when obtained from the relation shown in Fig. 3 of * Taddia et al. (2018A&A...609A.134F). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Optical photometry of the local sequences in the standard system (for 34 SN) integer*4 nr__1 parameter (nr__1=612) ! Number of records character*119 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*6 SN_1 (nr__1) ! SN name integer*4 ID (nr__1) ! ID sequence number integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000) integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000) real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000) character*1 DE__1 (nr__1) ! Declination sign (J2000) integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000) integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000) real*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000) real*4 umag (nr__1) ! (mag) ?=- u magnitude real*4 e_umag (nr__1) ! (mag) ?=- rms uncertainty on umag real*4 gmag (nr__1) ! (mag) ?=- g magnitude real*4 e_gmag (nr__1) ! (mag) ?=- rms uncertainty on gmag real*4 rmag (nr__1) ! (mag) ?=- r magnitude real*4 e_rmag (nr__1) ! (mag) ?=- rms uncertainty on rmag real*4 imag (nr__1) ! (mag) ?=- i magnitude real*4 e_imag (nr__1) ! (mag) ?=- rms uncertainty on imag real*4 Bmag (nr__1) ! (mag) ?=- B magnitude real*4 e_Bmag (nr__1) ! (mag) ?=- rms uncertainty on Bmag real*4 Vmag (nr__1) ! (mag) ?=- V magnitude real*4 e_Vmag (nr__1) ! (mag) ?=- rms uncertainty on Vmag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! NIR photometry of the local sequences in the standard system (for 25 SN) integer*4 nr__2 parameter (nr__2=548) ! Number of records character*78 ar__2 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension J2000 real*8 DEdeg_2 (nr__2) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*6 SN_2 (nr__2) ! SN name integer*4 ID_1 (nr__2) ! ID sequence number integer*4 RAh_2 (nr__2) ! (h) Right ascension (J2000) integer*4 RAm_2 (nr__2) ! (min) Right ascension (J2000) real*4 RAs_2 (nr__2) ! (s) Right ascension (J2000) character*1 DE__2 (nr__2) ! Declination sign (J2000) integer*4 DEd_2 (nr__2) ! (deg) Declination (J2000) integer*4 DEm_2 (nr__2) ! (arcmin) Declination (J2000) real*4 DEs_2 (nr__2) ! (arcsec) Declination (J2000) real*4 Ymag (nr__2) ! (mag) ?=- Y magnitude real*4 e_Ymag (nr__2) ! (mag) ?=- rms uncertainty on Ymag real*4 Jmag (nr__2) ! (mag) ?=- J magnitude real*4 e_Jmag (nr__2) ! (mag) ?=- rms uncertainty on Jmag real*4 Hmag (nr__2) ! (mag) ?=- H magnitude real*4 e_Hmag (nr__2) ! (mag) ?=- rms uncertainty on Hmag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Optical photometry of SE SNe in the Swope natural system (for 34 SN) integer*4 nr__3 parameter (nr__3=597) ! Number of records character*101 ar__3 ! Full-size record character*6 SN_3 (nr__3) ! SN name real*8 JD (nr__3) ! (d) Julian date real*4 umag_1 (nr__3) ! (mag) ?=- u magnitude real*4 e_umag_1 (nr__3) ! (mag) ?=- rms uncertainty on umag real*4 gmag_1 (nr__3) ! (mag) ?=- g magnitude real*4 e_gmag_1 (nr__3) ! (mag) ?=- rms uncertainty on gmag real*4 rmag_1 (nr__3) ! (mag) ?=- r magnitude real*4 e_rmag_1 (nr__3) ! (mag) ?=- rms uncertainty on rmag real*4 imag_1 (nr__3) ! (mag) ?=- i magnitude real*4 e_imag_1 (nr__3) ! (mag) ?=- rms uncertainty on imag real*4 Bmag_1 (nr__3) ! (mag) ?=- B magnitude real*4 e_Bmag_1 (nr__3) ! (mag) ?=- rms uncertainty on Bmag real*4 Vmag_1 (nr__3) ! (mag) ?=- V magnitude real*4 e_Vmag_1 (nr__3) ! (mag) ?=- rms uncertainty on Vmag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Optical photometry of SE SNe in the du Pont natural system (for 5 SN) integer*4 nr__4 parameter (nr__4=10) ! Number of records character*101 ar__4 ! Full-size record character*6 SN_4 (nr__4) ! SN name real*8 JD_1 (nr__4) ! (d) Julian date real*4 umag_2 (nr__4) ! (mag) ?=- u magnitude real*4 e_umag_2 (nr__4) ! (mag) ?=- rms uncertainty on umag real*4 gmag_2 (nr__4) ! (mag) ?=- g magnitude real*4 e_gmag_2 (nr__4) ! (mag) ?=- rms uncertainty on gmag real*4 rmag_2 (nr__4) ! (mag) ?=- r magnitude real*4 e_rmag_2 (nr__4) ! (mag) ?=- rms uncertainty on rmag real*4 imag_2 (nr__4) ! (mag) ?=- i magnitude real*4 e_imag_2 (nr__4) ! (mag) ?=- rms uncertainty on imag real*4 Bmag_2 (nr__4) ! (mag) ?=- B magnitude real*4 e_Bmag_2 (nr__4) ! (mag) ?=- rms uncertainty on Bmag real*4 Vmag_2 (nr__4) ! (mag) ?=- V magnitude real*4 e_Vmag_2 (nr__4) ! (mag) ?=- rms uncertainty on Vmag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! NIR photometry of SE SNe in the Swope natural system (for 13 SN) integer*4 nr__5 parameter (nr__5=178) ! Number of records character*59 ar__5 ! Full-size record character*6 SN_5 (nr__5) ! SN name real*8 JD_2 (nr__5) ! (d) Julian date real*4 Ymag_1 (nr__5) ! (mag) ?=- Y magnitude real*4 e_Ymag_1 (nr__5) ! (mag) ?=- rms uncertainty on Ymag real*4 Jmag_1 (nr__5) ! (mag) ?=- J(RC1) magnitude real*4 e_Jmag_1 (nr__5) ! (mag) ?=- rms uncertainty on Jmag real*4 Hmag_1 (nr__5) ! (mag) ?=- H magnitude real*4 e_Hmag_1 (nr__5) ! (mag) ?=- rms uncertainty on Hmag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! NIR photometry of SE SNe in the Swope natural system (for 5 SN) integer*4 nr__6 parameter (nr__6=68) ! Number of records character*59 ar__6 ! Full-size record character*6 SN_6 (nr__6) ! SN name real*8 JD_3 (nr__6) ! (d) Julian date real*4 Ymag_2 (nr__6) ! (mag) ?=- Y magnitude real*4 e_Ymag_2 (nr__6) ! (mag) ?=- rms uncertainty on Ymag real*4 Jmag_2 (nr__6) ! (mag) ?=- J(RC2) magnitude real*4 e_Jmag_2 (nr__6) ! (mag) ?=- rms uncertainty on Jmag real*4 Hmag_2 (nr__6) ! (mag) ?=- H magnitude real*4 e_Hmag_2 (nr__6) ! (mag) ?=- rms uncertainty on Hmag c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! NIR photometry of SE SNe in the du Pont natural system (for 23 SN) integer*4 nr__7 parameter (nr__7=164) ! Number of records character*59 ar__7 ! Full-size record character*6 SN_7 (nr__7) ! SN name real*8 JD_4 (nr__7) ! (d) Julian date real*4 Ymag_3 (nr__7) ! (mag) ?=- Y(dw) magnitude real*4 e_Ymag_3 (nr__7) ! (mag) ?=- rms uncertainty on Ymag real*4 Jmag_3 (nr__7) ! (mag) ?=- J magnitude real*4 e_Jmag_3 (nr__7) ! (mag) ?=- rms uncertainty on Jmag real*4 Hmag_3 (nr__7) ! (mag) ?=- H magnitude real*4 e_Hmag_3 (nr__7) ! (mag) ?=- rms uncertainty on Hmag C============================================================================= C Loading file 'table1.dat' ! General properties of 34 Type Ibc and Type IIb * supernovae C Format for file interpretation 1 format( + A6,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A23,1X,A19, + 1X,F6.4,A1,1X,A9,1X,A18,1X,A6,1X,I3,1X,I4,1X,F7.2,1X,F4.2,A1, + 1X,I1,1X,I1,1X,I1,1X,I1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,34 read(1,'(A149)')ar__ read(ar__,1) + SN(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__), + DEs(i__),Host(i__),Morph(i__),z(i__),n_z(i__),RDisc(i__), + Disc(i__),Type(i__),Phase1(i__),Phase2(i__),T_B_max(i__), + e_T_B_max(i__),n_T_B_max(i__),T5(i__),T6(i__),T7(i__),T8(i__) if(ar__(143:143) .EQ. '') T5(i__) = iNULL__ if(ar__(145:145) .EQ. '') T6(i__) = iNULL__ if(ar__(147:147) .EQ. '') T7(i__) = iNULL__ if(ar__(149:149) .EQ. '') T8(i__) = iNULL__ RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,1) + SN(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__), + DEs(i__),Host(i__),Morph(i__),z(i__),n_z(i__),RDisc(i__), + Disc(i__),Type(i__),Phase1(i__),Phase2(i__),T_B_max(i__), + e_T_B_max(i__),n_T_B_max(i__),T5(i__),T6(i__),T7(i__),T8(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Optical photometry of the local sequences in * the standard system (for 34 SN) C Format for file interpretation 2 format( + A6,1X,I2,2X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X,F6.3, + 1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X, + F6.3,1X,F5.3,2X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,612 read(1,'(A119)')ar__1 read(ar__1,2) + SN_1(i__),ID(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__), + DEd_1(i__),DEm_1(i__),DEs_1(i__),umag(i__),e_umag(i__), + gmag(i__),e_gmag(i__),rmag(i__),e_rmag(i__),imag(i__), + e_imag(i__),Bmag(i__),e_Bmag(i__),Vmag(i__),e_Vmag(i__) if (idig(ar__1(38:43)).EQ.0) umag(i__) = rNULL__ if (idig(ar__1(45:49)).EQ.0) e_umag(i__) = rNULL__ if (idig(ar__1(52:57)).EQ.0) gmag(i__) = rNULL__ if (idig(ar__1(59:63)).EQ.0) e_gmag(i__) = rNULL__ if (idig(ar__1(66:71)).EQ.0) rmag(i__) = rNULL__ if (idig(ar__1(73:77)).EQ.0) e_rmag(i__) = rNULL__ if (idig(ar__1(80:85)).EQ.0) imag(i__) = rNULL__ if (idig(ar__1(87:91)).EQ.0) e_imag(i__) = rNULL__ if (idig(ar__1(94:99)).EQ.0) Bmag(i__) = rNULL__ if (idig(ar__1(101:105)).EQ.0) e_Bmag(i__) = rNULL__ if (idig(ar__1(108:113)).EQ.0) Vmag(i__) = rNULL__ if (idig(ar__1(115:119)).EQ.0) e_Vmag(i__) = rNULL__ RAdeg_1(i__) = rNULL__ DEdeg_1(i__) = rNULL__ c Derive coordinates RAdeg_1 and DEdeg_1 from input data c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown) if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15. if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4. if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240. if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__) if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60. if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600. if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__) c ..............Just test output........... write(6,2) + SN_1(i__),ID(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__), + DEd_1(i__),DEm_1(i__),DEs_1(i__),umag(i__),e_umag(i__), + gmag(i__),e_gmag(i__),rmag(i__),e_rmag(i__),imag(i__), + e_imag(i__),Bmag(i__),e_Bmag(i__),Vmag(i__),e_Vmag(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! NIR photometry of the local sequences in the * standard system (for 25 SN) C Format for file interpretation 3 format( + A6,1X,I3,2X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X,F6.3, + 1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,548 read(1,'(A78)')ar__2 read(ar__2,3) + SN_2(i__),ID_1(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__), + DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),Ymag(i__), + e_Ymag(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__) if (idig(ar__2(39:44)).EQ.0) Ymag(i__) = rNULL__ if (idig(ar__2(46:50)).EQ.0) e_Ymag(i__) = rNULL__ if (idig(ar__2(53:58)).EQ.0) Jmag(i__) = rNULL__ if (idig(ar__2(60:64)).EQ.0) e_Jmag(i__) = rNULL__ if (idig(ar__2(67:72)).EQ.0) Hmag(i__) = rNULL__ if (idig(ar__2(74:78)).EQ.0) e_Hmag(i__) = rNULL__ RAdeg_2(i__) = rNULL__ DEdeg_2(i__) = rNULL__ c Derive coordinates RAdeg_2 and DEdeg_2 from input data c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown) if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15. if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4. if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240. if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__) if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60. if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600. if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__) c ..............Just test output........... write(6,3) + SN_2(i__),ID_1(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__), + DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),Ymag(i__), + e_Ymag(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Optical photometry of SE SNe in the * Swope natural system (for 34 SN) C Format for file interpretation 4 format( + A6,1X,F10.2,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3, + 2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,597 read(1,'(A101)')ar__3 read(ar__3,4) + SN_3(i__),JD(i__),umag_1(i__),e_umag_1(i__),gmag_1(i__), + e_gmag_1(i__),rmag_1(i__),e_rmag_1(i__),imag_1(i__), + e_imag_1(i__),Bmag_1(i__),e_Bmag_1(i__),Vmag_1(i__), + e_Vmag_1(i__) if (idig(ar__3(20:25)).EQ.0) umag_1(i__) = rNULL__ if (idig(ar__3(27:31)).EQ.0) e_umag_1(i__) = rNULL__ if (idig(ar__3(34:39)).EQ.0) gmag_1(i__) = rNULL__ if (idig(ar__3(41:45)).EQ.0) e_gmag_1(i__) = rNULL__ if (idig(ar__3(48:53)).EQ.0) rmag_1(i__) = rNULL__ if (idig(ar__3(55:59)).EQ.0) e_rmag_1(i__) = rNULL__ if (idig(ar__3(62:67)).EQ.0) imag_1(i__) = rNULL__ if (idig(ar__3(69:73)).EQ.0) e_imag_1(i__) = rNULL__ if (idig(ar__3(76:81)).EQ.0) Bmag_1(i__) = rNULL__ if (idig(ar__3(83:87)).EQ.0) e_Bmag_1(i__) = rNULL__ if (idig(ar__3(90:95)).EQ.0) Vmag_1(i__) = rNULL__ if (idig(ar__3(97:101)).EQ.0) e_Vmag_1(i__) = rNULL__ c ..............Just test output........... write(6,4) + SN_3(i__),JD(i__),umag_1(i__),e_umag_1(i__),gmag_1(i__), + e_gmag_1(i__),rmag_1(i__),e_rmag_1(i__),imag_1(i__), + e_imag_1(i__),Bmag_1(i__),e_Bmag_1(i__),Vmag_1(i__), + e_Vmag_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Optical photometry of SE SNe in the * du Pont natural system (for 5 SN) C Format for file interpretation 5 format( + A6,1X,F10.2,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3, + 2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,10 read(1,'(A101)')ar__4 read(ar__4,5) + SN_4(i__),JD_1(i__),umag_2(i__),e_umag_2(i__),gmag_2(i__), + e_gmag_2(i__),rmag_2(i__),e_rmag_2(i__),imag_2(i__), + e_imag_2(i__),Bmag_2(i__),e_Bmag_2(i__),Vmag_2(i__), + e_Vmag_2(i__) if (idig(ar__4(20:25)).EQ.0) umag_2(i__) = rNULL__ if (idig(ar__4(27:31)).EQ.0) e_umag_2(i__) = rNULL__ if (idig(ar__4(34:39)).EQ.0) gmag_2(i__) = rNULL__ if (idig(ar__4(41:45)).EQ.0) e_gmag_2(i__) = rNULL__ if (idig(ar__4(48:53)).EQ.0) rmag_2(i__) = rNULL__ if (idig(ar__4(55:59)).EQ.0) e_rmag_2(i__) = rNULL__ if (idig(ar__4(62:67)).EQ.0) imag_2(i__) = rNULL__ if (idig(ar__4(69:73)).EQ.0) e_imag_2(i__) = rNULL__ if (idig(ar__4(76:81)).EQ.0) Bmag_2(i__) = rNULL__ if (idig(ar__4(83:87)).EQ.0) e_Bmag_2(i__) = rNULL__ if (idig(ar__4(90:95)).EQ.0) Vmag_2(i__) = rNULL__ if (idig(ar__4(97:101)).EQ.0) e_Vmag_2(i__) = rNULL__ c ..............Just test output........... write(6,5) + SN_4(i__),JD_1(i__),umag_2(i__),e_umag_2(i__),gmag_2(i__), + e_gmag_2(i__),rmag_2(i__),e_rmag_2(i__),imag_2(i__), + e_imag_2(i__),Bmag_2(i__),e_Bmag_2(i__),Vmag_2(i__), + e_Vmag_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! NIR photometry of SE SNe in the * Swope natural system (for 13 SN) C Format for file interpretation 6 format( + A6,1X,F10.2,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,178 read(1,'(A59)')ar__5 read(ar__5,6) + SN_5(i__),JD_2(i__),Ymag_1(i__),e_Ymag_1(i__),Jmag_1(i__), + e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__) if (idig(ar__5(20:25)).EQ.0) Ymag_1(i__) = rNULL__ if (idig(ar__5(27:31)).EQ.0) e_Ymag_1(i__) = rNULL__ if (idig(ar__5(34:39)).EQ.0) Jmag_1(i__) = rNULL__ if (idig(ar__5(41:45)).EQ.0) e_Jmag_1(i__) = rNULL__ if (idig(ar__5(48:53)).EQ.0) Hmag_1(i__) = rNULL__ if (idig(ar__5(55:59)).EQ.0) e_Hmag_1(i__) = rNULL__ c ..............Just test output........... write(6,6) + SN_5(i__),JD_2(i__),Ymag_1(i__),e_Ymag_1(i__),Jmag_1(i__), + e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! NIR photometry of SE SNe in the * Swope natural system (for 5 SN) C Format for file interpretation 7 format( + A6,1X,F10.2,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,68 read(1,'(A59)')ar__6 read(ar__6,7) + SN_6(i__),JD_3(i__),Ymag_2(i__),e_Ymag_2(i__),Jmag_2(i__), + e_Jmag_2(i__),Hmag_2(i__),e_Hmag_2(i__) if (idig(ar__6(20:25)).EQ.0) Ymag_2(i__) = rNULL__ if (idig(ar__6(27:31)).EQ.0) e_Ymag_2(i__) = rNULL__ if (idig(ar__6(34:39)).EQ.0) Jmag_2(i__) = rNULL__ if (idig(ar__6(41:45)).EQ.0) e_Jmag_2(i__) = rNULL__ if (idig(ar__6(48:53)).EQ.0) Hmag_2(i__) = rNULL__ if (idig(ar__6(55:59)).EQ.0) e_Hmag_2(i__) = rNULL__ c ..............Just test output........... write(6,7) + SN_6(i__),JD_3(i__),Ymag_2(i__),e_Ymag_2(i__),Jmag_2(i__), + e_Jmag_2(i__),Hmag_2(i__),e_Hmag_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! NIR photometry of SE SNe in the * du Pont natural system (for 23 SN) C Format for file interpretation 8 format( + A6,1X,F10.2,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3,2X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,164 read(1,'(A59)')ar__7 read(ar__7,8) + SN_7(i__),JD_4(i__),Ymag_3(i__),e_Ymag_3(i__),Jmag_3(i__), + e_Jmag_3(i__),Hmag_3(i__),e_Hmag_3(i__) if (idig(ar__7(20:25)).EQ.0) Ymag_3(i__) = rNULL__ if (idig(ar__7(27:31)).EQ.0) e_Ymag_3(i__) = rNULL__ if (idig(ar__7(34:39)).EQ.0) Jmag_3(i__) = rNULL__ if (idig(ar__7(41:45)).EQ.0) e_Jmag_3(i__) = rNULL__ if (idig(ar__7(48:53)).EQ.0) Hmag_3(i__) = rNULL__ if (idig(ar__7(55:59)).EQ.0) e_Hmag_3(i__) = rNULL__ c ..............Just test output........... write(6,8) + SN_7(i__),JD_4(i__),Ymag_3(i__),e_Ymag_3(i__),Jmag_3(i__), + e_Jmag_3(i__),Hmag_3(i__),e_Hmag_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end C Locate position of first digit in string; or return 0 integer function idig(c) character*(*) c character*1 c1 integer lc,i lc=len(c) idig=0 do i=1,lc if(c(i:i).ne.' ') go to 1 end do 1 if(i.gt.lc) return c1=c(i:i) if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1 if(i.gt.lc) return c1=c(i:i) if(c1.ge.'0'.and.c1.le.'9') idig=i return end