I coded a Matlab script to evalue the analitic solution of the temperature distribution reported on the paper of Jackson et al.

I show the plot of T along z (depth) and r (radius)

Laser power = 1W

Surface:

absorption = 12ppm;

coating thickness = 10micron;

Bulk:

absorption = 20ppm/cm

thickness = 5mm;