Index: /issm/trunk/test/NightlyRun/test1601.m
===================================================================
--- /issm/trunk/test/NightlyRun/test1601.m	(revision 10429)
+++ /issm/trunk/test/NightlyRun/test1601.m	(revision 10430)
@@ -27,8 +27,20 @@
 disp(['Error between Cartesian and rotated CS: ' num2str(max(abs(vel0-vel1))/(max(abs(vel0))+eps)) ]);
 
+%Now, put CS back to normal except on the side where the spc are applied
+pos=find(x==0 | x==1000000);
+md.diagnostic.referential(:)=NaN;
+md.diagnostic.referential(pos,1:3)=repmat([cos(theta),sin(theta),0],size(pos,1),1);
+md.diagnostic.referential(pos,4:6)=repmat([0,0,1],size(pos,1),1);
+md=solve(md,DiagnosticSolutionEnum);
+vel2=md.results.DiagnosticSolution.Vel;
+
+plotmodel(md,'data',vel0,'data',vel2,'data',vel2-vel0,'title','Cartesian CS','title','Rotated CS','title','difference')
+disp(['Error between Cartesian and rotated CS: ' num2str(max(abs(vel0-vel2))/(max(abs(vel0))+eps)) ]);
+
 %Fields and tolerances to track changes
-field_names     ={'vel1'};
-field_tolerances={1e-11};
+field_names     ={'vel1','vel2'};
+field_tolerances={1e-11,1e-11};
 field_values={...
 	vel1, ...
+	vel2, ...
 	};