alm_i2t returns a real or complex array, containing the with and . The negative m are therefore ignored.
fits2alm, i1, a1, 'alm1.fits' | |
ac1 = alm_i2t(i1, a1, /complex, lmax=100, mmax=100) | |
fits2alm, i2, a2, 'alm2.fits' | |
ac2 = alm_i2t(i2, a2, /complex, lmax=100, mmax=100) | |
ac = 0.9*ac1 + 0.1*ac2 | |
alm_t2i, ac, i, a | |
alm2fits, i, a, 'almsum.fits' |
The example above reads 2 sets of from FITS files, puts the alm's with in tabular arrays, and then make a weighted sum of the alm's. The resulting alm or put back into a indexed list in order to be written to FITS.
Version 3.82, 2022-07-28