Index: /issm/trunk/src/m/model/basalstress.m
===================================================================
--- /issm/trunk/src/m/model/basalstress.m	(revision 6041)
+++ /issm/trunk/src/m/model/basalstress.m	(revision 6042)
@@ -9,6 +9,6 @@
 
 %compute exponents
-s=averaging(md,1./md.p,0);
-r=averaging(md,md.p./md.q,0);
+s=averaging(md,1./md.drag_p,0);
+r=averaging(md,md.drag_p./md.drag_q,0);
 
 %compute horizontal velocity
@@ -18,5 +18,5 @@
 
 %compute basal drag
-bx=(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag).^2.*ubx.^s;
-by=(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag).^2.*uby.^s;
-b=(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag).^2.*ub.^s;
+bx=(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag_coefficient).^2.*ubx.^s;
+by=(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag_coefficient).^2.*uby.^s;
+b=(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag_coefficient).^2.*ub.^s;