clear

load AIS_2010_monthly_dhdt.mat
load AIS_2011_monthly_dhdt.mat
load AIS_2012_monthly_dhdt.mat
load AIS_2013_monthly_dhdt.mat
load AIS_2014_monthly_dhdt.mat
load AIS_2015_monthly_dhdt.mat
load AIS_2016_monthly_dhdt.mat
load AIS_2017_monthly_dhdt.mat

data = [AIS_2010_monthly_dhdt AIS_2011_monthly_dhdt(:,3:14) AIS_2012_monthly_dhdt(:,3:14) AIS_2013_monthly_dhdt(:,3:14) AIS_2014_monthly_dhdt(:,3:14)  AIS_2015_monthly_dhdt(:,3:14) AIS_2016_monthly_dhdt(:,3:14) AIS_2017_monthly_dhdt(:,3:14)] ;

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

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



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

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

