energysystem.add(
        Sink(
            label="demand_th",
            inputs={bth: Flow(nominal_value=40, fix=data["demand_th"])},