uk2jysr

 sfs=uk2jysr(freqs)

 Calc scale factors from uK to Jy/Sr assuming RJ conversion formula
 I = 2*k*(nu/c)^2 * T

 For CMB maps even though the temp diffs are small this is presumably
 not that good an approximation at higher frequencies - what one
 cares about is the error with respect to the Planck formula for
 T+deltaT1 vs T+deltaT2

 Accurate replacement is:
 delI=planck_dIdT(nu,T)*delT*1e26*1e-6
 where 1e26 is W->Jy and 1e-6 is K->uK

0001 function sfs=uk2jysr(freqs)
0002 % sfs=uk2jysr(freqs)
0003 %
0004 % Calc scale factors from uK to Jy/Sr assuming RJ conversion formula
0005 % I = 2*k*(nu/c)^2 * T
0006 %
0007 % For CMB maps even though the temp diffs are small this is presumably
0008 % not that good an approximation at higher frequencies - what one
0009 % cares about is the error with respect to the Planck formula for
0010 % T+deltaT1 vs T+deltaT2
0011 %
0012 % Accurate replacement is:
0013 % delI=planck_dIdT(nu,T)*delT*1e26*1e-6
0014 % where 1e26 is W->Jy and 1e-6 is K->uK
0015 
0016 sfs=1e-6*1e26*2*1.38e-23*freqs.^2./3e8^2;