Index: /issm/trunk-jpl/test/NightlyRun/test101.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test101.m	(revision 22500)
+++ /issm/trunk-jpl/test/NightlyRun/test101.m	(revision 22501)
@@ -7,5 +7,5 @@
 
 %output
-md.stressbalance.requested_outputs={'default','DeviatoricStressxx','DeviatoricStressyy','DeviatoricStressxy','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6','DeviatoricStresseffective'};
+md.stressbalance.requested_outputs={'default','DeviatoricStressxx','DeviatoricStressyy','DeviatoricStressxy','MassFlux1','MassFlux2','MassFlux3','MassFlux4','MassFlux5','MassFlux6'};
 md.outputdefinition.definitions={...
 	massfluxatgate('name','MassFlux1','profilename',['../Exp/MassFlux1.exp'],'definitionstring','Outputdefinition1'),...
Index: /issm/trunk-jpl/test/NightlyRun/test1101.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1101.m	(revision 22500)
+++ /issm/trunk-jpl/test/NightlyRun/test1101.m	(revision 22501)
@@ -12,14 +12,13 @@
 for i=1:length(L_list),
 	L=L_list{i};
-	nx=10; %numberof nodes in x direction
-	ny=10;
+	nx=20; %numberof nodes in x direction
+	ny=20;
 	md=model();
 	md=squaremesh(md,L,L,nx,ny);
 	md=setmask(md,'',''); %ice sheet test
 	md=parameterize(md,'../Par/ISMIPA.par');
-	md=extrude(md,4,1.);
+	md=extrude(md,9,1.);
 
-	md=setflowequation(md,'FS','all');
-	md.flowequation.fe_FS='MINI';
+	md=setflowequation(md,'HO','all');
 
 	%Create dirichlet on the bed only
@@ -43,7 +42,5 @@
 	%Compute the stressbalance
 	md.cluster=generic('name',oshostname(),'np',8);
-	md.stressbalance.maxiter=1;
 	md=solve(md,'Stressbalance');
-	error('S');
 
 	%Plot the results and save them
Index: /issm/trunk-jpl/test/NightlyRun/test613.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test613.m	(revision 22500)
+++ /issm/trunk-jpl/test/NightlyRun/test613.m	(revision 22501)
@@ -15,20 +15,13 @@
 md.inversion.iscontrol=1;
 md.inversion.nsteps=2;
-%md.inversion.control_parameters={'Vx','Vy'};
-md.inversion.control_parameters={'BalancethicknessSpcthickness'};
+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.gradient_scaling=[100];
-%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.min_parameters=[md.geometry.thickness-100];
-md.inversion.max_parameters=[md.geometry.thickness+100];
+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;
 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1);
 md.inversion.maxiter_per_step=4*ones(md.inversion.nsteps,1);
 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
-md.inversion = m1qn3inversion(md.inversion);
-
-md.inversion.surface_obs = ones(md.mesh.numberofvertices,1);
 
 md.verbose.control=1;