Index: /issm/trunk/src/c/objects/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Tria.cpp	(revision 397)
+++ /issm/trunk/src/c/objects/Tria.cpp	(revision 398)
@@ -1674,19 +1674,4 @@
 		gauss_l1l2l3[2]=*(third_gauss_area_coord+ig);
 
-		
-		/*Compute velocities at gaussian point: */
-		GetParameterValue(&velocity_x, &vx_list[0],gauss_l1l2l3);
-		GetParameterValue(&velocity_y, &vy_list[0],gauss_l1l2l3);
-		#ifdef _DEBUG_ 
-			printf("Velocity: %lf %lf\n", velocity_x,velocity_y);
-		#endif
-	
-		/*Compute obs_velocities at gaussian point: */
-		GetParameterValue(&obs_velocity_x, &obs_vx_list[0],gauss_l1l2l3);
-		GetParameterValue(&obs_velocity_y, &obs_vy_list[0],gauss_l1l2l3);
-		#ifdef _DEBUG_ 
-			printf("Observed velocity: %g %g\n", obs_velocity_x,obs_velocity_y);
-		#endif
-
 		/* Get Jacobian determinant: */
 		GetJacobianDeterminant2d(&Jdet, &xyz_list[0][0],gauss_l1l2l3);
@@ -1948,5 +1933,5 @@
 		/*Build gradje_g_gaussian vector (actually -dJ/ddrag): */
 		for (i=0;i<numgrids;i++){
-			grade_g_gaussian[i]=2*drag*alpha_complement*( (lambda*vx+mu*vy))*Jdet*gauss_weight*l1l2l3[i]; 
+			grade_g_gaussian[i]=-2*drag*alpha_complement*((lambda*vx+mu*vy))*Jdet*gauss_weight*l1l2l3[i]; 
 		}