A high-order accuracy flow solver for high-Reynolds-number stratified flows