Index: /issm/trunk/src/m/solutions/cielo/GradJCompute.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/GradJCompute.m	(revision 171)
+++ /issm/trunk/src/m/solutions/cielo/GradJCompute.m	(revision 172)
@@ -2,7 +2,13 @@
 
 %Recover solution for this stiffness and right hand side: 
+if m.parameters.debug
+	disp('      computing velocities...')
+end
 [u_g K_ff0 K_fs0 ]=diagnostic_core_nonlinear(m,inputs);
 
 %Buid Du, difference between observed velocity and model velocity.
+if m.parameters.debug
+	disp('      computing Du...')
+end
 [Du_g]=Du(m.elements,m.nodes,m.loads,m.materials,m.parameters, u_g,u_g_obs,inputs);
 
@@ -11,4 +17,7 @@
 
 %Solve for adjoint vector: 
+if m.parameters.debug
+	disp('      computing adjoint state...')
+end
 lambda_f=Solver(K_ff0,Du_f,[],m.parameters);
 
Index: /issm/trunk/src/m/solutions/cielo/control.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/control.m	(revision 171)
+++ /issm/trunk/src/m/solutions/cielo/control.m	(revision 172)
@@ -31,5 +31,5 @@
 
 		%Update inputs in datasets
-		[m.elements,m.nodes, loads,m.materials]=UpdateFromInputs(m.elements,m.nodes, loads,m.materials,inputs);
+		[m.elements,m.nodes,m.loads,m.materials]=UpdateFromInputs(m.elements,m.nodes,m.loads,m.materials,inputs);
 
 		disp('      computing gradJ...');
@@ -85,5 +85,4 @@
 	disp('      done.');
 	
-	
 	t2=clock;
 	disp(['      overall time spent on control code ' num2str(etime(t2,t1)) ' seconds'])
Index: /issm/trunk/src/m/solutions/cielo/plot_direction.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/plot_direction.m	(revision 171)
+++ /issm/trunk/src/m/solutions/cielo/plot_direction.m	(revision 172)
@@ -1,4 +1,1 @@
-grad_g=c(n).grad_g;
-grad_gx=grad_g(1:2:end);grad_gy=grad_g(2:2:end); grad_g=sqrt(grad_gx.^2+grad_gy.^2);
-plotmodel(md,'data',grad_g,'title',['Normalized Direction for ' m.parameters.control_type],'figure',1,'colorbar#all','on'); drawnow;
-clear grad_g grad_gx grad_gy
+plotmodel(md,'data',sqrt(c(n).grad_g(1:2:end).^2+c(n).grad_g(2:2:end).^2),'title',['Normalized Direction for ' m.parameters.control_type],'figure',1,'colorbar#all','on'); drawnow; 
Index: /issm/trunk/src/m/solutions/cielo/plot_newdistribution.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/plot_newdistribution.m	(revision 171)
+++ /issm/trunk/src/m/solutions/cielo/plot_newdistribution.m	(revision 172)
@@ -1,3 +1,1 @@
-p_gx=p_g(1:2:end);p_gy=p_g(2:2:end); p_g=sqrt(p_gx.^2+p_gy.^2);
-plotmodel(md,'data',p_g,'title',['Distribution of ' m.parameters.control_type 'at iteration' num2str(n)],'figure',1,'colorbar#all','on'); drawnow;
-clear p_gx p_gy
+plotmodel(md,'data',p_g(1:2:end),'title',['Distribution of ' m.parameters.control_type ' at iteration ' num2str(n)],'figure',2,'colorbar#all','on'); pause(1);