Clean: refactor code from thermosoil_getdiff into its own subroutine Clean: properly define the argument soilc_total in the stomate part as out