Multiscale modelling of neutron star oceans