Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 9599)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 9600)
@@ -1512,5 +1512,5 @@
 	double     dux,duy;
 	double     epsvel=2.220446049250313e-16;
-	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
+	double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
 	double     scalex=0,scaley=0,scale=0,S=0;
 	double     vx,vy,vxobs,vyobs,weight;
@@ -1614,7 +1614,11 @@
 						velocity_mag    =sqrt(pow(vx,   2.)+pow(vy,   2.))+epsvel;
 						obs_velocity_mag=sqrt(pow(vxobs,2.)+pow(vyobs,2.))+epsvel;
+						printf("-------------- file: Tria.cpp line: %i\n",__LINE__); 
+						printf("%g %g\n",velocity_mag,obs_velocity_mag);
 						scale=-8*pow(meanvel,2.)/pow(velocity_mag,2.)*log(velocity_mag/obs_velocity_mag);
+						printf("%g\n",meanvel);
 						dux=scale*vx;
 						duy=scale*vy;
+						printf("%g %g %g\n",vx,vy,scale);
 						pe->values[i*NDOF2+0]+=dux*weight*Jdet*gauss->weight*basis[i]; 
 						pe->values[i*NDOF2+1]+=duy*weight*Jdet*gauss->weight*basis[i]; 
@@ -1688,5 +1692,5 @@
 	double     dux,duy;
 	double     epsvel=2.220446049250313e-16;
-	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
+	double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
 	double     scalex=0,scaley=0,scale=0,S=0;
 	double     vx,vy,vxobs,vyobs,weight;
@@ -4858,5 +4862,5 @@
 	double     misfit,Jdet;
 	double     epsvel=2.220446049250313e-16;
-	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
+	double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
 	double     velocity_mag,obs_velocity_mag;
 	double     xyz_list[NUMVERTICES][3];
@@ -4923,5 +4927,5 @@
 	double     Jelem=0, S=0;
 	double     epsvel=2.220446049250313e-16;
-	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
+	double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
 	double     misfit, Jdet;
 	double     vx,vy,vxobs,vyobs,weight;
@@ -5013,5 +5017,5 @@
 	double     misfit,Jdet;
 	double     epsvel=2.220446049250313e-16;
-	double     meanvel=1000/(365*24*3600); /*1000 m/yr*/
+	double     meanvel=3.170979198376458e-05; /*1000 m/yr*/
 	double     vx,vy,vxobs,vyobs,weight;
 	double     xyz_list[NUMVERTICES][3];