Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/956/108 into FORTRAN code for reading data files line by line.
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-08
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/ApJ/956/108 TESS 4yrs obs. of Type Ia supernovae (Fausnaugh+, 2023)
*================================================================================
*Four years of type Ia supernovae observed by TESS: early-time light-curve shapes
*and constraints on companion interaction models.
* Fausnaugh M.M., Vallely P.J., Tucker M.A., Kochanek C.S., Shappee B.J.,
* Stanek K.Z., Ricker G.R., Vanderspek R., Agarwal M., Daylan T.,
* Jayaraman R., Hounsell R., Muthukrishna D.
* <Astrophys. J., 956, 108 (2023)>
* =2023ApJ...956..108F
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Properties of the supernova sample
integer*4 nr__
parameter (nr__=307) ! Number of records
character*104 ar__ ! Full-size record
character*10 Name ! SN identifier
character*13 OName ! Other identifier
real*4 TmagP ! (mag) [14/21.3] Peak TESS band magnitude
* (T_peak_)
real*4 e_TmagP ! (mag) [0.05/0.5] Uncertainty in TmagP (1)
real*4 TmagL ! (mag) [17.9/22.4] Limiting TESS band magnitude
* (T_lim_)
real*4 delTmag ! (mag) [-0.94/5.4] Magnitude increase in TESS band
* ({Delta}T)
real*8 tPeak ! (d) [1352.4/2691.6] BJD of peak TESS band
* magnitude; tPeak-2457000.0
real*4 TMag ! (mag) [-21.51/-15.21] Absolute TESS band magnitude
* (2)
real*4 e_TMag ! (mag) [0.05/0.9] Uncertainty in TMag
real*4 nuLum ! (10+35W) [0.13/48.2] Luminosity; 1e+42erg/s
real*4 e_nuLum ! (10+35W) [0.01/9.6] Uncertainty in nuLum
real*4 z ! [0.006/0.2] Spectroscopic redshift of host
* galaxy
real*4 E_B_V ! (mag) [0.004/0.4] THe (B-V) color excess
integer*4 Reduc ! [0/110] Data reduction bitmask (3)
integer*4 Treduc ! [0/110] TESSreduce bitmask (4)
integer*4 Qual ! [0/1110] Quality bitmask (5)
*Note (1): Uncertainty includes the TESS instrument flux calibration (0.05mag,
* Appendix A), uncertainty in the fit to the baseline flux
* (Section 3.1), and any uncertainty in the calibration offset of
* the second sector, if applicable (Section 3.1).
*Note (2): Absolute TESS magnitude. Includes uncertainty in Tpeak (column 2)
* and an assumed uncertainty of 1000km/s in the redshift if the
* host galaxy is not reported to TNS.
*Note (3): Marks SN where we applied the "background model" correction (first
* bit, "100"), where we detrended variable star signals (second bit,
* "010"), or where the light curve is affected by asteroids (third
* bit, "001"). See Section A.1 for details.
*Note (4): Marks SN where we extracted a TESSreduce light curve (first bit,
* "100"), where we elected to use the TESSreduce light curve rather
* than the ISIS light curve (second bit, "010"), or where the two
* light curves show different shapes/amplitudes (third bit, "001").
* See Section A.2 for details.
*Note (5): Marks SN where TESS did not observe the peak flux (first bit,
* "1000", see Appendix C), the SN is part of the "No First Light"
* subsample (second bit, "0100", see Appendix C), the SN is part of
* the "High Quality" subsample (third bit, "0010", see Appendix C.1),
* or the SN is part of the "High Dynamic Range" subsample (fourth bit,
* "0001", see Appendix C.2).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! SN light curves
integer*4 nr__1
parameter (nr__1=680384) ! Number of records
character*81 ar__1 ! Full-size record
character*10 Name_1 ! SN identifier
real*8 BJD ! (d) [1325.3/2691.6] Barycentric JD - TDB;
* BJD-2457000
real*8 CRate ! (ct/s) [-84.2/974.8] Count rate
real*8 e_CRate ! (ct/s) [0.16/14.5] Uncertainty in CRate
real*8 Frac ! [-16.2/16] Fraction of Peak
real*4 e_Frac ! [0.007/4] Uncertainty in Frac
real*8 Tmag_1 ! (mag) [14/22.4]? Apparent TESS band magnitude (1)
real*8 e_Tmag_1 ! (mag) [0.008/28.5]? Uncertainty in Tmag (1)
real*8 CalOff ! [-0.7/363]? Calibration offset (2)
*Note (1): For fluxes below the 8 hour detection limits (including zero and
* negative flux due to noise), the magnitudes are represented by TmagL
* from Table 1 and the uncertainties are blank.
*Note (2): Gives the shift applied to flux calibrate the second sector
* of TESS observations for supernova observed near the ecliptic
* poles (see Appendix B).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Curved Power Law posterior medians
integer*4 nr__2
parameter (nr__2=307) ! Number of records
character*91 ar__2 ! Full-size record
character*10 Name_2 ! SN identifier
real*8 BJD_1 ! (d) [1330.74/2688.48] Time of first light;
* BJD-2470000.0
real*4 e_BJD ! (d) [0.01/3.6] Uncertainty in BJD
real*4 tRise ! (d) [3.49/28.81]? Rise time (1)
real*4 PLbeta1 ! [0.02/3.75] Mean power law index, t^{beta}1^
real*4 e_PLbeta1 ! [0.01/0.92] Uncertainty in PLbeta1
real*8 PLbeta2 ! [-0.095/-0.0001] Mean power law index,
* t^{beta}2^
real*4 e_PLbeta2 ! [0.0001/0.03] Uncertainty in PLbeta2
real*8 lnZ ! ([-]) [-7151.2/11349.7] Natural log of Bayesian
* Evidence Z from dynesty
real*4 e_lnZ ! ([-]) [0.04/0.2] Uncertainty in lnZ
real*8 BIC ! [-22731.4/14303.7] Bayesian Information
* Criterion (2)
real*4 Ratio ! [0.06/0.98] Noise metric,
* rms_Gauss_/rms_observed_ (3)
integer*4 Mask ! [0/1110] BitMask (4)
*Note (1): The uncertainty on the rise time is determined by the uncertainty
* on the time of first light (column 3).
*Note (2): BIC=klnN-lnL
*
* where k is the number of parameters, N is the number of data points in
* the fit, and L is the likelihood. See Appendix B.
*Note (3): Defined in Section 4, which quantifies departures of the light curves
* from random Gaussian noise. Values near zero indicate very little
* improvement from binning (systematic errors dominate over random
* noise), while 1 indicates perfect Gaussian noise scaling (random
* noise only).
*Note (4): The same Quality Bitmask from Table 1, which marks the subsamples of
* SNe: TESS did not observe the peak flux of this SN (first bit,
* "1000", see Appendix C), the SN is part of the "No First Light"
* subsample (second bit, "0100", see Appendix C), the SN is part of
* the "High Quality" subsample (third bit, "0010", see Appendix C.1),
* or the SN is part of the "High Dynamic Range" subsample (fourth bit,
* "0001", see Appendix C.2).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Curved Power Law with companion interaction
posterior medians
integer*4 nr__3
parameter (nr__3=307) ! Number of records
character*120 ar__3 ! Full-size record
character*10 Name_3 ! SN identifier
real*8 BJD_2 ! (d) [1331.15/2689.25] Time of first light;
* BJD-2470000.0
real*4 e_BJD_1 ! (d) [0.01/3.3] Uncertainty in BJD
real*4 tRise_1 ! (d) [2.8/28.01]? Rise time (1)
real*4 PLbeta1_1 ! [0.01/3.9] Mean power law index,
* t^{beta}_1_^
real*4 e_PLbeta1_1 ! [0.01/0.95] Uncertainty in PLbeta1
real*8 PLbeta2_1 ! [-0.094/-0.0001] Mean power law index,
* t^{beta}_2_^
real*4 e_PLbeta2_1 ! [0.0001/0.03] Uncertainty in PLbeta2
real*4 MedRad ! (Rsun) [1/223.2] Median (50th percentile)
* companion Roche lobe radius
real*4 e_MedRad ! (Rsun) [0/55.5] Uncertainty in MedRad
real*4 MaxRad ! (Rsun) [1/223.2] The 3{sigma} (99.9th
* percentile) companion Roche lobe radius
real*4 MedSep ! (10+12cm) [0.2/40.2] Median (50th percentile)
* companion Roche lobe radius
real*4 e_MedSep ! (10+12cm) [0/9.99] Uncertainty in MedSep
real*4 MaxSep ! (10+12cm) [0.2/40.2] The 3{sigma} (99.9th percentile)
* companion Roche lobe radius
real*8 lnZ_1 ! ([-]) [-7151.4/11354.3] Natural log of Bayesian
* Evidence Z from dynesty
real*4 e_lnZ_1 ! ([-]) [0.04/0.2] Uncertainty in lnZ
real*8 BIC_1 ! [-22733.8/14311.4] Bayesian Information
* Criterion
integer*4 Mask_1 ! [0/1110] BitMask (2)
*Note (1): The uncertainty on the rise time is determined by the uncertainty
* on the time of first light (column 3).
*Note (2): The same Quality Bitmask from Table 1, which marks the subsamples of
* SNe: TESS did not observe the peak flux of this SN (first bit,
* "1000", see Appendix C), the SN is part of the "No First Light"
* subsample (second bit, "0100", see Appendix C), the SN is part of
* the "High Quality" subsample (third bit, "0010", see Appendix C.1),
* or the SN is part of the "High Dynamic Range" subsample (fourth bit,
* "0001", see Appendix C.2).
C=============================================================================
C Loading file 'table1.dat' ! Properties of the supernova sample
C Format for file interpretation
1 format(
+ A10,1X,A13,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F7.2,1X,F6.2,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,F6.4,1X,F5.3,1X,I3,1X,I3,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,307
read(1,'(A104)')ar__
read(ar__,1)
+ Name,OName,TmagP,e_TmagP,TmagL,delTmag,tPeak,TMag,e_TMag,
+ nuLum,e_nuLum,z,E_B_V,Reduc,Treduc,Qual
c ..............Just test output...........
write(6,1)
+ Name,OName,TmagP,e_TmagP,TmagL,delTmag,tPeak,TMag,e_TMag,
+ nuLum,e_nuLum,z,E_B_V,Reduc,Treduc,Qual
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! SN light curves
C Format for file interpretation
2 format(
+ A10,1X,F10.5,1X,F8.4,1X,F7.4,1X,F8.4,1X,F6.4,1X,F7.4,1X,F7.4,
+ 1X,F10.6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,680384
read(1,'(A81)')ar__1
read(ar__1,2)
+ Name_1,BJD,CRate,e_CRate,Frac,e_Frac,Tmag_1,e_Tmag_1,CalOff
if(ar__1(56:62) .EQ. '') Tmag_1 = rNULL__
if(ar__1(64:70) .EQ. '') e_Tmag_1 = rNULL__
if(ar__1(72:81) .EQ. '') CalOff = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,BJD,CRate,e_CRate,Frac,e_Frac,Tmag_1,e_Tmag_1,CalOff
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Curved Power Law posterior medians
C Format for file interpretation
3 format(
+ A10,1X,F7.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F7.4,1X,F6.4,
+ 1X,F8.2,1X,F4.2,1X,F9.2,1X,F6.4,1X,I5)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,307
read(1,'(A91)')ar__2
read(ar__2,3)
+ Name_2,BJD_1,e_BJD,tRise,PLbeta1,e_PLbeta1,PLbeta2,e_PLbeta2,
+ lnZ,e_lnZ,BIC,Ratio,Mask
if(ar__2(25:29) .EQ. '') tRise = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,BJD_1,e_BJD,tRise,PLbeta1,e_PLbeta1,PLbeta2,e_PLbeta2,
+ lnZ,e_lnZ,BIC,Ratio,Mask
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Curved Power Law with companion interaction
* posterior medians
C Format for file interpretation
4 format(
+ A10,1X,F7.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F7.4,1X,F6.4,
+ 1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F8.2,1X,
+ F4.2,1X,F9.2,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,307
read(1,'(A120)')ar__3
read(ar__3,4)
+ Name_3,BJD_2,e_BJD_1,tRise_1,PLbeta1_1,e_PLbeta1_1,PLbeta2_1,
+ e_PLbeta2_1,MedRad,e_MedRad,MaxRad,MedSep,e_MedSep,MaxSep,
+ lnZ_1,e_lnZ_1,BIC_1,Mask_1
if(ar__3(25:29) .EQ. '') tRise_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3,BJD_2,e_BJD_1,tRise_1,PLbeta1_1,e_PLbeta1_1,PLbeta2_1,
+ e_PLbeta2_1,MedRad,e_MedRad,MaxRad,MedSep,e_MedSep,MaxSep,
+ lnZ_1,e_lnZ_1,BIC_1,Mask_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end