Index: /issm/trunk-jpl/src/m/plot/plot_basaldrag.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_basaldrag.m	(revision 27595)
+++ /issm/trunk-jpl/src/m/plot/plot_basaldrag.m	(revision 27596)
@@ -16,15 +16,12 @@
 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')