ax.plot(profile['x'][:-1], profile['y'][:-1], '.k', label='mass averaged path')