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

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-10-26 12:53:25 -0700 (Thu, 26 Oct 2006)
New Revision: 5102

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/read_bc.f
Log:
Fixed routine so BC arrays are zeroed out before returning even when
there are no BC entries.



Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/read_bc.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/read_bc.f	2006-10-26 19:43:52 UTC (rev 5101)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/read_bc.f	2006-10-26 19:53:25 UTC (rev 5102)
@@ -84,6 +84,8 @@
 c
       ierr=izero
       nconcforce=izero
+      call fill(bond,zero,ndof*numnp)
+      call ifill(ibond,izero,ndof*numnp)
 c
 c...  return if no BC defined for this partition
 c
@@ -111,8 +113,6 @@
 c
 c...  read BC, returning an error code if node is out of range
 c
-      call fill(bond,zero,ndof*numnp)
-      call ifill(ibond,izero,ndof*numnp)
       call pskip(kr)
       do i=1,numbc
         read(kr,*,end=30,err=30) n,(ibond(j,n),j=1,ndof),