Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16960)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 16961)
@@ -3507,6 +3507,6 @@
 	IssmDouble LprimeFSSSA[4][numdof2dm];
 	IssmDouble DLFSSSA[4][4]={0.0};
-	IssmDouble Ke_drag_gaussian[numdof2dm][numdofs];
-	IssmDouble Ke_drag_gaussian2[numdof2d][numdof2dm];
+	IssmDouble Ke_drag[numdof2dm][numdofs];
+	IssmDouble Ke_drag2[numdof2d][numdof2dm];
 	IssmDouble *xyz_list      = NULL;
 	IssmDouble *xyz_list_tria = NULL;
@@ -3586,14 +3586,15 @@
 					&DLSSAFS[0][0],8,8,0,
 					&LprimeSSAFS[0][0],8,numdofs,0,
-					&Ke_drag_gaussian[0][0],0);
+					&Ke_drag[0][0],1);
 
 		TripleMultiply( &LFSSSA[0][0],4,numdof2d,1,
 					&DLFSSSA[0][0],4,4,0,
 					&LprimeFSSSA[0][0],4,numdof2dm,0,
-					&Ke_drag_gaussian2[0][0],0);
-		for(i=0;i<numdof2dm;i++) for(j=0;j<numdofs;j++) Ke->values[i*numdoftot+j+numdofm]+=Ke_drag_gaussian[i][j];
-		for(i=0;i<numdof2d;i++) for(j=0;j<numdof2dm;j++) Ke->values[(i+numdofm)*numdoftot+j]+=Ke_drag_gaussian2[i][j];
-	}
-
+					&Ke_drag2[0][0],1);
+	}
+
+	for(i=0;i<numdof2dm;i++) for(j=0;j<numdofs;j++) Ke->values[i*numdoftot+j+numdofm]+=Ke_drag[i][j];
+	for(i=0;i<numdof2d;i++) for(j=0;j<numdof2dm;j++) Ke->values[(i+numdofm)*numdoftot+j]+=Ke_drag2[i][j];
+		
 	/*Transform Coordinate System*/
 	element->TransformStiffnessMatrixCoord(Ke,node_list,numnodes,cs_list);
@@ -3917,5 +3918,5 @@
 		vzHO_input->GetInputDerivativeValue(&dw[0],xyz_list,gauss);
 
-		for(i=0;i<vnumnodes;i++){
+		for(i=0;i<6;i++){
 			pe->values[i*3+0]+=-Jdet*gauss->weight*viscosity*dw[0]*dbasis[2][i];
 			pe->values[i*3+1]+=-Jdet*gauss->weight*viscosity*dw[1]*dbasis[2][i];