High-performance 2D solver for Kelvin-Helmholtz instability in incompressible flows, featuring Numba JIT compilation for 10-50x speedup and multi-core parallelization. Includes predefined scenarios, NetCDF output, automated visualization, and both CLI and Python API for the geophysical fluid dynamics studies.