Index: /issm/trunk/src/m/solutions/control_core.m
===================================================================
--- /issm/trunk/src/m/solutions/control_core.m	(revision 5465)
+++ /issm/trunk/src/m/solutions/control_core.m	(revision 5466)
@@ -27,4 +27,5 @@
 	%Initialize misfits with a vector of zeros
 	J=zeros(nsteps,1);
+	search_scalar=0;
 
 	%Get core from solution type
@@ -34,5 +35,4 @@
 	%Launch once a complete solution to set up all inputs
 	disp('      preparing initial solution');
-	isstokes
 	if isstokes,
 		eval(['femmodel=' solutioncore '(femmodel);']);
@@ -53,5 +53,5 @@
 
 		displaystring(verbose,'\n%s',['      computing gradJ...']);
-		femmodel=gradient_core(femmodel);
+		femmodel=gradient_core(femmodel,n,search_scalar);
 
 		%Return gradient if asked
Index: /issm/trunk/src/m/solutions/gradient_core.m
===================================================================
--- /issm/trunk/src/m/solutions/gradient_core.m	(revision 5465)
+++ /issm/trunk/src/m/solutions/gradient_core.m	(revision 5466)
@@ -35,5 +35,5 @@
 	end
 
-	if (step>0 && search_scalar==0),
+	if (step>1 && search_scalar==0),
 		displaystring(verbose,'\n%s',['      orthogonalization...']);
 		new_gradient=Orth(grad,old_gradient);