This implements a discretized bimolecular recombination operator
Parameters: I V Kr 
I is a field or function
V is a field or function
Kr is a variable
is a variable, field, or function
is a variable, field, or function
Note, if or is a field or a function, the other one
must also be a field or a function (see const for how to
convert a variable to a function).