Subroutine "disnrg" has been used in "nmrred.f
", k2 is the "rk2" and k3 is the "rk3", so in the
program line it seems you didn't think about "1/2",
and I'm not sure whether my understanding is correct.