Index: /issm/trunk-jpl/src/m/contrib/morlighem/sia.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/morlighem/sia.m	(revision 28062)
+++ /issm/trunk-jpl/src/m/contrib/morlighem/sia.m	(revision 28063)
@@ -30,8 +30,11 @@
 else
 
+	%Smooth surface slope a little bit first
 	[sx,sy,s]=slope(md);
-	sx = average(md, sx, 3);
-	sy = average(md, sy, 3);
+	sx = averaging(md, sx, 3);
+	sy = averaging(md, sy, 3);
 	surf = sqrt(sx.^2 + sy.^2);
+
+	%Get other variables
 	s  = md.geometry.surface;
 	z  = md.mesh.z;
@@ -44,6 +47,9 @@
 
 	chi = 1-(s-z)./H;
-	vz  = -md.smb.mass_balance.*((n+2)*chi + (1-chi).^(n+2) -1)/(n+1)
+	vz  = -md.smb.mass_balance.*((n+2)*chi + (1-chi).^(n+2) -1)/(n+1);
 end
 
+vx = vx*md.constants.yts;
+vy = vy*md.constants.yts;
+
 vel=sqrt(vx.^2+vy.^2+vz.^2);