ps2cmbspec

 [l,Cs_l]=ps2cmbspec(u,C_u)

 Convert output of powspec to Cscript_l form

 e.g.: ad=calc_ad(20,256);
       [l,Cs_l,C_l] = get_cmb_model(0.05,0.35,0.60,65,1,0);
       cmb=cmb_skysim(ad,l,C_l);
       [C_u,u]=powspec(ad,i2f(ad,cmb));
       [lo,Cs_lo]=ps2cmbspec(u,ad.del_u,C_u);
       plot(l,Cs_l,'r');
       hold on; plot(lo,Cs_lo,'b'); hold off

0001 function [l,Cs_l]=ps2cmbspec(u,del_u,C_u)
0002 % [l,Cs_l]=ps2cmbspec(u,C_u)
0003 %
0004 % Convert output of powspec to Cscript_l form
0005 %
0006 % e.g.: ad=calc_ad(20,256);
0007 %       [l,Cs_l,C_l] = get_cmb_model(0.05,0.35,0.60,65,1,0);
0008 %       cmb=cmb_skysim(ad,l,C_l);
0009 %       [C_u,u]=powspec(ad,i2f(ad,cmb));
0010 %       [lo,Cs_lo]=ps2cmbspec(u,ad.del_u,C_u);
0011 %       plot(l,Cs_l,'r');
0012 %       hold on; plot(lo,Cs_lo,'b'); hold off
0013 
0014 C_l=C_u*del_u^2;
0015 l=u*2*pi;
0016 Cs_l=C_l.*l.*(l+1)/(2*pi);
0017