Changeset 12848

Timestamp:

07/31/12 18:22:37 (13 years ago)

Author:

cborstad

Message:

reshaped x and y arrays to rank-1 for InterpFromMeshToMesh2d

File:

• issm/trunk-jpl/test/Par/SquareShelfConstrained.py (modified) (1 diff)

issm/trunk-jpl/test/Par/SquareShelfConstrained.py

@@ -18 +18 @@
 #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)