# model DDC25 from Blondin et al. (2017b), accepted to MNRAS
# (arXiv:1711.09107)

########################################################################
#								       #
#		PLEASE REFERENCE BLONDIN ET AL. (2017b) 	       #
#								       #
#	      IF YOU USE THESE RESULTS IN A PUBLICATION		       #
#								       #
########################################################################


### HYDRODYNAMICAL INPUT

The initial conditions are given at 0.50 days past explosion in the
file named DDC25_SN_HYDRO_DATA_0.50d. This file gives the usual hydro
variables at each grid point (radius, velocity, temperature, density)
+ the rosseland mean optical depth scale and mass absorption
coefficient, followed by the mass fractions of 26 elements and 105
isotopes.

#Variable	      	 unit	 	  NOTES
Radius		      	 10^10 cm	  assumes homologous expansion (R = V * t)
Velocity	      	 km/s
Sigma		      	 N/A		  velocity gradient = dlnV/dlnR - 1 = 0 for homologous expansion
Temperature	      	 10^4 K
Density		      	 g/cm^3
Atom density	      	 /cm^3
Electron density      	 /cm^3
Rosseland mean opacity	 /10^10 cm	  Rosseland-mean optical-depth scale
Kappa	       		 cm^2/g		  = (Rosseland mean opacity) / Density / 10^10

COMMENTS:

- the spatial grid starts from the outermost point (R=3.0239971e14 cm,
  V=6.9999933e4 km/s) and progresses inwards until the innermost point
  (R=3.6653746e12 cm, V=8.4846634e2 km/s).

- the elemental mass fractions are set to the sum of the corresponding
  isotope mass fractions when present (e.g. CHL(I) = CHL_35(I) +
  CHL_36(I) + CHL_37(I) at each spatial index I).

- ALL MASS FRACTIONS ARE GIVEN AT 0.50 DAYS PAST EXPLOSION (i.e. the
  isotope mass fractions take into account radioactive decays). See
  Tables 3-5 in Dessart et al. (2014), MNRAS, 441, 532 for the full
  list of 1-step and 2-step decay chains considered in these
  calculations.

- the following naming convention is used to designate elements:

  Hydrogen 	HYD
  Helium	HE
  Carbon 	CARB
  Nitrogen 	NIT
  Oxygen 	OXY
  Fluorine	FLU
  Neon 		NEON
  Sodium 	SOD
  Magnesium 	MAG
  Aluminum 	ALUM
  Silicon 	SIL
  Phosphorous 	PHOS
  Sulfur 	SUL
  Chlorine	CHL
  Argon 	ARG
  Potassium	POT
  Calcium 	CAL
  Scandium	SCAN
  Titanium 	TIT
  Vanadium 	VAN
  Chromium 	CHRO
  Manganese 	MAN
  Iron 		IRON
  Cobalt 	COB
  Nickel 	NICK
  Barium	BAR
  

### SYNTHETIC SPECTRA

Synthetic spectra are given at each time step between 1.30 d and
223.4 d since explosion (55 spectra in total). They were interpolated
on a 5-Angstrom scale between 2500 and 25500 Angstroms (i.e., 4601
wavelengths) based on the original observer-frame calculation, which
computes the flux at >80000 wavelengths in the range 2500-25500 A.

The filenames give the time in days since explosion, e.g.,
DDC25_spec_2500_25500_interp5_018.720d.dat is at 18.72 days since
explosion. The first column gives the wavelength in Angstroms, the
second column gives the corresponding flux at that wavelength
(i.e. not the average flux over a wavelength bin) in erg/s/cm^2/A at
10pc, such that direct integration of the spectra using a given filter
bandpass + corresponding zero point will yield the absolute magnitude
in that filter. We have checked that using interpolated over binned
fluxes has negligible impact (< 0.001 mag) on the inferred
magnitudes.

For further information, please contact:

Stéphane Blondin (stephane.blondin@lam.fr)
