


DRIVINGSTRESS - evaluates the driving stress
The driving stress is computed according to the following formula:
driving stress= rho_ice*g*H*slope
Usage:
[Fx,Fy,Fmag]=drivingstress(md)

0001 function [px,py,pmag]=drivingstress(md) 0002 %DRIVINGSTRESS - evaluates the driving stress 0003 % 0004 % The driving stress is computed according to the following formula: 0005 % driving stress= rho_ice*g*H*slope 0006 % 0007 % Usage: 0008 % [Fx,Fy,Fmag]=drivingstress(md) 0009 0010 %Get slope 0011 [sx,sy,s]=slope(md); 0012 0013 %Average thickness over elements 0014 thickness_bar=(md.thickness(md.elements(:,1))+md.thickness(md.elements(:,2))+md.thickness(md.elements(:,3)))/3; 0015 0016 px=md.rho_ice*md.g*thickness_bar.*sx; 0017 py=md.rho_ice*md.g*thickness_bar.*sy; 0018 pmag=sqrt(px.^2+py.^2);