Index: /issm/trunk-jpl/test/NightlyRun/test1104.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1104.m	(revision 11030)
+++ /issm/trunk-jpl/test/NightlyRun/test1104.m	(revision 11031)
@@ -15,5 +15,5 @@
 	md=extrude(md,10,1);
 
-	md=setflowequation(md,'pattyn','all','stokes','all');
+	md=setflowequation(md,'stokes','all');
 
 	%Create dirichlet on the bed only
Index: /issm/trunk-jpl/test/NightlyRun/test1106.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1106.m	(revision 11030)
+++ /issm/trunk-jpl/test/NightlyRun/test1106.m	(revision 11031)
@@ -15,5 +15,5 @@
 	md=extrude(md,6,1);
 
-	md=setflowequation(md,'pattyn','all','stokes','all'); 
+	md=setflowequation(md,'stokes','all'); 
 
 	%Create MPCs to have periodic boundary conditions
Index: /issm/trunk-jpl/test/NightlyRun/test1108.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1108.m	(revision 11030)
+++ /issm/trunk-jpl/test/NightlyRun/test1108.m	(revision 11031)
@@ -15,5 +15,5 @@
 	md=extrude(md,6,1);
 
-	md=setflowequation(md,'pattyn','all','stokes','all');
+	md=setflowequation(md,'stokes','all');
 
 	%We need one grd on dirichlet: the 4 corners are set to zero
Index: sm/trunk-jpl/test/NightlyRun/test1603.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1603.m	(revision 11030)
+++ 	(revision )
@@ -1,34 +1,0 @@
-md=triangle(model,'../Exp/Square.exp',150000);
-md=setmask(md,'../Exp/SquareShelf.exp','');
-md=parameterize(md,'../Par/SquareSheetShelf.par');
-md=extrude(md,5,1);
-md=setflowequation(md,'stokes','all');
-md.diagnostic.spcvx(find(md.mesh.y>0))=NaN;
-md.initialization.vx(:)=0;
-md.initialization.vy(:)=0;
-md.initialization.vel(:)=0;
-
-md.cluster=generic('name',oshostname,'np',3);
-md=solve(md,DiagnosticSolutionEnum);
-vel0=md.results.DiagnosticSolution.Vel;
-
-theta=30*pi/180;
-x=md.mesh.x;
-y=md.mesh.y;
-md.mesh.x=cos(theta)*x-sin(theta)*y;
-md.mesh.y=sin(theta)*x+cos(theta)*y;
-
-md.diagnostic.referential(:,1:3)=repmat([cos(theta),sin(theta),0],md.mesh.numberofvertices,1);
-md.diagnostic.referential(:,4:6)=repmat([0,0,1],md.mesh.numberofvertices,1);
-md=solve(md,DiagnosticSolutionEnum);
-vel1=md.results.DiagnosticSolution.Vel;
-
-plotmodel(md,'data',vel0,'data',vel1,'data',vel1-vel0,'title','Cartesian CS','title','Rotated CS','title','difference','view#all',2)
-disp(['Error between Cartesian and rotated CS: ' num2str(max(abs(vel0-vel1))/(max(abs(vel0))+eps)) ]);
-
-%Fields and tolerances to track changes
-field_names     ={'vel1'};
-field_tolerances={1e-9};
-field_values={...
-	vel1, ...
-	};