subroutine, public surface_perturbation::ppfbet	(	real(kind=kind_phys), intent(in)	pr,
		real(kind=kind_phys), intent(in)	p,
		real(kind=kind_phys), intent(in)	q,
		integer	iflag,
		real(kind=kind_phys), intent(out)	x
	)