Index: sm/trunk-jpl/src/m/psl/icelm.m
===================================================================
--- /issm/trunk-jpl/src/m/psl/icelm.m	(revision 19323)
+++ 	(revision )
@@ -1,24 +1,0 @@
-function ice_lm = icelm(lMax,nPix,ice,sh) 
-
-%---------------------------------------------------------------------
-% icelm :: a function to compute SH coefficients of "ice load function" 
-%---------------------------------------------------------------------
-% This code is written as a part of the ISSM-PSL project
-% (c) S. Adhikari 
-%     Jet Propulsion Laboratory, Caltech 
-%     November 3, 2014
-%---------------------------------------------------------------------
-
-p1 = 0; 
-
-ice_lm = zeros(1,(lMax+1)^2); 
-
-for l=0:lMax
-   for m=-l:l
-      ice_lm(1+p1) = sum(sh(:,1+p1).*ice');
-      p1 = p1+1;
-   end
-end
-
-ice_lm = ice_lm*4*pi/nPix; 
-
Index: sm/trunk-jpl/src/m/psl/ocelm.m
===================================================================
--- /issm/trunk-jpl/src/m/psl/ocelm.m	(revision 19323)
+++ 	(revision )
@@ -1,25 +1,0 @@
-function oce_lm = ocelm(sh,lMax,ocean,areaEle) 
-
-%OCElm :: a function to compute SH coefficients of ocean function 
-% 
-%USAGE: oce_lm = ocelm(sh,lMax,ocean,areaEle); 
-%
-%sh (spherical harmonics) 
-%lMax (maximum SH degree to be considered) 
-%ocean (ocean function) 
-%areaEle (area of elements) 
-%
-
-p1 = 0; 
-
-oce_lm = zeros(1,(lMax+1)^2); 
-
-% weighted area integration 
-for l=0:lMax
-   for m=-l:l
-      oce_lm(1+p1) = sum((sh(:,1+p1).*ocean').*areaEle');
-      p1 = p1+1;
-   end
-end
-oce_lm = oce_lm/sum(areaEle);
-