The non-homogeneous bivariate compound Poisson process with short term periodic intensity function is used for modeling the events with seasonal patterns or periodic trends. In this paper, this process is carefully introduced. In order to characterize the dependence structure between jumps, the Levy copula function is provided. For estimating the parameters of the model, the inference for margins method is used. As an application, this model is fitted to an automobile insurance dataset with inference for margins method and its accuracy is compared with the full maximum likelihood method. By using the goodness of fit test, it is confirmed that this model is appropriate for describing the data.