Index: /issm/trunk-jpl/src/m/materials/buddjacka.m
===================================================================
--- /issm/trunk-jpl/src/m/materials/buddjacka.m	(revision 20683)
+++ /issm/trunk-jpl/src/m/materials/buddjacka.m	(revision 20684)
@@ -43,4 +43,5 @@
 %semilogy(Temp,Bcall,'--k',Temp,Bbjall,'k')
 
+rigidity=zeros(size(T));
 pos=find(T<=-40);
 rigidity(pos)=1e8*(-0.000237326134296*(T(pos)+50).^3+ 0.017054655749852*(T(pos)+50).^2-0.496435984007500*(T(pos)+50)+7.670967258739796);