Index: /issm/trunk-jpl/test/Par/SquareShelfConstrained.py
===================================================================
--- /issm/trunk-jpl/test/Par/SquareShelfConstrained.py	(revision 12847)
+++ /issm/trunk-jpl/test/Par/SquareShelfConstrained.py	(revision 12848)
@@ -18,10 +18,15 @@
 #Initial velocity 
 mat=matio.loadmat('../Data/SquareShelfConstrained.data')
+#Reshape as Rank-1 arrays
+x=reshape(mat['x'],(size(mat['x']),))
+y=reshape(mat['y'],(size(mat['y']),))
+vx=mat['vx']
+vy=mat['vy']
 #deal with 'F' oriented matlab matrices!
 index=mat['index'].astype(float)
 index=reshape(index.T,(len(index),3),order='F')
 
-md.initialization.vx = im.InterpFromMeshToMesh2d(index, mat['x'], mat['y'], mat['vx'], md.mesh.x, md.mesh.y)
-md.initialization.vy = im.InterpFromMeshToMesh2d(index, mat['x'], mat['y'], mat['vy'], md.mesh.x, md.mesh.y)
+md.initialization.vx = im.InterpFromMeshToMesh2d(index, x, y, vx, md.mesh.x, md.mesh.y)
+md.initialization.vy = im.InterpFromMeshToMesh2d(index, x, y, vy, md.mesh.x, md.mesh.y)
 md.initialization.vz = zeros(md.mesh.numberofvertices)
 md.initialization.pressure = zeros(md.mesh.numberofvertices)