<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>