#define IM		2
#define NW		(((lmax + 1) * (lmax + 2)) / 2)