Index: /issm/trunk-jpl/src/m/solve/WriteData.py
===================================================================
--- /issm/trunk-jpl/src/m/solve/WriteData.py	(revision 13898)
+++ /issm/trunk-jpl/src/m/solve/WriteData.py	(revision 13899)
@@ -17,5 +17,5 @@
 	#process options
 	#  the import above the function level doesn't work?
-	from pairoptions import *
+	from pairoptions import pairoptions
 	options=pairoptions(*args)
 
@@ -101,6 +101,8 @@
 	elif strcmpi(format,'BooleanMat'):    # {{{
 
-		if not isinstance(data,numpy.ndarray):
+		if   isinstance(data,bool):
 			data=numpy.array([data])
+		elif isinstance(data,(list,tuple)):
+			data=numpy.array(data).reshape(-1,1)
 		if len(data.shape) == 1:
 			if numpy.size(data):
@@ -132,6 +134,8 @@
 	elif strcmpi(format,'IntMat'):    # {{{
 
-		if not isinstance(data,numpy.ndarray):
+		if   isinstance(data,(int,long)):
 			data=numpy.array([data])
+		elif isinstance(data,(list,tuple)):
+			data=numpy.array(data).reshape(-1,1)
 		if len(data.shape) == 1:
 			if numpy.size(data):
@@ -163,6 +167,8 @@
 	elif strcmpi(format,'DoubleMat'):    # {{{
 
-		if not isinstance(data,numpy.ndarray):
+		if   isinstance(data,(bool,int,long,float)):
 			data=numpy.array([data])
+		elif isinstance(data,(list,tuple)):
+			data=numpy.array(data).reshape(-1,1)
 		if len(data.shape) == 1:
 			if numpy.size(data):
@@ -197,6 +203,8 @@
 		recordlength=4+4    #number of records + code
 		for matrix in data:
-			if not isinstance(matrix,numpy.ndarray):
+			if   isinstance(matrix,(bool,int,long,float)):
 				matrix=numpy.array([matrix])
+			elif isinstance(matrix,(list,tuple)):
+				matrix=numpy.array(matrix).reshape(-1,1)
 			if len(matrix.shape) == 1:
 				if numpy.size(matrix):
@@ -219,6 +227,8 @@
 		#write each matrix: 
 		for matrix in data:
-			if not isinstance(matrix,numpy.ndarray):
+			if   isinstance(matrix,(bool,int,long,float)):
 				matrix=numpy.array([matrix])
+			elif isinstance(matrix,(list,tuple)):
+				matrix=numpy.array(matrix).reshape(-1,1)
 			if len(matrix.shape) == 1:
 				matrix=matrix.reshape(numpy.size(matrix),1)