Index: /issm/trunk-jpl/src/m/classes/surfaceload.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/surfaceload.m	(revision 25220)
+++ /issm/trunk-jpl/src/m/classes/surfaceload.m	(revision 25221)
@@ -26,5 +26,5 @@
 		
 		end % }}}
-		function md = checkconsistency(self,md,solution,analyses) % {{{
+			function md = checkconsistency(self,md,solution,analyses) % {{{
 
 			if ~ismember('SealevelriseAnalysis',analyses) | (strcmp(solution,'TransientSolution') & md.transient.isslr==0), 
@@ -32,5 +32,11 @@
 			end
 			if ~isempty(self.icethicknesschange),
-				md = checkfield(md,'fieldname','solidearth.surfaceload.icethicknesschange','timeseries',1,'NaN',1,'Inf',1);
+				if isa(self.icethicknesschange,'cell'),
+					for i=1:length(self.icethicknesschange),
+						md = checkfield(md,'field',self.icethicknesschange{i},'NaN',0,'Inf',1,'timeserieslength',1,'Inf',1);
+					end
+				else
+					md = checkfield(md,'field',self.icethicknesschange,'NaN',1,'Inf',1,'timeserieslength',1,'Inf',1);
+				end
 			end
 			if ~isempty(self.waterheightchange),
@@ -69,6 +75,12 @@
 				self.other=zeros(md.mesh.numberofelements+1,1);
 			end
-			WriteData(fid,prefix,'object',self,'fieldname','icethicknesschange','name','md.solidearth.surfaceload.icethicknesschange',...
+			if isa(self.icethicknesschange,'cell'),
+				WriteData(fid,prefix,'object',self,'fieldname','icethicknesschange','name','md.solidearth.surfaceload.icethicknesschange',...
+				'format','MatArray','timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts);
+			else
+				WriteData(fid,prefix,'object',self,'fieldname','icethicknesschange','name','md.solidearth.surfaceload.icethicknesschange',...
 				'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts);
+			end
+
 			WriteData(fid,prefix,'object',self,'fieldname','waterheightchange','name','md.solidearth.surfaceload.waterheightchange',...
 				'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts);
@@ -78,5 +90,5 @@
 		end % }}}
 		function savemodeljs(self,fid,modelname) % {{{
-		
+
 			writejs1Darray(fid,[modelname '.surfaceload.icethicknesschange'],self.icethicknesschange);
 			writejs1Darray(fid,[modelname '.surfaceload.waterheightchange'],self.waterheightchange);