Changeset 17879

Timestamp:

04/29/14 13:26:29 (11 years ago)

Author:

cborstad

Message:

CHG: fixed array indexing problem, nodes list passed to md.extract should be rank-1

File:

• issm/trunk-jpl/src/m/classes/model.py (modified) (1 diff)

issm/trunk-jpl/src/m/classes/model.py

@@ -268 +268 @@
                                                 #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