Changeset 16720 for issm/trunk-jpl/src/c/classes/Elements/Penta.cpp

Timestamp:

11/12/13 15:54:29 (11 years ago)

Author:

seroussi

Message:

NEW: adding mass transport in analyses

File:

• issm/trunk-jpl/src/c/classes/Elements/Penta.cpp (modified) (1 diff)

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

@@ -126 +126 @@
         _assert_(this->inputs);
         this->inputs->AddInput(new PentaInput(input_enum,values,interpolation_enum));
+}
+/*}}}*/
+/*FUNCTION Penta::AddBasalInput{{{*/
+void  Penta::AddBasalInput(int input_enum,IssmDouble* values, int interpolation_enum){
+
+        _assert_(this->inputs);
+        if(!IsOnBed()) return;
+        else{
+                if(interpolation_enum==P1Enum){
+                        int        i;
+                        IssmDouble extrudedvalues[NUMVERTICES];
+                        Penta*     penta=NULL;
+
+                        for(i=1;i<NUMVERTICES2D;i++){
+                                extrudedvalues[i]=values[i];
+                                extrudedvalues[i+NUMVERTICES2D]=values[i];
+                        }
+                        this->inputs->AddInput(new PentaInput(input_enum,&extrudedvalues[0],P1Enum));
+                        penta=this;
+                        for(;;){
+                                penta->inputs->AddInput(new PentaInput(input_enum,&extrudedvalues[0],P1Enum));
+                                if (penta->IsOnSurface()) break;
+                                penta=penta->GetUpperElement(); _assert_(penta->Id()!=this->id);
+                        }
+                }
+                else _error_("not implemented yet");
+        }
 }
 /*}}}*/