diff --git a/dune/phasefield/porenetwork/pn_2pflow.hh b/dune/phasefield/porenetwork/pn_2pflow.hh
index 9af50e0faef5a1a362423201548197589557f7d7..c1dfca88f5b46519d272ff18cb7de188bc4e3fc2 100644
--- a/dune/phasefield/porenetwork/pn_2pflow.hh
+++ b/dune/phasefield/porenetwork/pn_2pflow.hh
@@ -349,7 +349,7 @@ private:
   void findStationaryState()
   {
 
-    writeVTK(0);
+    //writeVTK(0);
     Dune::FieldVector<double,1> integral;
     param.reaction.reactionRate = 0;
     int iter = 0;
@@ -358,6 +358,7 @@ private:
 
     do
     {
+      std::cout << "Find Stationary State, Iteration " << iter << " of (maximum) " << param.PCS_Search_MaxCHNSIter << std::endl;
       iter++;
       applyChPressureSaturation(1);
       applyNsPressureSaturation(1);
@@ -440,14 +441,12 @@ public:
 
     initializeGridOperators();
 
-/*    applyNs();
-    updatePhiSolid();
+    findStationaryState();
+    calculateKf();
+    updateSaturation();
     initialPhiSolid = phiSolid;
-    updateKf();
-*/
-    //findStationaryState();
-    updateKf();
-    std::cout << kf_11 << ", "
+
+    std::cout << "Pn_2PFlow Initialized, initial Kf (viscosity ratio possibly not right)" << kf_11 << ", "
     << kf_12 <<  ", "
     << kf_21 << ", "
     << kf_22 << std::endl;
@@ -640,6 +639,7 @@ void calculateKf(RF ViscosityRatio = 0)
 
 void findSaturation(RF Curvature)
 {
+  std::cout << "Pn_2PFlow, function find Saturation not yet implemented" << std::endl;
   //TODO!!!!
 }