Changeset 21698


Ignore:
Timestamp:
05/03/17 11:58:22 (8 years ago)
Author:
jbondzio
Message:

BUG: signum error in directional components of driving stress

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/m/mech/drivingstress.m

    r13007 r21698  
    1414thickness_bar=(md.geometry.thickness(md.mesh.elements(:,1))+md.geometry.thickness(md.mesh.elements(:,2))+md.geometry.thickness(md.mesh.elements(:,3)))/3;
    1515
    16 px=md.materials.rho_ice*md.constants.g*thickness_bar.*sx;
    17 py=md.materials.rho_ice*md.constants.g*thickness_bar.*sy;
     16px=-md.materials.rho_ice*md.constants.g*thickness_bar.*sx;
     17py=-md.materials.rho_ice*md.constants.g*thickness_bar.*sy;
    1818pmag=sqrt(px.^2+py.^2);
Note: See TracChangeset for help on using the changeset viewer.