Changeset 28063

Timestamp:

01/23/24 13:35:32 (22 months ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing some stuff

File:

• issm/trunk-jpl/src/m/contrib/morlighem/sia.m (modified) (2 diffs)

issm/trunk-jpl/src/m/contrib/morlighem/sia.m

@@ -30 +30 @@
 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 +47 @@
 
         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);