fig, ax = plt.subplots() ax.set_xlabel(
'
Basin area [$km^2$]
'
)