1) generation map component
TQU_dust.pro
TQU_synchrotron.pro                                                                         the outputs are maps in MJy/sr
T_freefree.pro
(/data1/sara/ref/refsky/Galaxy)

2) coadding components and conversion MJy/sr ---> mK antenna temperature
cieli.pro                        (only diffuse components)                                                    
cieli_cmb.pro                (diffuse components + cmb)               
cieli_diff_sources.pro     (diffuse components+sources+cmb)                                      
(/data1/sara/ref/refsky/cieli/wmap_sky)

conversion MJy/sr ---> mK
nu=float(v)                          ; GHz
c=2.99792458e10                 ; cm / s
h=6.626176e-27                  ; erg * s
k=1.380662e-16                  ; erg / K
T_cmb=2.726                     ; K
conversion_K_MJy_sr=2.*k*(nu*1.e9/c)^2*1.e17
conversion_MJy_sr_mK=(1./conversion_K_MJy_sr)*1000.

Now the coadded maps are in the same units of WMAP maps (mK)