Changeset 17096

Timestamp:

01/10/14 14:10:18 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: keep solution on the same finite element (commented out for now)

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp ¶

@@ -2409 +2409 @@
         element->InputChangeName(PressureEnum,PressurePicardEnum);
 
-        /*Add vx and vy as inputs to the tria element: */
+        /*Add vx and vy as inputs to the element: */
+        //element->AddInput(VxEnum,  vx,element->GetElementType());
+        //element->AddInput(VyEnum,  vy,element->GetElementType());
+        //element->AddInput(VelEnum,vel,element->GetElementType());
         element->AddInput(VxEnum,vx,P1Enum);
         element->AddInput(VyEnum,vy,P1Enum);
-        element->AddBasalInput(VelEnum,vel,P1Enum);
+        element->AddInput(VelEnum,vel,P1Enum);
 
         /*Free ressources:*/
@@ -3374 +3377 @@
 
         /*Add vx and vy as inputs to the tria element: */
+        //element->AddInput(VxEnum,      vx,      element->VelocityInterpolation());
+        //element->AddInput(VyEnum,      vy,      element->VelocityInterpolation());
+        //element->AddInput(VelEnum,     vel,     element->VelocityInterpolation());
+        //element->AddInput(PressureEnum,pressure,element->PressureInterpolation());
+        //if(dim==3) element->AddInput(VzEnum,vz, element->VelocityInterpolation());
         element->AddInput(VxEnum,vx,P1Enum);
         element->AddInput(VyEnum,vy,P1Enum);