This will require dealing with non-linear materials. I guess a saturable absorber could be defined around a specific wavelength as a 3 level system.

But frankly, avoid it if you can. While FDTD can simulate non-linear materials, such computations are very complicated to setup properly, if not almost impossible if time constants are above a few femtoseconds.