Home > src > main > matlab > gen_geomodel.m

gen_geomodel

PURPOSE ^

SYNOPSIS ^

function f = gen_geomodel( L, P )

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function f = gen_geomodel( L, P )
0002 
0003 maxdepth = 3000;
0004 
0005 nodes = slag_sampling(P, maxdepth);
0006 
0007 for r = 0:P-1
0008     depth = maxdepth - nodes(r+1);
0009     cosimod = interpJRmodel(depth,'xpcube');
0010     [alm_trunc, L_trunc] = cosi2alm( cosimod );
0011     if r == 0
0012         f = zeros(P, L, 2*L-1);
0013     end
0014     alm = zeros(L^2, 1);
0015     if L_trunc < L
0016         alm(1:L_trunc^2, 1) = alm_trunc(:, 1);
0017     else
0018         alm(:, 1) = alm_trunc(1:L^2, 1);
0019     end
0020     temp = ssht_inverse(alm, L);
0021     f(r+1,:,:) = real(temp);
0022 end
0023 
0024 end

Generated on Mon 24-Sep-2012 12:26:33 by m2html © 2005