source: issm/oecreview/Archive/21724-22754/ISSM-22541-22542.diff@ 22755

Last change on this file since 22755 was 22755, checked in by Mathieu Morlighem, 7 years ago

CHG: added 21724-22754
File size: 1.1 KB
RevLine 
[22755]1Index: ../trunk-jpl/src/m/mech/basalstress.m
2===================================================================
3--- ../trunk-jpl/src/m/mech/basalstress.m (revision 22541)
4+++ ../trunk-jpl/src/m/mech/basalstress.m (revision 22542)
5@@ -18,6 +18,6 @@
6 uby=md.initialization.vy/md.constants.yts;
7
8 %compute basal drag
9-bx=(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base)).^r.*(md.friction.coefficient).^2.*ubx.^s;
10-by=(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base)).^r.*(md.friction.coefficient).^2.*uby.^s;
11-b=(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base)).^r.*(md.friction.coefficient).^2.*ub.^s;
12+bx = - (md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base)).^r.*(md.friction.coefficient).^2.*ubx.^s;
13+by = - (md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base)).^r.*(md.friction.coefficient).^2.*uby.^s;
14+b = sqrt(bx.^2+by.^2);
Note: See TracBrowser for help on using the repository browser.