Index: /issm/trunk/src/c/modules/ModelProcessorx/DistributeNumDofs.cpp
===================================================================
--- /issm/trunk/src/c/modules/ModelProcessorx/DistributeNumDofs.cpp	(revision 5546)
+++ /issm/trunk/src/c/modules/ModelProcessorx/DistributeNumDofs.cpp	(revision 5547)
@@ -26,4 +26,7 @@
 		else if (vertices_type[0]==HutterApproximationEnum){
 			numdofs=2;
+		}
+		else if (vertices_type[0]==StokesApproximationEnum){
+			numdofs=4;
 		}
 		else if (vertices_type[0]==MacAyealPattynApproximationEnum){
Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5546)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5547)
@@ -774,4 +774,7 @@
 			CreateKMatrixDiagnosticPattyn( Kgg);
 		}
+		else if(approximation==StokesApproximationEnum){
+			CreateKMatrixDiagnosticStokes( Kgg);
+		}
 		else if(approximation==HutterApproximationEnum){
 			return;
@@ -854,4 +857,7 @@
 		else if(approximation==HutterApproximationEnum){
 			return;
+		}
+		else if(approximation==StokesApproximationEnum){
+			CreatePVectorDiagnosticStokes( pg);
 		}
 		else if(approximation==MacAyealPattynApproximationEnum){
Index: /issm/trunk/src/c/objects/Loads/Icefront.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 5546)
+++ /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 5547)
@@ -385,14 +385,14 @@
 	/*Branck on the type of icefront: */
 	if (type==SegmentIcefrontEnum){
-		CreatePVectorDiagnosticHorizSegment(pg);
+		CreatePVectorDiagnosticMacAyeal(pg);
 	}
 	else if (type==QuadIceFrontEnum){
-		CreatePVectorDiagnosticHorizQuad(pg);
+		CreatePVectorDiagnosticPattyn(pg);
 	}
 	else ISSMERROR("type %i not supported",type);
 }	
 /*}}}*/
-/*FUNCTION Icefront::CreatePVectorDiagnosticHorizSegment{{{1*/
-void Icefront::CreatePVectorDiagnosticHorizSegment( Vec pg){
+/*FUNCTION Icefront::CreatePVectorDiagnosticMacAyeal{{{1*/
+void Icefront::CreatePVectorDiagnosticMacAyeal( Vec pg){
 
 	int       i;
@@ -520,6 +520,6 @@
 }
 /*}}}*/
-/*FUNCTION Icefront::CreatePVectorDiagnosticHorizQuad {{{1*/
-void Icefront::CreatePVectorDiagnosticHorizQuad( Vec pg){
+/*FUNCTION Icefront::CreatePVectorDiagnosticPattyn{{{1*/
+void Icefront::CreatePVectorDiagnosticPattyn( Vec pg){
 
 	int i,j;
Index: /issm/trunk/src/c/objects/Loads/Icefront.h
===================================================================
--- /issm/trunk/src/c/objects/Loads/Icefront.h	(revision 5546)
+++ /issm/trunk/src/c/objects/Loads/Icefront.h	(revision 5547)
@@ -72,6 +72,6 @@
 		/*Load management: {{{1*/
 		void  CreatePVectorDiagnosticHoriz( Vec pg);
-		void  CreatePVectorDiagnosticHorizSegment( Vec pg);
-		void  CreatePVectorDiagnosticHorizQuad( Vec pg);
+		void  CreatePVectorDiagnosticMacAyeal( Vec pg);
+		void  CreatePVectorDiagnosticPattyn( Vec pg);
 		void  CreatePVectorDiagnosticStokes( Vec pg);
 		void  GetDofList(int** pdoflist,int approximation_enum=0);