Changeset 21306

Timestamp:

10/23/16 13:11:06 (9 years ago)

Author:

adhikari

Message:

CHG: love numbers added for center-of-figure reference frame

Location:

issm/trunk-jpl

Files:

• src/m/boundaryconditions/love_numbers.m (modified) (2 diffs)
• test/NightlyRun/test2002.m (modified) (1 diff)
• test/NightlyRun/test2101.m (modified) (1 diff)
• test/NightlyRun/test2110.m (modified) (1 diff)

issm/trunk-jpl/src/m/boundaryconditions/love_numbers.m

@@ -1 @@
-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 +10036 @@
         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
+

issm/trunk-jpl/test/NightlyRun/test2002.m

@@ -20 +20 @@
 %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)=[];
 
 %}}}

issm/trunk-jpl/test/NightlyRun/test2101.m

@@ -12 +12 @@
 %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:  {{{

issm/trunk-jpl/test/NightlyRun/test2110.m

@@ -17 +17 @@
 %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:  {{{