ax2.hlines(y=z[indVelZero], xmin=sXY[indVelStart],
               xmax=sXY[indVelZero],
               color='grey', linestyle='--', linewidth=1.5, label=(r'$\Delta z = %.1f$$m$' % deltaz))