  Conventional radially-local neoclassical calculations become inadequate if
the radial gradient scale lengths of the H-mode pedestal become as small as the
poloidal ion gyroradius. Here, we describe a radially global $\delta f$
continuum code that generalizes neoclassical calculations to allow stronger
gradients. As with conventional neoclassical calculations, the formulation is
time-independent and requires only the solution of a single sparse linear
system. We demonstrate precise agreement with an asymptotic analytic solution
of the radially global kinetic equation in the appropriate limits of aspect
ratio and collisionality. This agreement depends crucially on accurate
treatment of finite orbit width effects.
