cellZone: example coded function object for generating a cellZone.

The example is provided in etc/caseDicts/functionTemplates, with other
example coded function objects.

test/Lagrangian/coupling: Switch to fan injection

This change means the only effect spreading the particles out of the
plane is the turbulent dispersion model. This makes it easier to see the
effect of that model and switch it on and off and otherwise confirm its
behaviour.

stackPlume: updated for two-way Lagrangian coupling

Created tag for OpenFOAM-13

zoneSet,zoneGenerators::lookup: Improved zone access checking and dia…

…gnostic messages

If the requested zone with the specified zoneType is not available lookup now
prints if it available with a different type, e.g. if in the ballValve
createZonesDict the lookup of the ball is specified to be of faceZone type:

ball
{
    type        point;
    ball
    {
        zoneType face;
    }
}

the error message

--> FOAM FATAL IO ERROR:
Cannot find faceZone ball
    Found cellZone ball
    Available faceZones: 0()

file: .../tutorials/incompressibleFluid/ballValve/system/createZonesDict/ball at line 21.

is printed.

etc/caseDicts/functions/surfaceFieldValue/patchValue.cfg: Removed red…

…undant patch entry

Merge branch 'master' of github.amrom.workers.dev-OpenFOAM:OpenFOAM/OpenFOAM-dev

tutorials/compressibleVoF/climbingRod/system/setFieldsDict: Corrected

BasicThermo: Removed BasicThermo::mixture method

This method is considered confusing as it does not provide public
interface to the mixture/composition in way that the old
basicSpecieMixture did. This access is now provided directly by the
multicomponentThermo and its derivations.

multiphaseEuler: nAlphaSubCycles is now a Function1 of the alpha Cour…

…ant number

The explicit MULES sub-cycling control nAlphaSubCycles is now a Function1 of the
alpha Courant number which allows a very convenient way to specify the number of
sub-cycles independent of the time-step the case is run with while maintaining
stability of the explicit solution of the phase-fraction.  For example with the
following entry in fvSolution of the pipeBend case:

solvers
{
    "alpha.*"
    {
        nAlphaSubCycles polynomial (0 3);
    }
.
.
.

nAlphaSubCycles is evaluated as three times the current maximum alpha Courant
number rounded upwards to an integer, i.e. given the maximum Courant number is
set to 1 for this case 3 sub-cycles will be executed.