Changeset 19310

Timestamp:

04/23/15 15:54:25 (10 years ago)

Author:

adhikari

Message:

weighted area integration over the surface of a sphere

File:

• issm/trunk-jpl/src/m/psl/ocelm.m (modified) (1 diff)

issm/trunk-jpl/src/m/psl/ocelm.m

@@ -1 @@
-function oce_lm = ocelm(lMax,nPix,ocean,sh) 
+function oce_lm = ocelm(sh,lMax,ocean,areaEle) 
 
-%---------------------------------------------------------------------
-% ocelm :: a function to compute SH coefficients of "ocean function" 
-%---------------------------------------------------------------------
-% This code is written as a part of the ISSM-PSL project
-% (c) S. Adhikari 
-%     Jet Propulsion Laboratory, Caltech 
-%     November 3, 2014
-%---------------------------------------------------------------------
+%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');
+      oce_lm(1+p1) = sum((sh(:,1+p1).*ocean').*areaEle');
       p1 = p1+1;
    end
 end
+oce_lm = oce_lm/sum(areaEle);
 
-oce_lm = oce_lm*4*pi/nPix;
-