
ncdata2 = nc_getall('d:\checkouts\kustlijnkaartenboek\scripts\BKL_TKL_TND.nc');
ncdata = nc_getall('d:\checkouts\kustlijnkaartenboek\scripts\BKL_TKL_TND.nc');

ncdata2 = nc_getall('d:\checkouts\kustlijnkaartenboek\scripts\BKL_TKL_TND.nc');
ncdata = nc_getall('d:\checkouts\kustlijnkaartenboek\scripts\BKL_TKL_TND.nc');

%%
f1 = clf;
hold on
names = {'trend','basal_coastline','testing_coastline'};
for k = 1:length(names)
    for i = 1:length(ncdata.id.data)
        j = find(ncdata2.id.data== ncdata.id.data(i));
        plot(ncdata.(names{k}).data(:,i),ncdata2.(names{k}).data(1:end-1,j))
    end
    xlabel('old')
    ylabel('new')
    title(names{k})
    box on
    saveas(f1,[names{k} '.png'])
end
