clear all
% make a simple plot of relative sea level change using matlab

% plot Jakobshavn Isbræ data
ncdisp('JI_rsl.nc');
x  = ncread('JI_rsl.nc','X');
y  = ncread('JI_rsl.nc','Y');
RSL  = ncread('JI_rsl.nc','RSL');
subplot(3,1,1)
[C,h]=contourf(x,y,RSL,'LevelStep',0.20,'LineWidth',0.5,'TextStep',0.20);
ylabel('UTM zone 22 - northing in km')
xlabel('UTM zone 22 - easting in km')
colorbar

% plot Kangerlussuaq Glacier data
ncdisp('KG_rsl.nc');
x  = ncread('KG_rsl.nc','X');
y  = ncread('KG_rsl.nc','Y');
RSL  = ncread('KG_rsl.nc','RSL');
subplot(3,1,2)
[C,h]=contourf(x,y,RSL,'LevelStep',0.20,'LineWidth',0.5,'TextStep',0.20);
ylabel('UTM zone 24 - northing in km')
xlabel('UTM zone 24 - easting in km')
colorbar

% plot Helheim Glacier data
ncdisp('HG_rsl.nc');
x  = ncread('HG_rsl.nc','X');
y  = ncread('HG_rsl.nc','Y');
RSL  = ncread('HG_rsl.nc','RSL');
subplot(3,1,3)
[C,h]=contourf(x,y,RSL,'LevelStep',0.10,'LineWidth',0.5,'TextStep',0.20);
ylabel('UTM zone 24 - northing in km')
xlabel('UTM zone 24 - easting in km')
colorbar