Changeset 17604

Timestamp:

03/31/14 13:36:06 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: for FS, do not reinterpolate results on P1 inputs, keep the same interpolation

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp (modified) (1 diff)

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

@@ -4219 +4219 @@
 
         /*Add vx and vy as inputs to the tria element: */
-        if(results_on_nodes){
-                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());
-        }
-        else{
-                element->AddInput(VxEnum,vx,P1Enum);
-                element->AddInput(VyEnum,vy,P1Enum);
-                element->AddInput(VelEnum,vel,P1Enum);
-                element->AddInput(PressureEnum,pressure,P1Enum);
-                if(dim==3) element->AddInput(VzEnum,vz,P1Enum);
-        }
+        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());
 
         /*Free ressources:*/