Simplification : removed shumdiag calculated in the initialization phase. It is not needed. Only shumdiag_perma is needed in thermosoil_initialized.