Multimodal and asymmetric bivariate circular data arise in several different disciplines and fitting appropriate distribution plays an important role in the analysis of such data. In this paper, we propose a new bivariate circular distribution which can be used to model both asymmetric and multimodal bivariate circular data simultaneously. In fact the proposed density covers unimodality as well as multimodality, symmetry as well as asymmetry of circular bivariate data. A number of properties of the proposed density are presented. A Bayesian approach with MCMC scheme is employed for statistical inference. Three real datasets and a simulation study are provided to illustrate the performance of the proposed model in comparison with alternative models such as finite mixture Cosine model.