cmap = [ ...
    1.0000    1.0000    1.0000
    0.6667    1.0000    0.6667
    0.3333    1.0000    0.3333
         0    0.9192         0
    1.0000    1.0000         0
    1.0000    0.5000         0
    0.8000         0         0
    1            0           1];
colorbar;

NR = length(cmap); %length of original cmap
NF=100; %number of rows in the fine colormap matrix

xr = 0:1/(NR-1):1; %vector of length NR spanning 0:1
xf= 0:1/(NF-1):1; %vector of length NF spanning 0:1;

%interpolate each column to the finer grid
cmapf = [ (interp1(xr,cmap(:,1), xf))' (interp1(xr,cmap(:,2), xf))'  (interp1(xr,cmap(:,3), xf))'];


