diff --git a/apps/cpu/JetBreakup/JetBreakup.cpp b/apps/cpu/JetBreakup/JetBreakup.cpp
index aabc3f4d03541213fc7dabdd4bf1b4e9af28fc96..03c1d34c4ea67a2a7a5440d1fce4fe3affcaddf5 100644
--- a/apps/cpu/JetBreakup/JetBreakup.cpp
+++ b/apps/cpu/JetBreakup/JetBreakup.cpp
@@ -81,7 +81,7 @@ void run(string configname)
 
         // Sleep(30000);
 
-        double rho_h, rho_l, r_rho, mu_h, mu_l, Uo, D, sigma;
+        double rho_h=0, rho_l=0, r_rho=0, mu_h=0, /*mu_l,*/ Uo=0, D=0, sigma=0;
 
         switch (caseN) {
             case 1: 
@@ -94,7 +94,7 @@ void run(string configname)
                 //dynamic viscosity of heavy fluid (Pa � s)
                 mu_h = 2.87e-3;
                 //dynamic viscosity of light fluid (Pa � s)
-                mu_l = 1.97e-5;
+                //mu_l = 1.97e-5;
                 //velocity (m/s)
                 Uo = 100;
                 //diameter of jet (m)
@@ -112,7 +112,7 @@ void run(string configname)
                 // dynamic viscosity of heavy fluid (Pa � s)
                 mu_h = 2.87e-3;
                 // dynamic viscosity of light fluid (Pa � s)
-                mu_l = 1.84e-5;
+                //mu_l = 1.84e-5;
                 // velocity (m/s)
                 Uo = 200;
                 // diameter of jet (m)
@@ -130,7 +130,7 @@ void run(string configname)
                 // dynamic viscosity of heavy fluid (Pa � s)
                 mu_h = 2.87e-3;
                 // dynamic viscosity of light fluid (Pa � s)
-                mu_l = 1.84e-5;
+                //mu_l = 1.84e-5;
                 // velocity (m/s)
                 Uo = 200;
                 // diameter of jet (m)
diff --git a/apps/cpu/RisingBubble2D/RisingBubble2D.cpp b/apps/cpu/RisingBubble2D/RisingBubble2D.cpp
index f5e215fe763b90b0e753c340a00e737d05c3ca97..4483533409cffc74605cf737211401c9465231d5 100644
--- a/apps/cpu/RisingBubble2D/RisingBubble2D.cpp
+++ b/apps/cpu/RisingBubble2D/RisingBubble2D.cpp
@@ -46,7 +46,7 @@ void run(string configname)
         double cpStart     = config.getValue<double>("cpStart");
         double cpStep      = config.getValue<double>("cpStep");
         bool newStart      = config.getValue<bool>("newStart");
-        double rStep = config.getValue<double>("rStep");
+        //double rStep = config.getValue<double>("rStep");
 
         std::shared_ptr<vf::mpi::Communicator> comm = vf::mpi::MPICommunicator::getInstance();
         int myid                = comm->getProcessID();
@@ -100,13 +100,13 @@ void run(string configname)
         //density of heavy fluid
         LBMReal rho_h = 1.0;
         //density of light fluid
-        LBMReal rho_l = rho_h / r_rho;
+        //LBMReal rho_l = rho_h / r_rho;
 
         //kinimatic viscosity
         LBMReal nu_h = nuL;
         //LBMReal nu_l = nuG;
         //#dynamic viscosity
-        LBMReal mu_h = rho_h * nu_h;
+        //LBMReal mu_h = rho_h * nu_h;
         
         //gravity
         LBMReal g_y = Re * Re * nu_h * nu_h / (D*D*D);