Index: /issm/trunk-jpl/src/m/classes/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.py	(revision 17878)
+++ /issm/trunk-jpl/src/m/classes/model.py	(revision 17879)
@@ -268,20 +268,20 @@
 						#size = number of nodes * n
 						if   fieldsize[0]==numberofvertices1:
-							setattr(getattr(md2,fieldi),fieldj,field[pos_node,:])
+							setattr(getattr(md2,fieldi),fieldj,field[pos_node])
 						elif fieldsize[0]==numberofvertices1+1:
-							setattr(getattr(md2,fieldi),fieldj,numpy.vstack((field[pos_node,:],field[-1,:])))
+							setattr(getattr(md2,fieldi),fieldj,numpy.vstack((field[pos_node],field[-1,:])))
 						#size = number of elements * n
 						elif fieldsize[0]==numberofelements1:
-							setattr(getattr(md2,fieldi),fieldj,field[pos_elem,:])
+							setattr(getattr(md2,fieldi),fieldj,field[pos_elem])
 			else:
 				if len(fieldsize):
 					#size = number of nodes * n
 					if   fieldsize[0]==numberofvertices1:
-						setattr(md2,fieldi,field[pos_node,:])
+						setattr(md2,fieldi,field[pos_node])
 					elif fieldsize[0]==numberofvertices1+1:
-						setattr(md2,fieldi,numpy.hstack((field[pos_node,:],field[-1,:])))
+						setattr(md2,fieldi,numpy.hstack((field[pos_node],field[-1,:])))
 					#size = number of elements * n
 					elif fieldsize[0]==numberofelements1:
-						setattr(md2,fieldi,field[pos_elem,:])
+						setattr(md2,fieldi,field[pos_elem])
 
 		#modify some specific fields
