Application of copula functions to construct confidence intervals of bivariate drought frequency curve