Index: /issm/trunk/src/c/objects/Loads/Icefront.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 4683)
+++ /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 4684)
@@ -654,4 +654,5 @@
 	Node**   nodes=NULL;
 	Element* element=NULL;
+	Penta*   penta=NULL;
 
 	/*quad grids: */
@@ -668,4 +669,5 @@
 	double v35[3];
 	double v45[3];
+	bool isstokes;
 
 	/*Recover hook objects: */
@@ -673,7 +675,12 @@
 	element=(Element*)helement->delivers();
 	nodes=(Node**)hnodes->deliverp();
-	
+
 	/*check icefront is associated to a pentaelem: */
 	if(element->Enum()!=PentaEnum) ISSMERROR("Only Penta supported yet");
+	penta=(Penta*)element;
+
+	/*If not Stokes, return*/
+	penta->inputs->GetParameterValue(&isstokes,IsStokesEnum);
+	if(!isstokes) return;
 
 	/* Get dof list and node coordinates: */
Index: /issm/trunk/src/c/objects/Loads/Pengrid.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Pengrid.cpp	(revision 4683)
+++ /issm/trunk/src/c/objects/Loads/Pengrid.cpp	(revision 4684)
@@ -558,4 +558,5 @@
 	double Ke[4][4]={0.0};
 	double penalty_offset;
+	bool isstokes;
 
 	/*pointers: */
@@ -569,4 +570,8 @@
 	node=(Node*)hnode->delivers();
 	penta=(Penta*)helement->delivers();
+
+	/*If not Stokes, return*/
+	penta->inputs->GetParameterValue(&isstokes,IsStokesEnum);
+	if(!isstokes) return;
 
 	//recover slope: */