Index: /issm/trunk/src/m/utils/Mesh/ComputeHessian.m
===================================================================
--- /issm/trunk/src/m/utils/Mesh/ComputeHessian.m	(revision 2921)
+++ /issm/trunk/src/m/utils/Mesh/ComputeHessian.m	(revision 2922)
@@ -43,4 +43,6 @@
 grad_elx=sum(field(index).*alpha,2); 
 grad_ely=sum(field(index).*beta,2);
+grad_elx(1)
+grad_ely(1)
 
 %Compute gradient for each grid (average of the elements around)
Index: /issm/trunk/src/m/utils/Mesh/MergeMetrics.m
===================================================================
--- /issm/trunk/src/m/utils/Mesh/MergeMetrics.m	(revision 2922)
+++ /issm/trunk/src/m/utils/Mesh/MergeMetrics.m	(revision 2922)
@@ -0,0 +1,15 @@
+function metric=MergeMetrics(metric1,metric2)
+
+M1xx=metric1(:,1); M1xy=metric1(:,2); M1yy=metric1(:,3);
+M1lambda1=0.5*((M1xx+M1yy)+sqrt(4*M1xy.^2+(M1xx-M1yy).^2));
+M1lambda2=0.5*((M1xx+M1yy)-sqrt(4*M1xy.^2+(M1xx-M1yy).^2));
+lambda1=min(M1lambda1,M1lambda2);
+
+M2xx=metric1(:,1); M2xy=metric1(:,2); M2yy=metric1(:,3);
+M2lambda1=0.5*((M2xx+M2yy)+sqrt(4*M2xy.^2+(M2xx-M2yy).^2));
+M2lambda2=0.5*((M2xx+M2yy)-sqrt(4*M2xy.^2+(M2xx-M2yy).^2));
+lambda2=min(M2lambda2,M2lambda2);
+
+metric=metric1;
+pos=find(lambda2<lambda1);
+metric(pos,:)=metric2(pos,:);