v4.1.0

New Features

• 5th-order upwind advection (Wicker & Skamarock 2002) for horizontal fluxes in 2D and 3D, configurable via advection_scheme: "upwind5" — provides implicit numerical dissipation, eliminating the need for explicit horizontal diffusion
• Beljaars (1994) convective velocity scale in MOST surface layer — prevents flux collapse in free-convective (low wind) conditions

Bug Fixes

• Fix MOST surface flux collapse when mean wind is weak by incorporating M_eff = sqrt(M² + (1.2·w*)²)
• Fix TKE over-initialization: use tke_min cold-start instead of convective scaling profile that caused excessive initial dissipation
• Increase theta perturbation amplitude to 0.1 K (standard LES practice)
• Set theta_surface to 307 K for consistency with target heat flux (~0.24 K m/s)

Other

• Add minimal 2D config (cbl_2d_minimal.yaml) for isolated physics testing