Index: /issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp	(revision 18053)
+++ /issm/trunk-jpl/src/c/analyses/AdjointBalancethickness2Analysis.cpp	(revision 18054)
@@ -49,5 +49,5 @@
 	int         num_responses,i;
 	IssmDouble  hobs,hu2,weight,NUMx,NUMy,DEN,Jdet;
-	IssmDouble  vx,vy,vbar2,nux,nuy,dphi[2];
+	IssmDouble  vx,vy,vbar2,nux,nuy,phi,dphi[2];
 	int        *responses = NULL;
 	IssmDouble *xyz_list  = NULL;
@@ -86,4 +86,5 @@
 		nux_input->GetInputValue(&nux,gauss);
 		nuy_input->GetInputValue(&nuy,gauss);
+		potential_input->GetInputValue(&phi,gauss);
 		potential_input->GetInputDerivativeValue(&dphi[0],xyz_list,gauss);
 		thicknessobs_input->GetInputValue(&hobs,gauss);
@@ -102,5 +103,6 @@
 			switch(responses[resp]){
 				case Balancethickness2MisfitEnum:
-					for(i=0;i<numnodes;i++) pe->values[i]+=(NUMx+NUMy)/DEN *weight*Jdet*gauss->weight;
+					//for(i=0;i<numnodes;i++) pe->values[i]+=(NUMx+NUMy)/DEN *weight*Jdet*gauss->weight;
+					for(i=0;i<numnodes;i++) pe->values[i]+= phi*basis[i]*weight*Jdet*gauss->weight;
 					break;
 				default: