Index: /issm/trunk-jpl/src/m/classes/lovenumbers.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/lovenumbers.m	(revision 25148)
+++ /issm/trunk-jpl/src/m/classes/lovenumbers.m	(revision 25149)
@@ -24,15 +24,16 @@
 			options=pairoptions(varargin{:});
 			maxdeg=getfieldvalue(options,'maxdeg',1000);
-			self=setdefaultparameters(self,maxdeg);
+			referenceframe=getfieldvalue(options,'referenceframe','CM');
+			self=setdefaultparameters(self,maxdeg,referenceframe);
 		end % }}}
-		function self = setdefaultparameters(self,maxdeg) % {{{
+		function self = setdefaultparameters(self,maxdeg,referenceframe) % {{{
 		
 			%initialize love numbers:
-			self.h=getlovenumbers('type','loadingverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
-			self.k=getlovenumbers('type','loadinggravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
-			self.l=getlovenumbers('type','loadinghorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
-			self.th=getlovenumbers('type','tidalverticaldisplacement','referenceframe','CM','maxdeg',maxdeg);
-			self.tk=getlovenumbers('type','tidalgravitationalpotential','referenceframe','CM','maxdeg',maxdeg);
-			self.tl=getlovenumbers('type','tidalhorizontaldisplacement','referenceframe','CM','maxdeg',maxdeg);
+			self.h=getlovenumbers('type','loadingverticaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
+			self.k=getlovenumbers('type','loadinggravitationalpotential','referenceframe',referenceframe,'maxdeg',maxdeg);
+			self.l=getlovenumbers('type','loadinghorizontaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
+			self.th=getlovenumbers('type','tidalverticaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
+			self.tk=getlovenumbers('type','tidalgravitationalpotential','referenceframe',referenceframe,'maxdeg',maxdeg);
+			self.tl=getlovenumbers('type','tidalhorizontaldisplacement','referenceframe',referenceframe,'maxdeg',maxdeg);
 
 			%secular fluid love number: