fvOptions/.../codedSource: Corrected application of phase sources

src/fvOptions/sources/general/codedSource/CodedSource.C

@@ -210,6 +210,26 @@ void Foam::fv::CodedSource<Type>::addSup
 }
 
 
+template<class Type>
+void Foam::fv::CodedSource<Type>::addSup
+(
+    const volScalarField& alpha,
+    const volScalarField& rho,
+    fvMatrix<Type>& eqn,
+    const label fieldi
+) const
+{
+    if (debug)
+    {
+        Info<< "CodedSource<"<< pTraits<Type>::typeName
+            << ">::addSup for source " << name_ << endl;
+    }
+
+    updateLibrary();
+    redirectFvOption().addSup(alpha, rho, eqn, fieldi);
+}
+
+
 template<class Type>
 void Foam::fv::CodedSource<Type>::constrain
 (

src/fvOptions/sources/general/codedSource/CodedSource.H

@@ -202,6 +202,16 @@ public:
                 const label fieldi
             ) const;
 
+            //- Explicit and implicit matrix contributions
+            //  to phase equation
+            virtual void addSup
+            (
+                const volScalarField& alpha,
+                const volScalarField& rho,
+                fvMatrix<Type>& eqn,
+                const label fieldi
+            ) const;
+
             //- Set value
             virtual void constrain
             (