Hello ewang,
I after all think it is totally normal that you get a much higher computation time when you couple Telemac with Tomawac. Indeed, the propagation step is solved in TOMAWAC by means of the method of characteristics. When there is no current and the depth is constant over time, the characteristics can be traced back only once, at the beginning of the simulation. But when the current is not stationary, as in your case, the characteristics should be traced back at each time step when the current is updated. This takes some of the computation time...
Regards,
Elodie