Index: ../trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 20964)
+++ ../trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 20965)
@@ -62,12 +62,6 @@
 	/*Allocate arrays*/
 	int numvertices = this->GetNumberOfVertices();
 	IssmDouble* lambdas = xNew<IssmDouble>(numvertices);
-	IssmDouble* vorticityx = xNew<IssmDouble>(numvertices);
-	IssmDouble* vorticityy = xNew<IssmDouble>(numvertices);
-	IssmDouble* vorticityz = xNew<IssmDouble>(numvertices);
-	IssmDouble* omega_corrx = xNew<IssmDouble>(numvertices);
-	IssmDouble* omega_corry = xNew<IssmDouble>(numvertices);
-	IssmDouble* omega_corrz = xNew<IssmDouble>(numvertices);
 
 	/* Start looping on the number of vertices: */
 	Gauss* gauss=this->NewGauss();
@@ -106,22 +100,10 @@
 		}
 		EstarStrainrateQuantities(&epso,&epsprime,vx,vy,vz,vmag,&dvx[0],&dvy[0],&dvz[0],&dvmag[0]);
 		lambdas[iv]=EstarLambdaS(epso,epsprime);
-		vorticityx[iv]=dvz[1]-dvy[2];
-		vorticityy[iv]=dvx[2]-dvz[0];
-		vorticityz[iv]=dvy[0]-dvx[1];
-		omega_corrx[iv] = -vorticityx[iv] + 2*(dvmag[1]*vz - dvmag[2]*vy)/vmag;
-		omega_corry[iv] = -vorticityy[iv] + 2*(dvmag[2]*vx - dvmag[0]*vz)/vmag;
-		omega_corrz[iv] = -vorticityz[iv] + 2*(dvmag[0]*vy - dvmag[1]*vx)/vmag;
 	}
 
 	/*Add Stress tensor components into inputs*/
 	this->AddInput(LambdaSEnum,lambdas,P1Enum);
-	this->AddInput(Outputdefinition21Enum,vorticityx,P1Enum);
-	this->AddInput(Outputdefinition22Enum,vorticityy,P1Enum);
-	this->AddInput(Outputdefinition23Enum,vorticityz,P1Enum);
-	this->AddInput(Outputdefinition31Enum,omega_corrx,P1Enum);
-	this->AddInput(Outputdefinition32Enum,omega_corry,P1Enum);
-	this->AddInput(Outputdefinition33Enum,omega_corrz,P1Enum);
 	
 	/*Clean up and return*/
 	delete gauss;