Index: /issm/trunk/src/c/objects/Pengrid.cpp
===================================================================
--- /issm/trunk/src/c/objects/Pengrid.cpp	(revision 460)
+++ /issm/trunk/src/c/objects/Pengrid.cpp	(revision 461)
@@ -197,6 +197,6 @@
 	
 	const int numgrids=1;
-	const int NDOF3=3;
-	const int numdof=numgrids*NDOF3;
+	const int NDOF4=4;
+	const int numdof=numgrids*NDOF4;
 	int       doflist[numdof];
 	int       numberofdofspernode;
@@ -223,6 +223,6 @@
 
 	//Create elementary matrix: add penalty to contrain wb (wb=ub*db/dx+vb*db/dy)
-	Ke[0][0]=-slope[0]*kmax*pow(10.0,penalty_offset);
-	Ke[1][1]=-slope[1]*kmax*pow(10.0,penalty_offset);
+	Ke[2][0]=-slope[0]*kmax*pow(10.0,penalty_offset);
+	Ke[2][1]=-slope[1]*kmax*pow(10.0,penalty_offset);
 	Ke[2][2]=kmax*pow(10,penalty_offset);