drivingstress

PURPOSE ^

DRIVINGSTRESS - evaluates the driving stress

SYNOPSIS ^

function [px,py,pmag]=drivingstress(md)

DESCRIPTION ^

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)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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);

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003