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