source: issm/oecreview/Archive/20545-21336/ISSM-21306-21307.diff@ 21337

Index: ../trunk-jpl/src/m/boundaryconditions/love_numbers.m
===================================================================
--- ../trunk-jpl/src/m/boundaryconditions/love_numbers.m        (revision 21306)
+++ ../trunk-jpl/src/m/boundaryconditions/love_numbers.m        (revision 21307)
@@ -3,26 +3,37 @@
 %      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'.
+%      Usage: 
+%      series = love_numbers(type) 
+%      series = love_numbers(type,reference_frame) 
+% 
+%      - type = one of 'h','k','l','gamma' and 'lambda'. 
+%      - reference_frame = one of 'CM' (default) and 'CF'.
 %
 %      Example:  
-%      love_k=love_numbers('k','CF');
+%      love_k = love_numbers('k');
+%      love_k = love_numbers('k','CF');
 % 
 
-       %some checks:
-       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
-
+% some checks:
+if nargin==1, 
+       frame='CM'; 
+       disp('Info: computation is done in Center of Mass (CM) reference frame'); 
+elseif nargin==2,
+       frame=reference_frame;
        if ~( strcmpi(reference_frame,'CM') | strcmpi(reference_frame,'CF')), 
                error('reference_frame should be one of ''CM'' or ''CF''');
        end
+else
+       help love_numbers
+       error('bad usage'); 
+end
 
-       love_numbers=[...
+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
+
+               love_numbers=[...
      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;
@@ -10040,9 +10051,9 @@
        end
 
        % choose degree 1 term for CF reference system 
-       if reference_frame=='CM', 
+       if frame=='CM', 
                return; 
-       elseif reference_frame=='CF', % from Blewitt, 2003, JGR 
+       elseif frame=='CF', % from Blewitt, 2003, JGR 
                if type=='h',
                        series(2,1) = -0.269; 
                elseif type=='k',
@@ -10051,6 +10062,6 @@
                        series(2,1) = 0.134; 
                end
        else 
-               error(['love_numbers error message: unknow reference frame:' reference_frame]);
+               error(['love_numbers error message: unknow reference frame:' frame]);
        end