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

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;