Index: /issm/trunk-jpl/src/m/materials/cuffey.m
===================================================================
--- /issm/trunk-jpl/src/m/materials/cuffey.m	(revision 20627)
+++ /issm/trunk-jpl/src/m/materials/cuffey.m	(revision 20628)
@@ -25,5 +25,5 @@
 % rigidity=fittedmodel(temperature);
 
-rigidity=zeros(length(T),1);
+rigidity=zeros(size(T));
 pos=find(T<=-45);         rigidity(pos)=10^8*(-0.000396645116301*(T(pos)+50).^3+ 0.013345579471334*(T(pos)+50).^2  -0.356868703259105*(T(pos)+50)+7.272363035371383);
 pos=find(-45<=T & T<-40); rigidity(pos)=10^8*(-0.000396645116301*(T(pos)+45).^3+ 0.007395902726819*(T(pos)+45).^2  -0.253161292268336*(T(pos)+45)+5.772078366321591);
Index: /issm/trunk-jpl/src/m/materials/cuffey.py
===================================================================
--- /issm/trunk-jpl/src/m/materials/cuffey.py	(revision 20627)
+++ /issm/trunk-jpl/src/m/materials/cuffey.py	(revision 20628)
@@ -16,5 +16,5 @@
 		raise RuntimeError("input temperature should be in Kelvin (positive)")
 	
-	T = temperature.reshape(-1,)-273.15
+	T = temperature-273.15
 
 	rigidity=numpy.zeros_like(T)
Index: /issm/trunk-jpl/src/m/materials/cuffeytemperate.m
===================================================================
--- /issm/trunk-jpl/src/m/materials/cuffeytemperate.m	(revision 20627)
+++ /issm/trunk-jpl/src/m/materials/cuffeytemperate.m	(revision 20628)
@@ -19,3 +19,3 @@
 end
 
-rigidity=reshape(cuffey(temperature), size(temperature)).*(1+181.25*max(0., min(0.01, waterfraction))).^(-1/stressexp);
+rigidity=cuffey(temperature).*(1+181.25*max(0., min(0.01, waterfraction))).^(-1/stressexp);
Index: /issm/trunk-jpl/src/m/materials/cuffeytemperate.py
===================================================================
--- /issm/trunk-jpl/src/m/materials/cuffeytemperate.py	(revision 20627)
+++ /issm/trunk-jpl/src/m/materials/cuffeytemperate.py	(revision 20628)
@@ -25,5 +25,5 @@
 	end
 
-	rigidity=numpy.multiply(cuffey(temperature).reshape(temperature.shape), (1*numpy.ones(waterfraction.shape)+181.25*numpy.maximum(numpy.zeros(waterfraction.shape), numpy.minimum(0.01*numpy.ones(waterfraction.shape), waterfraction)))**(-1/stressexp));
+	rigidity=numpy.multiply(cuffey(temperature), (1*numpy.ones(waterfraction.shape)+181.25*numpy.maximum(numpy.zeros(waterfraction.shape), numpy.minimum(0.01*numpy.ones(waterfraction.shape), waterfraction)))**(-1/stressexp));
 
 	return rigidity