How the math works
1. Wake = unit cone parented to craft, tip at origin, axis along craft-local +Z (behind).
2. Cone is scaled each frame: radius = a + b·speed, length = c + d·speed.
3. Each ocean vertex is transformed into the cone's unscaled local space via wake.matrixWorld-1.
4. Inside test: 0 < z < 1 AND √(x²+y²) < z (the unit cone).
5. While inside, inject a Mexican-hat force into hv: hv -= F·speed·(1 - 1.5·r/z)·(1 - z)·sat(h)·dt. Down on axis, UP at rim. sat = 900/(900+h²) throttles force in already-deep cells.
6. Klein-Gordon wave eq on ALL inner cells: hv = hv·e-λeffdt + (c²·∇²h/dx² - m²·h)·dt; h += hv·dt. -m²·h = per-cell spring (uniform ω=√m²). λeff = λ + 0.025·|h| = self-limiting damping.
7. Final height = ambient sines + h[i]. Stiffness sets bounce speed; wave speed sets propagation; damping sets persistence; saturation + |h|-damping prevent runaway at extreme settings.