Numerical modelling of the thermosphere-ionosphere-protonosphere system