Changeset 25237

Timestamp:

07/09/20 09:30:37 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: we don't need to use VxAverage under ice shelves for GL distance

File:

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

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

@@ -95 +95 @@
         /*Retrieve all inputs and parameters*/
         basalelement->GetVerticesCoordinates(&xyz_list);
-        switch(domaintype){
-                case Domain2DhorizontalEnum:
-                        vx_input=basalelement->GetInput2(VxEnum); _assert_(vx_input);
-                        vy_input=basalelement->GetInput2(VyEnum); _assert_(vy_input);
-                   bc_input=basalelement->GetInput2(MeshVertexonboundaryEnum); _assert_(bc_input);
-                break;
-                case Domain3DEnum:
-                        vx_input=basalelement->GetInput2(VxAverageEnum); _assert_(vx_input);
-                        vy_input=basalelement->GetInput2(VyAverageEnum); _assert_(vy_input);
-                        bc_input=basalelement->GetInput2(MeshVertexonboundaryEnum); _assert_(bc_input);
-                break;
-                default: _error_("mesh "<<EnumToStringx(domaintype)<<" not supported yet");
-        }
+        vx_input=basalelement->GetInput2(VxEnum); _assert_(vx_input);
+        vy_input=basalelement->GetInput2(VyEnum); _assert_(vy_input);
+        bc_input=basalelement->GetInput2(MeshVertexonboundaryEnum); _assert_(bc_input);
+
         IssmDouble h = basalelement->CharacteristicLength();