Dear Stephen,
Running the script that you posted,the program returns this error, "slope of the log calculation is outside the boundaries of natural diffusion"
Can you please provide some explanation or a reference for the boundary values that you listed in the script (0.9 to 1.10).

Thanks in advance.

Phani