Index: /issm/trunk/src/m/utils/Analysis/plot_basins.m
===================================================================
--- /issm/trunk/src/m/utils/Analysis/plot_basins.m	(revision 2995)
+++ /issm/trunk/src/m/utils/Analysis/plot_basins.m	(revision 2995)
@@ -0,0 +1,14 @@
+%run in top level Agu09.
+
+figure(1),clf,hold on;
+
+%domain=expread('Exp/DomainOutline50km.exp',1);
+%plot(domain.x,domain.y,'k-');
+
+list=listfiles;
+for i=1:length(list),
+	name=list{i};
+	contour=expread(name,1);
+	plot(contour.x,contour.y,'r-');
+	text(mean(contour.x),mean(contour.y),num2str(name(1:end-4)),'FontSize',18);
+end
Index: /issm/trunk/src/m/utils/Analysis/resetspcs.m
===================================================================
--- /issm/trunk/src/m/utils/Analysis/resetspcs.m	(revision 2995)
+++ /issm/trunk/src/m/utils/Analysis/resetspcs.m	(revision 2995)
@@ -0,0 +1,14 @@
+function vel=resetspcs(md);
+
+	resolution=20000; x_ellsworth=-1.24*10^6; y_ellsworth=1*10^5;
+	pos=find((md.x>(x_ellsworth-resolution)) & (md.x<(x_ellsworth+resolution)) & (md.y>(y_ellsworth-resolution)) & (md.y<(y_ellsworth+resolution)));
+	md.spcvelocity=zeros(md.numberofgrids,6);
+	md.spcvelocity(pos,1:3)=1;
+
+	%recompute with spcs reinitialized
+	md.control_analysis=0;
+	md=solve(md,'analysis_type','diagnostic');
+	md=tres(md,'diagnostic');
+
+	vel=md.vel;
+