Index: /issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp	(revision 26474)
+++ /issm/trunk-jpl/src/c/analyses/AdjointHorizAnalysis.cpp	(revision 26475)
@@ -510,5 +510,5 @@
 					for(i=0;i<vnumnodes;i++){
 						if (domaintype!=Domain2DverticalEnum){
-							scale=1./(S*2*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
+							scale=1./(S*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
 							dux=scale*(vxobs-vx);
 							duy=scale*(vyobs-vy);
@@ -517,5 +517,5 @@
 						}
 						else{
-							scale=1./(S*2*fabs(vx-vxobs)+epsvel);
+							scale=1./(S*fabs(vx-vxobs)+epsvel);
 							dux=scale*(vxobs-vx);
 							pe->values[i*2+0]+=dux*weight*Jdet*gauss->weight*vbasis[i]; 
@@ -746,5 +746,5 @@
 					for(i=0;i<numnodes;i++){
 						if(domaintype!=Domain2DverticalEnum){
-							scale=1./(S*2*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
+							scale=1./(S*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
 							dux=scale*(vxobs-vx);
 							duy=scale*(vyobs-vy);
@@ -998,5 +998,5 @@
 					for(i=0;i<numnodes;i++){
 						if(domaintype!=Domain2DverticalEnum){
-							scale=1./(S*2*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
+							scale=1./(S*sqrt(pow(vx-vxobs,2)+pow(vy-vyobs,2))+epsvel);
 							dux=scale*(vxobs-vx);
 							duy=scale*(vyobs-vy);
@@ -1005,5 +1005,5 @@
 						}
 						else{
-							scale=1./(S*2*fabs(vx-vxobs)+epsvel);
+							scale=1./(S*fabs(vx-vxobs)+epsvel);
 							dux=scale*(vxobs-vx);
 							pe->values[i]+=dux*weight*Jdet*gauss->weight*basis[i];