Index: /issm/trunk-jpl/test/NightlyRun/test1401.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1401.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test1401.m	(revision 13669)
@@ -72,6 +72,6 @@
 
 %Fields and tolerances to track changes
-field_names     ={'xyams' 'yyams' 'xbamg' 'ybamg'};
-field_tolerances={1e-13 1e-13 1e-13 1e-13};
+field_names     ={'xyams','yyams','xbamg','ybamg'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
 field_values={...
 	x1, y1,...
Index: /issm/trunk-jpl/test/NightlyRun/test1402.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1402.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test1402.m	(revision 13669)
@@ -102,6 +102,6 @@
 
 %Fields and tolerances to track changes
-field_names     ={'xyams' 'yyams' 'xbamg' 'ybamg'};
-field_tolerances={1e-13 1e-13 1e-13 1e-13};
+field_names     ={'xyams','yyams','xbamg','ybamg'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
 field_values={...
 	x1, y1,...
Index: /issm/trunk-jpl/test/NightlyRun/test272.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test272.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test272.m	(revision 13669)
@@ -26,5 +26,5 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Gradient' 'Misfits' 'MaterialsRheologyZbar' 'Pressure' 'Vel' 'Vx' 'Vy'};
+field_names     ={'Gradient','Misfits','MaterialsRheologyZbar','Pressure','Vel','Vx','Vy'};
 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
 field_values={...
Index: /issm/trunk-jpl/test/NightlyRun/test415.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test415.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test415.m	(revision 13669)
@@ -23,5 +23,5 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Gradient' 'Misfits' 'FrictionCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalforcingsMeltingRate'};
+field_names     ={'Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy','Vz','Temperature','BasalforcingsMeltingRate'};
 field_tolerances={1e-10,1e-10,1e-11,1e-13,1e-10,1e-09,1e-10,1e-8,1e-09,1e-6};
 field_values={...
Index: /issm/trunk-jpl/test/NightlyRun/test416.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test416.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test416.m	(revision 13669)
@@ -23,5 +23,5 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Gradient' 'Misfits' 'FrictionCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalforcingsMeltingRate'};
+field_names     ={'Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy','Vz','Temperature','BasalforcingsMeltingRate'};
 field_tolerances={1e-08,1e-07,1e-08,1e-08,1e-08,1e-08,1e-08,1e-07,1e-08,1e-05};
 field_values={...
Index: /issm/trunk-jpl/test/NightlyRun/test511.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test511.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test511.m	(revision 13669)
@@ -27,5 +27,5 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Gradient' 'Misfits' 'MaterialsRheologyB' 'Pressure' 'Vel' 'Vx' 'Vy'};
+field_names     ={'Gradient','Misfits','MaterialsRheologyB','Pressure','Vel','Vx','Vy'};
 field_tolerances={1e-11,1e-12,1e-12,1e-09,1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
 field_values={...
Index: /issm/trunk-jpl/test/NightlyRun/test610.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test610.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test610.m	(revision 13669)
@@ -1,6 +1,6 @@
-md=triangle(model(),'../Exp/79North.exp',10000);
+md=triangle(model(),'../Exp/79North.exp',10000.);
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
-md=extrude(md,3,1);
+md=extrude(md,3,1.);
 md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
Index: /issm/trunk-jpl/test/NightlyRun/test610.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test610.py	(revision 13669)
+++ /issm/trunk-jpl/test/NightlyRun/test610.py	(revision 13669)
@@ -0,0 +1,24 @@
+import numpy
+from model import *
+from EnumDefinitions import *
+from MatlabFuncs import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from solve import *
+
+md=triangle(model(),'../Exp/79North.exp',10000.)
+md=setmask(md,'../Exp/79NorthShelf.exp','')
+md=parameterize(md,'../Par/79North.py')
+md.extrude(3,1.)
+md=setflowequation(md,'macayeal','all')
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,BalancethicknessSolutionEnum())
+
+#Fields and tolerances to track changes
+field_names     =['Thickness']
+field_tolerances=[1e-13]
+field_values=[\
+	md.results['BalancethicknessSolution'][1]['Thickness'],\
+	]
Index: /issm/trunk-jpl/test/NightlyRun/test611.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test611.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test611.m	(revision 13669)
@@ -1,3 +1,3 @@
-md=triangle(model(),'../Exp/79North.exp',10000);
+md=triangle(model(),'../Exp/79North.exp',10000.);
 md=setmask(md,'../Exp/79NorthShelf.exp','');
 md=parameterize(md,'../Par/79North.par');
@@ -10,9 +10,9 @@
 md.inversion.control_parameters={'BalancethicknessThickeningRate'};
 md.inversion.thickness_obs=md.geometry.thickness;
-md.inversion.min_parameters=-50*ones(md.mesh.numberofvertices,1);
-md.inversion.max_parameters=50*ones(md.mesh.numberofvertices,1);
+md.inversion.min_parameters=-50.*ones(md.mesh.numberofvertices,1);
+md.inversion.max_parameters=50.*ones(md.mesh.numberofvertices,1);
 md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
-md.inversion.gradient_scaling=10/md.constants.yts*ones(md.inversion.nsteps,1);
+md.inversion.gradient_scaling=10./md.constants.yts*ones(md.inversion.nsteps,1);
 md.inversion.maxiter_per_step=4*ones(md.inversion.nsteps,1);
 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
@@ -22,5 +22,5 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Gradient' 'Misfits' 'BalancethicknessThickeningRate' 'Thickness'};
+field_names     ={'Gradient','Misfits','BalancethicknessThickeningRate','Thickness'};
 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
 field_values={...
Index: /issm/trunk-jpl/test/NightlyRun/test611.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test611.py	(revision 13669)
+++ /issm/trunk-jpl/test/NightlyRun/test611.py	(revision 13669)
@@ -0,0 +1,41 @@
+import numpy
+from model import *
+from EnumDefinitions import *
+from MatlabFuncs import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from solve import *
+
+md=triangle(model(),'../Exp/79North.exp',10000.)
+md=setmask(md,'../Exp/79NorthShelf.exp','')
+md=parameterize(md,'../Par/79North.py')
+md=setflowequation(md,'macayeal','all')
+
+#control parameters
+md.inversion.nsteps=2
+md.prognostic.stabilization=1
+md.inversion.iscontrol=1
+md.inversion.control_parameters=['BalancethicknessThickeningRate']
+md.inversion.thickness_obs=md.geometry.thickness
+md.inversion.min_parameters=-50.*numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.max_parameters=50.*numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.cost_functions=201*numpy.ones((md.inversion.nsteps,1))
+md.inversion.cost_functions_coefficients=numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.gradient_scaling=10./md.constants.yts*numpy.ones((md.inversion.nsteps,1))
+md.inversion.maxiter_per_step=4*numpy.ones((md.inversion.nsteps,1))
+md.inversion.step_threshold=0.99*numpy.ones((md.inversion.nsteps,1))
+
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,BalancethicknessSolutionEnum())
+
+#Fields and tolerances to track changes
+field_names     =['Gradient','Misfits','BalancethicknessThickeningRate','Thickness']
+field_tolerances=[1e-12,1e-12,1e-12,1e-12,1e-12,1e-12]
+field_values=[\
+	md.results['BalancethicknessSolution'][1]['Gradient1'],\
+	md.results['BalancethicknessSolution'][1]['J'],\
+	md.results['BalancethicknessSolution'][1]['BalancethicknessThickeningRate'],\
+	md.results['BalancethicknessSolution'][1]['Thickness']
+]
Index: /issm/trunk-jpl/test/NightlyRun/test612.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test612.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test612.m	(revision 13669)
@@ -1,3 +1,3 @@
-md=triangle(model(),'../Exp/79North.exp',10000);
+md=triangle(model(),'../Exp/79North.exp',10000.);
 md=meshconvert(md);
 md=setmask(md,'../Exp/79NorthShelf.exp','');
@@ -11,9 +11,9 @@
 md.inversion.control_parameters={'BalancethicknessThickeningRate'};
 md.inversion.thickness_obs=md.geometry.thickness;
-md.inversion.min_parameters=-50*ones(md.mesh.numberofvertices,1);
-md.inversion.max_parameters=50*ones(md.mesh.numberofvertices,1);
+md.inversion.min_parameters=-50.*ones(md.mesh.numberofvertices,1);
+md.inversion.max_parameters=50.*ones(md.mesh.numberofvertices,1);
 md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
-md.inversion.gradient_scaling=10/md.constants.yts*ones(md.inversion.nsteps,1);
+md.inversion.gradient_scaling=10./md.constants.yts*ones(md.inversion.nsteps,1);
 md.inversion.maxiter_per_step=4*ones(md.inversion.nsteps,1);
 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
@@ -23,5 +23,5 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Gradient' 'Misfits' 'BalancethicknessThickeningRate' 'Thickness'};
+field_names     ={'Gradient','Misfits','BalancethicknessThickeningRate','Thickness'};
 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
 field_values={...
Index: /issm/trunk-jpl/test/NightlyRun/test612.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test612.py	(revision 13669)
+++ /issm/trunk-jpl/test/NightlyRun/test612.py	(revision 13669)
@@ -0,0 +1,43 @@
+import numpy
+from model import *
+from EnumDefinitions import *
+from MatlabFuncs import *
+from triangle import *
+from meshconvert import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from solve import *
+
+md=triangle(model(),'../Exp/79North.exp',10000.)
+md=meshconvert(md)
+md=setmask(md,'../Exp/79NorthShelf.exp','')
+md=parameterize(md,'../Par/79North.py')
+md=setflowequation(md,'macayeal','all')
+
+#control parameters
+md.inversion.nsteps=2
+md.balancethickness.stabilization=3
+md.inversion.iscontrol=1
+md.inversion.control_parameters=['BalancethicknessThickeningRate']
+md.inversion.thickness_obs=md.geometry.thickness
+md.inversion.min_parameters=-50.*numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.max_parameters=50.*numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.cost_functions=201*numpy.ones((md.inversion.nsteps,1))
+md.inversion.cost_functions_coefficients=numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.gradient_scaling=10./md.constants.yts*numpy.ones((md.inversion.nsteps,1))
+md.inversion.maxiter_per_step=4*numpy.ones((md.inversion.nsteps,1))
+md.inversion.step_threshold=0.99*numpy.ones((md.inversion.nsteps,1))
+
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,BalancethicknessSolutionEnum())
+
+#Fields and tolerances to track changes
+field_names     =['Gradient','Misfits','BalancethicknessThickeningRate','Thickness']
+field_tolerances=[1e-12,1e-12,1e-12,1e-12,1e-12,1e-12]
+field_values=[\
+	md.results['BalancethicknessSolution'][1]['Gradient1'],\
+	md.results['BalancethicknessSolution'][1]['J'],\
+	md.results['BalancethicknessSolution'][1]['BalancethicknessThickeningRate'],\
+	md.results['BalancethicknessSolution'][1]['Thickness']
+]
Index: /issm/trunk-jpl/test/NightlyRun/test613.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test613.m	(revision 13668)
+++ /issm/trunk-jpl/test/NightlyRun/test613.m	(revision 13669)
@@ -1,3 +1,3 @@
-md=triangle(model(),'../Exp/79North.exp',10000);
+md=triangle(model(),'../Exp/79North.exp',10000.);
 md=meshconvert(md);
 md=setmask(md,'../Exp/79NorthShelf.exp','');
@@ -15,9 +15,9 @@
 md.inversion.thickness_obs=md.geometry.thickness;
 md.inversion.nsteps=2;
-md.inversion.control_parameters={'Vx' 'Vy'};
+md.inversion.control_parameters={'Vx','Vy'};
 md.balancethickness.stabilization=1;
-md.inversion.gradient_scaling=[10/md.constants.yts*ones(md.inversion.nsteps,1) 10/md.constants.yts*ones(md.inversion.nsteps,1)];
-md.inversion.min_parameters=[-2000*ones(md.mesh.numberofvertices,1) -2000*ones(md.mesh.numberofvertices,1)];
-md.inversion.max_parameters=[+2000*ones(md.mesh.numberofvertices,1) +2000*ones(md.mesh.numberofvertices,1)];
+md.inversion.gradient_scaling=[10./md.constants.yts*ones(md.inversion.nsteps,1) 10./md.constants.yts*ones(md.inversion.nsteps,1)];
+md.inversion.min_parameters=[-2000.*ones(md.mesh.numberofvertices,1) -2000.*ones(md.mesh.numberofvertices,1)];
+md.inversion.max_parameters=[+2000.*ones(md.mesh.numberofvertices,1) +2000.*ones(md.mesh.numberofvertices,1)];
 md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
@@ -30,5 +30,5 @@
 
 %Fields and tolerances to track changes
-field_names     ={'Gradient1' 'Gradient2' 'Misfits' 'Vx' 'Vy' 'Thickness'};
+field_names     ={'Gradient1','Gradient2','Misfits','Vx','Vy','Thickness'};
 field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
 field_values={...