Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 20964)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 20965)
@@ -63,10 +63,4 @@
 	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: */
@@ -107,20 +101,8 @@
 		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*/