qgutils provides a couple of functions useful in the QG world and elsewhere.
python setup.py install --user- Suppose we have 2 layers of thickness 500 and 2500 m
- The Brunt-Vaisala frequency (squared) N^2 = 10^-5 s^-2
- the coriolis parameter is f0 = 10^-4 s^-1
import numpy as np
import qgutils as qg
dh = np.array([500,2500]) # thickness of the layers
N2 = np.array([1e-5])
Rd = qg.comp_modes(dh,N2,f0=1e-4)
print("First deformation radius = {0:.0f} m".format(Rd[1]))