[cig-commits] r6531 - short/3D/PyLith/branches/pylith-0.8/pylith3d/utils

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-04-09 11:03:28 -0700 (Mon, 09 Apr 2007)
New Revision: 6531

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/utils/prescomp.f
Log:
Fixed direction of normal (it was backwards).


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/utils/prescomp.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/utils/prescomp.f	2007-04-09 15:56:24 UTC (rev 6530)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/utils/prescomp.f	2007-04-09 18:03:28 UTC (rev 6531)
@@ -123,9 +123,9 @@
 c...  compute cross product to get normal.
 c     If nvert == 4, the average of 2 vectors is used.
 c
-      call cross(x(1,ivert(1)),x(1,ivert(2)),x(1,ivert(3)),pvec)
+      call cross(x(1,ivert(3)),x(1,ivert(2)),x(1,ivert(1)),pvec)
       if(nvert.eq.4) then
-        call cross(x(1,ivert(3)),x(1,ivert(4)),x(1,ivert(1)),v)
+        call cross(x(1,ivert(1)),x(1,ivert(4)),x(1,ivert(3)),v)
         do i=1,nsd
           pvec(i)=0.5d0*(pvec(i)+v(i))
         end do