Index: /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp	(revision 17179)
+++ /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp	(revision 17180)
@@ -173,11 +173,4 @@
 		D[0] = D_scalar*vx;
 		if(dim==2) D[1*dim+1] = D_scalar*vy;
-		TripleMultiply(B,dim,numnodes,1,
-					D,dim,dim,0,
-					B,dim,numnodes,0,
-					&Ke->values[0],1);
-
-		D[0*dim+0]=D_scalar*vx;
-		if(dim==2) D[1*dim+1]=D_scalar*vy;
 		TripleMultiply(B,dim,numnodes,1,
 					D,dim,dim,0,
Index: /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp	(revision 17179)
+++ /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp	(revision 17180)
@@ -181,11 +181,5 @@
 		for(int i=0;i<dim*dim;i++) D[i]=0.;
 		D[0] = D_scalar*vx;
-		if(dim==2) D[1*dim+1] = D_scalar*vy;
-		TripleMultiply(B,dim,numnodes,1,
-					D,dim,dim,0,
-					B,dim,numnodes,0,
-					&Ke->values[0],1);
-
-		D[0*dim+0]=D_scalar*vx;
+
 		if(dim==2) D[1*dim+1]=D_scalar*vy;
 		TripleMultiply(B,dim,numnodes,1,