combineColormaps

function cmap = combineColormaps(low, high, cutoff)
% Creates a new colormap by combining two colormaps. Low and high specificy
% the minimum and maximum of your data (clim) and cutoff the value at which
% to switch to the other colormap
%
% Example:
%
% [X,Y] = meshgrid(-10:0.1:10,-10:0.1:10);
% f = sinc(sqrt((X/pi).^2+(Y/pi).^2));
% clim = [min(f(:)),max(f(:))];
% imagesc(f,clim)
% cmap = combineColormaps(clim(1), clim(2), 0.2);
% colormap(cmap)
% colorbar