For turbulence, we use the 1.5-order TKE (Turbulent Kinetic Energy) closure, solving a prognostic equation for TKE and diagnosing eddy diffusivity ( K_m = C_k l \sqrte ). The constant (C_k) is tuned —not derived. This empirical calibration is not a weakness; it is analogous to the Moody chart for pipe friction.
For engineers, the "magic" of climate modeling lies in numerical methods. The Earth is divided into a three-dimensional grid. The challenge is balancing with computational cost .
Managing petabytes of output data (NetCDF being the industry standard format).
(2014) by John B. Drake, a former researcher at Oak Ridge National Laboratory.