Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15475)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15476)
@@ -7029,7 +7029,4 @@
 	/*Transform Coordinate System*/
 	TransformStiffnessMatrixCoord(Ke,nodes,NUMVERTICES,XYZPEnum);
-	if(id==24){
-		//printarray(Ke->values,24,24);
-	}
 
 	/*Clean up and return*/
@@ -7985,7 +7982,7 @@
 			for(p=0;p<6;p++){
 				for(j=0;j<4;j++){
-					pe->values[p*4+j] += gauss->weight*Jdet*tau*rho_ice*forcex*SW[p][j][0];
-					pe->values[p*4+j] += gauss->weight*Jdet*tau*rho_ice*forcey*SW[p][j][1];
-					pe->values[p*4+j] += gauss->weight*Jdet*tau*rho_ice*(forcez-gravity)*SW[p][j][2];
+					pe->values[p*4+j] += gauss->weight*Jdet*tau*forcex*SW[p][j][0];
+					pe->values[p*4+j] += gauss->weight*Jdet*tau*forcey*SW[p][j][1];
+					pe->values[p*4+j] += gauss->weight*Jdet*tau*(forcez-rho_ice*gravity)*SW[p][j][2];
 				}
 			}
Index: /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp	(revision 15475)
+++ /issm/trunk-jpl/src/c/classes/Elements/PentaRef.cpp	(revision 15476)
@@ -509,5 +509,5 @@
 		B_prime[(NDOF4*NUMNODESP1)*5+NDOF4*i+3] = 0.;
 		B_prime[(NDOF4*NUMNODESP1)*6+NDOF4*i+3] = 0.;
-		B_prime[(NDOF4*NUMNODESP1)*7+NDOF4*i+3] = basis[i];
+		B_prime[(NDOF4*NUMNODESP1)*7+NDOF4*i+3] = - basis[i];
 	}