In the Colby sand load method down in SEDTRN, the log10 function returns NaN in this code, when using the Anaconda python 3.9.7. I've verified that it returns a good value in python 3.7, and it even returns a good number when I comment out the ‘njit’ declaration. It’s just the numba compiled version that returns the NaN. To make it even more perplexing, when I add the print statement below (which I did as I was trying to debug), it works just fine in all versions, even with numba.
In the Colby sand load method down in SEDTRN, the log10 function returns NaN in this code, when using the Anaconda python 3.9.7. I've verified that it returns a good value in python 3.7, and it even returns a good number when I comment out the ‘njit’ declaration. It’s just the numba compiled version that returns the NaN. To make it even more perplexing, when I add the print statement below (which I did as I was trying to debug), it works just fine in all versions, even with numba.