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

Timestamp:

05/10/10 15:29:06 (15 years ago)

Author:

Eric.Larour

Message:

Added vx,vy,vz, vx_obs,vy_obs,vz_obs,weights in elements

File:

• issm/trunk/src/c/objects/Elements/Penta.cpp (modified) (1 diff)

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

@@ -118 +118 @@
                 this->inputs->AddInput(new PentaVertexInput(GeothermalFluxEnum,nodeinputs));
         }       
+
+        /*vx,vy and vz: */
+        if (iomodel->vx) {
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->vx[penta_node_ids[i]-1];
+                this->inputs->AddInput(new PentaVertexInput(VxEnum,nodeinputs));
+        }
+        if (iomodel->vy) {
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->vy[penta_node_ids[i]-1];
+                this->inputs->AddInput(new PentaVertexInput(VyEnum,nodeinputs));
+        }
+        if (iomodel->vz) {
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->vz[penta_node_ids[i]-1];
+                this->inputs->AddInput(new PentaVertexInput(VzEnum,nodeinputs));
+        }
+        if (iomodel->vx_obs) {
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->vx_obs[penta_node_ids[i]-1];
+                this->inputs->AddInput(new PentaVertexInput(VxObsEnum,nodeinputs));
+        }
+        if (iomodel->vy_obs) {
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->vy_obs[penta_node_ids[i]-1];
+                this->inputs->AddInput(new PentaVertexInput(VyObsEnum,nodeinputs));
+        }
+        if (iomodel->vz_obs) {
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->vz_obs[penta_node_ids[i]-1];
+                this->inputs->AddInput(new PentaVertexInput(VzObsEnum,nodeinputs));
+        }
+        if (iomodel->weights) {
+                for(i=0;i<6;i++)nodeinputs[i]=iomodel->weights[penta_node_ids[i]-1];
+                this->inputs->AddInput(new PentaVertexInput(WeightsEnum,nodeinputs));
+        }
 
         if (iomodel->elementoniceshelf) this->inputs->AddInput(new BoolInput(ElementOnIceShelfEnum,(IssmBool)iomodel->elementoniceshelf[index]));