Index: /issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp	(revision 25236)
+++ /issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp	(revision 25237)
@@ -95,17 +95,8 @@
 	/*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();