def yearly_footprint(embodied_kg, watts, hours_per_day, years, grid_g_per_kwh):
kwh_per_year = watts * hours_per_day * 365 / 1000
operational_kg = kwh_per_year * grid_g_per_kwh / 1000
return embodied_kg / years + operational_kg
scenarios = [("Ontario, 3 years", 3, 53), ("Ontario, 6 years", 6, 53),
("500 g/kWh grid, 3 years", 3, 500), ("500 g/kWh grid, 6 years", 6, 500)]
for name, years, grid in scenarios:
kg = yearly_footprint(300, 20, 6, years, grid)
print(f"{name:26}{kg:7.1f} kg CO2e per year")