[l,Cs_l]=load_cmbfast(filename) Load cmbfast spectrum and normalize to uk^2 Note Cs_l is "C script l" To get C_l do C_l=(Cscript_l*2*pi)./(l.*(l+1));
0001 function [l,Cs_l]=load_cmbfast(filename) 0002 % [l,Cs_l]=load_cmbfast(filename) 0003 % 0004 % Load cmbfast spectrum and normalize to uk^2 0005 % Note Cs_l is "C script l" 0006 % To get C_l do C_l=(Cscript_l*2*pi)./(l.*(l+1)); 0007 0008 temp=load(filename); 0009 l=temp(:,1); 0010 Cs_l=(2.73^2)*(1e6^2)*[temp(:,2)]; 0011 0012 return