Index: /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp	(revision 17180)
+++ /issm/trunk-jpl/src/c/analyses/FreeSurfaceBaseAnalysis.cpp	(revision 17181)
@@ -169,8 +169,8 @@
 
 		D_scalar=dt*gauss->weight*Jdet;
-
 		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,
Index: /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp	(revision 17180)
+++ /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp	(revision 17181)
@@ -178,9 +178,8 @@
 
 		D_scalar=dt*gauss->weight*Jdet;
-
 		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,