diff --git a/appl/conjugateheattransfer/iterative/precicewrapper.cc b/appl/conjugateheattransfer/iterative/precicewrapper.cc
index ae13e3817880fd5c59c29a9cc543e039d06bb8a1..844c43daa5e1fa4f2edb1da0b7fbd1db6a888bc0 100644
--- a/appl/conjugateheattransfer/iterative/precicewrapper.cc
+++ b/appl/conjugateheattransfer/iterative/precicewrapper.cc
@@ -184,6 +184,16 @@ bool PreciceWrapper::isCoupledEntity(const int faceID) const
   return indexMapper_.isDumuxIdMapped( faceID );
 }
 
+std::vector<double>& PreciceWrapper::getHeatFluxToWrite()
+{
+  assert( wasCreated_ );
+  assert( heatFluxToWrite_ != HeatFluxType::UNDEFINED);
+  if ( writeHeatFluxType_ == HeatFluxType::FreeFlow )
+    return freeFlowHeatFlux_;
+  else
+    return solidHeatFlux_;
+}
+
 //void PreciceWrapper::readScalarQuantitiy(const int dataID, std::vector<double> &data)
 //{
 //  assert( wasCreated_ );
diff --git a/appl/conjugateheattransfer/iterative/precicewrapper.hh b/appl/conjugateheattransfer/iterative/precicewrapper.hh
index f9f18c452a2fe5a505d8fad9d7b7f9906afc66c6..98887c0930c3d5b9cf617900027e3ccdb75906ed 100644
--- a/appl/conjugateheattransfer/iterative/precicewrapper.hh
+++ b/appl/conjugateheattransfer/iterative/precicewrapper.hh
@@ -105,6 +105,8 @@ public:
   bool isCoupledEntity( const int faceID ) const;
 
 
+  std::vector<double>& getHeatFluxToWrite();
+
 //  static void readScalarQuantitiy( const int dataID, std::vector<double>& data );
 //  static void writeScalarQuantitiy( const int dataID, std::vector<double>& data );