Ah.  I must have looked at that line dozens of times and never noticed the typo.

Adding `IMPLICIT NONE` would have caught it.  Apparently `IMPLICIT NONE` is not standard Fortran 77, but we're already using it in special/specfun/specfun.f.  Any objections to adding it to this subroutine?