  Integrable PDEs on the line can be analyzed by the so-called Inverse
Scattering Transform (IST) method. A particularly powerful aspect of the IST is
its ability to predict the large $t$ behavior of the solution. Namely, starting
with initial data $u(x,0)$, IST implies that the solution $u(x,t)$ asymptotes
to a collection of solitons as $t \to \infty$, $x/t = O(1)$; moreover the
shapes and speeds of these solitons can be computed from $u(x,0)$ using only
{\it linear} operations. One of the most important developments in this area
has been the generalization of the IST from initial to initial-boundary value
(IBV) problems formulated on the half-line. It can be shown that again $u(x,t)$
asymptotes into a collection of solitons, where now the shapes and the speeds
of these solitons depend both on $u(x,0)$ and on the given boundary conditions
at $x = 0$. A major complication of IBV problems is that the computation of the
shapes and speeds of the solitons involves the solution of a {\it nonlinear}
Volterra integral equation. However, for a certain class of boundary
conditions, called linearizable, this complication can be bypassed and the
relevant computation is as effective as in the case of the problem on the line.
Here, after reviewing the general theory for KdV, we analyze three different
types of linearizable boundary conditions. For these cases, the initial
conditions are: (a) restrictions of one and two soliton solutions at $t = 0$;
(b) profiles of certain exponential type; (c) box-shaped profiles. For each of
these cases, by computing explicitly the shapes and the speeds of the
asymptotic solitons, we elucidate the influence of the boundary.
