The software constitues a two-layer quasi-gesotrophic basin ocean model solving for the flow in each layer described through the layer streamfunctions. The model solves the full non-linear (non-linearity can be turned off) quasi-geostrophic potential vorticity equation in each layer with no-normal-flow boundary conditions. There is a bottom topography (can be specified by user) and Rayleigh friction on the lower layer relative vorticity. There is also support for upper layer forcing (e.g. wind stress). The model uses third order Adam-Bashforth time-stepping for the potential vorticity and sine transforms for inverting PV to streamfunctions. A 2-D Quadratic Upwind Interpolation for Convective Kinematics (QUICK) scheme is used for advection. The flow is solved in terms of the barotropic an baroclinic streamfunctions, but outputted as the layer streamfunctions. Model originally written by Joseph H. LaCasce in Fortran 77. Model updated to Fortran 2003 and implemented topography by Jostein Brændshøi as part of a master project supervised by LaCasce in 2017/2018. See README.txt for more details on the model.