Index: /issm/trunk-jpl/src/c/python/io/FetchPythonData.cpp
===================================================================
--- /issm/trunk-jpl/src/c/python/io/FetchPythonData.cpp	(revision 13483)
+++ /issm/trunk-jpl/src/c/python/io/FetchPythonData.cpp	(revision 13484)
@@ -75,10 +75,14 @@
 	M=dims[0]; N=dims[1];
 	
-	/*retrieve internal value: */
-	dmatrix=(double*)PyArray_DATA((PyArrayObject*)py_matrix);
-
-	/*copy matrix: */
-	matrix=xNew<double>(M*N);
-	memcpy(matrix,dmatrix,(M*N)*sizeof(double));
+	if (M && N) {
+		/*retrieve internal value: */
+		dmatrix=(double*)PyArray_DATA((PyArrayObject*)py_matrix);
+
+		/*copy matrix: */
+		matrix=xNew<double>(M*N);
+		memcpy(matrix,dmatrix,(M*N)*sizeof(double));
+	}
+	else
+		matrix=NULL;
 
 	/*output: */
@@ -107,10 +111,14 @@
 	M=dims[0]; N=dims[1];
 	
-	/*retrieve internal value: */
-	dmatrix=(double*)PyArray_DATA((PyArrayObject*)py_matrix);
-
-	/*transform into integer matrix: */
-	matrix=xNew<int>(M*N);
-	for(i=0;i<M*N;i++)matrix[i]=(int)dmatrix[i];
+	if (M && N) {
+		/*retrieve internal value: */
+		dmatrix=(double*)PyArray_DATA((PyArrayObject*)py_matrix);
+
+		/*transform into integer matrix: */
+		matrix=xNew<int>(M*N);
+		for(i=0;i<M*N;i++)matrix[i]=(int)dmatrix[i];
+	}
+	else
+		matrix=NULL;
 
 	/*output: */
@@ -136,10 +144,14 @@
 	M=dims[0]; 
 	
-	/*retrieve internal value: */
-	dvector=(double*)PyArray_DATA((PyArrayObject*)py_vector);
-
-	/*copy vector: */
-	vector=xNew<double>(M);
-	memcpy(vector,dvector,(M)*sizeof(double));
+	if (M) {
+		/*retrieve internal value: */
+		dvector=(double*)PyArray_DATA((PyArrayObject*)py_vector);
+
+		/*copy vector: */
+		vector=xNew<double>(M);
+		memcpy(vector,dvector,(M)*sizeof(double));
+	}
+	else
+		vector=NULL;
 
 	/*output: */
Index: /issm/trunk-jpl/src/m/classes/bamggeom.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/bamggeom.py	(revision 13483)
+++ /issm/trunk-jpl/src/m/classes/bamggeom.py	(revision 13484)
@@ -12,5 +12,5 @@
 		self.Vertices=numpy.empty((0,3))
 		self.Edges=numpy.empty((0,3))
-		self.TangentAtEdges=numpy.empty((0,0))
+		self.TangentAtEdges=numpy.empty((0,4))
 		self.Corners=numpy.empty((0,1))
 		self.RequiredVertices=numpy.empty((0,1))