A parallel domain decomposition algorithm for coastal ocean circulation models based on integer linear programming