added subroutine SetBoundaryValuesST_w_s(w_s) for density part of the solid phase, have to add(define)  Win_s, code and executing