source: issm/oecreview/Archive/20545-21336/ISSM-20964-20965.diff@ 21337

../trunk-jpl/src/c/classes/Elements/Element.cpp

@@ -62 +62 @@
         /*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 +100 @@
                 }
                 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;