Changeset 27596 for issm/trunk-jpl/src/m/plot/plot_basaldrag.m

Timestamp:

02/20/23 08:45:18 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: call basalstress instead of recalculating manually

File:

• issm/trunk-jpl/src/m/plot/plot_basaldrag.m (modified) (1 diff)

issm/trunk-jpl/src/m/plot/plot_basaldrag.m

@@ -16 +16 @@
 end
 
-%compute exponents
-s=averaging(md,1./md.friction.p,0);
-r=averaging(md,md.friction.q./md.friction.p,0);
-
-ub_mag=sqrt(md.initialization.vx.^2+md.initialization.vy.^2)/md.constants.yts;
-drag_mag=(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_mag.^s/1000;
+tau_b    = basalstress(md);
+drag_mag = tau_b/1000;
+ub_mag   = sqrt(md.initialization.vx.^2+md.initialization.vy.^2)/md.constants.yts;
 sig=1;
 
 %compute horizontal velocity 
 if strcmpi(type,'basal_drag')
-        ub=sqrt(md.initialization.vx.^2+md.initialization.vy.^2)/md.constants.yts;
+        ub = ub_mag;
         title_str='Basal drag [kPa]';
 elseif strcmpi(type,'basal_dragx')