Simulating realistic disc galaxies with a novel sub-resolution ISM model