Index: /issm/trunk-jpl/src/m/mech/calcbackstress.m
===================================================================
--- /issm/trunk-jpl/src/m/mech/calcbackstress.m	(revision 22949)
+++ /issm/trunk-jpl/src/m/mech/calcbackstress.m	(revision 22950)
@@ -62,4 +62,5 @@
 
 % analytical backstress solution
-backstress=T-(1.-D).*B.*sign(ex0).*(2+a0).*abs(ex0).^(1./n)./((1+a0+a0.^2+b0.^2).^((n-1)/2./n));
+%backstress=T-(1.-D).*B.*sign(ex0).*(2+a0).*abs(ex0).^(1./n)./((1+a0+a0.^2+b0.^2).^((n-1)/2./n));
+backstress=1-((1.-D).*B.*sign(ex0).*(2+a0).*abs(ex0).^(1./n)./((1+a0+a0.^2+b0.^2).^((n-1)/2./n)))./T; %fix
 backstress(find(backstress<0))=0;
Index: /issm/trunk-jpl/src/m/plot/plot_tensor_components.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_tensor_components.m	(revision 22949)
+++ /issm/trunk-jpl/src/m/plot/plot_tensor_components.m	(revision 22950)
@@ -74,5 +74,5 @@
 		string=type;
 	end
-	options=addfielddefault(options,'title',[upper(string(1)) string(2:end) ' ' component]);
+	options=changefieldvalue(options,'title',[upper(string(1)) string(2:end) ' ' component]);
 	applyoptions(md,[],options);
 end
Index: /issm/trunk-jpl/src/m/plot/plot_tensor_principal.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_tensor_principal.m	(revision 22949)
+++ /issm/trunk-jpl/src/m/plot/plot_tensor_principal.m	(revision 22950)
@@ -72,5 +72,5 @@
 	string=type;
 end
-options=addfielddefault(options,'title',[upper(string(1)) string(2:end) ' ' component]);
+options=changefieldvalue(options,'title',[upper(string(1)) string(2:end) ' ' component]);
 options=changefieldvalue(options,'colorbar',2);
 applyoptions(md,[],options);
Index: /issm/trunk-jpl/src/m/solve/loadresultsfromcluster.m
===================================================================
--- /issm/trunk-jpl/src/m/solve/loadresultsfromcluster.m	(revision 22949)
+++ /issm/trunk-jpl/src/m/solve/loadresultsfromcluster.m	(revision 22950)
@@ -11,5 +11,5 @@
 options=pairoptions(varargin{:});
 nolog=getfieldvalue(options,'nolog',0);
-md.private.runtimename = getfieldvalue(options,'runtimename',md.private.runtimename)
+md.private.runtimename = getfieldvalue(options,'runtimename',md.private.runtimename);
 
 %retrieve cluster, to be able to call its methods