The following FORTRAN subroutine implements the approximate formula for the grain-assisted ion removal rate coefficient in Weingartner & Draine (2001, ApJ, 563, 842), equation 8: get_alpha_g.f