Important Note:

in April 2008 G. Mellor found a bug in profq that has been 
fixed in pom2k.f (pom2k_old.f is previous version).

The change concerns boundary conditions in the loops after "...= -q2lb"
(correct loop is k=3,kbm1, not k=2,kbm1)