Changeset 17876 for issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py

Timestamp:

04/29/14 10:36:31 (12 years ago)

Author:

cborstad

Message:

CHG: treat stressbalance spc's as rank-1

File:

• issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py (modified) (2 diffs)

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,))
-        md.stressbalance.spcvy=float('nan')*numpy.ones((md.mesh.numberofvertices,))
-        md.stressbalance.spcvz=float('nan')*numpy.ones((md.mesh.numberofvertices,))
+        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))
@@ -68 +68 @@
         else:
                 pos=numpy.nonzero(md.mesh.vertexonboundary)[0]
-        md.stressbalance.spcvx[pos[:]]=0
-        md.stressbalance.spcvy[pos[:]]=0
-        md.stressbalance.spcvz[pos[:]]=0
+        md.stressbalance.spcvx[pos]=0
+        md.stressbalance.spcvy[pos]=0
+        md.stressbalance.spcvz[pos]=0
 
         #Dirichlet Values