TorusOcean

TorusOcean is a regional 2D/3D ocean model implementing nonlinear shallow
water and primitive equations in sigma/S coordinates with ROMS-style numerics:
Arakawa C-grid, mode-split barotropic/baroclinic solver, Adams-Bashforth 3rd
order time stepping, and optional GPU compute via wgpu (Vulkan/Metal/D3D12).

Features:
- Nonlinear shallow water equations (2D barotropic)
- Primitive equations in sigma/S coordinates (3D baroclinic)
- Arakawa C-grid staggering
- Mode-split time stepping (Adams-Bashforth 3rd order)
- Multithreaded via Rayon
- JSON-based configuration
- CSV and NetCDF output (NetCDF requires system libnetcdf)