Index: /issm/trunk-jpl/src/m/plot/applyoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 26218)
+++ /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 26219)
@@ -289,4 +289,41 @@
 		shpdisp(filenamei,1,stylei,linewidthi,getfieldvalue(options,'unit',1));
 	end
+end
+if exist(options,'contours'),
+
+	hold on;
+	contours=getfieldvalue(options,'contours');
+	style=getfieldvalue(options,'contourstyle',{'-'});
+	linewidth=getfieldvalue(options,'linewidth',{1});
+	color=getfieldvalue(options,'contourcolor',{'r'});
+	contourheight=getfieldvalue(options,'contourheight',1); 
+
+	radius=md.solidearth.planetradius;
+	ratio=1+(contourheight*1000/radius);
+
+	
+	if ~isa(contours,'cell'),
+		contours={contours};
+	end
+	nc=length(contours);
+	if ~isa(style,'cell'), error('contour style should be a cell array'); end
+	if ~isa(linewidth,'cell'), error('contour line width should be a cell array'); end
+	if ~isa(color,'cell'), error('contour color should be a cell array'); end
+
+	for i=1:length(contours),
+		ci=contours{i};
+		if length(style)==1, sti=style{1}; else sti=style{i}; end
+		if length(color)==1, coli=color{1}; else coli=color{i}; end
+		if length(linewidth)==1, li=linewidth{1}; else li=linewidth{i}; end
+
+		for j=1:length(ci),
+			cijx=ci(j).x*ratio;
+			cijy=ci(j).y*ratio;
+			cijz=ci(j).z*ratio;
+
+			plot3(cijx,cijy,cijz,'LineWidth',li,'LineStyle',sti,'Color',coli);
+		end
+	end
+
 end