def total_energy_ejected(t):
    """
    Thermal and kinetic energy released by each type of SN up to the time t after the explosion
    where tc is the characteristic cooling time of the shell sorrounding the remnant (53000 yrs)
    from Ferrini & Poggiantti, 1993, ApJ, 410, 44F

    """
    if t <= 0:
        return 0.0
    tc = 5.3e-05
    if t > tc:
        rt = (tc / t) ** 0.4
        return 1 - 0.44 * rt ** 2 * (1 - 0.41 * rt) - 0.22 * rt ** 2
    else:
        return 9811.32 * t
