pp1 = solph.components.Source(
        label="power_plant1",
        outputs={bel: solph.Flow(nominal_value=10, variable_costs=10.25)},