Simulation of saturated and unsaturated flow in karst systems at catchment scale using a double continuum approach