Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 26441)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 26442)
@@ -2915,6 +2915,6 @@
 	Input* thickness_input= element->GetInput(ThicknessEnum); _assert_(thickness_input);
 	Input* surface_input  = element->GetInput(SurfaceEnum);   _assert_(surface_input);
-	Input* vx_input       = element->GetInput(VxEnum);        _assert_(vx_input); //vertically integrated vx
-	Input* vy_input       = element->GetInput(VyEnum);        _assert_(vy_input); //vertically integrated vy
+	//Input* vx_input       = element->GetInput(VxEnum);        _assert_(vx_input); //vertically integrated vx
+	//Input* vy_input       = element->GetInput(VyEnum);        _assert_(vy_input); //vertically integrated vy
 	Input* vxbase_input   = element->GetInput(VxBaseEnum);    _assert_(vxbase_input); 
 	Input* vybase_input   = element->GetInput(VyBaseEnum);    _assert_(vybase_input);
@@ -3144,4 +3144,8 @@
 		ice_pressure_sh   = gravity*rho_ice*thickness*thickness*(n+1)/(2*(n+3));
 		pressure_sh = ice_pressure_sh + water_pressure_sh;
+		if (b>0) {
+			pressure = ice_pressure;
+			pressure_sh = ice_pressure_sh;
+		}
 
 		for (int i=0;i<numnodes;i++){
@@ -3276,4 +3280,7 @@
 		vx[i]=vbx[i]+vshx[i]*(n[i]+1)/(n[i]+2);
 		vy[i]=vby[i]+vshy[i]*(n[i]+1)/(n[i]+2);
+		// HOTFIX: set surface velocity to VxEnum and VyEnum
+	//	vx[i]=vbx[i]+vshx[i];
+	//	vy[i]=vby[i]+vshy[i];
 	}