Index: /issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp	(revision 23036)
+++ /issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp	(revision 23037)
@@ -87,5 +87,5 @@
 	//Input* vx_input=element->GetInput(BaseSlopeXEnum); _assert_(vx_input);
 	//Input* vy_input=element->GetInput(BaseSlopeYEnum); _assert_(vy_input);
-	Input* vx_input=element->GetInput(VxEnum); _assert_(vx_input);
+   Input* vx_input=element->GetInput(VxEnum); _assert_(vx_input);
 	Input* vy_input=element->GetInput(VyEnum); _assert_(vy_input);
 
@@ -105,6 +105,6 @@
 
 		/*Diffusion */
-		if(sqrt(vx*vx+vy*vy)<10./31536000.){
-			IssmPDouble kappa = 1.e-4;
+		if(sqrt(vx*vx+vy*vy)<1000./31536000.){
+			IssmPDouble kappa = -10;
 			for(int i=0;i<numnodes;i++){
 				for(int j=0;j<numnodes;j++){
@@ -113,5 +113,4 @@
 			}
 		}
-
 		/*Advection: */
 		vx_input->GetInputValue(&vx,gauss);
@@ -119,5 +118,5 @@
 		for(int i=0;i<numnodes;i++){
 			for(int j=0;j<numnodes;j++){
-				Ke->values[i*numnodes+j] += D_scalar*(vx*dbasis[0*numnodes+j]*basis[i] + vy*dbasis[1*numnodes+j]*basis[i]);
+				Ke->values[i*numnodes+j] += (D_scalar*(vx*dbasis[0*numnodes+j]*basis[i] + vy*dbasis[1*numnodes+j]*basis[i]))*1e-2;
 			}
 		}