Index: /issm/trunk-jpl/src/m/mech/analyticaldamage.py
===================================================================
--- /issm/trunk-jpl/src/m/mech/analyticaldamage.py	(revision 17902)
+++ /issm/trunk-jpl/src/m/mech/analyticaldamage.py	(revision 17903)
@@ -53,5 +53,5 @@
 	assert len(kwargs)==0, 'error, unexpected or misspelled kwargs'
 
-	if len(sigmab)==1:
+	if isinstance(sigmab,(int,float)):
 		sigmab=sigmab*npy.ones((md.mesh.numberofvertices,))
 
@@ -83,5 +83,5 @@
 	
 	backstress=npy.zeros((md.mesh.numberofvertices,))
-	
+
 	# backstress to bring D down to one 
 	backstress[pos]=T[pos]-(1.-D[pos])*B[pos]*npy.sign(ex[pos])*(2.+a[pos])*npy.abs(ex[pos])**(1./n[pos])/(1.+a[pos]+a[pos]**2)**((n[pos]-1.)/2./n[pos])
Index: /issm/trunk-jpl/src/m/mech/thomasparams.py
===================================================================
--- /issm/trunk-jpl/src/m/mech/thomasparams.py	(revision 17902)
+++ /issm/trunk-jpl/src/m/mech/thomasparams.py	(revision 17903)
@@ -134,5 +134,5 @@
 	elif eq=='Weertman2D':
 		theta=1./9
-		a=npy.ones((md.mesh.numberofvertices,1))
+		a=npy.ones((md.mesh.numberofvertices,))
 	elif eq=='Thomas':
 		theta=((1.+a+a**2+b**2)**((n-1.)/2.))/(npy.abs(2.+a)**n)