The results of 2D fluid plasma simulations of the ion flow in the boundary plasma of DIII-D, JET, and JT-60U are reported. The model includes the effects of drifts and of impurity radiation using a multi-species model of intrinsic carbon impurities. Drift effects are important in determining the primary and impurity ion fluxes in the SOL and the private region, with ExB drifts dominant. Simulated parallel velocities are consistent with experimental measurement on the high field side of JT-60U, where the parallel flow is determined by- ion sources. Simulated parallel velocities are significantly less than seer. in experiment at the outer midplane of JT-60U, and at the top of JET where the Bow is sensitive to poloidal variations of the turbulence driven transport, suggesting the velocity in these regions is determined by transport. Parallel flows are reversed by changing the direction of the ion VB drift relative to the X-point.