Dear mfb,
I think I found the solution. [itex]\partial{P}/\partial{t}[/itex] is some kind of current (bound charges are moving). So if we write total current density as [itex]J=J_f+\nabla\times{M}+J_p[/itex] where [itex]J_p[/itex] is polarization current density, we can solve the problem.
