Changeset 27595
- Timestamp:
- 02/20/23 08:45:01 (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/mech/basalstress.m
r24943 r27595 13 13 r=averaging(md,md.friction.q./md.friction.p,0); 14 14 15 %compute horizontal velocity 15 %Compute effective pressure 16 switch(md.friction.coupling) 17 case 0 18 N = max(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base),0); 19 case 3 20 N = max(md.friction.effective_pressure, 0); 21 otherwise 22 error('not supported yet'); 23 end 24 25 %compute sliding velocity 16 26 ub=sqrt(md.initialization.vx.^2+md.initialization.vy.^2)/md.constants.yts; 17 27 ubx=md.initialization.vx/md.constants.yts; … … 19 29 20 30 %compute basal drag (S.I.) 21 alpha2 =(max(md.constants.g*(md.materials.rho_ice*md.geometry.thickness+md.materials.rho_water*md.geometry.base),0).^r).*(md.friction.coefficient.^2).*(ub.^(s-1));31 alpha2 = (N.^r).*(md.friction.coefficient.^2).*(ub.^(s-1)); 22 32 b = alpha2.*ub; 23 33 bx = -alpha2.*ubx; 24 34 by = -alpha2.*uby; 25 35 26 clear alpha2 27 36 %reeeturn magnitude of only one output is requested 37 if nargout==1 38 bx = b; 28 39 end
Note:
See TracChangeset
for help on using the changeset viewer.