Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5522)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5523)
@@ -2206,13 +2206,4 @@
 	}
 
-	//if diagnostic Full-Stokes: check that the element is Stokes
-	if(analysis_type==DiagnosticStokesAnalysisEnum){
-		ISSMASSERT(iomodel->elements_type);
-		if(iomodel->elements_type[index]==StokesApproximationEnum)
-		 this->inputs->AddInput(new BoolInput(IsStokesEnum,true));
-		else
-		 this->inputs->AddInput(new BoolInput(IsStokesEnum,false));
-	}
-
 }
 /*}}}*/
@@ -2944,5 +2935,4 @@
 	bool onbed;
 	bool shelf;
-	bool isstokes;
 
 	/*inputs: */
@@ -2958,9 +2948,8 @@
 	inputs->GetParameterValue(&onbed,ElementOnBedEnum);
 	inputs->GetParameterValue(&shelf,ElementOnIceShelfEnum);
-	inputs->GetParameterValue(&isstokes,IsStokesEnum);
 	inputs->GetParameterValue(&approximation,ApproximationEnum);
 
 	/*If on water or not Stokes, skip stiffness: */
-	if(onwater || !isstokes || approximation!=StokesApproximationEnum) return;
+	if(onwater || approximation!=StokesApproximationEnum) return;
 
 	/*recovre material parameters: */
@@ -4131,5 +4120,4 @@
 	bool onbed;
 	bool shelf;
-	bool isstokes;
 	int  approximation;
 	Input* vx_input=NULL;
@@ -4142,5 +4130,4 @@
 	inputs->GetParameterValue(&onbed,ElementOnBedEnum);
 	inputs->GetParameterValue(&shelf,ElementOnIceShelfEnum);
-	inputs->GetParameterValue(&isstokes,IsStokesEnum);
 	inputs->GetParameterValue(&approximation,ApproximationEnum);
 
@@ -4149,5 +4136,5 @@
 
 	/*If on water or not Stokes, skip load: */
-	if(onwater || !isstokes || approximation!=StokesApproximationEnum) return;
+	if(onwater || approximation!=StokesApproximationEnum) return;
 
 	/*If on water, skip load: */
Index: /issm/trunk/src/c/objects/Loads/Icefront.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 5522)
+++ /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 5523)
@@ -697,5 +697,5 @@
 	double v35[3];
 	double v45[3];
-	bool isstokes;
+	int approximation;
 
 	/*Recover hook objects: */
@@ -709,6 +709,6 @@
 
 	/*If not Stokes, return*/
-	penta->inputs->GetParameterValue(&isstokes,IsStokesEnum);
-	if(!isstokes) return;
+	penta->inputs->GetParameterValue(&approximation,ApproximationEnum);
+	if(approximation!=StokesApproximationEnum) 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 5522)
+++ /issm/trunk/src/c/objects/Loads/Pengrid.cpp	(revision 5523)
@@ -574,5 +574,5 @@
 	double Ke[4][4]={0.0};
 	double penalty_offset;
-	bool isstokes;
+	int approximation;
 
 	/*pointers: */
@@ -588,6 +588,6 @@
 
 	/*If not Stokes, return*/
-	penta->inputs->GetParameterValue(&isstokes,IsStokesEnum);
-	if(!isstokes) return;
+	penta->inputs->GetParameterValue(&approximation,ApproximationEnum);
+	if(approximation!=StokesApproximationEnum) return;
 
 	//recover slope: */