[cig-commits] r6010 - short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-02-12 12:06:58 -0800 (Mon, 12 Feb 2007)
New Revision: 6010

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/getjac2d.f
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/traction_cmp_ss.f
Log:
Fixed 2D Jacobian determinant (needed square root).
Fixed problem where element contributions were replacing traction
vector contents rather than updating them.



Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/getjac2d.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/getjac2d.f	2007-02-12 20:05:19 UTC (rev 6009)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/getjac2d.f	2007-02-12 20:06:58 UTC (rev 6010)
@@ -80,6 +80,7 @@
         ierr=113
         write(errstrng,700) ientry,det
       end if
+      det=sqrt(det)
 c
  700  format("getjac2d:  entry # ",i7,2x,1pe15.8)
       return

Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/traction_cmp_ss.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/traction_cmp_ss.f	2007-02-12 20:05:19 UTC (rev 6009)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/traction_cmp_ss.f	2007-02-12 20:06:58 UTC (rev 6010)
@@ -105,7 +105,7 @@
       do i=1,nsnodes
         do j=1,ndof
           k=id(j,tractionverts(i))
-          btraction(k)=p(j,i)
+          if(k.gt.0) btraction(k)=btraction(k)+p(j,i)
         end do
       end do
       return