subroutine surface_perturbation::cdfbet	(	real(kind=kind_phys), intent(in)	x,
		real(kind=kind_phys), intent(in)	p,
		real(kind=kind_phys), intent(in)	q,
		real(kind=kind_phys), intent(in)	eps,
		integer	iflag,
		real(kind=kind_phys), intent(out)	cdfx
	)