Changeset 24886

Timestamp:

05/22/20 09:30:47 (5 years ago)

Author:

Eric.Larour

Message:

CHG: added functionality to support mdanalysis routine, viewer.

File:

• issm/trunk-jpl/src/m/classes/sealevelmodel.m (modified) (7 diffs)

issm/trunk-jpl/src/m/classes/sealevelmodel.m ¶

@@ -19 +19 @@
                 settings         = 0;
                 private          = 0;
-                range            = 0;
                 mergedcaps       = 0;
                 transitions      = {};
@@ -94 +93 @@
                         slm.private           = private();
                         slm.cluster           = generic();
-                        slm.range             = {};
                         slm.transitions       = {};
                         slm.eltransitions     = {};
@@ -106 +104 @@
                         disp(sprintf('%19s: %-22s -- %s','cluster'         ,['[1x1 ' class(self.cluster) ']'],'cluster parameters (number of cpus...)'));
                         disp(sprintf('%19s: %-22s -- %s','miscellaneous'   ,['[1x1 ' class(self.miscellaneous) ']'],'miscellaneous fields'));
-                        disp(sprintf('%19s: %-22s -- %s','range'   ,['[1x1 ' class(self.range) ']'],'ranges'));
                 end % }}}
                 function self=mergeresults(self) % {{{
@@ -112 +109 @@
                         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});
+                                %icecaps=self.icecaps(self.range{2*(i-1)+2});
                                 for j=1:length(self.icecaps{1}.results.TransientSolution),
                                         for k=1:length(champs),
@@ -149 +146 @@
                         end
                 end % }}}
+                function list=continents(self) % {{{
+                        list={};
+                        for  i=1:length(self.basins),
+                                list{end+1}=self.basins{i}.continent;
+                        end
+                        list=unique(list);
+                end % }}}
+                function list=basinsfromcontinent(self,continent) % {{{
+                        list={};
+                        for  i=1:length(self.icecaps),
+                                if strcmpi(self.basins{i}.continent,continent),
+                                        list{end+1}=self.basins{i}.name;
+                                end
+                        end
+                        list=unique(list);
+                end % }}}
                 function addbasin(self,bas) % {{{
                 if ~strcmpi(class(bas),'basin')
@@ -341 +354 @@
                         self.earth=md;
 
+                        %Create mesh radius:
+                        self.earth.mesh.r=planetradius('earth')*ones(md.mesh.numberofvertices,1);
+
                 end % }}}
                 function viscousiterations(self) % {{{
@@ -433 +449 @@
                         self.earth=ic;
                 end % }}}
+                function self=initializemodels(self) % {{{
+
+                        for i=1:length(self.basins),
+                                md=model();
+                                md.miscellaneous.name=self.basins{i}.name;
+                                self.addicecap(md);
+                        end
+                end % }}}
         end
 end