  We describe a fast method to simulate in a semi-analytical way consistent
maps of the thermal, kinetic and polarised Sunyaev-Zel'dovich effect, featuring
both cluster spatial correlations and large scale velocity flows.
