hypsometric_variable_g_pref

exo_skryer.vert_alt.hypsometric_variable_g_pref(p_lev: Array, T_lay: Array, mu_lay: Array, params)

Compute an altitude profile with altitude-dependent gravity anchored at p_ref.

This scheme defines z = 0 at a reference pressure p_ref and integrates both upward and downward to fill the full level grid.

Like hypsometric_variable_g(), this implementation uses an exact (per-layer) analytic update for spherical gravity assuming layer-wise constant temperature and mean molecular weight:

g(z) = g_ref * (R0 / (R0 + z))**2

Notes

• Levels at pressures greater than p_ref end up with negative altitudes.

• p_ref is clipped to lie within the provided pressure grid.