Subgrid flux model for scalars

Currently, for Large Eddy Simulations, subgrid flux of scalars are obtained using a constant subgrid Prandtl number (default value is 1). It would be interesting to have more options as scalar fluctuations are not always similar to velocity fluctuations.

(transferred from Mantis bug tracker issue 215, opened by Mathtrack).