fixed pb0r - now gives results the same as the IDL version down to the 7th decimal place approximately.