Changeset 4797 for issm/trunk/src/c/objects/Elements/Penta.cpp

Timestamp:

07/26/10 07:48:41 (15 years ago)

Author:

Mathieu Morlighem

Message:

new vel must be computed after Stokes

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (4 diffs)

issm/trunk/src/c/objects/Elements/Penta.cpp

@@ -2375 +2375 @@
         double       vy[numvertices];
         double       vz[numvertices];
+        double       vel[numvertices];
         double       pressure[numvertices];
         int          dummy;
@@ -2387 +2388 @@
         }
 
-        /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
+        /*Ok, we have vx and vy in values, fill in all arrays: */
         for(i=0;i<numvertices;i++){
                 vx[i]=values[i*numdofpervertex+0];
@@ -2400 +2401 @@
                 pressure[i]=pressure[i]*stokesreconditioning;
         }
+
+        /*Now Compute vel*/
+        for(i=0;i<numvertices;i++) vel[i]=pow( pow(vx[i],2.0) + pow(vy[i],2.0) + pow(vz[i],2.0) , 0.5);
         
         /*Now, we have to move the previous inputs  to old 
@@ -2412 +2416 @@
         this->inputs->AddInput(new PentaVertexInput(VyEnum,vy));
         this->inputs->AddInput(new PentaVertexInput(VzEnum,vz));
+        this->inputs->AddInput(new PentaVertexInput(VelEnum,vel));
         this->inputs->AddInput(new PentaVertexInput(PressureEnum,pressure));
 }