Index: /issm/trunk-jpl/test/NightlyRun/test270.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test270.m	(revision 13624)
+++ /issm/trunk-jpl/test/NightlyRun/test270.m	(revision 13625)
@@ -1,12 +1,12 @@
-md=triangle(model,'../Exp/Square.exp',150000);
+md=triangle(model(),'../Exp/Square.exp',150000);
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
 md.materials=matdamageice();
 md.materials.rheology_B=paterson(md.initialization.temperature);
-md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 md.materials.rheology_Z=0.5*ones(md.mesh.numberofvertices,1);
 md=setflowequation(md,'macayeal','all');
 md.cluster=generic('name',oshostname(),'np',3);
-md=solve(md,DiagnosticSolutionEnum);
+md=solve(md,DiagnosticSolutionEnum());
 
 %Fields and tolerances to track changes
Index: /issm/trunk-jpl/test/NightlyRun/test270.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test270.py	(revision 13625)
+++ /issm/trunk-jpl/test/NightlyRun/test270.py	(revision 13625)
@@ -0,0 +1,32 @@
+import numpy
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from matdamageice import *
+from paterson import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+
+md=triangle(model(),'../Exp/Square.exp',150000.)
+md=setmask(md,'all','')
+md=parameterize(md,'../Par/SquareShelf.py')
+md.materials=matdamageice()
+md.materials.rheology_B=paterson(md.initialization.temperature)
+md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
+md.materials.rheology_Z=0.5*numpy.ones((md.mesh.numberofvertices,1))
+md=setflowequation(md,'macayeal','all')
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,DiagnosticSolutionEnum())
+
+#Fields and tolerances to track changes
+field_names     =['Vx','Vy','Vel','Pressure']
+field_tolerances=[1e-13,1e-13,1e-13,1e-13]
+field_values=[\
+	md.results['DiagnosticSolution'][1]['Vx'],\
+	md.results['DiagnosticSolution'][1]['Vy'],\
+	md.results['DiagnosticSolution'][1]['Vel'],\
+	md.results['DiagnosticSolution'][1]['Pressure'],\
+	]
Index: /issm/trunk-jpl/test/NightlyRun/test272.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test272.m	(revision 13624)
+++ /issm/trunk-jpl/test/NightlyRun/test272.m	(revision 13625)
@@ -1,8 +1,8 @@
-md=triangle(model,'../Exp/Square.exp',150000);
+md=triangle(model(),'../Exp/Square.exp',150000);
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
 md.materials=matdamageice();
 md.materials.rheology_B=paterson(md.initialization.temperature);
-md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 md.materials.rheology_Z=0.5*ones(md.mesh.numberofvertices,1);
 md=setflowequation(md,'macayeal','all');
@@ -14,8 +14,8 @@
 md.inversion.max_parameters=ones(md.mesh.numberofvertices,1);
 md.inversion.nsteps=2;
-md.inversion.cost_functions=101*ones(md.inversion.nsteps,1);
+md.inversion.cost_functions=101.*ones(md.inversion.nsteps,1);
 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
 md.inversion.gradient_scaling=0.9*ones(md.inversion.nsteps,1);
-md.inversion.maxiter_per_step=2*ones(md.inversion.nsteps,1);
+md.inversion.maxiter_per_step=2.*ones(md.inversion.nsteps,1);
 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
 md.inversion.vx_obs=md.initialization.vx; 
@@ -23,5 +23,5 @@
 
 md.cluster=generic('name',oshostname(),'np',3);
-md=solve(md,DiagnosticSolutionEnum);
+md=solve(md,DiagnosticSolutionEnum());
 
 %Fields and tolerances to track changes
@@ -30,5 +30,5 @@
 field_values={...
    (md.results.DiagnosticSolution.Gradient1),...
-   md.results.DiagnosticSolution.J,...
+   (md.results.DiagnosticSolution.J),...
    (md.results.DiagnosticSolution.MaterialsRheologyZbar),...
    (md.results.DiagnosticSolution.Pressure),...
Index: /issm/trunk-jpl/test/NightlyRun/test272.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test272.py	(revision 13625)
+++ /issm/trunk-jpl/test/NightlyRun/test272.py	(revision 13625)
@@ -0,0 +1,50 @@
+import numpy
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from matdamageice import *
+from paterson import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+
+md=triangle(model(),'../Exp/Square.exp',150000.)
+md=setmask(md,'all','')
+md=parameterize(md,'../Par/SquareShelf.py')
+md.materials=matdamageice()
+md.materials.rheology_B=paterson(md.initialization.temperature)
+md.materials.rheology_n=3.*numpy.ones((md.mesh.numberofelements,1))
+md.materials.rheology_Z=0.5*numpy.ones((md.mesh.numberofvertices,1))
+md=setflowequation(md,'macayeal','all')
+
+#control parameters
+md.inversion.iscontrol=1
+md.inversion.control_parameters=['MaterialsRheologyZbar']
+md.inversion.min_parameters=10**-13*numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.max_parameters=numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.nsteps=2
+md.inversion.cost_functions=101.*numpy.ones((md.inversion.nsteps,1))
+md.inversion.cost_functions_coefficients=numpy.ones((md.mesh.numberofvertices,1))
+md.inversion.gradient_scaling=0.9*numpy.ones((md.inversion.nsteps,1))
+md.inversion.maxiter_per_step=2.*numpy.ones((md.inversion.nsteps,1))
+md.inversion.step_threshold=0.99*numpy.ones((md.inversion.nsteps,1))
+md.inversion.vx_obs=md.initialization.vx 
+md.inversion.vy_obs=md.initialization.vy
+
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,DiagnosticSolutionEnum())
+
+#Fields and tolerances to track changes
+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=[\
+   md.results['DiagnosticSolution'][1]['Gradient1'],\
+   md.results['DiagnosticSolution'][1]['J'],\
+   md.results['DiagnosticSolution'][1]['MaterialsRheologyZbar'],\
+   md.results['DiagnosticSolution'][1]['Pressure'],\
+   md.results['DiagnosticSolution'][1]['Vel'],\
+   md.results['DiagnosticSolution'][1]['Vx'],\
+   md.results['DiagnosticSolution'][1]['Vy'],\
+]
Index: /issm/trunk-jpl/test/NightlyRun/test274.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test274.m	(revision 13624)
+++ /issm/trunk-jpl/test/NightlyRun/test274.m	(revision 13625)
@@ -1,3 +1,3 @@
-md=triangle(model,'../Exp/SquareHole.exp','../Exp/Rifts.exp',50000);
+md=triangle(model(),'../Exp/SquareHole.exp','../Exp/Rifts.exp',50000);
 md=meshprocessrifts(md,'../Exp/Square.exp');
 md=setmask(md,'all','');
@@ -5,10 +5,10 @@
 md.materials=matdamageice();
 md.materials.rheology_B=paterson(md.initialization.temperature);
-md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.materials.rheology_n=3.*ones(md.mesh.numberofelements,1);
 md.materials.rheology_Z=0.5*ones(md.mesh.numberofvertices,1);
 md=setflowequation(md,'macayeal','all');
 
 md.cluster=generic('name',oshostname(),'np',3);
-md=solve(md,DiagnosticSolutionEnum);
+md=solve(md,DiagnosticSolutionEnum());
 
 %Fields and tolerances to track changes