gas_resource = solph.components.Source(
        label="rgas",
        outputs={
            bgas: solph.Flow(
                nominal_value=fossil_share