d=DATA{:,1};D=DATA{:,2};
min_year=year(min(d));
max_year=year(max(d));
A=zeros(12,max_year-min_year+1);
for k=1:length(d)
m=month(d(k));y=year(d(k))-min_year+1;
A(12-m+1,y)=D(k);
end
in=[min_year:max_year];%% added
contourf(in,1:12,A)%% edited
map = [0.8 0.1 0.1
0.9 0.9 0.1
1 1 1
0 0.9 0.1
0.1 0.1 1];
colormap(map)
colorbar
you can change the color by changing values in map
