Refactor: move thermosoil_coef into its own module thermosoil_coefficients. Add all missing arguments to the subroutinte header