Index: /issm/trunk-jpl/src/m/classes/inversion.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/inversion.py	(revision 13641)
+++ /issm/trunk-jpl/src/m/classes/inversion.py	(revision 13642)
@@ -96,8 +96,8 @@
 		#inversed parameter (10^8 for B, 50 for drag) and can be decreased
 		#after the first iterations
-		self.gradient_scaling=50*numpy.ones(self.nsteps)
+		self.gradient_scaling=50*numpy.ones((self.nsteps,1))
 
 		#several responses can be used:
-		self.cost_functions=101*numpy.ones(self.nsteps)
+		self.cost_functions=101*numpy.ones((self.nsteps,1))
 
 		#step_threshold is used to speed up control method. When
@@ -124,5 +124,5 @@
 
 		num_controls=numpy.size(md.inversion.control_parameters)
-		num_costfunc=numpy.size(md.inversion.cost_functions,1)
+		num_costfunc=numpy.size(md.inversion.cost_functions,axis=1)
 
 		md = checkfield(md,'inversion.iscontrol','values',[0,1])
@@ -177,5 +177,5 @@
 
 		#process cost functions
-		num_cost_functions=numpy.size(self.cost_functions,1)
+		num_cost_functions=numpy.size(self.cost_functions,axis=1)
 		data=self.cost_functions
 		data[numpy.nonzero(data==101)]=SurfaceAbsVelMisfitEnum()