Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/956/108 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-12
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 (nr__) ! SN identifier
character*13 OName (nr__) ! Other identifier
real*4 TmagP (nr__) ! (mag) [14/21.3] Peak TESS band magnitude
* (T_peak_)
real*4 e_TmagP (nr__) ! (mag) [0.05/0.5] Uncertainty in TmagP (1)
real*4 TmagL (nr__) ! (mag) [17.9/22.4] Limiting TESS band magnitude
* (T_lim_)
real*4 delTmag (nr__) ! (mag) [-0.94/5.4] Magnitude increase in TESS band
* ({Delta}T)
real*8 tPeak (nr__) ! (d) [1352.4/2691.6] BJD of peak TESS band
* magnitude; tPeak-2457000.0
real*4 TMag (nr__) ! (mag) [-21.51/-15.21] Absolute TESS band magnitude
* (2)
real*4 e_TMag (nr__) ! (mag) [0.05/0.9] Uncertainty in TMag
real*4 nuLum (nr__) ! (10+35W) [0.13/48.2] Luminosity; 1e+42erg/s
real*4 e_nuLum (nr__) ! (10+35W) [0.01/9.6] Uncertainty in nuLum
real*4 z (nr__) ! [0.006/0.2] Spectroscopic redshift of host
* galaxy
real*4 E_B_V (nr__) ! (mag) [0.004/0.4] THe (B-V) color excess
integer*4 Reduc (nr__) ! [0/110] Data reduction bitmask (3)
integer*4 Treduc (nr__) ! [0/110] TESSreduce bitmask (4)
integer*4 Qual (nr__) ! [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 (nr__1) ! SN identifier
real*8 BJD (nr__1) ! (d) [1325.3/2691.6] Barycentric JD - TDB;
* BJD-2457000
real*8 CRate (nr__1) ! (ct/s) [-84.2/974.8] Count rate
real*8 e_CRate (nr__1) ! (ct/s) [0.16/14.5] Uncertainty in CRate
real*8 Frac (nr__1) ! [-16.2/16] Fraction of Peak
real*4 e_Frac (nr__1) ! [0.007/4] Uncertainty in Frac
real*8 Tmag_1 (nr__1) ! (mag) [14/22.4]? Apparent TESS band magnitude (1)
real*8 e_Tmag_1 (nr__1) ! (mag) [0.008/28.5]? Uncertainty in Tmag (1)
real*8 CalOff (nr__1) ! [-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 (nr__2) ! SN identifier
real*8 BJD_1 (nr__2) ! (d) [1330.74/2688.48] Time of first light;
* BJD-2470000.0
real*4 e_BJD (nr__2) ! (d) [0.01/3.6] Uncertainty in BJD
real*4 tRise (nr__2) ! (d) [3.49/28.81]? Rise time (1)
real*4 PLbeta1 (nr__2) ! [0.02/3.75] Mean power law index, t^{beta}1^
real*4 e_PLbeta1 (nr__2) ! [0.01/0.92] Uncertainty in PLbeta1
real*8 PLbeta2 (nr__2) ! [-0.095/-0.0001] Mean power law index,
* t^{beta}2^
real*4 e_PLbeta2 (nr__2) ! [0.0001/0.03] Uncertainty in PLbeta2
real*8 lnZ (nr__2) ! ([-]) [-7151.2/11349.7] Natural log of Bayesian
* Evidence Z from dynesty
real*4 e_lnZ (nr__2) ! ([-]) [0.04/0.2] Uncertainty in lnZ
real*8 BIC (nr__2) ! [-22731.4/14303.7] Bayesian Information
* Criterion (2)
real*4 Ratio (nr__2) ! [0.06/0.98] Noise metric,
* rms_Gauss_/rms_observed_ (3)
integer*4 Mask (nr__2) ! [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 (nr__3) ! SN identifier
real*8 BJD_2 (nr__3) ! (d) [1331.15/2689.25] Time of first light;
* BJD-2470000.0
real*4 e_BJD_1 (nr__3) ! (d) [0.01/3.3] Uncertainty in BJD
real*4 tRise_1 (nr__3) ! (d) [2.8/28.01]? Rise time (1)
real*4 PLbeta1_1 (nr__3) ! [0.01/3.9] Mean power law index,
* t^{beta}_1_^
real*4 e_PLbeta1_1(nr__3) ! [0.01/0.95] Uncertainty in PLbeta1
real*8 PLbeta2_1 (nr__3) ! [-0.094/-0.0001] Mean power law index,
* t^{beta}_2_^
real*4 e_PLbeta2_1(nr__3) ! [0.0001/0.03] Uncertainty in PLbeta2
real*4 MedRad (nr__3) ! (Rsun) [1/223.2] Median (50th percentile)
* companion Roche lobe radius
real*4 e_MedRad (nr__3) ! (Rsun) [0/55.5] Uncertainty in MedRad
real*4 MaxRad (nr__3) ! (Rsun) [1/223.2] The 3{sigma} (99.9th
* percentile) companion Roche lobe radius
real*4 MedSep (nr__3) ! (10+12cm) [0.2/40.2] Median (50th percentile)
* companion Roche lobe radius
real*4 e_MedSep (nr__3) ! (10+12cm) [0/9.99] Uncertainty in MedSep
real*4 MaxSep (nr__3) ! (10+12cm) [0.2/40.2] The 3{sigma} (99.9th percentile)
* companion Roche lobe radius
real*8 lnZ_1 (nr__3) ! ([-]) [-7151.4/11354.3] Natural log of Bayesian
* Evidence Z from dynesty
real*4 e_lnZ_1 (nr__3) ! ([-]) [0.04/0.2] Uncertainty in lnZ
real*8 BIC_1 (nr__3) ! [-22733.8/14311.4] Bayesian Information
* Criterion
integer*4 Mask_1 (nr__3) ! [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(i__),OName(i__),TmagP(i__),e_TmagP(i__),TmagL(i__),
+ delTmag(i__),tPeak(i__),TMag(i__),e_TMag(i__),nuLum(i__),
+ e_nuLum(i__),z(i__),E_B_V(i__),Reduc(i__),Treduc(i__),
+ Qual(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),OName(i__),TmagP(i__),e_TmagP(i__),TmagL(i__),
+ delTmag(i__),tPeak(i__),TMag(i__),e_TMag(i__),nuLum(i__),
+ e_nuLum(i__),z(i__),E_B_V(i__),Reduc(i__),Treduc(i__),
+ Qual(i__)
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(i__),BJD(i__),CRate(i__),e_CRate(i__),Frac(i__),
+ e_Frac(i__),Tmag_1(i__),e_Tmag_1(i__),CalOff(i__)
if(ar__1(56:62) .EQ. '') Tmag_1(i__) = rNULL__
if(ar__1(64:70) .EQ. '') e_Tmag_1(i__) = rNULL__
if(ar__1(72:81) .EQ. '') CalOff(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),BJD(i__),CRate(i__),e_CRate(i__),Frac(i__),
+ e_Frac(i__),Tmag_1(i__),e_Tmag_1(i__),CalOff(i__)
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(i__),BJD_1(i__),e_BJD(i__),tRise(i__),PLbeta1(i__),
+ e_PLbeta1(i__),PLbeta2(i__),e_PLbeta2(i__),lnZ(i__),
+ e_lnZ(i__),BIC(i__),Ratio(i__),Mask(i__)
if(ar__2(25:29) .EQ. '') tRise(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),BJD_1(i__),e_BJD(i__),tRise(i__),PLbeta1(i__),
+ e_PLbeta1(i__),PLbeta2(i__),e_PLbeta2(i__),lnZ(i__),
+ e_lnZ(i__),BIC(i__),Ratio(i__),Mask(i__)
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(i__),BJD_2(i__),e_BJD_1(i__),tRise_1(i__),
+ PLbeta1_1(i__),e_PLbeta1_1(i__),PLbeta2_1(i__),
+ e_PLbeta2_1(i__),MedRad(i__),e_MedRad(i__),MaxRad(i__),
+ MedSep(i__),e_MedSep(i__),MaxSep(i__),lnZ_1(i__),e_lnZ_1(i__),
+ BIC_1(i__),Mask_1(i__)
if(ar__3(25:29) .EQ. '') tRise_1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),BJD_2(i__),e_BJD_1(i__),tRise_1(i__),
+ PLbeta1_1(i__),e_PLbeta1_1(i__),PLbeta2_1(i__),
+ e_PLbeta2_1(i__),MedRad(i__),e_MedRad(i__),MaxRad(i__),
+ MedSep(i__),e_MedSep(i__),MaxSep(i__),lnZ_1(i__),e_lnZ_1(i__),
+ BIC_1(i__),Mask_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end