A three-dimensional primitive equation, baroclinic numerical model incorporating a range of turbulence closure schemes is used to investigate the effects of vertical diffusion of momentum and density upon the spread of a freshwater plume, with particular reference to the Ebro plume. Initial calculations show that there are some differences in the horizontal spread and vertical mixing of the plume when diffusion coefficients are computed from a two-equation turbulence energy model compared with a one-equation model. To understand results from the turbulence energy models, the sensitivity of the plume dynamics to variations in the coefficient of vertical eddy viscosity and diffusivity is also considered, with increases in these parameters having a significant effect upon the cross-shore and along-shore spread of the plume. Also, increasing these parameters changes the plume characteristics from supercritical to subcritical and reduces the occurrence of meandering and baroclinic instability along the plume's off-shore edge. However, differences in the southerly spread (the direction of Kelvin wave propagation in the model) of the plume (although not its northerly spread) produced by changes in diffusion coefficients are small compared with the influence of changes in the bottom slope, upon the along-shore southerly spread of the plume, which moves in the direction of Kelvin wave propagation in the near coastal region. Results from the series of calculations are used as a guide in experimental design, with reference to a planned experiment in the Ebro region involving a coastal HF Radar deployment, as well as off-shore measurements. Calculations suggest that surface current measurements from a coastal HF Radar, together with a detailed survey of the density field associated with the plume, may be an appropriate, although indirect, means of determining suitable mixing coefficients to use in plume discharge problems. Detailed measurements of water depth variation will also be required.