Index: sm/trunk/src/m/utils/Analysis/resetspcs.m
===================================================================
--- /issm/trunk/src/m/utils/Analysis/resetspcs.m	(revision 9367)
+++ 	(revision )
@@ -1,18 +1,0 @@
-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.spcvx=NaN*ones(md.numberofnodes);
-	md.spcvy=NaN*ones(md.numberofnodes);
-	md.spcvz=NaN*ones(md.numberofnodes);
-	md.spcvx(pos)=0;
-	md.spcvy(pos)=0;
-	md.spcvz(pos)=0;
-
-	%recompute with spcs reinitialized
-	md.control_analysis=0;
-	md=solve(md,'analysis_type','diagnostic');
-	md=tres(md,'diagnostic');
-
-	vel=md.vel;
-
Index: sm/trunk/src/m/utils/Analysis/runmemessage.m
===================================================================
--- /issm/trunk/src/m/utils/Analysis/runmemessage.m	(revision 9367)
+++ 	(revision )
@@ -1,12 +1,0 @@
-function runmemessage(messages,num,type)
-
-found=0;
-for i=1:size(messages,1),
-	if messages{i,1}==num,
-		disp(sprintf('\n%s: %g\n',[messages{i,2} ' in ' type],num));
-		found=1;
-	end
-end
-if ~found, 
-	error(['runmemessage error message: could not find run ' num]);
-end
Index: /issm/trunk/src/m/utils/Cluster/QueueRequirements.m
===================================================================
--- /issm/trunk/src/m/utils/Cluster/QueueRequirements.m	(revision 9367)
+++ /issm/trunk/src/m/utils/Cluster/QueueRequirements.m	(revision 9368)
@@ -1,8 +1,7 @@
 function QueueRequirements(available_queues,queue_requirements_time,queue_requirements_np,queue,np,time)
-%QUEUEREQUIREMENTS queue requirements in time, number of cpus, by name of queue.
+%QUEUEREQUIREMENTS - queue requirements in time, number of cpus, by name of queue.
 %
-%
-% Usage: QueueRequirements(available_queues,queue_requirements_time,queue_requirements_np,np,time)
-%
+%   Usage: 
+%      QueueRequirements(available_queues,queue_requirements_time,queue_requirements_np,np,time)
 
 %Ok, go through requirements for current queue:
Index: sm/trunk/src/m/utils/Cluster/clusterreq.m
===================================================================
--- /issm/trunk/src/m/utils/Cluster/clusterreq.m	(revision 9367)
+++ 	(revision )
@@ -1,21 +1,0 @@
-function clusterreq(nel),
-%CLUSTERREQ - display memory requisites for cluster
-%
-%   Usage:
-%      clusterreq(nel)
-
-elements=[120000;440000;1700000;5000000];
-memory  =[6898;8825;17130;45013];
-
-
-%from elements,memory, figure out the entire memory needed for nel elements.
-fittedmodel=fit(elements,memory,'linearinterp');
-mem=fittedmodel(nel)
-
-%display memory requisites for several cluster ring sizes: 
-
-disp('Cluster size, Memory per CPU');
-
-for i=1:5:256,
-	disp(sprintf('%i        %f',i,mem/i));
-end
Index: /issm/trunk/src/m/utils/Cluster/parallelrange.m
===================================================================
--- /issm/trunk/src/m/utils/Cluster/parallelrange.m	(revision 9367)
+++ /issm/trunk/src/m/utils/Cluster/parallelrange.m	(revision 9368)
@@ -1,15 +1,13 @@
 function [i1,i2]=parallelrange(rank,numprocs,globalsize)
-%PARALLELRANGE from a rank, and a number of processors, figure out a range, for parallel tasks.
+%PARALLELRANGE - from a rank, and a number of processors, figure out a range, for parallel tasks.
 %
-% usage: [i1,i1]=parallelrange(rank,numprocs,globalsize)
-%
+%   Usage: 
+%      [i1,i1]=parallelrange(rank,numprocs,globalsize)
 
 num_local_rows=zeros(numprocs,1);
 
 for i=1:numprocs,
-
 	%we use floor. we under distribute rows. The rows left  are then redistributed, therefore resulting in a more even distribution.
 	num_local_rows(i)=floor(globalsize/numprocs);
-
 end
 
Index: sm/trunk/src/m/utils/Meca/tau_principal.m
===================================================================
--- /issm/trunk/src/m/utils/Meca/tau_principal.m	(revision 9367)
+++ 	(revision )
@@ -1,37 +1,0 @@
-function [f,g,h,i]=tau_principal(tau_xx,tau_yy,tau_xy)
-%TAU_PRINCIPAL - compute principal component of stress tensor
-%
-%   Usage:
-%      [f,g,h,i]=tau_principal(tau_xx,tau_yy,tau_xy)
-
-nel=length(tau_xx);
-tau_xx=tau_xx/(10^5);
-tau_yy=tau_yy/(10^5);
-tau_xy=tau_xy/(10^5);
-
-%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%determination of principal components of stress tensor%%%%%%%%%%%%
-teta1=zeros(nel,1);
-teta2=zeros(nel,1);
-ind=find((tau_xx-tau_yy)~=0);
-teta1(ind)=1/2*atan(tau_xy(ind)./(tau_xx(ind)-tau_yy(ind)));
-ind=find((tau_xx-tau_yy)==0);
-teta1(ind)=pi/2;
-teta2=teta1+pi/2;
-
-tmp1=tau_xx.*cos(teta1).^2+tau_yy.*sin(teta1).^2+2*tau_xy.*sin(teta1).*cos(teta1);
-tmp2=tau_xx.*cos(teta2).^2+tau_yy.*sin(teta2).^2+2*tau_xy.*sin(teta2).*cos(teta2);
-
-a=tmp1<=tmp2;
-b=tmp1>=tmp2;
-tau_1=tmp1.*b+a.*tmp2;
-tau_2=tmp1.*a+b.*tmp2;
-
-teta1_pp=teta1.*b+a.*teta2;
-teta2_pp=teta1.*a+b.*teta2;
-%teta1_pp=mod(teta1_pp,2*pi);
-%teta2_pp=mod(teta2_pp,2*pi);
-f=tau_1;
-g=tau_2;
-h=teta1_pp;
-i=teta2_pp;
-