Index: /issm/trunk-jpl/src/m/boundaryconditions/love_numbers.m
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/love_numbers.m	(revision 21305)
+++ /issm/trunk-jpl/src/m/boundaryconditions/love_numbers.m	(revision 21306)
@@ -1,17 +1,24 @@
-function series=love_numbers(type)
-%LOVE_NUMBERS: provide love numbers (type 'h','k','l','gamma' and 'lambda'
-%			   retrieved from: http://www.srosat.com/iag-jsg/loveNb.php
-%    Usage:   series=love_numbers(type) 
-%             where type is one of 'h','k','l','gamma' and 'lambda'. 
+function series=love_numbers(type,reference_frame)
+%LOVE_NUMBERS: provide love numbers 
+%	retrieved from: http://www.srosat.com/iag-jsg/loveNb.php 
+%	in a chosen reference frame 
+%  
+%	Usage: series=love_numbers(type,reference_frame) 
+%	type is one of 'h','k','l','gamma' and 'lambda'. 
+%	reference_frame is one of 'CM' and 'CF'.
 %
-%    Example:  
-%          love_k=love_numbers('k');
+%	Example:  
+%	love_k=love_numbers('k','CF');
 % 
 
 	%some checks:
-	if nargin~=1, error('love_numbers error message: wrong usage'); end
+	if nargin~=2, error('love_numbers error message: wrong usage'); end
 
 	if ~( strcmpi(type,'h') | strcmpi(type,'k') | strcmpi(type,'l') | strcmpi(type,'gamma') | strcmpi(type,'lambda') ), 
 		error('type should be one of ''h'',''k'',''l'',''gamma'' and ''lambda''');
+	end
+
+	if ~( strcmpi(reference_frame,'CM') | strcmpi(reference_frame,'CF')), 
+		error('reference_frame should be one of ''CM'' or ''CF''');
 	end
 
@@ -10029,4 +10036,21 @@
 	elseif type=='lambda',
 		series=love_numbers(:,5);
-	else error(['love_numbers error message: unknow type:' type]);
+	else 
+		error(['love_numbers error message: unknow type:' type]);
 	end
+
+	% choose degree 1 term for CF reference system 
+	if reference_frame=='CM', 
+		return; 
+	elseif reference_frame=='CF', % from Blewitt, 2003, JGR 
+		if type=='h',
+			series(2,1) = -0.269; 
+		elseif type=='k',
+			series(2,1) = 0.021;  
+		elseif type=='l',
+			series(2,1) = 0.134; 
+		end
+	else 
+		error(['love_numbers error message: unknow reference frame:' reference_frame]);
+	end
+
Index: /issm/trunk-jpl/test/NightlyRun/test2002.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 21305)
+++ /issm/trunk-jpl/test/NightlyRun/test2002.m	(revision 21306)
@@ -20,7 +20,7 @@
 %elastic loading from love numbers: 
 nlov=101;
-md.slr.love_h = love_numbers('h'); md.slr.love_h(nlov+1:end)=[];
-md.slr.love_k = love_numbers('k'); md.slr.love_k(nlov+1:end)=[];
-md.slr.love_l = love_numbers('l'); md.slr.love_l(nlov+1:end)=[];
+md.slr.love_h = love_numbers('h','CM'); md.slr.love_h(nlov+1:end)=[];
+md.slr.love_k = love_numbers('k','CM'); md.slr.love_k(nlov+1:end)=[];
+md.slr.love_l = love_numbers('l','CM'); md.slr.love_l(nlov+1:end)=[];
 
 %}}}
Index: /issm/trunk-jpl/test/NightlyRun/test2101.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2101.m	(revision 21305)
+++ /issm/trunk-jpl/test/NightlyRun/test2101.m	(revision 21306)
@@ -12,6 +12,6 @@
 %love numbers:
 nlov=10001;
-md.esa.love_h = love_numbers('h'); md.esa.love_h(nlov+1:end)=[];
-md.esa.love_l = love_numbers('l'); md.esa.love_l(nlov+1:end)=[];
+md.esa.love_h = love_numbers('h','CM'); md.esa.love_h(nlov+1:end)=[];
+md.esa.love_l = love_numbers('l','CM'); md.esa.love_l(nlov+1:end)=[];
 
 %mask:  {{{
Index: /issm/trunk-jpl/test/NightlyRun/test2110.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2110.m	(revision 21305)
+++ /issm/trunk-jpl/test/NightlyRun/test2110.m	(revision 21306)
@@ -17,6 +17,6 @@
 %love numbers:
 nlov=10001;	% horizontal displacements do not work for low degree truncation, e.g., 101
-md.esa.love_h = love_numbers('h'); md.esa.love_h(nlov+1:end)=[];
-md.esa.love_l = love_numbers('l'); md.esa.love_l(nlov+1:end)=[];
+md.esa.love_h = love_numbers('h','CM'); md.esa.love_h(nlov+1:end)=[];
+md.esa.love_l = love_numbers('l','CM'); md.esa.love_l(nlov+1:end)=[];
 
 %mask:  {{{