A mathematical formulation for robust targeting in heat integrated water allocation network. (1st February 2023)