Index: /issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/sealevelmodel.m
===================================================================
--- /issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/sealevelmodel.m	(revision 22243)
+++ /issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/sealevelmodel.m	(revision 22244)
@@ -108,21 +108,25 @@
 							if length(icecaps{1}.results.TransientSolution(j).(champs{k}))==icecaps{1}.mesh.numberofvertices,
 								md.results.TransientSolution(j).(champs{k})=zeros(md.mesh.numberofvertices,1);
-							else
-								continue;
-							end
-						elseif strcmpi(class(icecaps{1}.results.TransientSolution(j).(champs{k})),'char'),
-							switch champs{k},
-							case 'IceVolume'
-								error('IceVolume mergeresults not implemented yet!');
-							otherwise continue;
+								for l=1:length(trans),
+									resultcap=icecaps{l}.results.TransientSolution(j).(champs{k});
+									md.results.TransientSolution(j).(champs{k})(trans{l})=resultcap;
+								end
+							else 
+								if strcmpi(champs{k},'IceVolume') | strcmpi(champs{k},'IceVolumeAboveFloatation') ,
+									md.results.TransientSolution(j).(champs{k})=0;
+									for l=1:length(trans),
+										resultcap=icecaps{l}.results.TransientSolution(j).(champs{k});
+										md.results.TransientSolution(j).(champs{k})= md.results.TransientSolution(j).(champs{k})+resultcap;
+									end
+								elseif strcmpi(champs{k},'time'),
+									md.results.TransientSolution(j).(champs{k})= icecaps{1}.results.TransientSolution(j).(champs{k});
+								else
+									continue;
+								end
 							end
 						else
 							continue;
 						end
-						for l=1:length(trans),
-							resultcap=icecaps{l}.results.TransientSolution(j).(champs{k});
-							md.results.TransientSolution(j).(champs{k})(trans{l})=resultcap;
-						end
-					end
+											end
 				end
 				self.mergedcaps{2*(i-1)+1}=md;