Changeset 21306
- Timestamp:
- 10/23/16 13:11:06 (8 years ago)
- Location:
- issm/trunk-jpl
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/boundaryconditions/love_numbers.m
r20181 r21306 1 function series=love_numbers(type) 2 %LOVE_NUMBERS: provide love numbers (type 'h','k','l','gamma' and 'lambda' 3 % retrieved from: http://www.srosat.com/iag-jsg/loveNb.php 4 % Usage: series=love_numbers(type) 5 % where type is one of 'h','k','l','gamma' and 'lambda'. 1 function series=love_numbers(type,reference_frame) 2 %LOVE_NUMBERS: provide love numbers 3 % retrieved from: http://www.srosat.com/iag-jsg/loveNb.php 4 % in a chosen reference frame 5 % 6 % Usage: series=love_numbers(type,reference_frame) 7 % type is one of 'h','k','l','gamma' and 'lambda'. 8 % reference_frame is one of 'CM' and 'CF'. 6 9 % 7 % 8 % love_k=love_numbers('k');10 % Example: 11 % love_k=love_numbers('k','CF'); 9 12 % 10 13 11 14 %some checks: 12 if nargin~= 1, error('love_numbers error message: wrong usage'); end15 if nargin~=2, error('love_numbers error message: wrong usage'); end 13 16 14 17 if ~( strcmpi(type,'h') | strcmpi(type,'k') | strcmpi(type,'l') | strcmpi(type,'gamma') | strcmpi(type,'lambda') ), 15 18 error('type should be one of ''h'',''k'',''l'',''gamma'' and ''lambda'''); 19 end 20 21 if ~( strcmpi(reference_frame,'CM') | strcmpi(reference_frame,'CF')), 22 error('reference_frame should be one of ''CM'' or ''CF'''); 16 23 end 17 24 … … 10029 10036 elseif type=='lambda', 10030 10037 series=love_numbers(:,5); 10031 else error(['love_numbers error message: unknow type:' type]); 10038 else 10039 error(['love_numbers error message: unknow type:' type]); 10032 10040 end 10041 10042 % choose degree 1 term for CF reference system 10043 if reference_frame=='CM', 10044 return; 10045 elseif reference_frame=='CF', % from Blewitt, 2003, JGR 10046 if type=='h', 10047 series(2,1) = -0.269; 10048 elseif type=='k', 10049 series(2,1) = 0.021; 10050 elseif type=='l', 10051 series(2,1) = 0.134; 10052 end 10053 else 10054 error(['love_numbers error message: unknow reference frame:' reference_frame]); 10055 end 10056 -
issm/trunk-jpl/test/NightlyRun/test2002.m
r21296 r21306 20 20 %elastic loading from love numbers: 21 21 nlov=101; 22 md.slr.love_h = love_numbers('h' ); md.slr.love_h(nlov+1:end)=[];23 md.slr.love_k = love_numbers('k' ); md.slr.love_k(nlov+1:end)=[];24 md.slr.love_l = love_numbers('l' ); md.slr.love_l(nlov+1:end)=[];22 md.slr.love_h = love_numbers('h','CM'); md.slr.love_h(nlov+1:end)=[]; 23 md.slr.love_k = love_numbers('k','CM'); md.slr.love_k(nlov+1:end)=[]; 24 md.slr.love_l = love_numbers('l','CM'); md.slr.love_l(nlov+1:end)=[]; 25 25 26 26 %}}} -
issm/trunk-jpl/test/NightlyRun/test2101.m
r21273 r21306 12 12 %love numbers: 13 13 nlov=10001; 14 md.esa.love_h = love_numbers('h' ); md.esa.love_h(nlov+1:end)=[];15 md.esa.love_l = love_numbers('l' ); md.esa.love_l(nlov+1:end)=[];14 md.esa.love_h = love_numbers('h','CM'); md.esa.love_h(nlov+1:end)=[]; 15 md.esa.love_l = love_numbers('l','CM'); md.esa.love_l(nlov+1:end)=[]; 16 16 17 17 %mask: {{{ -
issm/trunk-jpl/test/NightlyRun/test2110.m
r21273 r21306 17 17 %love numbers: 18 18 nlov=10001; % horizontal displacements do not work for low degree truncation, e.g., 101 19 md.esa.love_h = love_numbers('h' ); md.esa.love_h(nlov+1:end)=[];20 md.esa.love_l = love_numbers('l' ); md.esa.love_l(nlov+1:end)=[];19 md.esa.love_h = love_numbers('h','CM'); md.esa.love_h(nlov+1:end)=[]; 20 md.esa.love_l = love_numbers('l','CM'); md.esa.love_l(nlov+1:end)=[]; 21 21 22 22 %mask: {{{
Note:
See TracChangeset
for help on using the changeset viewer.