EnKF-C provides a compact generic framework for off-line data assimilation (DA) into large-scale layered geophysical models with the ensemble Kalman filter (EnKF). Following are its other main features:
-
coded in C for GNU/Linux platform;
-
model-agnostic;
-
can conduct DA in EnKF, ensemble optimal interpolation (EnOI), or hybrid EnKF/EnOI modes;
-
permits multiple model grids;
-
can handle rectangular, curvilinear, or unstructured horizontal grids, z, sigma or hybrid vertical grids.
EnKF-C is coded for simplicity, scalability and robustness. To handle as large systems as possible it uses shared memory capabilities of MPI-3. Here is a snapshot of ensemble spread of sea surface temperature from the 96-member EnKF ocean forecasting system with MOM5 based OFAM3 model (51 x 1500 x 3600 grid), assimilating about 14M super-observations at each 3-day cycle.
For more information see README and user guide. (An older version of the user guide is also available from arXiv.) Have a feel for how the code works by running the included example.
Checkout EnKF-C by running git clone https://github.com/sakov/enkf-c.