function DM_plot_dune(handles)

global data

%Plot data on axes
data.handles.profile_line.dune = [];
hold(handles.dune,'off') %release any hold
h = pcolor(handles.dune,data.x,data.y,data.znew);
data.handles.pcolor.dune = h;
shading(handles.dune,'flat')
caxis(handles.dune,[-3 3])
%set(gca,'dataaspectratio',[1 1 1])
set(handles.dune,'fontsize',8)
set(handles.dune,'xtick',[])
set(handles.dune,'ytick',[])
%zoom(handles.dune,'on')
hold(handles.dune,'on')
[c,h] = contour(handles.dune,data.x,data.y,data.znew,[0:0.5:3],'k');
%clabel(h);
data.handles.contour.dune = h;
xlim(handles.dune,get(handles.nodune,'XLim'))
ylim(handles.dune,get(handles.nodune,'YLim'))