source: issm/oecreview/Archive/19101-20495/ISSM-20330-20331.diff@ 20498

../trunk-jpl/src/m/psl/shlm_element.m

@@ -24 +24 @@
    end
 end 
 sh_lm = sh_lm*4*pi/sum(areas);
-
-

../trunk-jpl/src/m/psl/shlm.m

@@ -14 +14 @@
 p = 0; 
 
 sh_lm = zeros((lMax+1)^2,1); 
-%sh_lm = zeros(1,(lMax+1)^2); 
 
 % weighted area integration over the unit sphere 
 for l=0:lMax
    for m=-l:l
-      %func0 = sh(:,1+p).*(func<1); 
       func0 = sh(:,1+p).*func; 
       sh_lm(1+p) = sum(mean(func0(index),2).*areas); 
       p = p+1;
    end
 end 
-%sh_lm = sh_lm/sum(areas);
 sh_lm = sh_lm*4*pi/sum(areas);