%% Simple example to plot a topview of the data
clear
close all
savefigs=0; %Toggle

% Load

%Mpath='F:\check_outs\tudelft\vlugtenburg\'
Mpath='C:\Madeschipper\check_out\OpenEarthRawData\tudelft\vlugtenburg\'


load([Mpath 'AllVbeachsurveys.mat'])

load('cmapVBmin10to5.mat')

%% plot gridded (from transects) data
 figure(227)
       set(gcf,'PaperPosition',[0  0  40 20])
        set(gcf,'Position',[ 8         179        1444         487])
 
 i=0;
    for i_date=[1 16 length(Vsurvey)]
       i=i+1;
       subplot(1,3,i)
  
        
     
        pk=pcolor(Vsurvey(i_date).transects.y,Vsurvey(i_date).transects.cr_sh,Vsurvey(i_date).transects.elevation)
        hold all
        
        [C,h] = contour(Vsurvey(i_date).transects.y,Vsurvey(i_date).transects.cr_sh,Vsurvey(i_date).transects.elevation,[-9:1:4],'--k');
        clabel(C,h,'LabelSpacing',72,'FontSize',10,'Color','k','Rotation',0) ;
        
        shading flat
        hold all
        
     %   set(gca,'Xtick',[flipud(Vsurvey(i_date).transects.y)])
     %   set(gca,'Xticklabel',num2str([flipud(Vsurvey(i_date).transects.y)],'%4.0f'))
        
        xlim([0 1750])
        
        xlabel('Alongshore distance [m]')
        
        ylim([0 1200])
        if i==1
            ylabel('Cross shore distance [m]')
        end
        
        caxis([-10 5])
        %colorbar
        colormap(cmap)
       % set(gca,'position',[0.1300    0.1100    0.7750    0.8150])
       % legend([pk  h], 'Elevation [m]',[ 'contour ' Vsurvey(i_date).date([4:6 1:3 7:10])] )
        title([ 'Topography survey  ' Vsurvey(i_date).date([1:3 7:10]) ', gridded (based on transects)'])
        
        
        x0=30;           y0=1050;  theta0 =-131.65534968173412+90;
        h=100; w=100 ;
        [x_ds3,y_ds3]=northarrow(h,w,x0,y0,theta0);
        
  
    end
      set(gcf,'renderer','zbuffer')
        hold off
        
        if savefigs==1
            print ([ 'EvolutionTopo.jpg'],'-djpeg')
        end