Refactor: new subroutine thermosoil_apparent_ground_heat_flux upgraded from thermosoil_coef
There are two implementatio from thermosoil_apparent_ground_heat_flux: