Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5669)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5670)
@@ -2897,6 +2897,6 @@
 
 			/*Get L matrix if viscous basal drag present: */
-			GetLStokes(&LStokes[0][0],  gauss_coord_tria);
-			GetLprimeStokes(&LprimeStokes[0][0], &xyz_list[0][0], gauss_coord_tria, gauss_coord);
+			GetLStokes(&LStokes[0][0],  gauss_coord);
+			GetLprimeStokes(&LprimeStokes[0][0], &xyz_list[0][0], gauss_coord);
 
 			/*Compute strain rate: */
Index: /issm/trunk/src/c/objects/Elements/PentaRef.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/PentaRef.cpp	(revision 5669)
+++ /issm/trunk/src/c/objects/Elements/PentaRef.cpp	(revision 5670)
@@ -477,5 +477,5 @@
 /*}}}*/
 /*FUNCTION PentaRef::GetLStokes {{{1*/
-void PentaRef::GetLStokes(double* LStokes, double* gauss_tria){
+void PentaRef::GetLStokes(double* LStokes, double* gauss){
 	/*
 	 * Compute L  matrix. L=[L1 L2 L3] where Li is square and of size numdof. 
@@ -507,7 +507,7 @@
 
 	/*Get l1l2l3 in actual coordinate system: */
-	l1l2l3[0]=gauss_tria[0];
-	l1l2l3[1]=gauss_tria[1];
-	l1l2l3[2]=gauss_tria[2];
+	l1l2l3[0]=gauss[0];
+	l1l2l3[1]=gauss[1];
+	l1l2l3[2]=gauss[2];
 
 	/*Build LStokes: */
@@ -574,5 +574,5 @@
 /*}}}*/
 /*FUNCTION PentaRef::GetLprimeStokes {{{1*/
-void PentaRef::GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss_tria, double* gauss){
+void PentaRef::GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss){
 
 	/*
@@ -604,7 +604,7 @@
 
 	/*Get l1l2l3 in actual coordinate system: */
-	l1l2l3[0]=gauss_tria[0];
-	l1l2l3[1]=gauss_tria[1];
-	l1l2l3[2]=gauss_tria[2];
+	l1l2l3[0]=gauss[0];
+	l1l2l3[1]=gauss[1];
+	l1l2l3[2]=gauss[2];
 
 	GetNodalFunctionsP1Derivatives(&dh1dh6[0][0],xyz_list,gauss);
Index: /issm/trunk/src/c/objects/Elements/PentaRef.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/PentaRef.h	(revision 5669)
+++ /issm/trunk/src/c/objects/Elements/PentaRef.h	(revision 5670)
@@ -43,6 +43,6 @@
 		void GetBVert(double* B, double* xyz_list, double* gauss);
 		void GetBprimeAdvec(double* Bprime_advec, double* xyz_list, double* gauss);
-		void GetLStokes(double* LStokes, double* gauss_tria);
-		void GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss_tria, double* gauss);
+		void GetLStokes(double* LStokes, double* gauss);
+		void GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss);
 		void GetParameterValue(double* pvalue,double* plist,double* gauss);
 		void GetParameterValue(double* pvalue,double* plist,GaussTria* gauss){ISSMERROR("only PentaGauss are supported");};