Index: /issm/trunk-jpl/src/m/shp/shp2exp.m
===================================================================
--- /issm/trunk-jpl/src/m/shp/shp2exp.m	(revision 21348)
+++ /issm/trunk-jpl/src/m/shp/shp2exp.m	(revision 21349)
@@ -33,14 +33,13 @@
 			expstruct(end).nods=length(x);
 			expstruct(end).density=1;
-			expstruct(end).closed=1;
+			expstruct(end).closed=0;
 			%exp(end).name=num2str(shp(i).id);
 		elseif strcmpi(shp(i).Geometry,'Line'),
 			x=shp(i).X; y=shp(i).Y;
-			x(end)=x(1); y(end)=y(1);
 			expstruct(end+1).x=x;
 			expstruct(end).y=y;
 			expstruct(end).nods=length(x);
 			expstruct(end).density=1;
-			expstruct(end).closed=1;
+			expstruct(end).closed=0;
 		end
 	end
Index: /issm/trunk-jpl/src/m/shp/shp2exp.py
===================================================================
--- /issm/trunk-jpl/src/m/shp/shp2exp.py	(revision 21348)
+++ /issm/trunk-jpl/src/m/shp/shp2exp.py	(revision 21349)
@@ -26,5 +26,5 @@
 
 	shp=shapefile.Reader(shapefilename)
-	expdict=dict(closed=1,density=1)
+	expdict=dict(density=1)
 
 	x=[]
@@ -33,4 +33,5 @@
 		geom=shp.shapes()[i].shapeType
 		if geom==5: # polygon
+         expdict['closed']=1
 			tmpx=[p[0] for p in shp.shapes()[i].points]
 			tmpy=[q[1] for q in shp.shapes()[i].points]
@@ -38,4 +39,5 @@
 			y.append(tmpy)
 		elif geom==3: # line
+         expdict['closed']=0
 			tmpx=[p[0] for p in shp.shapes()[i].points]
 			tmpy=[q[1] for q in shp.shapes()[i].points]
@@ -43,4 +45,5 @@
 			y.append(tmpy)
 		elif geom==1: # point
+         expdict['closed']=0
 			x.append(shp.shapes()[i].points[0][0])
 			y.append(shp.shapes()[i].points[0][1])
Index: /issm/trunk-jpl/src/py3/shp/shp2exp.py
===================================================================
--- /issm/trunk-jpl/src/py3/shp/shp2exp.py	(revision 21348)
+++ /issm/trunk-jpl/src/py3/shp/shp2exp.py	(revision 21349)
@@ -26,5 +26,5 @@
 
 	shp=shapefile.Reader(shapefilename)
-	expdict=dict(closed=1,density=1)
+	expdict=dict(density=1)
 
 	x=[]
@@ -33,4 +33,5 @@
 		geom=shp.shapes()[i].shapeType
 		if geom==5: # polygon
+         expdict['closed']=1
 			tmpx=[p[0] for p in shp.shapes()[i].points]
 			tmpy=[q[1] for q in shp.shapes()[i].points]
@@ -38,4 +39,5 @@
 			y.append(tmpy)
 		elif geom==3: # line
+         expdict['closed']=0
 			tmpx=[p[0] for p in shp.shapes()[i].points]
 			tmpy=[q[1] for q in shp.shapes()[i].points]
@@ -43,4 +45,5 @@
 			y.append(tmpy)
 		elif geom==1: # point
+         expdict['closed']=0
 			x.append(shp.shapes()[i].points[0][0])
 			y.append(shp.shapes()[i].points[0][1])