Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15476)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 15477)
@@ -7982,7 +7982,7 @@
 			for(p=0;p<6;p++){
 				for(j=0;j<4;j++){
-					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];
+					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];
 				}
 			}