diff --git a/src/derive/incflo_derive.cpp b/src/derive/incflo_derive.cpp
index 935d227b..2311125f 100644
--- a/src/derive/incflo_derive.cpp
+++ b/src/derive/incflo_derive.cpp
@@ -265,10 +265,15 @@ void incflo::compute_nodal_inertial_num_at_level (int lev,
         amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
         {
              // Regularized Pressure
-             Real p_reg = std::sqrt(p_nd_arr(i,j,k)*p_nd_arr(i,j,k)
-                                    + eps*eps);
-             p_reg += p_nd_arr(i,j,k);
-             p_reg *= Real(0.5);
+             //Real p_reg = std::sqrt(p_nd_arr(i,j,k)*p_nd_arr(i,j,k)
+             //                       + eps*eps);
+             //p_reg += p_nd_arr(i,j,k);
+             //p_reg *= Real(0.5);
+
+             // Note: This version of Regularized Pressure only works for
+             // Static pressure, i.e., p_s >= 0
+             Real p_reg = p_nd_arr(i,j,k) + eps;
+
              // Strainrate in incflo is two-times the actual value
              inrt_num_arr(i,j,k) = std::sqrt(ro_scnd/p_reg)*
                                    diam_scnd*Real(0.5)*sr_arr(i,j,k);
diff --git a/src/rheology/incflo_rheology.cpp b/src/rheology/incflo_rheology.cpp
index 477809b5..9324a6ba 100644
--- a/src/rheology/incflo_rheology.cpp
+++ b/src/rheology/incflo_rheology.cpp
@@ -263,8 +263,7 @@ void incflo::compute_nodal_viscosity_at_level (int lev,
                    amrex::ParallelFor(bx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
                    {
                         // Regularized strain rate
-                        Real sr_reg = std::sqrt(Real(0.25)*sr_arr(i,j,k)*sr_arr(i,j,k)
-                                                + eps*eps);
+                        Real sr_reg = Real(0.5)*sr_arr(i,j,k) + eps;
                         vel_eta_snd_arr(i,j,k) *= p_static_arr(i,j,k);
                         vel_eta_snd_arr(i,j,k) /= (Real(2.0)*sr_reg);
                    });
@@ -307,8 +306,7 @@ void incflo::compute_nodal_viscosity_at_level (int lev,
                         vel_eta_snd_arr(i,j,k) += mu_1_scnd;
                         // The above value is stress ratio
                         // Regularized strain rate
-                        Real sr_reg = std::sqrt(Real(0.25)*sr_arr(i,j,k)*sr_arr(i,j,k)
-                                                + eps*eps);
+                        Real sr_reg = Real(0.5)*sr_arr(i,j,k) + eps;
                         vel_eta_snd_arr(i,j,k) *= p_static_arr(i,j,k);
                         vel_eta_snd_arr(i,j,k) /= (Real(2.0)*sr_reg);
                    });