<p><b>mperego@fsu.edu</b> 2011-09-13 10:49:32 -0600 (Tue, 13 Sep 2011)</p><p>in error_checking routine, vectors p, q, r, s were normalized (using rutine unit_vector_in_3space)  before computing t=q-p and s=s-r, which led to wrong computation of angles between mesh entities. Those calls to unit_vector_in_3space have been removed.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice/icesheet/src/icesheet.F
===================================================================
--- branches/land_ice/icesheet/src/icesheet.F        2011-09-09 16:24:02 UTC (rev 988)
+++ branches/land_ice/icesheet/src/icesheet.F        2011-09-13 16:49:32 UTC (rev 989)
@@ -285,9 +285,6 @@
   p(1)=xCellNew(iCell1); p(2)=yCellNew(iCell1); p(3)=zCellNew(iCell1)
   q(1)=xCellNew(iCell2); q(2)=yCellNew(iCell2); q(3)=zCellNew(iCell2)
   r(1)=xEdgeNew(iEdge);  r(2)=yEdgeNew(iEdge);  r(3)=zEdgeNew(iEdge)
-  call unit_vector_in_3space(p)
-  call unit_vector_in_3space(q)
-  call unit_vector_in_3space(r)
   t = q - p
   s = r - p
   call unit_vector_in_3space(t)
@@ -312,10 +309,6 @@
   q(1)=xCellNew(iCell2); q(2)=yCellNew(iCell2); q(3)=zCellNew(iCell2)
   r(1)=xVertexNew(iVertex1); r(2)=yVertexNew(iVertex1); r(3)=zVertexNew(iVertex1)
   s(1)=xVertexNew(iVertex2); s(2)=yVertexNew(iVertex2); s(3)=zVertexNew(iVertex2)
-  call unit_vector_in_3space(p)
-  call unit_vector_in_3space(q)
-  call unit_vector_in_3space(r)
-  call unit_vector_in_3space(s)
   t = q - p
   s = s - r
   call unit_vector_in_3space(t)

</font>
</pre>