Index: ../trunk-jpl/src/m/materials/buddjacka.m
===================================================================
--- ../trunk-jpl/src/m/materials/buddjacka.m	(revision 20683)
+++ ../trunk-jpl/src/m/materials/buddjacka.m	(revision 20684)
@@ -42,6 +42,7 @@
 %semilogy(Temp,Acall,'--k',Temp,Abjall,'k')
 %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);
 pos=find(-40<T & T<=-35);