added parameter for regularization around delta_h=0 in thermostatic3WayValve