Index: /issm/trunk/test/Validation/ControlMethods/Square.par
===================================================================
--- /issm/trunk/test/Validation/ControlMethods/Square.par	(revision 279)
+++ /issm/trunk/test/Validation/ControlMethods/Square.par	(revision 280)
@@ -37,5 +37,5 @@
 
 %Parallel options
-md.np=8;
+md.np=3;
 md.time=50;
 md.waitonlock=1;
Index: /issm/trunk/test/Validation/ControlMethods/runme.m
===================================================================
--- /issm/trunk/test/Validation/ControlMethods/runme.m	(revision 279)
+++ /issm/trunk/test/Validation/ControlMethods/runme.m	(revision 280)
@@ -1,5 +1,6 @@
 step=[2];
 division=2;
-density=80000;
+density=50000;
+startup
 
 %step1: initial velocity -> observed velocity
@@ -30,4 +31,5 @@
 	%diagnostic with this B
 	md=solve(md,'diagnostic','cielo');
+	%md.cluster='wilkes'; md=solve(md,'diagnostic_horiz','cielo');
 
 	%save observed velocities and exact B
@@ -55,5 +57,5 @@
 
 	%control method
-	md.nsteps=3;
+	md.nsteps=50;
 	md.control_type='B';
 	md.mincontrolconstraint=0;
@@ -61,12 +63,12 @@
 	md.plot=1;
 	md.debug=0;
-	md.optscal=10^6*ones(md.nsteps,1);
+	md.optscal=10^8*ones(md.nsteps,1);
 	md.fit=0*ones(md.nsteps,1);
-	md.maxiter=2*ones(md.nsteps,1);
-	md.eps_rel=0.1;
+	md.maxiter=10*ones(md.nsteps,1);
+	md.eps_rel=1;
 	md.eps_abs=NaN;
 
-	%md.cluster='wilkes';
-	md.eps_abs=NaN;
+%	md.cluster='wilkes';
+%	md=solve(md,'control','ice');
 	md=solve(md,'control','cielo');
 end