diff --git a/applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H
new file mode 100644
index 0000000000..68cd06ed56
--- /dev/null
+++ b/applications/solvers/multiphase/compressibleMultiphaseInterFoam/contErr.H
@@ -0,0 +1,15 @@
+volScalarField::Internal contErr
+(
+    (fvc::ddt(rho) + fvc::div(mixture.rhoPhi()))()
+);
+
+forAllConstIter
+(
+    PtrDictionary<phaseModel>,
+    mixture.phases(),
+    phase
+)
+{
+    const volScalarField& rho = phase().thermo().rho()();
+    contErr -= (fvModels.source(phase(), rho)&rho);
+}