Bimodal behavior in the zonal mean flow of a baroclinic beta-channel model