ax1.plot(row['geom_path_ln3d_%s_resampled' % projstr].xy[0], row['geom_path_ln3d_%s_resampled' % projstr].xy[1],
            'b', linestyle='dashed', label='thalweg ($S_{xy}$)')