For the second case you want to divide by the force, not the velocity.
You are basically computing what fraction of the time you spend at a particular point in phase space.
However what you have is a probability density. So $P(z,p)$ is something you multiply by a volume in phase space to get a probability.
A correct way to get it would be to consider a whole period (e.g. time from top to top or from bottom to bottom or from any point in phase space to itself).
A Dirac delta has units. It has units of one over the argument so since $\int \delta (x) dx=1$ then $\delta (x)$ has units of inverse distance and $\int \delta (p) dp=1$ then $\delta (p)$ has units of inverse momentum.
When you integrate over momentum you get a probability density for position. This has units that you multiply by a distance to get a probability of being in that range of distance.
These probability densities over positions are the ones where it is inversely proportional to velocity. So indeed after you integrate over momentum you have a probability per unit length and probabilities per unit length are inversely proportional to velocity. So it should work.
However if you integrate over position then you have a probability density for momentum. This is something that you multiply by a range of momentums to get a probability of having momentum in that range. There is no rule that this kind of a density is inversely proportional to velocity. And it isn't, it is inversely proportional to force.
This is similar to how energy density is different for frequency and for wave length because each is something you multiply by a range of things to get the actual amount of stuff. See my answer at https://physics.stackexchange.com/a/165256 if you need more details about densities in different united domains.
Let's review where that inversely proportional to velocity rule came from. As a simplification during a time interval $\Delta t$ you might travel a distance $\Delta z$ = $v\Delta t$ And really the probability you are computing is $\Delta t/T$ where $T$ is that period I mentioned.
So the probability you are in a region of width $\Delta z$ is proportional to the amount of time $\Delta t$ you spend in there. Thus inversely proportional to velocity, in fact it is $\Delta z/(vT).$ So the spatial density is $1/(vT).$ And this rule doesn't apply if you are asking how much time something spends in a region of size $\Delta p.$ So don't divide by velocity if you want a probability density of being within a range of momentums of width $\Delta p.$
To do the probability density for momentum intervals, use $\Delta p$ = $F\Delta t$ then the probability is $\Delta t/T$ = $\Delta p/(FT)$ so the probability density for momentum intervals is $1/(FT).$
So a density is just like the Dirac delta it is something you integrate to get the probability but it needs to be a per interval of distance or a per interval of momentum for it to be something you can integrate.
And the momentum version should give you pause since you generally need to correctly account for the force when it turns around at the ground. So over counting is an issue to consider.