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