Index: /issm/trunk-jpl/src/m/shp/shpread.m
===================================================================
--- /issm/trunk-jpl/src/m/shp/shpread.m	(revision 16338)
+++ /issm/trunk-jpl/src/m/shp/shpread.m	(revision 16339)
@@ -40,5 +40,9 @@
 		Struct(end).density=1;
 		Struct(end).closed=1;
-		Struct(end).name=num2str(shp(i).id);
+		if isfield(shp,'id'),
+			Struct(end).name=num2str(shp(i).id);
+		else
+			Struct(end).name='';
+		end
 	end
 end
Index: /issm/trunk-jpl/src/m/shp/shpwrite.m
===================================================================
--- /issm/trunk-jpl/src/m/shp/shpwrite.m	(revision 16339)
+++ /issm/trunk-jpl/src/m/shp/shpwrite.m	(revision 16339)
@@ -0,0 +1,26 @@
+function shpwrite(shp,filename)
+%SHPWRITE - write a shape file from a contour structure
+%
+%   Usage:
+%      shpwrite(shp,filename)
+%
+%   Example:
+%      shpwrite(shp,'domainoutline.shp')
+%
+%   See also SHPREAD
+
+
+%initialize number of profile
+count=0;
+
+
+contours=struct([]);
+for i=1:length(shp),
+	contours(i).Geometry='Polygon';
+	contours(i).id=i;
+	contours(i).X=shp(i).x;
+	contours(i).Y=shp(i).y;
+end
+	
+shapewrite(contours,filename);
+end