; 
zi = interp1(d,z,di); 

figure
plot(d,z,'ko')
hold on
plot(di,zi,'b.-')

%% 
% We can even compare that to other datasets:  

plot(di,bedmap2_interp(lati,loni,'surfw'),'r')
plot(di,cryosat2_interp(lati,loni),'g')

legend('original flight data','100 m flight data','bedmap2','cryosat2')

%% 
% Note: The CryoSat-2 elevations seem silly compared to the other data, but there's a good chance
% they're right.  I originally made the sample dataset by adding some noise to Bedmap2.  

%% Citing AMT
% If this funct