clear

load AIS_2003_monthly_dhdt.mat
load AIS_2004_monthly_dhdt.mat
load AIS_2005_monthly_dhdt.mat
load AIS_2006_monthly_dhdt.mat
load AIS_2007_monthly_dhdt.mat
load AIS_2008_monthly_dhdt.mat
load AIS_2009_monthly_dhdt.mat

data = [AIS_2003_monthly_dhdt AIS_2004_monthly_dhdt(:,3:14) AIS_2005_monthly_dhdt(:,3:14) AIS_2006_monthly_dhdt(:,3:14) AIS_2007_monthly_dhdt(:,3:14) AIS_2008_monthly_dhdt(:,3:14) AIS_2009_monthly_dhdt(:,3:14) ] ;


RowSums = sum(data(:, 3:86),2);

figure()
scatter(AIS_2003_monthly_dhdt(:,1),AIS_2003_monthly_dhdt(:,2),10,(12/84)*RowSums ,'.')
caxis([-1 1])
colormap(flipud(bluewhitered(100)))
colorbar
hold
A1=load('AIS_Grounded.xyz');
plot(A1(:,1),A1(:,2),'k.')


%%
load index_land_ice.mat
figure()
scatter(AIS_2003_monthly_dhdt(index_land_ice,1),AIS_2003_monthly_dhdt(index_land_ice,2),10,(12/84)*RowSums(index_land_ice) ,'.')
caxis([-1 1])
colormap(flipud(bluewhitered(100)))
colorbar


sum((12/84)*RowSums(index_land_ice))*0.001

