Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5805)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5806)
@@ -1866,4 +1866,18 @@
 				}
 			}
+			if(*(iomodel->elements_type+index)==PattynStokesApproximationEnum){
+				/*Create VzPattyn and VzStokes Enums*/
+				if(iomodel->vz && iomodel->gridonstokes){
+					for(i=0;i<6;i++) nodeinputs[i]=iomodel->vz[penta_vertex_ids[i]-1]/iomodel->yts*iomodel->gridonstokes[penta_vertex_ids[i]-1];
+					this->inputs->AddInput(new PentaVertexInput(VzStokesEnum,nodeinputs));
+					for(i=0;i<6;i++) nodeinputs[i]=iomodel->vz[penta_vertex_ids[i]-1]/iomodel->yts*(1-iomodel->gridonstokes[penta_vertex_ids[i]-1]);
+					this->inputs->AddInput(new PentaVertexInput(VzPattynEnum,nodeinputs));
+				}
+				else{
+					for(i=0;i<6;i++)nodeinputs[i]=0;
+					this->inputs->AddInput(new PentaVertexInput(VzStokesEnum,nodeinputs));
+					this->inputs->AddInput(new PentaVertexInput(VzPattynEnum,nodeinputs));
+				}
+			}
 			break;
 
@@ -1890,4 +1904,7 @@
 		else if (*(iomodel->elements_type+index)==StokesApproximationEnum){
 			this->inputs->AddInput(new IntInput(ApproximationEnum,StokesApproximationEnum));
+		}
+		else if (*(iomodel->elements_type+index)==PattynStokesApproximationEnum){
+			this->inputs->AddInput(new IntInput(ApproximationEnum,PattynStokesApproximationEnum));
 		}
 		else if (*(iomodel->elements_type+index)==NoneApproximationEnum){
@@ -4921,4 +4938,10 @@
 	int         approximation;
 
+	/*Get the approximation and do nothing if the element in Stokes or None*/
+	inputs->GetParameterValue(&approximation,ApproximationEnum);
+	if(approximation==StokesApproximationEnum || approximation==NoneApproximationEnum){
+		return;
+	}
+
 	/*Get dof list: */
 	GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
@@ -4926,5 +4949,4 @@
 	/*Get node data: */
 	GetVerticesCoordinates(&xyz_list[0][0], nodes, NUMVERTICES);
-	inputs->GetParameterValue(&approximation,ApproximationEnum);
 
 	/*Use the dof list to index into the solution vector: */