Transport consumption reporting 

### Checklist

- [X] I am using the current [`master`](https://github.com/PyPSA/pypsa-usa/tree/master) branch
- [X] I am running on an up-to-date [`pypsa-usa` environment](https://github.com/PyPSA/pypsa-usa/blob/master/workflow/envs/environment.yaml). Update via `conda env update -f envs/environment.yaml`

### The Issue

I believe the plotter for transport consumption is double counting the non-road transport. It is grabbing both the individual links and the aggregators. The Sankey diagram PR has a potential solution.

 https://github.com/PyPSA/pypsa-usa/blob/e1d057a06ac607d37182fd46969f0e7b9dc361dc/workflow/scripts/summary_sector.py#L460-L487

### Steps To Reproduce

_No response_

### Expected Behavior

_No response_

### Error Message

_No response_

### Anything else?

_No response_

	def get_transport_consumption(
	n: pypsa.Network,
	state: Optional[str] = None,
	) -> pd.DataFrame:
	"""
	Takes load from p0 link as loads are in kVMT or similar.
	"""
	loads = n.links[
	(n.links.carrier.str.startswith("trn-"))
	& ~(n.links.carrier.str.endswith("-veh"))
	& ~(n.links.carrier == ("trn-air"))
	& ~(n.links.carrier == ("trn-rail"))
	& ~(n.links.carrier == ("trn-boat"))
	]
	if state:
	l = _get_links_in_state(n, state)
	loads = loads[loads.index.isin(l)]
	df = n.links_t.p0[loads.index].mul(n.snapshot_weightings["objective"], axis=0).T
	df.index = df.index.map(lambda x: x.split("trn-")[1])
	return df.groupby(level=0).sum().T

	match sector:
	case "res" \| "com" \| "ind":
	return get_service_consumption(n, sector, state)
	case "trn":
	return get_transport_consumption(n, state)
	case _:
	raise NotImplementedError

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Transport consumption reporting #485

Checklist

The Issue

Steps To Reproduce

Expected Behavior

Error Message

Anything else?

Metadata

Assignees

Labels

Type

Projects

Milestone

Relationships

Development