Index: /issm/trunk/src/m/model/parameterization/parametercontrolB.m
===================================================================
--- /issm/trunk/src/m/model/parameterization/parametercontrolB.m	(revision 8453)
+++ /issm/trunk/src/m/model/parameterization/parametercontrolB.m	(revision 8454)
@@ -39,8 +39,11 @@
 end
 
+
 %cm_min
-cm_min=getfieldvalue(options,'cm_min',1);
-if (length(cm_min)~=1)
-	md.cm_min=1;
+cm_min=getfieldvalue(options,'cm_min',paterson(273.15+5)*ones(md.numberofnodes,1));
+if (length(cm_min)==1)
+	md.cm_min=cm_min*ones(md.numberofnodes,1);
+elseif (length(cm_min)==md.numberofnodes)
+	md.cm_min=cm_min;
 else
 	md.cm_min=cm_min;
@@ -48,7 +51,9 @@
 
 %cm_max
-cm_max=getfieldvalue(options,'cm_max',10^10);
-if (length(cm_max)~=1)
-	md.cm_max=10^10;
+cm_max=getfieldvalue(options,'cm_max',paterson(273.15-70)*ones(md.numberofnodes,1));
+if (length(cm_max)==1)
+	md.cm_max=cm_max*ones(md.numberofnodes,1);
+elseif (length(cm_max)==md.numberofnodes)
+	md.cm_max=cm_max;
 else
 	md.cm_max=cm_max;
Index: /issm/trunk/src/m/model/plot/plot_basaldrag.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_basaldrag.m	(revision 8453)
+++ /issm/trunk/src/m/model/plot/plot_basaldrag.m	(revision 8454)
@@ -30,5 +30,5 @@
 
 %compute basal drag
-drag=(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed)).^r.*(md.drag_coefficient).^2.*ub.^s/1000;
+drag=(max(md.g*(md.rho_ice*md.thickness+md.rho_water*md.bed),0)).^r.*(md.drag_coefficient).^2.*ub.^s/1000;
 
 %Figure out if this is a Section plot