Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16350)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16351)
@@ -2063,4 +2063,5 @@
 				name==EplHeadEnum ||
 				name==HydrologydcMaskEplactiveEnum ||
+				name==MeshVertexonbedEnum ||
 				name==WaterTransferEnum ||
 				name==QmuVxEnum ||
@@ -2113,22 +2114,19 @@
 void Tria::EdgeOnBedIndices(int* pindex1,int* pindex2){
 
-	bool found=false;
 	IssmDouble values[NUMVERTICES];
+	int        indices[3][2] = {{1,2},{2,0},{0,1}};
 
 	/*Retrieve all inputs and parameters*/
 	GetInputListOnVertices(&values[0],MeshVertexonbedEnum);
 
-	for(int i=0;i<NUMVERTICES;i++){
-		if(values[i]==1.){
-			if(found){
-				*pindex2 = i;
-				return;
-			}
-			else{
-				*pindex1 = i;
-			}
-		}
-	}
-
+	for(int i=0;i<3;i++){
+		if(values[indices[i][0]] == 1. && values[indices[i][1]] == 1.){
+			*pindex1 = indices[i][0];
+			*pindex2 = indices[i][1];
+			return;
+		}
+	}
+
+	_printf_("list of vertices on bed: "<<values[0]<<" "<<values[1]<<" "<<values[2]);
 	_error_("Could not find 2 vertices on bed");
 }
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp	(revision 16350)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp	(revision 16351)
@@ -114,4 +114,5 @@
 	}
 	if(isFS){
+		iomodel->FetchDataToInput(elements,MeshVertexonbedEnum);
 		iomodel->FetchDataToInput(elements,PressureEnum,0.);
 		if(dakota_analysis)elements->InputDuplicate(PressureEnum,QmuPressureEnum);