Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/135/1766 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 2026-May-15
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/135/1766 Sloan Digital Sky Survey-II supernova survey (Zheng+, 2008)
*================================================================================
*First-year spectroscopy for the Sloan Digital Sky Survey-II supernova survey.
* Zheng C., Romani R.W., Sako M., Marriner J., Bassett B., Becker A.,
* Choi C., Cinabro D., Dejongh F., Depoy D.L., Dilday B., Doi M.,
* Frieman J.A., Garnavich P.M., Hogan C.J., Holtzman J., Im M., Jha S.,
* Kessler R., Konishi K., Lampeitl H., Marshall J.L., McGinnis D.,
* Miknaitis G., Nichol R.C., Prieto J.L., Riess A.G., Richmond M.W.,
* Schneider D.P., Smith M., Takanashi N., Tokita K., Van Der Heyden K.,
* Yasuda N., Assef R.J., Barentine J., Bender R., Blandford R.D., Bremer M.,
* Brewington H., Collins C.A., Crotts A., Dembicky J., Eastman J., Edge A.,
* Elson E., Eyler M.E., Filippenko A.V., Foley R.J., Frank S., Goobar A.,
* Harvanek M., Hopp U., Ihara Y., Kahn S., Ketzeback W., Kleinman S.J.,
* Kollatschny W., Krzesinski J., Leloudas G., Long D.C., Lucey J.,
* Malanushenko E., Malanushenko V., Mcmillan R.J., Morgan C.W., Morokuma T.,
* Nitta A., Ostman L., Pan K., Romer A.K., Saurage G., Schlesinger K.,
* Snedden S.A., Sollerman J., Stritzinger M., Watson L.C., Watters S.,
* Wheeler J.C., York D.
* <Astron. J., 135, 1766-1784 (2008)>
* =2008AJ....135.1766Z
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! SDSS-II SN spectroscopic follow-up observations
integer*4 nr__
parameter (nr__=260) ! Number of records
character*58 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__) ! Official IAU name (YYYYAA), when assigned
character*1 n_SN (nr__) ! Note on SN (G1)
integer*4 SNID (nr__) ! SDSS Internal SN candidate ID,
* <SDSS-II SN NNNNN> in Simbad
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*8 Tel (nr__) ! Telescope
character*10 ObsDate (nr__) ! ("YYYY/MM/DD") Observation date (UT)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! SDSS-II SN redshift cross-correlation and type
analysis
integer*4 nr__1
parameter (nr__1=260) ! Number of records
character*69 ar__1 ! Full-size record
character*6 SN_1 (nr__1) ! Official IAU name (YYYYAA), when assigned
character*1 n_SN_1 (nr__1) ! Note on SN (G1)
integer*4 SNID_1 (nr__1) ! SDSS Internal SN candidate ID
integer*4 Ep_p (nr__1) ! (d) ? Spectrum day (in observer frame) relative to
* epoch of maximum brightness from the
* photometric typing
integer*4 Ep_s (nr__1) ! (d) ? Spectrum day (in observer frame) relative to
* epoch of maximum brightness from the
* spectroscopic cross-correlation analysis
character*7 SNtype (nr__1) ! Best-fitted type of the SN (2)
real*4 zSN (nr__1) ! ? SN redshift (3)
character*6 SNtemp (nr__1) ! Template spectrum (4)
real*4 rSN (nr__1) ! ? r value of the rvsao.xcsao comparison with
* SN templates (5)
real*4 zG (nr__1) ! ? Redshift from cross-correlation against
* galaxy templates (6)
character*6 GType (nr__1) ! Best-fit morphological type of the host galaxy
real*4 rG (nr__1) ! ? r value for rvsao.xcsao correlation against
* host-galaxy templates (7)
*Note (2): peculiar types:
* * "Ia?" indicates a best match to an SN Ia template, although some
* features are not detected.
* * "II?" indicates a general match to the SN II template, although
* expected SN II features are not all detected.
* * "Ia-pec" indicates peculiar SN Ia.
*Note (3): Redshift measured from cross-correlation against SN templates.
* This value has been corrected for the systematic offset 0.003.
* The expected redshift error is 0.005.
*Note (4): Template spectrum providing the best match to the observed spectrum.
* If blank, then the best-fitted template belongs to the Nugent SN Ia set.
* SN templates from Pun et al. (1995ApJS...99..223P),
* Turatto et al. (1996MNRAS.283....1T),
* Hamuy et al. (2002AJ....124..417H),
* Gal-Yam et al. (2002MNRAS.332L..73G),
* Leonard et al. (2002PASP..114...35L),
* Elmhamdi et al. (2004A&A...426..963E),
* and Nugent et al. (2002PASP..114..803N)
*Note (5): Ratio of peak height h to the root-mean square (rms).
* r_{SN}_>3 is a criterion for a "good" fit. See text for details.
*Note (6): The expected redshift error is 0.0005.
*Note (7): Ratio of peak height h to the root-mean square (rms).
* r_g_>3 is a criterion for a "good" fit. See text for details.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! SDSS-II SNe spectral decomposition and host
analysis results
integer*4 nr__2
parameter (nr__2=96) ! Number of records
character*69 ar__2 ! Full-size record
character*1 Flag (nr__2) ! [AB] Flag about extinction (1)
integer*4 SNID_2 (nr__2) ! SDSS Internal SN candidate ID
character*7 Tel_1 (nr__2) ! Spectroscopic follow-up telescope
integer*4 Ep_p_1 (nr__2) ! (d) Observation date of the SN relative to B-band
* maximum brightness estimated from light-curve fits
character*15 SNs (nr__2) ! The best and second-best SN type from
* cross-correlation analysis (2)
integer*4 Ep_sh (nr__2) ! (d) Age of the SN relative to B-band maximum
* brightness estimated from spectrum
* decomposition
character*7 Gal_sh (nr__2) ! Host-galaxy morphological type from spectral
* decomposition (3)
real*4 SNfrac (nr__2) ! Estimate of the SN light fraction in the
* observed spectrum
real*4 e_SNfrac (nr__2) ! Statistical error on SNFrac (the maximum value
* of the asymmetric errors obtained from the
* PCA code)
real*4 chi2 (nr__2) ! Reduced {chi}^2^ value
real*4 Q (nr__2) ! Quality factor (4)
*Note (1): Flag as follows:
* A = SNe have extinction estimates from MCLS2k2 analysis (Jha et al.
* 2007, Cat. J/ApJ/659/122)
* B = SNe have A_V from light-curve fitting using the frame difference
* photometry
*Note (2): "Ia?" indicates best fit to an SN Ia template, but a reasonable
* fit to the next-best type.
*Note (3): Unknown indicates too little host light for type ID.
*Note (4): Quality factor Q={Delta}({chi}^2^_{nu}_)sqrt({nu})/{chi}^2^_{nu}_,
* which estimates the significance of the {chi}^2^ differentiation between
* the best two SN types. Negative values indicate a selected type
* differing from that determined by the cross-correlation analysis.
C=============================================================================
C Loading file 'table1.dat' ! SDSS-II SN spectroscopic follow-up observations
C Format for file interpretation
1 format(
+ A6,1X,A1,1X,I5,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,
+ A8,1X,A10)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,260
read(1,'(A58)')ar__
read(ar__,1)
+ SN(i__),n_SN(i__),SNID(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Tel(i__),ObsDate(i__)
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__),n_SN(i__),SNID(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Tel(i__),ObsDate(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 'table3.dat' ! SDSS-II SN redshift cross-correlation and type
* analysis
C Format for file interpretation
2 format(
+ A6,1X,A1,1X,I5,1X,I3,1X,I3,1X,A7,1X,F6.4,2X,A6,1X,F4.1,1X,
+ F6.4,1X,A6,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,260
read(1,'(A69)')ar__1
read(ar__1,2)
+ SN_1(i__),n_SN_1(i__),SNID_1(i__),Ep_p(i__),Ep_s(i__),
+ SNtype(i__),zSN(i__),SNtemp(i__),rSN(i__),zG(i__),GType(i__),
+ rG(i__)
if(ar__1(16:18) .EQ. '') Ep_p(i__) = iNULL__
if(ar__1(20:22) .EQ. '') Ep_s(i__) = iNULL__
if(ar__1(32:37) .EQ. '') zSN(i__) = rNULL__
if(ar__1(47:50) .EQ. '') rSN(i__) = rNULL__
if(ar__1(52:57) .EQ. '') zG(i__) = rNULL__
if(ar__1(66:69) .EQ. '') rG(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ SN_1(i__),n_SN_1(i__),SNID_1(i__),Ep_p(i__),Ep_s(i__),
+ SNtype(i__),zSN(i__),SNtemp(i__),rSN(i__),zG(i__),GType(i__),
+ rG(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! SDSS-II SNe spectral decomposition and host
* analysis results
C Format for file interpretation
3 format(
+ A1,1X,I5,1X,A7,I3,1X,A15,I3,1X,A7,1X,F5.3,1X,F5.3,1X,F5.2,1X,
+ F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,96
read(1,'(A69)')ar__2
read(ar__2,3)
+ Flag(i__),SNID_2(i__),Tel_1(i__),Ep_p_1(i__),SNs(i__),
+ Ep_sh(i__),Gal_sh(i__),SNfrac(i__),e_SNfrac(i__),chi2(i__),
+ Q(i__)
c ..............Just test output...........
write(6,3)
+ Flag(i__),SNID_2(i__),Tel_1(i__),Ep_p_1(i__),SNs(i__),
+ Ep_sh(i__),Gal_sh(i__),SNfrac(i__),e_SNfrac(i__),chi2(i__),
+ Q(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end