Index: /issm/trunk-jpl/src/m/classes/stressbalance.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/stressbalance.m	(revision 20090)
+++ /issm/trunk-jpl/src/m/classes/stressbalance.m	(revision 20091)
@@ -89,4 +89,5 @@
 			%Early return
 			if ~ismember(StressbalanceAnalysisEnum(),analyses), return; end
+			if (solution==TransientSolutionEnum & md.transient.isstressbalance == 0), return; end
 
 			md = checkfield(md,'fieldname','stressbalance.spcvx','Inf',1,'timeseries',1);
@@ -195,7 +196,9 @@
 			WriteData(fid,'object',self,'class','stressbalance','fieldname','referential','format','DoubleMat','mattype',1);
 
-			WriteData(fid,'data',self.loadingforce(:,1),'format','DoubleMat','mattype',1,'enum',LoadingforceXEnum);
-			WriteData(fid,'data',self.loadingforce(:,2),'format','DoubleMat','mattype',1,'enum',LoadingforceYEnum);
-			WriteData(fid,'data',self.loadingforce(:,3),'format','DoubleMat','mattype',1,'enum',LoadingforceZEnum);
+			if size(self.loadingforce,2)==3,
+				WriteData(fid,'data',self.loadingforce(:,1),'format','DoubleMat','mattype',1,'enum',LoadingforceXEnum);
+				WriteData(fid,'data',self.loadingforce(:,2),'format','DoubleMat','mattype',1,'enum',LoadingforceYEnum);
+				WriteData(fid,'data',self.loadingforce(:,3),'format','DoubleMat','mattype',1,'enum',LoadingforceZEnum);
+			end
 
 			%process requested outputs