Index: /issm/trunk/test/Validation/ControlMethods/runme.m
===================================================================
--- /issm/trunk/test/Validation/ControlMethods/runme.m	(revision 173)
+++ /issm/trunk/test/Validation/ControlMethods/runme.m	(revision 174)
@@ -1,4 +1,5 @@
 step=[2];
 division=2;
+density=80000;
 
 %step1: initial velocity -> observed velocity
@@ -8,5 +9,5 @@
 	%create model
 	md=model;
-	md=mesh(md,'DomainOutline.exp',100000);
+	md=mesh(md,'DomainOutline.exp',density);
 	md=geography(md,'all','');
 	md=parameterize(md,'Square.par');
@@ -44,5 +45,5 @@
 	%create model
 	md=model;
-	md=mesh(md,'DomainOutline.exp',100000);
+	md=mesh(md,'DomainOutline.exp',density);
 	md=geography(md,'all','');
 	md=parameterize(md,'Square.par');
@@ -54,11 +55,15 @@
 
 	%control method
-	md.nsteps=10;
+	md.nsteps=3;
 	md.control_type='B';
+	md.mincontrolconstraint=0;
+	md.maxcontrolconstraint=1.2*10^9;
 	md.plot=1;
-	md.debug=1;
-	md.optscal=10^7*ones(md.nsteps,1);
+	md.debug=0;
+	md.optscal=10^6*ones(md.nsteps,1);
 	md.fit=0*ones(md.nsteps,1);
-	md.maxiter=30*ones(md.nsteps,1);
+	md.maxiter=2*ones(md.nsteps,1);
+	md.eps_rel=0.1;
+	md.eps_abs=NaN;
 
 	%md.cluster='wilkes';