diff --git a/src/finiteVolume/finiteVolume/divSchemes/gaussDivScheme/vectorGaussDivScheme.H b/src/finiteVolume/finiteVolume/divSchemes/gaussDivScheme/vectorGaussDivScheme.H
new file mode 100644
index 000000000..1d1a9e9ea
--- /dev/null
+++ b/src/finiteVolume/finiteVolume/divSchemes/gaussDivScheme/vectorGaussDivScheme.H
@@ -0,0 +1,79 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | foam-extend: Open Source CFD
+ \\ / O peration | Version: 3.2
+ \\ / A nd | Web: http://www.foam-extend.org
+ \\/ M anipulation | For copyright notice see file Copyright
+-------------------------------------------------------------------------------
+License
+ This file is part of foam-extend.
+
+ foam-extend is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by the
+ Free Software Foundation, either version 3 of the License, or (at your
+ option) any later version.
+
+ foam-extend is distributed in the hope that it will be useful, but
+ WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with foam-extend. If not, see .
+
+Typedef
+ Foam::fv::vectorGaussDivScheme
+
+Description
+ Specialisation of gaussDivScheme for vectors. Needed for implicit fvmDiv
+ operator for block coupled systems.
+
+SourceFiles
+ vectorGaussDivScheme.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef vectorGaussDivScheme_H
+#define vectorGaussDivScheme_H
+
+#include "gaussDivScheme.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace fv
+{
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+template<>
+tmp > gaussDivScheme::fvmUDiv
+(
+ const GeometricField&
+) const;
+
+template<>
+tmp > gaussDivScheme::fvmUDiv
+(
+ const surfaceScalarField& flux,
+ const GeometricField&
+) const;
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace fv
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //