Index: /issm/trunk/src/m/utils/Meca/paterson.m
===================================================================
--- /issm/trunk/src/m/utils/Meca/paterson.m	(revision 6967)
+++ /issm/trunk/src/m/utils/Meca/paterson.m	(revision 6968)
@@ -38,2 +38,5 @@
 pos10=find(-5<=T & T<-2);    rigidity(pos10)=10^8*(-0.003748937622487*(T(pos10)+5).^3+0.015290593619213*(T(pos10)+5).^2  -0.048160403003748*(T(pos10)+5)+  0.854987973338348);
 pos11=find(-2<=T);           rigidity(pos11)=10^8*(-0.003748937622488*(T(pos11)+2).^3-0.018449844983174*(T(pos11)+2).^2  -0.057638157095631*(T(pos11)+2)+  0.746900791092860);
+
+%Now make sure that rigidity is positive
+pos=find(rigidity<0);        rigidity(pos)=10^6;