def electricity_marginal_cost(m, z, tp):
    """Return marginal cost of production per MWh in load_zone z during timepoint tp."""
    return m.dual[m.Energy_Balance[z, tp]]/m.bring_timepoint_costs_to_base_year[tp]