Index: ../trunk-jpl/src/m/psl/shlm_element.m
===================================================================
--- ../trunk-jpl/src/m/psl/shlm_element.m	(revision 20330)
+++ ../trunk-jpl/src/m/psl/shlm_element.m	(revision 20331)
@@ -24,5 +24,3 @@
    end
 end 
 sh_lm = sh_lm*4*pi/sum(areas);
-
-
Index: ../trunk-jpl/src/m/psl/shlm.m
===================================================================
--- ../trunk-jpl/src/m/psl/shlm.m	(revision 20330)
+++ ../trunk-jpl/src/m/psl/shlm.m	(revision 20331)
@@ -14,18 +14,15 @@
 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);