Index: /issm/trunk-jpl/src/m/plot/colormaps/demmap.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/colormaps/demmap.m	(revision 25620)
+++ /issm/trunk-jpl/src/m/plot/colormaps/demmap.m	(revision 25621)
@@ -12,8 +12,8 @@
 %Input checks
 if nargin<3,
-	help landseacolor
+	help demmap
 	error('3 or 4 arguments necessary');
 elseif nargin>4,
-	help landseacolor
+	help demmap
 	error('3 or 4 arguments necessary');
 end
Index: /issm/trunk-jpl/src/m/plot/processdata.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/processdata.m	(revision 25620)
+++ /issm/trunk-jpl/src/m/plot/processdata.m	(revision 25621)
@@ -205,6 +205,9 @@
 
 %layer projection? 
-if getfieldvalue(options,'layer',0)>=1,
+if getfieldvalue(options,'layer',0)>=1
 	data=project2d(md,data,getfieldvalue(options,'layer')); %project onto 2d mesh
+end
+if getfieldvalue(options,'depthaverage',0)
+	data=depthaverage(md,data)); %project onto 2d mesh
 end
 
Index: /issm/trunk-jpl/src/m/plot/processmesh.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/processmesh.m	(revision 25620)
+++ /issm/trunk-jpl/src/m/plot/processmesh.m	(revision 25621)
@@ -64,5 +64,5 @@
 
 %is it a 2d plot?
-if md.mesh.dimension()==2 | getfieldvalue(options,'layer',0)>=1,
+if md.mesh.dimension()==2 || getfieldvalue(options,'layer',0)>=1 || getfieldvalue(options,'depthaverage',0)
 	is2d=1;
 else
@@ -71,5 +71,5 @@
 
 %layer projection? 
-if getfieldvalue(options,'layer',0)>=1,
+if getfieldvalue(options,'layer',0)>=1 || getfieldvalue(options,'depthaverage',0)
 	if strcmpi(getfieldvalue(options,'coord','xy'),'latlon'),
 		error('processmesh error message: cannot work with 3D meshes for now');