Index: sm/trunk-jpl/src/m/utils/Meca/effective_value.m
===================================================================
--- /issm/trunk-jpl/src/m/utils/Meca/effective_value.m	(revision 11505)
+++ 	(revision )
@@ -1,18 +1,0 @@
-function effective_value=effective_value(A)
-%EFFECTIVE_VALUE - compute the effective value of any matrix
-%
-%   effective(A)= 1/sqrt(2)*sqrt(sum(Aij^2))
-%
-%   Usage:
-%      effective_value=effective_value(A)
-
-[row,col]=size(A);
-
-SUM=0;
-for i=1:row
-    for j=1:col
-        SUM=SUM+A(i,j)^2;
-    end
-end
-
-effective_value=1/sqrt(2)*sqrt(SUM);
Index: sm/trunk-jpl/src/m/utils/Meca/glen_coefficients.m
===================================================================
--- /issm/trunk-jpl/src/m/utils/Meca/glen_coefficients.m	(revision 11505)
+++ 	(revision )
@@ -1,26 +1,0 @@
-function [A ,B]=glen_coefficients(Ao,T);
-%GLEN_COEFFICIENTS - compute A and B, given the temperature and Ao
-%
-%   Ao,T according to patterson,Ao in s-1kPa-3.
-%   values range between 3.9 and 5.3*10^-16;
-%   Return: A in s-1kPa-3. B in Pa.a1/3
-%
-%   Usage:
-%      [A ,B]=glen_coefficients(Ao,T)
-
-To=263;
-R=8.314;
-E1=139000;
-E2=60000;
-
-%For T>263,
-Ac1=Ao*exp(E1/R/To);
-%For T<263,
-Ac2=Ao*exp(E2/R/To);
-
-A=Ac1*exp(-E1/R./T);
-pos=find(T<263);
-A(pos)=Ac2*exp(-E2/R./T(pos));
-
-B=(A*365*24*3600).^(-1/3);
-