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