Modelling the turbulent magnetized ISM