Changeset 20965

Timestamp:

07/20/16 16:29:14 (9 years ago)

Author:

felicity

Message:

DEL: removing old code not needed anymore

File:

• issm/trunk-jpl/src/c/classes/Elements/Element.cpp (modified) (2 diffs)

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

@@ -63 +63 @@
         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 +101 @@
                 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*/