Eddy parameterisation for primitive equation ocean models