Changeset 17873

Timestamp:

04/29/14 09:26:55 (11 years ago)

Author:

cborstad

Message:

BUG: reshape no longer needed if spc's are rank-1

File:

• issm/trunk-jpl/src/m/parameterization/setflowequation.py (modified) (1 diff)

issm/trunk-jpl/src/m/parameterization/setflowequation.py

@@ -89 +89 @@
                                               numpy.logical_not(numpy.isnan(md.stressbalance.spcvy)).astype(int)+ \
                                               numpy.logical_not(numpy.isnan(md.stressbalance.spcvz)).astype(int)==3, \
-                                              numpy.logical_and(nodeonHO,nodeonFS).reshape(-1,1)).astype(int)    #find all the nodes on the boundary of the domain without icefront
+                                              numpy.logical_and(nodeonHO,nodeonFS)).astype(int)    #find all the nodes on the boundary of the domain without icefront
 #               fullspcelems=double(sum(fullspcnodes(md.mesh.elements),2)==6);         %find all the nodes on the boundary of the domain without icefront
                 fullspcelems=(numpy.sum(fullspcnodes[md.mesh.elements-1],axis=1)==6).astype(int)    #find all the nodes on the boundary of the domain without icefront