Changeset 17870

Timestamp:

04/29/14 07:53:30 (11 years ago)

Author:

cborstad

Message:

CHG: starting to work on the rank-1 versus rank-2 issue in python

File:

• issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py (modified) (1 diff)

issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py

@@ -42 +42 @@
                 print "SetMarineIceSheetBC warning: ice front all around the glacier, no dirichlet found. Dirichlet must be added manually."
 
-        md.stressbalance.spcvx=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
-        md.stressbalance.spcvy=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
-        md.stressbalance.spcvz=float('nan')*numpy.ones((md.mesh.numberofvertices,1))
+        md.stressbalance.spcvx=float('nan')*numpy.ones((md.mesh.numberofvertices,))
+        md.stressbalance.spcvy=float('nan')*numpy.ones((md.mesh.numberofvertices,))
+        md.stressbalance.spcvz=float('nan')*numpy.ones((md.mesh.numberofvertices,))
         md.stressbalance.referential=float('nan')*numpy.ones((md.mesh.numberofvertices,6))
         md.stressbalance.loadingforce=0*numpy.ones((md.mesh.numberofvertices,3))