Changeset 25355

Timestamp:

08/06/20 17:17:11 (5 years ago)

Author:

adhikari

Message:

CHG: adding a switch that allows us to validate our solutions against the Spada benchmarks

File:

• issm/trunk-jpl/test/NightlyRun/test2084.m (modified) (3 diffs)

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

@@ -8 +8 @@
 md=model();
 md.cluster=generic('name',oshostname(),'np',1);
+
+% set validation=1 for comparing against the Spada benchark.   
+validation=0; 
 
 md.materials=materials('litho');
@@ -45 +48 @@
         };
 
+% validate elastic loading solutions against the Spada benchmark. {{{ 
+if validation 
+        spada_solutions = load('spada_elastic_loading_deg_h_l_k'); 
+        spada_d = spada_solutions(:,1); 
+        spada_h = spada_solutions(:,2); 
+        spada_l = spada_solutions(:,3); 
+        spada_k = spada_solutions(:,4); 
+
+        %rename ISSM solutions.  
+        issm_d = [1:md.love.sh_nmax];  
+        issm_h = md.results.LoveSolution.LoveHr(2:end,1); 
+        issm_l = md.results.LoveSolution.LoveLr(2:end,1); 
+        issm_k = md.results.LoveSolution.LoveKr(2:end,1); 
+
+        % relative difference for each degree, except for zero.  
+        diff_h = 1 - issm_h./spada_h;  
+        diff_l = 1 - issm_l./spada_l;  
+        diff_k = 1 - issm_k./spada_k;  
+
+        figure 
+        plot(spada_d,[diff_h diff_l diff_k]); grid on; 
+        legend('h','l','k'); title('loading love'); 
+
+else
+        % 
+end 
+% }}} 
 
 md.love.frequencies=([1e-3 1e-2 1e-1 1 -1e-3 -1e-2 -1e-1 -1]*2*pi)'/cst;
@@ -71 +101 @@
 md.love.nfreq=length(md.love.frequencies);
 md=solve(md,'lv');
+
+% validate elastic tidal solutions against the Spada benchmark. {{{ 
+if validation 
+        spada_solutions = load('spada_elastic_tidal_deg_h_l_k'); 
+        spada_d = spada_solutions(:,1); 
+        spada_h = spada_solutions(:,2); 
+        spada_l = spada_solutions(:,3); 
+        spada_k = spada_solutions(:,4); 
+
+        %rename ISSM solutions.  
+        issm_d = [2:md.love.sh_nmax];  
+        issm_h = md.results.LoveSolution.LoveHr(3:end,1); 
+        issm_l = md.results.LoveSolution.LoveLr(3:end,1); 
+        issm_k = md.results.LoveSolution.LoveKr(3:end,1); 
+
+        % relative difference for each degree, except for zero and one. 
+        diff_h = 1 - issm_h./spada_h;  
+        diff_l = 1 - issm_l./spada_l;  
+        diff_k = 1 - issm_k./spada_k;  
+
+        % k seems to have larger relative difference at higher degrees, although values approach zero. 
+        figure 
+        plot(spada_d,issm_k./spada_k);
+
+        figure 
+        plot(spada_d,[diff_h diff_l diff_k]); grid on; 
+        legend('h','l','k'); title('tidal love'); 
+else
+        % 
+end 
+% }}} 
 
 %tidal love numbers