Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3869)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3870)
@@ -1523,5 +1523,4 @@
 	GetVerticesCoordinates(&xyz_list[0][0], nodes, numgrids);
 	GetDofList(&doflist[0],&numberofdofspernode);
-
 
 	/*Get gaussian points and weights. Penta is an extrusion of a Tria, we therefore 
@@ -2422,5 +2421,4 @@
 	int i;
 
-
 	/* node data: */
 	const int    numgrids=6;
@@ -2473,5 +2471,4 @@
 	inputs->GetParameterValue(&onbed,ElementOnBedEnum);
 
-	
 	/*If on water, skip: */
 	if(onwater)return;
@@ -2513,8 +2510,7 @@
 
 			inputs->GetParameterDerivativeValue(&du[0],&xyz_list[0][0],gauss_coord,VxEnum);
-			inputs->GetParameterDerivativeValue(&du[0],&xyz_list[0][0],gauss_coord,VyEnum);
+			inputs->GetParameterDerivativeValue(&dv[0],&xyz_list[0][0],gauss_coord,VyEnum);
 			dudx=du[0];
 			dvdy=dv[1];
-
 
 			/* Get Jacobian determinant: */
@@ -2609,5 +2605,4 @@
 /*FUNCTION Penta::CreatePVectorThermal {{{1*/
 void Penta::CreatePVectorThermal( Vec pg, int analysis_type,int sub_analysis_type){
-
 
 	/*indexing: */