function [cmap] = cejulien(nc);
% Function [cmap] = cejulien(nc);
%
%    A very cute colormap with nc= 4*q+1 levels, q a small integer
%       (this is best for centered colorscales, with zero as white) 
 %  
%  INPUT  :   nc : small odd integer (13 or 17 recommended)
%   
%  OUTPUT :	cmap(nc,3)	- colormap RGB-value array
%     (c)    Julien Emile-Geay, LDEO, Aug 2006  
%            inspired from Gerd Krahman's ce(STUFF) series of colormaps  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
% set 5 points
C = [ 0.1324    0.1011    0.8799,
      0.3163    0.7411    1.0000,
      1.0    1.0000    1.0 ,
      1.0000    0.7783    0.1665,
      0.8734    0.1333    0.1333];
 %     0.6471    0.1647    0.1647];
%      0.6275    0.3216    0.1765 ];
%      0.9358    0.0877    0.0419 ]; 
%
nk=round((nc-1)/4);
cmap=zeros(nc,3);

for ki=1:4
   xo=[ki ki+1]; % old axis
   xn=linspace(ki,ki+1,nk+1); % new axis
   bounds=(ki-1)*nk+[1:nk+1];
   cmap(bounds,1)=interp1(xo,C(ki:ki+1,1),xn);
   cmap(bounds,2)=interp1(xo,C(ki:ki+1,2),xn);
   cmap(bounds,3)=interp1(xo,C(ki:ki+1,3),xn);
end