Changeset 6968

Timestamp:

01/06/11 10:56:16 (14 years ago)

Author:

seroussi

Message:

force paterson to be positive (in matlab)

File:

• issm/trunk/src/m/utils/Meca/paterson.m (modified) (1 diff)

issm/trunk/src/m/utils/Meca/paterson.m ¶

@@ -38 +38 @@
 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;