Nonlinear hydromagnetic convection in a horizontal layer of fluid rotating about the vertical axis is investigated using the mean field approximation. The boundary layer method is used assuming large Rayleigh number R for different ranges of the Chandrasekhar number Q and the Taylor number Ta. The heat flux F is determined for the value of the horizontal wave number α which maximizes F. It is shown that, for certain regions of the parameter space (R, Q, Ta), F and α change discontinuously for Ta greater than some critical value (given R and Q). Thus, for Ta about this critical value, wave numbers and heat fluxes of two different values will be predicted simultaneously. Also, for certian regions of the parameter space, the field can facilitate convection, but rotation can facilitate convection only for sufficiently large Ta.