Index: /issm/trunk-jpl/test/NightlyRun/test2084.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2084.m	(revision 25354)
+++ /issm/trunk-jpl/test/NightlyRun/test2084.m	(revision 25355)
@@ -8,4 +8,7 @@
 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,4 +48,31 @@
 	};
 
+% 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,4 +101,35 @@
 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