Index: /issm/trunk/src/m/solutions/cielo/GradJCompute.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/GradJCompute.m	(revision 253)
+++ /issm/trunk/src/m/solutions/cielo/GradJCompute.m	(revision 254)
@@ -3,5 +3,5 @@
 %Recover solution for this stiffness and right hand side: 
 if m.parameters.debug
-	disp('      computing velocities...')
+	disp('         computing velocities...')
 end
 [u_g K_ff0 K_fs0 ]=diagnostic_core_nonlinear(m,inputs);
@@ -9,5 +9,5 @@
 %Buid Du, difference between observed velocity and model velocity.
 if m.parameters.debug
-	disp('      computing Du...')
+	disp('         computing Du...')
 end
 [Du_g]=Du(m.elements,m.nodes,m.loads,m.materials,m.parameters, u_g,u_g_obs,inputs);
@@ -18,5 +18,5 @@
 %Solve for adjoint vector: 
 if m.parameters.debug
-	disp('      computing adjoint state...')
+	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 253)
+++ /issm/trunk/src/m/solutions/cielo/control.m	(revision 254)
@@ -30,5 +30,5 @@
 		inputs=add(inputs,m.parameters.control_type,p_g,'doublevec',2,m.parameters.numberofnodes);
 		inputs=add(inputs,'velocity',m.parameters.u_g,'doublevec',3,m.parameters.numberofnodes);
-		inputs=add(inputs,'fit',m.parameters.fit(n),'scalar');
+		inputs=add(inputs,'fit',m.parameters.fit(n),'double');
 
 		%Update inputs in datasets
Index: /issm/trunk/src/m/solutions/cielo/controlfinalsol.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/controlfinalsol.m	(revision 253)
+++ /issm/trunk/src/m/solutions/cielo/controlfinalsol.m	(revision 254)
@@ -2,4 +2,6 @@
 
 %From parameters, build inputs for icediagnostic_core, using the final parameters
+inputs=parameterinputs;
+inputs=add(inputs,'velocity',m.parameters.u_g,'doublevec',3,m.parameters.numberofnodes);
 inputs=add(inputs,m.parameters.control_type,p_g,'doublevec',2,m.parameters.numberofnodes);
 u_g=diagnostic_core_nonlinear(m,inputs);
Index: /issm/trunk/src/m/solutions/cielo/objectivefunctionC.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/objectivefunctionC.m	(revision 253)
+++ /issm/trunk/src/m/solutions/cielo/objectivefunctionC.m	(revision 254)
@@ -10,4 +10,6 @@
 
 %Plug parameter into inputs
+inputs=parameterinputs;
+inputs=add(inputs,'velocity',m.parameters.u_g,'doublevec',3,m.parameters.numberofnodes);
 inputs=add(inputs,m.parameters.control_type,parameter,'doublevec',2,m.parameters.numberofnodes);
 
@@ -16,4 +18,4 @@
 
 %Compute misfit for this velocity field. 
-inputs=add(inputs,'fit',m.parameters.fit(n),,'scalar');
+inputs=add(inputs,'fit',m.parameters.fit(n),'double');
 J=Misfit(m.elements,m.nodes,m.loads,m.materials,m.parameters, u_g, u_g_obs,inputs);