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

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-09-25 12:22:38 -0700 (Mon, 25 Sep 2006)
New Revision: 4610

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/preshape2d.f
Log:
Fixed dimensions of 2D shape function arrays.


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/preshape2d.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/preshape2d.f	2006-09-25 19:21:06 UTC (rev 4609)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/preshape2d.f	2006-09-25 19:22:38 UTC (rev 4610)
@@ -46,8 +46,8 @@
 c       nsgauss                 = number of element face gauss points
 c
 c     Output:
-c       sh(nsd+1,nsnodes,nsgauss)   = shape functions and their derivatives
-c       gauss(nsd+1,nsgauss)        = Gauss point coordinates and weights
+c       sh(nsd,nsnodes,nsgauss)   = shape functions and their derivatives
+c       gauss(nsd,nsgauss)        = Gauss point coordinates and weights
 c
       include "implicit.inc"
 c
@@ -60,8 +60,8 @@
 c...  subroutine arguments
 c
       integer intord,ietype,nsnodes,nsgauss,ierr
-      double precision sh(nsd+1,nsnodes,nsgauss)
-      double precision gauss(nsd+1,nsgauss)
+      double precision sh(nsd,nsnodes,nsgauss)
+      double precision gauss(nsd,nsgauss)
       character*(*) errstrng
 c
 c...  local constants
@@ -77,7 +77,6 @@
 c...  local variables
 c
       integer i,l,k,n,ind,nshsize,ngssize
-      double precision shtmp(2160)
       double precision rr,ss,drr,dss
       integer io(3)
 c
@@ -86,14 +85,13 @@
 c...  definitions
 c
       ierr=izero
-      ngssize=(nsd+1)*nsgauss
+      ngssize=nsd*nsgauss
       nshsize=ngssize*nsnodes
 c
 c...  initialize arrays
 c
       call fill(gauss,zero,ngssize)
       call fill(sh,zero,nshsize)
-      call fill(shtmp,zero,2160)
 c
 c... First type:  linear hex
 c