An interactive Geophysical Fluid Dynamics Textbook Notation

Ch 4. Scaling & the dimensionless numbers

Overview

Ch. 2 computed five dimensionless numbers from raw physical inputs. Two of them, RoRo and BuBu, do almost all the work of predicting which of this book’s reduced models applies at a given scale. This chapter turns that observation into a map: place any flow on the (Ro,Bu)(Ro,Bu) plane and read off directly whether it’s unbalanced, baroclinically active, or columnar/barotropic — and see exactly where every earlier chapter’s example flow actually lands.

The model

Notebook: notebooks/ch04_scaling.py → exported to /nb/ch04_scaling.html. Governing equations use the symbols in notation.

Ro=Uf0L ≷ 1,Bu=(LRL)2 ≷ 1,LR=NHf0.Ro=\frac{U}{f_0L}\ \gtrless\ 1,\qquad Bu=\left(\frac{L_R}{L}\right)^2\ \gtrless\ 1,\qquad L_R=\frac{NH}{f_0}.

gfdlib.scaling.classify_regime labels every (Ro,Bu)(Ro,Bu) point unbalanced / QG-baroclinic / QG-barotropic; the notebook renders it as a filled region map with five reference points from earlier chapters’ example flows placed directly on it.

Prefer to run it yourself? Download as a Jupyter notebook — it installs its own dependencies in the first cell. The sliders are static in Jupyter; edit value= and re-run to change a parameter.

Both fluids

Exercises

  1. (analytic) Show algebraically that Bu=(Ro/Fr)2Bu=(Ro/Fr)^2 (both are ratios built from the same U,L,H,N,f0U,L,H,N,f_0; note Fr/Ro=f0L/(NH)=L/LRFr/Ro=f_0L/(NH)=L/L_R), so the (Ro,Bu)(Ro,Bu) map is equivalent to an (Ro,Fr)(Ro,Fr) map — confirm this against the definitions in gfdlib.scaling.
  2. (computational) Modify ch04_scaling.py to add a sixth reference point for Ch. 12’s internal-wave beams, using that chapter’s typical N^,f^,ω^\hat N,\hat f,\hat\omega to estimate an effective Ro,BuRo,Bu (this requires some judgment — internal waves aren’t really a “balanced flow” scale in the same sense as the other five points, which is itself worth noting in your answer).
  3. (exploratory) Sweep the marker along the Bu=1Bu=1 line from very small to very large RoRo. At what RoRo does the marker cross into “unbalanced”? Does that crossing point depend on where along the Bu=1Bu=1 line you started? (It shouldn’t — re-read classify_regime’s logic to see why the RoRo boundary is checked independently of BuBu.)

Further reading

Vallis, AOFD 2nd ed., §2.7-2.8 (regime diagrams); Pedlosky, Geophysical Fluid Dynamics 2nd ed., §6.1-6.3; Cushman-Roisin & Beckers, Introduction to Geophysical Fluid Dynamics, ch. 4-5.