new function `pth_color_sfcsat()`

[ijlyttle]

get a color using a surface, a luminance, and a saturation:

pth_color_sfcsat <- function(surface, luminance, saturation = 1)