“…Derivativesarechosenfrom appropriate subdomainsso that wave componentsare "upwinded". The other approachuses a correction procedure( [25], [26], [1], [4]). To implement the correction method, the interior point approximationis integratedeverywhere, including at the boundaries.As a result, multiple solution valuesare available at eachinterfacepoint.…”