Changeset 14800

Timestamp:

04/30/13 07:41:52 (12 years ago)

Author:

Mathieu Morlighem

Message:

NEW: surface can be an argument of the function

File:

• issm/trunk-jpl/src/m/geometry/slope.m (modified) (2 diffs)

issm/trunk-jpl/src/m/geometry/slope.m

@@ -1 @@
-function [sx,sy,s]=slope(md)
+function [sx,sy,s]=slope(md,surf)
 %SLOPE - compute the surface slope
 %
 %   Usage:
 %      [sx,sy,s]=slope(md)
+%      [sx,sy,s]=slope(md,md.results.TransientSolution(1).Surface)
 
 %load some variables (it is much faster if the variab;es are loaded from md once for all) 
@@ -18 +19 @@
 end
 
+if nargin==1,
+        surf=md.geometry.surface;
+end
 %compute nodal functions coefficients N(x,y)=alpha x + beta y + gamma
 [alpha beta]=GetNodalFunctionsCoeff(index,x,y);
 
 summation=[1;1;1];
-sx=(md.geometry.surface(index).*alpha)*summation;
-sy=(md.geometry.surface(index).*beta)*summation;
+sx=(surf(index).*alpha)*summation;
+sy=(surf(index).*beta)*summation;
 s=sqrt(sx.^2+sy.^2);