Index: /issm/trunk-jpl/src/m/plot/colormaps/demmap.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/colormaps/demmap.m	(revision 14516)
+++ /issm/trunk-jpl/src/m/plot/colormaps/demmap.m	(revision 14517)
@@ -1,15 +1,27 @@
-function cmap = demmap(ncolors,minZ,maxZ);
+function cmap = demmap(ncolors,minZ,maxZ,varargin);
 %DEMMAP - concatenate sea and land color depending on zmin and zmax
 %
 %   Usage:
-%      cmap = demmap(n,zmin,zmax)
+%      cmap = demmap(n,zmin,zmax,varargin)
+%
+%   Example:
+%      cmap = demmap(50,-300,1200);
+%      cmap = demmap(50,-300,1200,'dem');
+%      cmap = demmap(50,-300,1200,'ibcao');
 
+%Input checks
 if nargin<3,
 	help landseacolor
-	error('3 arguments necessary');
+	error('3 or 4 arguments necessary');
+elseif nargin>4,
+	help landseacolor
+	error('3 or 4 arguments necessary');
 end
 
-if minZ == maxZ;
-	maxZ = minZ+1;
+if nargin==4,
+	colorscheme = varargin{1};
+	if ~ischar(colorscheme), error('color scheme should be a string'); end
+else
+	colorscheme = 'dem';
 end
 
@@ -60,3 +72,8 @@
 
 clim = [cmn cmx];
-cmap = [seacolor(nsea);landcolor(nland).^1.3];
+
+if strcmpi(colorscheme,'dem'),
+	cmap = [seacolor(nsea);landcolor(nland).^1.3];
+elseif strcmpi(colorscheme,'ibcao');
+	cmap = ibcao(nsea,nland);
+end
Index: /issm/trunk-jpl/src/m/plot/colormaps/ibcao.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/colormaps/ibcao.m	(revision 14517)
+++ /issm/trunk-jpl/src/m/plot/colormaps/ibcao.m	(revision 14517)
@@ -0,0 +1,69 @@
+function map = ibcap(nsea,nland)
+%IBCAO - IBCAO color map
+%
+%   Usage:
+%      map = ibcap(nsea,nland)
+%
+%Downloaded from http://soliton.vm.bytemark.co.uk/pub/cpt-city/ibcao
+
+Jsea = [ ...
+0.18039 0.29020 0.57255
+0.18039 0.29020 0.57255
+0.05882 0.44314 0.65490
+0.05882 0.44314 0.65490
+0.02745 0.49804 0.73725
+0.02745 0.49804 0.73725
+0.01176 0.54510 0.78824
+0.01176 0.54510 0.78824
+0.00784 0.63529 0.83922
+0.00784 0.63529 0.83922
+0.06667 0.71765 0.86667
+0.06667 0.71765 0.86667
+0.17647 0.75294 1.00000
+0.17647 0.75294 1.00000
+0.23529 0.76471 0.85882
+0.23529 0.76471 0.85882
+0.24314 0.76471 0.83922
+0.24314 0.76471 0.83922
+0.25882 0.76078 0.81176
+0.25882 0.76078 0.81176
+0.27451 0.76078 0.76078
+0.27451 0.76078 0.76078
+0.41961 0.78431 0.74902
+0.41961 0.78431 0.74902
+0.60000 0.83137 0.74902
+0.60000 0.83137 0.74902
+];
+
+Jland = [ ...
+0.85098 0.84314 0.30588
+0.85098 0.84314 0.30588
+0.93333 0.89020 0.41961
+0.93333 0.89020 0.41961
+0.93725 0.80784 0.35686
+0.93725 0.80784 0.35686
+0.89804 0.74510 0.31765
+0.89804 0.74510 0.31765
+0.85098 0.63922 0.21961
+0.85098 0.63922 0.21961
+0.75686 0.55294 0.22353
+0.75686 0.55294 0.22353
+0.71765 0.50980 0.22353
+0.71765 0.50980 0.22353
+0.68627 0.48235 0.21961
+0.68627 0.48235 0.21961
+0.65490 0.45882 0.21569
+0.65490 0.45882 0.21569
+0.58824 0.39608 0.20392
+0.58824 0.39608 0.20392
+1.00000 1.00000 1.00000
+];
+
+lsea = length(Jsea);
+ysea = interp1(1:lsea,Jsea,linspace(1,lsea,nsea),'*linear');
+
+lland = length(Jland);
+yland = interp1(1:lland,Jland,linspace(1,lland,nland),'*linear');
+
+
+map=[ysea;yland];