A numerical model of the zonal mean circulation of the middle atmosphere