Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4796)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 4797)
@@ -2375,4 +2375,5 @@
 	double       vy[numvertices];
 	double       vz[numvertices];
+	double       vel[numvertices];
 	double       pressure[numvertices];
 	int          dummy;
@@ -2387,5 +2388,5 @@
 	}
 
-	/*Ok, we have vx and vy in values, fill in vx and vy arrays: */
+	/*Ok, we have vx and vy in values, fill in all arrays: */
 	for(i=0;i<numvertices;i++){
 		vx[i]=values[i*numdofpervertex+0];
@@ -2400,4 +2401,7 @@
 		pressure[i]=pressure[i]*stokesreconditioning;
 	}
+
+	/*Now Compute vel*/
+	for(i=0;i<numvertices;i++) vel[i]=pow( pow(vx[i],2.0) + pow(vy[i],2.0) + pow(vz[i],2.0) , 0.5);
 	
 	/*Now, we have to move the previous inputs  to old 
@@ -2412,4 +2416,5 @@
 	this->inputs->AddInput(new PentaVertexInput(VyEnum,vy));
 	this->inputs->AddInput(new PentaVertexInput(VzEnum,vz));
+	this->inputs->AddInput(new PentaVertexInput(VelEnum,vel));
 	this->inputs->AddInput(new PentaVertexInput(PressureEnum,pressure));
 }