Computes surface heat flux from solar geometry, terrain slope, land
cover (ESA WorldCover), and atmospheric transmittance. The ground
temperature anomaly is resolved per-cell and converted to a thermal
strength index — fast enough to run a full day of hourly snapshots
in seconds.
Outputs: surface temperature anomaly · thermal
index · shadow mask · sun elevation
Full 3D incompressible Navier-Stokes solver under the Boussinesq
approximation: solar heating warms the ground, which injects
buoyancy into the boundary layer, driving convective thermal plumes.
Smagorinsky LES turbulence closure and a Red-Black SOR pressure
solver. Runs in a Web Worker; takes minutes for a coarse grid.
Outputs: 3D temperature anomaly θ' · vertical wind
w · convergence · cross-section slices
Solves the 2D shallow water equations (SWE) with reduced gravity
over real terrain. The atmosphere is modelled as a single layer of
depth h flowing over the orographic floor. Terrain-blocking,
orographic lifting, Froude number transitions, hydraulic jumps, and
vorticity wakes emerge naturally. Runs to quasi-steady state in
seconds to minutes.
3D Navier-Stokes · log profile · land-cover roughness
Full 3D incompressible Navier-Stokes solver for mechanical wind
flow. Terrain cells are blocked (no-slip); inflow uses a logarithmic
wind profile anchored to a per-cell roughness length z₀ derived from
ESA WorldCover. Semi-Lagrangian advection and a Jacobi pressure
Poisson solver. Resolves lee-side recirculation, flow separation,
and surface-layer structure.