Changeset 22230

Timestamp:

11/03/17 19:37:33 (7 years ago)

Author:

Eric.Larour

Message:

CHG: diverse

File:

• issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/sealevelmodel.m (modified) (2 diffs)

issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/sealevelmodel.m

@@ -18 +18 @@
                 settings         = 0;
                 private          = 0;
-                range           = 0;
+                range            = 0;
+                mergedcaps     = 0;
                 %}}}
         end
@@ -96 +97 @@
                         disp(sprintf('%19s: %-22s -- %s','range'   ,['[1x1 ' class(self.range) ']'],'ranges'));
                 end % }}}
+                function self=mergeresults(self) % {{{
+                        champs=fields(self.icecaps{1}.results.TransientSolution);
+                        for i=1:length(self.mergedcaps)/2,
+                                md=self.mergedcaps{2*(i-1)+1}; trans=self.mergedcaps{2*(i-1)+2};
+                                icecaps=self.icecaps(self.range{2*(i-1)+2});
+                                for j=1:length(self.icecaps{1}.results.TransientSolution),
+                                        for k=1:length(champs),
+                                                if strcmpi(class(icecaps{1}.results.TransientSolution(j).(champs{k})),'double'),
+                                                        %vertex or element? 
+                                                        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;
+                                                        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
+                                self.mergedcaps{2*(i-1)+1}=md;
+                        end
+                end % }}}
         end
 end