In atmospheric physics we use often a variety of vertical coordinates, $z$ the height above some reference surface, or $r$ as distance from the planetary center.
However sometimes it simplifies the equations we want to solve by using other coordinates, Pressure $P$, LogPressure, Potential temperature $\theta$,...
One vertical coordinate I want to specifically ask about is the shallow water height. It is defined by looking at the hydrostatic stability approximation $$ \frac{\partial P}{\partial z} = \rho_0 g$$ then integrating over a variable height $h(x,y)$ so that we arrive at constant pressure $P(h)$ we get easily $$P(z)=\rho_0 g (h-z)+P(h) $$
So $h(x,y)$ is now the variable height between two Iso-Pressure surfaces in some planetary fluid.
Now we can express the navier-stokes-equation, mass equation, thermodynamic equation in terms of $h$. I guess this will be very familiar to those who could answer me.
My question is now:
- Assuming no vertical mixing (=geostrophic movement on isobaric surfaces), can I now build a vertical model of the planet consisting of a multitude of shallow fluid-layers? Namely I'd take $h1(x,y)$ as between $p1$ and $p2$, $h2(x,y)$ between $p2$ and $p3$,... and thus build up and atmosphere of a gas giant between $\sim 10^{-1} bar$ and several $10^2$ bars.
- And more important: when could this approach be totally wrong?
Naively I'd expect this to be OK in my case, as I'm still effectively integrating over the whole planet, but I hoped for some experts opinion on this.