Index: /issm/trunk-jpl/src/m/classes/sealevelmodel.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/sealevelmodel.m	(revision 24885)
+++ /issm/trunk-jpl/src/m/classes/sealevelmodel.m	(revision 24886)
@@ -19,5 +19,4 @@
 		settings         = 0;
 		private          = 0;
-		range            = 0;
 		mergedcaps       = 0;
 		transitions      = {};
@@ -94,5 +93,4 @@
 			slm.private           = private();
 			slm.cluster           = generic();
-			slm.range             = {};
 			slm.transitions       = {};
 			slm.eltransitions     = {};
@@ -106,5 +104,4 @@
 			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,5 +109,5 @@
 			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,4 +146,20 @@
 			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,4 +354,7 @@
 			self.earth=md;
 
+			%Create mesh radius:
+			self.earth.mesh.r=planetradius('earth')*ones(md.mesh.numberofvertices,1);
+
 		end % }}}
 		function viscousiterations(self) % {{{
@@ -433,4 +449,12 @@
 			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