Index: /issm/trunk-jpl/src/m/boundaryconditions/love_numbers.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/love_numbers.py	(revision 22352)
+++ /issm/trunk-jpl/src/m/boundaryconditions/love_numbers.py	(revision 22353)
@@ -7,18 +7,26 @@
 #			   retrieved from: http://www.srosat.com/iag-jsg/loveNb.php
 #    Usage:   series=love_numbers(value) 
+#             series=love_numbers(value,reference_frame) 
+# 
 #             where value is one of 'h','k','l','gamma' and 'lambda'. 
+#	      reference_frame = one of 'CM' (default) and 'CF'.
 #
 #    Example:  
 #          love_k=love_numbers('k');
+#          love_k=love_numbers('k','CF');
 # 
 
-	#some checks:
-	if len(varargin)>0:
-		raise RuntimeError('love_numbers error message: wrong usage')
-
-	if value not in ['h','k','l','gamma','lambda']:
-		raise RuntimeError('value should be one of ''h'',''k'',''l'',''gamma'' and ''lambda''')
-
-	love_numbers=np.array([[    0         , 0          ,0          ,0          ,0          ,0          ,0          ],
+        # some checks:
+        if varargin==1:
+            frame='CM'; 
+        elif varargin==2: 
+            frame=reference_frame;
+        else:
+            raise RuntimeError('love_numbers error message: bad usage') 
+        
+        if value not in ['h','k','l','gamma','lambda']:
+            raise RuntimeError('value should be one of ''h'',''k'',''l'',''gamma'' and ''lambda''') 
+        
+        love_numbers=np.array([[    0         , 0          ,0          ,0          ,0          ,0          ,0          ],
 												 [	-1.28740059,-1.00000000,-0.89858519,1.28740059, 0.42519882  ,0.89858519 ,0.00000000 ],
 												 [	-1.00025365, -0.30922675, 0.02060926, 1.69102690, 0.46358648, 0.67016399, 0.61829668],
@@ -10021,7 +10029,6 @@
 												 [	-6.27342771, -0.00030948, 0.00018957, 7.27311822, 0.99905470, 0.99950094, 0.49327205],
 												 [	-6.27342778, -0.00030945, 0.00018956, 7.27311833, 0.99905480, 0.99950099, 0.49327194]]);
-
-
-	if value=='h':
+                                                                                                 
+        if value=='h':
 		series=love_numbers[:,0];
 	elif value=='k':
@@ -10035,4 +10042,14 @@
 	else:
 		raise RuntimeError(['love_numbers error message: unknow value:', value])
-		
-	return series
+	
+	% choose degree 1 term for CF reference system 
+        if frame=='CF': # from Blewitt, 2003, JGR 
+            if value=='h':
+                    series[2,1] = -0.269; 
+            elif value=='k':
+                    series[2,1] = 0.021;  
+            elif value=='l':
+                    series[2,1] = 0.134; 
+	
+        return series
+