The numerical implementation of an immersed boundary method in the atmospheric solver Meso-NH is presented. This technique models fluid–solid interaction and allows for the simulation of urban flows by considering buildings to be part of the resolved scales. This study constitutes a first robust step towards a better understanding of the interactions between weather and cities and better predictions of such interactions.